wxvdriver.py 30 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903
  1. """!
  2. @package wxvdriver.py
  3. @brief wxGUI vector digitizer (display driver)
  4. Code based on wxVdigit C++ component from GRASS 6.4.0
  5. (gui/wxpython/vdigit). Converted to Python in 2010/12-2011/01.
  6. List of classes:
  7. - DisplayDriver
  8. (C) 2007-2011 by the GRASS Development Team
  9. This program is free software under the GNU General Public License
  10. (>=v2). Read the file COPYING that comes with GRASS for details.
  11. @author Martin Landa <landa.martin gmail.com>
  12. """
  13. import math
  14. import wx
  15. from debug import Debug
  16. from preferences import globalSettings as UserSettings
  17. from grass.lib.grass import *
  18. from grass.lib.vector import *
  19. from grass.lib.vedit import *
  20. log = None
  21. progress = None
  22. def print_error(msg, type):
  23. """!Redirect stderr"""
  24. global log
  25. if log:
  26. log.write(msg)
  27. else:
  28. print msg
  29. return 0
  30. def print_progress(value):
  31. """!Redirect progress info"""
  32. global progress
  33. if progress:
  34. progress.SetValue(value)
  35. else:
  36. print value
  37. return 0
  38. errtype = CFUNCTYPE(UNCHECKED(c_int), String, c_int)
  39. errfunc = errtype(print_error)
  40. pertype = CFUNCTYPE(UNCHECKED(c_int), c_int)
  41. perfunc = pertype(print_progress)
  42. class DisplayDriver:
  43. def __init__(self, device, deviceTmp, mapObj, window, glog, gprogress):
  44. """!Display driver used by vector digitizer
  45. @param device wx.PseudoDC device where to draw vector objects
  46. @param deviceTmp wx.PseudoDC device where to draw temporary vector objects
  47. @param mapOng Map Object (render.Map)
  48. @param windiow parent window for dialogs
  49. @param glog logging device (None to discard messages)
  50. @param gprogress progress bar device (None to discard message)
  51. """
  52. global errfunc, perfunc, log, progress
  53. log = glog
  54. progress = gprogress
  55. G_gisinit('') # initialize GRASS libs
  56. G_set_error_routine(errfunc)
  57. G_set_percent_routine(perfunc)
  58. self.mapInfo = None # open vector map (Map_Info structure)
  59. self.poMapInfo = None # pointer to self.mapInfo
  60. self.is3D = False # is open vector map 3D
  61. self.dc = device # PseudoDC devices
  62. self.dcTmp = deviceTmp
  63. self.mapObj = mapObj
  64. self.region = mapObj.GetCurrentRegion()
  65. self.window = window
  66. self.log = log # log device
  67. # GRASS lib
  68. self.poPoints = Vect_new_line_struct()
  69. self.poCats = Vect_new_cats_struct()
  70. # selected objects
  71. self.selected = {
  72. 'field' : -1, # field number
  73. 'cats' : list(), # list of cats
  74. 'ids' : list(), # list of ids
  75. 'idsDupl' : list(), # list of duplicated features
  76. }
  77. # digitizer settings
  78. self.settings = {
  79. 'highlight' : None,
  80. 'highlightDupl' : { 'enabled' : False,
  81. 'color' : None },
  82. 'point' : { 'enabled' : False,
  83. 'color' : None },
  84. 'line' : { 'enabled' : False,
  85. 'color' : None },
  86. 'boundaryNo' : { 'enabled' : False,
  87. 'color' : None },
  88. 'boundaryOne' : { 'enabled' : False,
  89. 'color' : None },
  90. 'boundaryTwo' : { 'enabled' : False,
  91. 'color' : None },
  92. 'centroidIn' : { 'enabled' : False,
  93. 'color' : None },
  94. 'centroidOut' : { 'enabled' : False,
  95. 'color' : None },
  96. 'centroidDup' : { 'enabled' : False,
  97. 'color' : None },
  98. 'nodeOne' : { 'enabled' : False,
  99. 'color' : None },
  100. 'nodeTwo' : { 'enabled' : False,
  101. 'color' : None },
  102. 'vertex' : { 'enabled' : False,
  103. 'color' : None },
  104. 'area' : { 'enabled' : False,
  105. 'color' : None },
  106. 'direction' : { 'enabled' : False,
  107. 'color' : None },
  108. 'lineWidth' : -1, # screen units
  109. }
  110. # topology
  111. self._resetTopology()
  112. self._drawSelected = False
  113. self._drawSegments = False
  114. self.UpdateSettings()
  115. def __del__(self):
  116. """!Close currently open vector map"""
  117. G_unset_error_routine()
  118. G_unset_percent_routine()
  119. if self.poMapInfo:
  120. self.CloseMap()
  121. Vect_destroy_line_struct(self.poPoints)
  122. Vect_destroy_cats_struct(self.poCats)
  123. def _resetTopology(self):
  124. """!Reset topology dict
  125. """
  126. self.topology = {
  127. 'highlight' : 0,
  128. 'point' : 0,
  129. 'line' : 0,
  130. 'boundaryNo' : 0,
  131. 'boundaryOne' : 0,
  132. 'boundaryTwo' : 0,
  133. 'centroidIn' : 0,
  134. 'centroidOut' : 0,
  135. 'centroidDup' : 0,
  136. 'nodeOne' : 0,
  137. 'nodeTwo' : 0,
  138. 'vertex' : 0,
  139. }
  140. def _cell2Pixel(self, east, north, elev):
  141. """!Conversion from geographic coordinates (east, north)
  142. to screen (x, y)
  143. @todo 3D stuff...
  144. @param east, north, elev geographical coordinates
  145. @return x, y screen coordinates (integer)
  146. """
  147. map_res = max(self.region['ewres'], self.region['nsres'])
  148. w = self.region['center_easting'] - (self.mapObj.width / 2) * map_res
  149. n = self.region['center_northing'] + (self.mapObj.height / 2) * map_res
  150. return int((east - w) / map_res), int((n - north) / map_res)
  151. def _drawCross(self, pdc, point, size = 5):
  152. """!Draw cross symbol of given size to device content
  153. Used for points, nodes, vertices
  154. @param[in,out] PseudoDC where to draw
  155. @param point coordinates of center
  156. @param size size of the cross symbol
  157. @return 0 on success
  158. @return -1 on failure
  159. """
  160. if not pdc or not point:
  161. return -1
  162. pdc.DrawLine(point.x - size, point.y, point.x + size, point.y)
  163. pdc.DrawLine(point.x, point.y - size, point.x, point.y + size)
  164. return 0
  165. def _drawObject(self, robj):
  166. """!Draw given object to the device
  167. The object is defined as robject() from vedit.h.
  168. @param robj object to be rendered
  169. @return 1 on success
  170. @return -1 on failure (vector feature marked as dead, etc.)
  171. """
  172. if not self.dc or not self.dcTmp:
  173. return -1
  174. Debug.msg(3, "_drawObject(): type=%d npoints=%d", robj.type, robj.npoints)
  175. brush = None
  176. if self._isSelected(robj.fid):
  177. pdc = self.dcTmp
  178. if self.settings['highlightDupl']['enabled'] and self._isDuplicated(robj.fid):
  179. pen = wx.Pen(self.settings['highlightDupl'], self.settings['lineWidth'], wx.SOLID)
  180. else:
  181. pen = wx.Pen(self.settings['highlight'], self.settings['lineWidth'], wx.SOLID)
  182. dcId = 1
  183. self.topology['highlight'] += 1
  184. if not self._drawSelected:
  185. return
  186. else:
  187. pdc = self.dc
  188. pen, brush = self._definePen(robj.type)
  189. dcId = 0
  190. pdc.SetPen(pen)
  191. if brush:
  192. pdc.SetBrush(brush)
  193. if robj.type & (TYPE_POINT | TYPE_CENTROIDIN | TYPE_CENTROIDOUT | TYPE_CENTROIDDUP |
  194. TYPE_NODEONE | TYPE_NODETWO | TYPE_VERTEX): # -> point
  195. for i in range(robj.npoints):
  196. p = robj.point[i]
  197. self._drawCross(pdc, p)
  198. else:
  199. if dcId > 0 and self._drawSegments:
  200. dcId = 2 # first segment
  201. i = 0
  202. while i < robj.npoints - 1:
  203. point_beg = wx.Point(robj.point[i].x, robj.point[i].y)
  204. point_end = wx.Point(robj.point[i+1].x, robj.point[i+1].y)
  205. pdc.SetId(dcId) # set unique id & set bbox for each segment
  206. pdc.SetPen(pen)
  207. pdc.SetIdBounds(dcId - 1, wx.Rect(point_beg.x, point_beg.y, 0, 0))
  208. pdc.SetIdBounds(dcId, wx.RectPP(point_beg, point_end))
  209. pdc.DrawLine(point_beg.x, point_beg.y,
  210. point_end.x, point_end.y)
  211. i += 1
  212. dcId += 2
  213. pdc.SetIdBounds(dcId - 1, wx.Rect(robj.point[robj.npoints - 1].x,
  214. robj.point[robj.npoints - 1].y,
  215. 0, 0))
  216. else:
  217. points = list()
  218. for i in range(robj.npoints):
  219. p = robj.point[i]
  220. points.append(wx.Point(p.x, p.y))
  221. if robj.type == TYPE_AREA:
  222. pdc.DrawPolygon(points)
  223. else:
  224. pdc.DrawLines(points)
  225. def _definePen(self, rtype):
  226. """!Define pen/brush based on rendered object)
  227. Updates also self.topology dict
  228. @return pen, brush
  229. """
  230. if rtype == TYPE_POINT:
  231. key = 'point'
  232. elif rtype == TYPE_LINE:
  233. key = 'line'
  234. elif rtype == TYPE_BOUNDARYNO:
  235. key = 'boundaryNo'
  236. elif rtype == TYPE_BOUNDARYTWO:
  237. key = 'boundaryTwo'
  238. elif rtype == TYPE_BOUNDARYONE:
  239. key = 'boundaryOne'
  240. elif rtype == TYPE_CENTROIDIN:
  241. key = 'centroidIn'
  242. elif rtype == TYPE_CENTROIDOUT:
  243. key = 'centroidOut'
  244. elif rtype == TYPE_CENTROIDDUP:
  245. key = 'centroidDup'
  246. elif rtype == TYPE_NODEONE:
  247. key = 'nodeOne'
  248. elif rtype == TYPE_NODETWO:
  249. key = 'nodeTwo'
  250. elif rtype == TYPE_VERTEX:
  251. key = 'vertex'
  252. elif rtype == TYPE_AREA:
  253. key = 'area'
  254. elif rtype == TYPE_ISLE:
  255. key = 'isle'
  256. elif rtype == TYPE_DIRECTION:
  257. key = 'direction'
  258. if key not in ('direction', 'area', 'isle'):
  259. self.topology[key] += 1
  260. if key in ('area', 'isle'):
  261. pen = wx.TRANSPARENT_PEN
  262. if key == 'area':
  263. brush = wx.Brush(self.settings[key]['color'], wx.SOLID)
  264. else:
  265. brush = wx.TRANSPARENT_BRUSH
  266. else:
  267. pen = wx.Pen(self.settings[key]['color'], self.settings['lineWidth'], wx.SOLID)
  268. brush = None
  269. return pen, brush
  270. def _getDrawFlag(self):
  271. """!Get draw flag from the settings
  272. See vedit.h for list of draw flags.
  273. @return draw flag (int)
  274. """
  275. ret = 0
  276. if self.settings['point']['enabled']:
  277. ret |= DRAW_POINT
  278. if self.settings['line']['enabled']:
  279. ret |= DRAW_LINE
  280. if self.settings['boundaryNo']['enabled']:
  281. ret |= DRAW_BOUNDARYNO
  282. if self.settings['boundaryTwo']['enabled']:
  283. ret |= DRAW_BOUNDARYTWO
  284. if self.settings['boundaryOne']['enabled']:
  285. ret |= DRAW_BOUNDARYONE
  286. if self.settings['centroidIn']['enabled']:
  287. ret |= DRAW_CENTROIDIN
  288. if self.settings['centroidOut']['enabled']:
  289. ret |= DRAW_CENTROIDOUT
  290. if self.settings['centroidDup']['enabled']:
  291. ret |= DRAW_CENTROIDDUP
  292. if self.settings['nodeOne']['enabled']:
  293. ret |= DRAW_NODEONE
  294. if self.settings['nodeTwo']['enabled']:
  295. ret |= DRAW_NODETWO
  296. if self.settings['vertex']['enabled']:
  297. ret |= DRAW_VERTEX
  298. if self.settings['area']['enabled']:
  299. ret |= DRAW_AREA
  300. if self.settings['direction']['enabled']:
  301. ret |= DRAW_DIRECTION
  302. return ret
  303. def _isSelected(self, line, force = False):
  304. """!Check if vector object selected?
  305. @param line feature id
  306. @return True if vector object is selected
  307. @return False if vector object is not selected
  308. """
  309. if len(self.selected['cats']) < 1 or force:
  310. # select by id
  311. if line in self.selected['ids']:
  312. return True
  313. else:
  314. # select by cat
  315. cats = self.poCats.contents
  316. for i in range(cats.n_cats):
  317. if cats.field[i] == self.selected['field'] and \
  318. cats.cat[i] in self.selected['cats']:
  319. # remember id
  320. # -> after drawing all features selected.cats is reseted */
  321. self.selected['ids'].append(line)
  322. return True
  323. return False
  324. def _isDuplicated(self, featId):
  325. return False
  326. def _getRegionBox(self):
  327. """!Get bound_box() from current region
  328. @return bound_box
  329. """
  330. box = bound_box()
  331. box.N = self.region['n']
  332. box.S = self.region['s']
  333. box.E = self.region['e']
  334. box.W = self.region['w']
  335. box.T = PORT_DOUBLE_MAX
  336. box.B = -PORT_DOUBLE_MAX
  337. return box
  338. def DrawMap(self, force = False):
  339. """!Draw content of the vector map to the device
  340. @param force force drawing
  341. @return number of drawn features
  342. @return -1 on error
  343. """
  344. Debug.msg(1, "DisplayDriver.DrawMap(): force=%d", force)
  345. if not self.poMapInfo or not self.dc or not self.dcTmp:
  346. return -1
  347. rlist = Vedit_render_map(self.poMapInfo, byref(self._getRegionBox()), self._getDrawFlag(),
  348. self.region['center_easting'], self.region['center_northing'],
  349. self.mapObj.width, self.mapObj.height,
  350. max(self.region['nsres'], self.region['ewres'])).contents
  351. self._resetTopology()
  352. self.dc.BeginDrawing()
  353. self.dcTmp.BeginDrawing()
  354. # draw objects
  355. for i in range(rlist.nitems):
  356. robj = rlist.item[i].contents
  357. self._drawObject(robj)
  358. self.dc.EndDrawing()
  359. self.dcTmp.EndDrawing()
  360. # reset list of selected features by cat
  361. # list of ids - see IsSelected()
  362. self.selected['field'] = -1
  363. self.selected['cats'] = list()
  364. def _getSelectType(self):
  365. """!Get type(s) to be selected
  366. Used by SelectLinesByBox() and SelectLinesByPoint()
  367. """
  368. ftype = 0
  369. for feature in (('point', GV_POINT),
  370. ('line', GV_LINE),
  371. ('centroid', GV_CENTROID),
  372. ('boundary', GV_BOUNDARY)):
  373. if UserSettings.Get(group = 'vdigit', key = 'selectType',
  374. subkey = [feature[0], 'enabled']):
  375. ftype |= feature[1]
  376. return ftype
  377. def _validLine(self, line):
  378. """!Check if feature id is valid
  379. @param line feature id
  380. @return True valid feature id
  381. @return False invalid
  382. """
  383. if line > 0 and line <= Vect_get_num_lines(self.poMapInfo):
  384. return True
  385. return False
  386. def SelectLinesByBox(self, bbox, drawSeg = False, poMapInfo = None):
  387. """!Select vector objects by given bounding box
  388. If line id is already in the list of selected lines, then it will
  389. be excluded from this list.
  390. @param bbox bounding box definition
  391. @param drawSeg True to draw segments of line
  392. @param poMapInfo use external Map_info, None for self.poMapInfo
  393. @return number of selected features
  394. @return None on error
  395. """
  396. thisMapInfo = poMapInfo is None
  397. if not poMapInfo:
  398. poMapInfo = self.poMapInfo
  399. if not poMapInfo:
  400. return None
  401. if thisMapInfo:
  402. self._drawSegments = drawSeg
  403. self._drawSelected = True
  404. # select by ids
  405. self.selected['cats'] = list()
  406. if thisMapInfo:
  407. selected = self.selected['ids']
  408. else:
  409. selected = list()
  410. poList = Vect_new_list()
  411. x1, y1 = bbox[0]
  412. x2, y2 = bbox[1]
  413. poBbox = Vect_new_line_struct()
  414. Vect_append_point(poBbox, x1, y1, 0.0)
  415. Vect_append_point(poBbox, x2, y1, 0.0)
  416. Vect_append_point(poBbox, x2, y2, 0.0)
  417. Vect_append_point(poBbox, x1, y2, 0.0)
  418. Vect_append_point(poBbox, x1, y1, 0.0)
  419. Vect_select_lines_by_polygon(poMapInfo, poBbox,
  420. 0, None, # isles
  421. self._getSelectType(), poList)
  422. flist = poList.contents
  423. nlines = flist.n_values
  424. for i in range(nlines):
  425. line = flist.value[i]
  426. if UserSettings.Get(group = 'vdigit', key = 'selectInside',
  427. subkey = 'enabled'):
  428. inside = True
  429. if not self._validLine(line):
  430. return None
  431. Vect_read_line(poMapInfo, self.poPoints, None, line)
  432. points = self.poPoints.contents
  433. for p in range(points.n_points):
  434. if not Vect_point_in_poly(points.x[p], points.y[p],
  435. poBbox):
  436. inside = False
  437. break
  438. if not inside:
  439. continue # skip lines just overlapping bbox
  440. if not self._isSelected(line):
  441. selected.append(line)
  442. else:
  443. del selected[line]
  444. Vect_destroy_line_struct(poBbox)
  445. Vect_destroy_list(poList)
  446. return len(selected)
  447. def SelectLineByPoint(self, point, poMapInfo = None):
  448. """!Select vector feature by given point in given
  449. threshold
  450. Only one vector object can be selected. Bounding boxes of
  451. all segments are stores.
  452. @param point points coordinates (x, y)
  453. @param poMapInfo use external Map_info, None for self.poMapInfo
  454. @return dict {'line' : feature id, 'point' : point on line}
  455. @return None nothing found
  456. """
  457. thisMapInfo = poMapInfo is None
  458. if not poMapInfo:
  459. poMapInfo = self.poMapInfo
  460. if not poMapInfo:
  461. return None
  462. if thisMapInfo:
  463. self._drawSelected = True
  464. # select by ids
  465. self.selected['cats'] = list()
  466. if thisMapInfo:
  467. selected = self.selected['ids']
  468. else:
  469. selected = list()
  470. poFound = Vect_new_list()
  471. line_nearest = Vect_find_line_list(poMapInfo, point[0], point[1], 0,
  472. self._getSelectType(), self.GetThreshold(), self.is3D,
  473. None, poFound)
  474. if line_nearest > 0:
  475. if not self._isSelected(line_nearest):
  476. selected.append(line_nearest)
  477. else:
  478. del selected[line_nearest]
  479. px = c_double()
  480. py = c_double()
  481. pz = c_double()
  482. if not self._validLine(line_nearest):
  483. return None
  484. ftype = Vect_read_line(poMapInfo, self.poPoints, self.poCats, line_nearest)
  485. Vect_line_distance (self.poPoints, point[0], point[1], 0.0, self.is3D,
  486. byref(px), byref(py), byref(pz),
  487. None, None, None)
  488. # check for duplicates
  489. if self.settings['highlightDupl']['enabled']:
  490. found = poFound.contents
  491. for i in range(found.n_values):
  492. line = found.value[i]
  493. if line != line_nearest:
  494. selected.append(line)
  495. self.getDuplicates()
  496. for i in range(found.n_values):
  497. line = found.value[i]
  498. if line != line_nearest and not self._isDuplicated(line):
  499. del selected[line]
  500. Vect_destroy_list(poFound)
  501. if thisMapInfo:
  502. # drawing segments can be very expensive
  503. # only one features selected
  504. self._drawSegments = True
  505. return { 'line' : line_nearest,
  506. 'point' : (px.value, py.value, pz.value) }
  507. def _listToIList(self, plist):
  508. """!Generate from list struct_ilist
  509. """
  510. ilist = Vect_new_list()
  511. for val in plist:
  512. Vect_list_append(ilist, val)
  513. return ilist
  514. def GetSelectedIList(self, ilist = None):
  515. """!Get list of selected objects as struct_ilist
  516. Returned IList must be freed by Vect_destroy_list().
  517. @return struct_ilist
  518. """
  519. if ilist:
  520. return self._listToIList(ilist)
  521. return self._listToIList(self.selected['ids'])
  522. def GetSelected(self, grassId = True):
  523. """!Get ids of selected objects
  524. @param grassId True for feature id, False for PseudoDC id
  525. @return list of ids of selected vector objects
  526. """
  527. if grassId:
  528. return self.selected['ids']
  529. dc_ids = list()
  530. if not self._drawSegments:
  531. dc_ids.append(1)
  532. elif len(self.selected['ids']) > 0:
  533. # only first selected feature
  534. Vect_read_line(self.poMapInfo, self.poPoints, None,
  535. self.selected['ids'][0])
  536. points = self.poPoints.contents
  537. # node - segment - vertex - segment - node
  538. for i in range(1, 2 * points.n_points):
  539. dc_ids.append(i)
  540. return dc_ids
  541. def SetSelected(self, ids, layer = -1):
  542. """!Set selected vector objects
  543. @param list of ids (None to unselect features)
  544. @param layer layer number for features selected based on category number
  545. """
  546. if ids:
  547. self._drawSelected = True
  548. else:
  549. self._drawSelected = False
  550. if layer > 0:
  551. selected.field = layer
  552. self.selected['cats'] = ids
  553. else:
  554. field = -1
  555. self.selected['ids'] = ids
  556. def GetSelectedVertex(self, pos):
  557. """!Get PseudoDC vertex id of selected line
  558. Set bounding box for vertices of line.
  559. @param pos position
  560. @return id of center, left and right vertex
  561. @return 0 no line found
  562. @return -1 on error
  563. """
  564. returnId = list()
  565. # only one object can be selected
  566. if len(self.selected['ids']) != 1 or not self._drawSegments:
  567. return returnId
  568. startId = 1
  569. line = self.selected['ids'][0]
  570. if not self._validLine(line):
  571. return -1
  572. ftype = Vect_read_line(self.poMapInfo, self.poPoints, self.poCats, line)
  573. minDist = 0.0
  574. Gid = -1
  575. # find the closest vertex (x, y)
  576. DCid = 1
  577. points = self.poPoints.contents
  578. for idx in range(points.n_points):
  579. dist = Vect_points_distance(pos[0], pos[1], 0.0,
  580. points.x[idx], points.y[idx], points.z[idx], 0)
  581. if idx == 0:
  582. minDist = dist
  583. Gid = idx
  584. else:
  585. if minDist > dist:
  586. minDist = dist
  587. Gid = idx
  588. vx, vy = self._cell2Pixel(points.x[idx], points.y[idx], points.z[idx])
  589. rect = wx.Rect(vx, vy, 0, 0)
  590. self.dc.SetIdBounds(DCid, rect)
  591. DCid += 2
  592. if minDist > self.GetThreshold():
  593. return returnId
  594. # translate id
  595. DCid = Gid * 2 + 1
  596. # add selected vertex
  597. returnId.append(DCid)
  598. # left vertex
  599. if DCid == startId:
  600. returnId.append(-1)
  601. else:
  602. returnId.append(DCid - 2)
  603. # right vertex
  604. if DCid == (points.n_points - 1) * 2 + startId:
  605. returnId.append(-1)
  606. else:
  607. returnId.append(DCid + 2)
  608. return returnId
  609. def DrawSelected(self, flag):
  610. """!Draw selected features
  611. @param flag True to draw selected features
  612. """
  613. self._drawSelected = bool(flag)
  614. def CloseMap(self):
  615. """!Close vector map
  616. @return 0 on success
  617. @return non-zero on error
  618. """
  619. ret = 0
  620. if self.poMapInfo:
  621. # rebuild topology
  622. Vect_build_partial(self.poMapInfo, GV_BUILD_NONE)
  623. Vect_build(self.poMapInfo)
  624. # close map and store topo/cidx
  625. ret = Vect_close(self.poMapInfo)
  626. del self.mapInfo
  627. self.poMapInfo = self.mapInfo = None
  628. return ret
  629. def OpenMap(self, name, mapset, update = True):
  630. """!Open vector map by the driver
  631. @param name name of vector map to be open
  632. @param mapset name of mapset where the vector map lives
  633. @return map_info
  634. @return None on error
  635. """
  636. Debug.msg("DisplayDriver.OpenMap(): name=%s mapset=%s updated=%d",
  637. name, mapset, update)
  638. if not self.mapInfo:
  639. self.mapInfo = Map_info()
  640. self.poMapInfo = pointer(self.mapInfo)
  641. # open existing map
  642. if update:
  643. ret = Vect_open_update(self.poMapInfo, name, mapset)
  644. else:
  645. ret = Vect_open_old(self.poMapInfo, name, mapset)
  646. self.is3D = Vect_is_3d(self.poMapInfo)
  647. if ret == -1: # error
  648. del self.mapInfo
  649. self.poMapInfo = self.mapInfo = None
  650. elif ret < 2:
  651. dlg = wx.MessageDialog(parent = self.window,
  652. message = _("Topology for vector map <%s> is not available. "
  653. "Topology is required by digitizer. Do you want to "
  654. "rebuild topology (takes some time) and open the vector map "
  655. "for editing?") % name,
  656. caption=_("Topology missing"),
  657. style = wx.YES_NO | wx.YES_DEFAULT | wx.ICON_QUESTION | wx.CENTRE)
  658. ret = dlg.ShowModal()
  659. if ret != wx.ID_YES:
  660. del self.mapInfo
  661. self.poMapInfo = self.mapInfo = None
  662. else:
  663. Vect_build(self.poMapInfo)
  664. return self.poMapInfo
  665. def GetMapBoundingBox(self):
  666. """!Get bounding box of (opened) vector map layer
  667. @return (w,s,b,e,n,t)
  668. """
  669. if not self.poMapInfo:
  670. return None
  671. bbox = bound_box()
  672. Vect_get_map_box(self.poMapInfo, byref(bbox))
  673. return bbox.W, bbox.S, bbox.B, \
  674. bbox.E, bbox.N, bbox.T
  675. def UpdateSettings(self, alpha = 255):
  676. """!Update display driver settings
  677. @todo map units
  678. @alpha color value for aplha channel
  679. """
  680. color = dict()
  681. for key in self.settings.keys():
  682. if key == 'lineWidth':
  683. self.settings[key] = int(UserSettings.Get(group = 'vdigit', key = 'lineWidth',
  684. subkey = 'value'))
  685. continue
  686. color = wx.Color(UserSettings.Get(group = 'vdigit', key = 'symbol',
  687. subkey = [key, 'color'])[0],
  688. UserSettings.Get(group = 'vdigit', key = 'symbol',
  689. subkey = [key, 'color'])[1],
  690. UserSettings.Get(group = 'vdigit', key = 'symbol',
  691. subkey = [key, 'color'])[2])
  692. if key == 'highlight':
  693. self.settings[key] = color
  694. continue
  695. if key == 'highlightDupl':
  696. self.settings[key]['enabled'] = bool(UserSettings.Get(group = 'vdigit', key = 'checkForDupl',
  697. subkey = 'enabled'))
  698. else:
  699. self.settings[key]['enabled'] = bool(UserSettings.Get(group = 'vdigit', key = 'symbol',
  700. subkey = [key, 'enabled']))
  701. self.settings[key]['color'] = color
  702. def UpdateRegion(self):
  703. """!Update geographical region used by display driver
  704. """
  705. self.region = self.mapObj.GetCurrentRegion()
  706. def GetThreshold(self, type = 'snapping', value = None, units = None):
  707. """!Return threshold value in map units
  708. @param type snapping mode (node, vertex)
  709. @param value threshold to be set up
  710. @param units units (map, screen)
  711. @return threshold value
  712. """
  713. if value is None:
  714. value = UserSettings.Get(group = 'vdigit', key = type, subkey = 'value')
  715. if units is None:
  716. units = UserSettings.Get(group = 'vdigit', key = type, subkey = 'units')
  717. if value < 0:
  718. value = (self.region['nsres'] + self.region['ewres']) / 2.0
  719. if units == "screen pixels":
  720. # pixel -> cell
  721. res = max(self.region['nsres'], self.region['ewres'])
  722. return value * res
  723. return value