mapdisp_window.py 100 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697989910010110210310410510610710810911011111211311411511611711811912012112212312412512612712812913013113213313413513613713813914014114214314414514614714814915015115215315415515615715815916016116216316416516616716816917017117217317417517617717817918018118218318418518618718818919019119219319419519619719819920020120220320420520620720820921021121221321421521621721821922022122222322422522622722822923023123223323423523623723823924024124224324424524624724824925025125225325425525625725825926026126226326426526626726826927027127227327427527627727827928028128228328428528628728828929029129229329429529629729829930030130230330430530630730830931031131231331431531631731831932032132232332432532632732832933033133233333433533633733833934034134234334434534634734834935035135235335435535635735835936036136236336436536636736836937037137237337437537637737837938038138238338438538638738838939039139239339439539639739839940040140240340440540640740840941041141241341441541641741841942042142242342442542642742842943043143243343443543643743843944044144244344444544644744844945045145245345445545645745845946046146246346446546646746846947047147247347447547647747847948048148248348448548648748848949049149249349449549649749849950050150250350450550650750850951051151251351451551651751851952052152252352452552652752852953053153253353453553653753853954054154254354454554654754854955055155255355455555655755855956056156256356456556656756856957057157257357457557657757857958058158258358458558658758858959059159259359459559659759859960060160260360460560660760860961061161261361461561661761861962062162262362462562662762862963063163263363463563663763863964064164264364464564664764864965065165265365465565665765865966066166266366466566666766866967067167267367467567667767867968068168268368468568668768868969069169269369469569669769869970070170270370470570670770870971071171271371471571671771871972072172272372472572672772872973073173273373473573673773873974074174274374474574674774874975075175275375475575675775875976076176276376476576676776876977077177277377477577677777877978078178278378478578678778878979079179279379479579679779879980080180280380480580680780880981081181281381481581681781881982082182282382482582682782882983083183283383483583683783883984084184284384484584684784884985085185285385485585685785885986086186286386486586686786886987087187287387487587687787887988088188288388488588688788888989089189289389489589689789889990090190290390490590690790890991091191291391491591691791891992092192292392492592692792892993093193293393493593693793893994094194294394494594694794894995095195295395495595695795895996096196296396496596696796896997097197297397497597697797897998098198298398498598698798898999099199299399499599699799899910001001100210031004100510061007100810091010101110121013101410151016101710181019102010211022102310241025102610271028102910301031103210331034103510361037103810391040104110421043104410451046104710481049105010511052105310541055105610571058105910601061106210631064106510661067106810691070107110721073107410751076107710781079108010811082108310841085108610871088108910901091109210931094109510961097109810991100110111021103110411051106110711081109111011111112111311141115111611171118111911201121112211231124112511261127112811291130113111321133113411351136113711381139114011411142114311441145114611471148114911501151115211531154115511561157115811591160116111621163116411651166116711681169117011711172117311741175117611771178117911801181118211831184118511861187118811891190119111921193119411951196119711981199120012011202120312041205120612071208120912101211121212131214121512161217121812191220122112221223122412251226122712281229123012311232123312341235123612371238123912401241124212431244124512461247124812491250125112521253125412551256125712581259126012611262126312641265126612671268126912701271127212731274127512761277127812791280128112821283128412851286128712881289129012911292129312941295129612971298129913001301130213031304130513061307130813091310131113121313131413151316131713181319132013211322132313241325132613271328132913301331133213331334133513361337133813391340134113421343134413451346134713481349135013511352135313541355135613571358135913601361136213631364136513661367136813691370137113721373137413751376137713781379138013811382138313841385138613871388138913901391139213931394139513961397139813991400140114021403140414051406140714081409141014111412141314141415141614171418141914201421142214231424142514261427142814291430143114321433143414351436143714381439144014411442144314441445144614471448144914501451145214531454145514561457145814591460146114621463146414651466146714681469147014711472147314741475147614771478147914801481148214831484148514861487148814891490149114921493149414951496149714981499150015011502150315041505150615071508150915101511151215131514151515161517151815191520152115221523152415251526152715281529153015311532153315341535153615371538153915401541154215431544154515461547154815491550155115521553155415551556155715581559156015611562156315641565156615671568156915701571157215731574157515761577157815791580158115821583158415851586158715881589159015911592159315941595159615971598159916001601160216031604160516061607160816091610161116121613161416151616161716181619162016211622162316241625162616271628162916301631163216331634163516361637163816391640164116421643164416451646164716481649165016511652165316541655165616571658165916601661166216631664166516661667166816691670167116721673167416751676167716781679168016811682168316841685168616871688168916901691169216931694169516961697169816991700170117021703170417051706170717081709171017111712171317141715171617171718171917201721172217231724172517261727172817291730173117321733173417351736173717381739174017411742174317441745174617471748174917501751175217531754175517561757175817591760176117621763176417651766176717681769177017711772177317741775177617771778177917801781178217831784178517861787178817891790179117921793179417951796179717981799180018011802180318041805180618071808180918101811181218131814181518161817181818191820182118221823182418251826182718281829183018311832183318341835183618371838183918401841184218431844184518461847184818491850185118521853185418551856185718581859186018611862186318641865186618671868186918701871187218731874187518761877187818791880188118821883188418851886188718881889189018911892189318941895189618971898189919001901190219031904190519061907190819091910191119121913191419151916191719181919192019211922192319241925192619271928192919301931193219331934193519361937193819391940194119421943194419451946194719481949195019511952195319541955195619571958195919601961196219631964196519661967196819691970197119721973197419751976197719781979198019811982198319841985198619871988198919901991199219931994199519961997199819992000200120022003200420052006200720082009201020112012201320142015201620172018201920202021202220232024202520262027202820292030203120322033203420352036203720382039204020412042204320442045204620472048204920502051205220532054205520562057205820592060206120622063206420652066206720682069207020712072207320742075207620772078207920802081208220832084208520862087208820892090209120922093209420952096209720982099210021012102210321042105210621072108210921102111211221132114211521162117211821192120212121222123212421252126212721282129213021312132213321342135213621372138213921402141214221432144214521462147214821492150215121522153215421552156215721582159216021612162216321642165216621672168216921702171217221732174217521762177217821792180218121822183218421852186218721882189219021912192219321942195219621972198219922002201220222032204220522062207220822092210221122122213221422152216221722182219222022212222222322242225222622272228222922302231223222332234223522362237223822392240224122422243224422452246224722482249225022512252225322542255225622572258225922602261226222632264226522662267226822692270227122722273227422752276227722782279228022812282228322842285228622872288228922902291229222932294229522962297229822992300230123022303230423052306230723082309231023112312231323142315231623172318231923202321232223232324232523262327232823292330233123322333233423352336233723382339234023412342234323442345234623472348234923502351235223532354235523562357235823592360236123622363236423652366236723682369237023712372237323742375237623772378237923802381238223832384238523862387238823892390239123922393239423952396239723982399240024012402240324042405240624072408240924102411241224132414241524162417241824192420242124222423242424252426242724282429243024312432243324342435243624372438243924402441244224432444244524462447244824492450245124522453245424552456245724582459246024612462246324642465246624672468246924702471247224732474247524762477247824792480248124822483248424852486248724882489249024912492249324942495249624972498249925002501250225032504250525062507250825092510251125122513251425152516251725182519252025212522252325242525252625272528252925302531253225332534253525362537253825392540254125422543254425452546254725482549255025512552255325542555255625572558255925602561256225632564256525662567256825692570257125722573257425752576257725782579258025812582258325842585258625872588258925902591259225932594259525962597259825992600260126022603260426052606260726082609261026112612261326142615261626172618261926202621262226232624262526262627262826292630263126322633263426352636263726382639264026412642264326442645264626472648264926502651265226532654265526562657265826592660266126622663266426652666266726682669267026712672267326742675267626772678267926802681268226832684268526862687268826892690269126922693269426952696269726982699270027012702270327042705270627072708270927102711271227132714271527162717271827192720272127222723
  1. """!
  2. @package mapdisp_window.py
  3. @brief GIS map display canvas, buffered window.
  4. Classes:
  5. - MapWindow
  6. - BufferedWindow
  7. (C) 2006-2009 by the GRASS Development Team
  8. This program is free software under the GNU General Public
  9. License (>=v2). Read the file COPYING that comes with GRASS
  10. for details.
  11. @author Michael Barton
  12. @author Jachym Cepicky
  13. @author Martin Landa <landa.martin gmail.com>
  14. """
  15. import os
  16. import time
  17. import math
  18. import sys
  19. import tempfile
  20. import wx
  21. import grass.script as grass
  22. import dbm
  23. import dbm_dialogs
  24. import gdialogs
  25. import gcmd
  26. import utils
  27. import globalvar
  28. import gselect
  29. from debug import Debug
  30. from preferences import globalSettings as UserSettings
  31. from units import ConvertValue as UnitsConvertValue
  32. from vdigit import GV_LINES as VDigit_Lines_Type
  33. from vdigit import VDigitCategoryDialog
  34. from vdigit import VDigitZBulkDialog
  35. from vdigit import VDigitDuplicatesDialog
  36. from vdigit import PseudoDC
  37. class MapWindow(object):
  38. """!
  39. Abstract map window class
  40. Parent for BufferedWindow class (2D display mode) and
  41. GLWindow (3D display mode)
  42. """
  43. def __init__(self, parent, id=wx.ID_ANY,
  44. pos=wx.DefaultPosition,
  45. size=wx.DefaultSize,
  46. style=wx.NO_FULL_REPAINT_ON_RESIZE,
  47. Map=None, tree=None, lmgr=None):
  48. self.parent = parent # MapFrame
  49. #
  50. # mouse attributes -- position on the screen, begin and end of
  51. # dragging, and type of drawing
  52. #
  53. self.mouse = {
  54. 'begin': [0, 0], # screen coordinates
  55. 'end' : [0, 0],
  56. 'use' : "pointer",
  57. 'box' : "point"
  58. }
  59. def EraseMap(self):
  60. """!
  61. Erase the canvas (virtual method)
  62. """
  63. pass
  64. def UpdateMap(self):
  65. """!
  66. Updates the canvas anytime there is a change to the
  67. underlaying images or to the geometry of the canvas.
  68. """
  69. pass
  70. def OnLeftDown(self, event):
  71. pass
  72. def OnLeftUp(self, event):
  73. pass
  74. def OnMouseMotion(self, event):
  75. pass
  76. def OnZoomToMap(self, event):
  77. pass
  78. def OnZoomToRaster(self, event):
  79. pass
  80. def GetSelectedLayer(self, type = 'layer', multi = False):
  81. """!
  82. Get selected layer from layer tree
  83. @param type 'item' / 'layer' / 'nviz'
  84. @param multi return first selected layer or all
  85. @return layer / map layer properties / nviz properties
  86. @return None / [] on failure
  87. """
  88. ret = []
  89. if not self.tree or \
  90. not self.tree.GetSelection():
  91. if multi:
  92. return []
  93. else:
  94. return None
  95. if multi and \
  96. type == 'item':
  97. return self.tree.GetSelections()
  98. for item in self.tree.GetSelections():
  99. if not item.IsChecked():
  100. if multi:
  101. continue
  102. else:
  103. return None
  104. if type == 'item': # -> multi = False
  105. return item
  106. try:
  107. if type == 'nviz':
  108. layer = self.tree.GetPyData(item)[0]['nviz']
  109. else:
  110. layer = self.tree.GetPyData(item)[0]['maplayer']
  111. except:
  112. layer = None
  113. if multi:
  114. ret.append(layer)
  115. else:
  116. return layer
  117. return ret
  118. class BufferedWindow(MapWindow, wx.Window):
  119. """!
  120. A Buffered window class.
  121. When the drawing needs to change, you app needs to call the
  122. UpdateMap() method. Since the drawing is stored in a bitmap, you
  123. can also save the drawing to file by calling the
  124. SaveToFile(self,file_name,file_type) method.
  125. """
  126. def __init__(self, parent, id,
  127. pos = wx.DefaultPosition,
  128. size = wx.DefaultSize,
  129. style=wx.NO_FULL_REPAINT_ON_RESIZE,
  130. Map=None, tree=None, lmgr=None):
  131. MapWindow.__init__(self, parent, id, pos, size, style,
  132. Map, tree, lmgr)
  133. wx.Window.__init__(self, parent, id, pos, size, style)
  134. self.Map = Map
  135. self.tree = tree
  136. self.lmgr = lmgr # Layer Manager
  137. #
  138. # Flags
  139. #
  140. self.resize = False # indicates whether or not a resize event has taken place
  141. self.dragimg = None # initialize variable for map panning
  142. #
  143. # Variable for drawing on DC
  144. #
  145. self.pen = None # pen for drawing zoom boxes, etc.
  146. self.polypen = None # pen for drawing polylines (measurements, profiles, etc)
  147. # List of wx.Point tuples defining a polyline (geographical coordinates)
  148. self.polycoords = []
  149. # ID of rubber band line
  150. self.lineid = None
  151. # ID of poly line resulting from cumulative rubber band lines (e.g. measurement)
  152. self.plineid = None
  153. #
  154. # Event bindings
  155. #
  156. self.Bind(wx.EVT_PAINT, self.OnPaint)
  157. self.Bind(wx.EVT_SIZE, self.OnSize)
  158. self.Bind(wx.EVT_IDLE, self.OnIdle)
  159. self.Bind(wx.EVT_MOTION, self.MouseActions)
  160. self.Bind(wx.EVT_MOUSE_EVENTS, self.MouseActions)
  161. self.processMouse = True
  162. #
  163. # Render output objects
  164. #
  165. self.mapfile = None # image file to be rendered
  166. self.img = "" # wx.Image object (self.mapfile)
  167. # used in digitization tool (do not redraw vector map)
  168. self.imgVectorMap = None
  169. # decoration overlays
  170. self.overlays = {}
  171. # images and their PseudoDC ID's for painting and dragging
  172. self.imagedict = {}
  173. self.select = {} # selecting/unselecting decorations for dragging
  174. self.textdict = {} # text, font, and color indexed by id
  175. self.currtxtid = None # PseudoDC id for currently selected text
  176. #
  177. # Zoom objects
  178. #
  179. self.zoomhistory = [] # list of past zoom extents
  180. self.currzoom = 0 # current set of extents in zoom history being used
  181. self.zoomtype = 1 # 1 zoom in, 0 no zoom, -1 zoom out
  182. self.hitradius = 10 # distance for selecting map decorations
  183. self.dialogOffset = 5 # offset for dialog (e.g. DisplayAttributesDialog)
  184. # OnSize called to make sure the buffer is initialized.
  185. # This might result in OnSize getting called twice on some
  186. # platforms at initialization, but little harm done.
  187. ### self.OnSize(None)
  188. # create PseudoDC used for background map, map decorations like scales and legends
  189. self.pdc = PseudoDC()
  190. # used for digitization tool
  191. self.pdcVector = None
  192. # decorations (region box, etc.)
  193. self.pdcDec = PseudoDC()
  194. # pseudoDC for temporal objects (select box, measurement tool, etc.)
  195. self.pdcTmp = PseudoDC()
  196. # redraw all pdc's, pdcTmp layer is redrawn always (speed issue)
  197. self.redrawAll = True
  198. # will store an off screen empty bitmap for saving to file
  199. self._buffer = ''
  200. self.Bind(wx.EVT_ERASE_BACKGROUND, lambda x:None)
  201. # vars for handling mouse clicks
  202. self.dragid = -1
  203. self.lastpos = (0, 0)
  204. def Draw(self, pdc, img=None, drawid=None, pdctype='image', coords=[0, 0, 0, 0]):
  205. """!
  206. Draws map and overlay decorations
  207. """
  208. if drawid == None:
  209. if pdctype == 'image' and img:
  210. drawid = self.imagedict[img]
  211. elif pdctype == 'clear':
  212. drawid == None
  213. else:
  214. drawid = wx.NewId()
  215. if img and pdctype == 'image':
  216. # self.imagedict[img]['coords'] = coords
  217. self.select[self.imagedict[img]['id']] = False # ?
  218. pdc.BeginDrawing()
  219. if drawid != 99:
  220. bg = wx.TRANSPARENT_BRUSH
  221. else:
  222. bg = wx.Brush(self.GetBackgroundColour())
  223. pdc.SetBackground(bg)
  224. ### pdc.Clear()
  225. Debug.msg (5, "BufferedWindow.Draw(): id=%s, pdctype=%s, coord=%s" % \
  226. (drawid, pdctype, coords))
  227. # set PseudoDC id
  228. if drawid is not None:
  229. pdc.SetId(drawid)
  230. if pdctype == 'clear': # erase the display
  231. bg = wx.WHITE_BRUSH
  232. # bg = wx.Brush(self.GetBackgroundColour())
  233. pdc.SetBackground(bg)
  234. pdc.RemoveAll()
  235. pdc.Clear()
  236. pdc.EndDrawing()
  237. self.Refresh()
  238. return
  239. if pdctype == 'image': # draw selected image
  240. bitmap = wx.BitmapFromImage(img)
  241. w,h = bitmap.GetSize()
  242. pdc.DrawBitmap(bitmap, coords[0], coords[1], True) # draw the composite map
  243. pdc.SetIdBounds(drawid, wx.Rect(coords[0],coords[1], w, h))
  244. elif pdctype == 'box': # draw a box on top of the map
  245. if self.pen:
  246. pdc.SetBrush(wx.Brush(wx.CYAN, wx.TRANSPARENT))
  247. pdc.SetPen(self.pen)
  248. x2 = max(coords[0],coords[2])
  249. x1 = min(coords[0],coords[2])
  250. y2 = max(coords[1],coords[3])
  251. y1 = min(coords[1],coords[3])
  252. rwidth = x2-x1
  253. rheight = y2-y1
  254. rect = wx.Rect(x1, y1, rwidth, rheight)
  255. pdc.DrawRectangleRect(rect)
  256. pdc.SetIdBounds(drawid, rect)
  257. elif pdctype == 'line': # draw a line on top of the map
  258. if self.pen:
  259. pdc.SetBrush(wx.Brush(wx.CYAN, wx.TRANSPARENT))
  260. pdc.SetPen(self.pen)
  261. pdc.DrawLine(wx.Point(coords[0], coords[1]),
  262. wx.Point(coords[2], coords[3]))
  263. pdc.SetIdBounds(drawid, wx.Rect(coords[0], coords[1], coords[2], coords[3]))
  264. # self.ovlcoords[drawid] = coords
  265. elif pdctype == 'polyline': # draw a polyline on top of the map
  266. if self.polypen:
  267. pdc.SetBrush(wx.Brush(wx.CYAN, wx.TRANSPARENT))
  268. pdc.SetPen(self.polypen)
  269. ### pdc.DrawLines(coords)
  270. if (len(coords) < 2):
  271. return
  272. i = 1
  273. while i < len(coords):
  274. pdc.DrawLine(wx.Point(coords[i-1][0], coords[i-1][1]),
  275. wx.Point(coords[i][0], coords[i][1]))
  276. i += 1
  277. # get bounding rectangle for polyline
  278. xlist = []
  279. ylist = []
  280. if len(coords) > 0:
  281. for point in coords:
  282. x,y = point
  283. xlist.append(x)
  284. ylist.append(y)
  285. x1=min(xlist)
  286. x2=max(xlist)
  287. y1=min(ylist)
  288. y2=max(ylist)
  289. pdc.SetIdBounds(drawid, wx.Rect(x1,y1,x2,y2))
  290. # self.ovlcoords[drawid] = [x1,y1,x2,y2]
  291. elif pdctype == 'point': # draw point
  292. if self.pen:
  293. pdc.SetPen(self.pen)
  294. pdc.DrawPoint(coords[0], coords[1])
  295. coordsBound = (coords[0] - 5,
  296. coords[1] - 5,
  297. coords[0] + 5,
  298. coords[1] + 5)
  299. pdc.SetIdBounds(drawid, wx.Rect(coordsBound))
  300. # self.ovlcoords[drawid] = coords
  301. elif pdctype == 'text': # draw text on top of map
  302. if not img['active']:
  303. return #only draw active text
  304. if img.has_key('rotation'):
  305. rotation = float(img['rotation'])
  306. else:
  307. rotation = 0.0
  308. w, h = self.GetFullTextExtent(img['text'])[0:2]
  309. pdc.SetFont(img['font'])
  310. pdc.SetTextForeground(img['color'])
  311. coords, w, h = self.TextBounds(img)
  312. if rotation == 0:
  313. pdc.DrawText(img['text'], coords[0], coords[1])
  314. else:
  315. pdc.DrawRotatedText(img['text'], coords[0], coords[1], rotation)
  316. pdc.SetIdBounds(drawid, wx.Rect(coords[0], coords[1], w, h))
  317. pdc.EndDrawing()
  318. self.Refresh()
  319. return drawid
  320. def TextBounds(self, textinfo):
  321. """!
  322. Return text boundary data
  323. @param textinfo text metadata (text, font, color, rotation)
  324. @param coords reference point
  325. """
  326. if textinfo.has_key('rotation'):
  327. rotation = float(textinfo['rotation'])
  328. else:
  329. rotation = 0.0
  330. coords = textinfo['coords']
  331. Debug.msg (4, "BufferedWindow.TextBounds(): text=%s, rotation=%f" % \
  332. (textinfo['text'], rotation))
  333. self.Update()
  334. ### self.Refresh()
  335. self.SetFont(textinfo['font'])
  336. w, h = self.GetTextExtent(textinfo['text'])
  337. if rotation == 0:
  338. coords[2], coords[3] = coords[0] + w, coords[1] + h
  339. return coords, w, h
  340. boxh = math.fabs(math.sin(math.radians(rotation)) * w) + h
  341. boxw = math.fabs(math.cos(math.radians(rotation)) * w) + h
  342. coords[2] = coords[0] + boxw
  343. coords[3] = coords[1] + boxh
  344. return coords, boxw, boxh
  345. def OnPaint(self, event):
  346. """!
  347. Draw PseudoDC's to buffered paint DC
  348. self.pdc for background and decorations
  349. self.pdcVector for vector map which is edited
  350. self.pdcTmp for temporaly drawn objects (self.polycoords)
  351. If self.redrawAll is False on self.pdcTmp content is re-drawn
  352. """
  353. Debug.msg(4, "BufferedWindow.OnPaint(): redrawAll=%s" % self.redrawAll)
  354. dc = wx.BufferedPaintDC(self, self.buffer)
  355. ### dc.SetBackground(wx.Brush("White"))
  356. dc.Clear()
  357. # use PrepareDC to set position correctly
  358. self.PrepareDC(dc)
  359. # create a clipping rect from our position and size
  360. # and update region
  361. rgn = self.GetUpdateRegion().GetBox()
  362. dc.SetClippingRect(rgn)
  363. switchDraw = False
  364. if self.redrawAll is None:
  365. self.redrawAll = True
  366. switchDraw = True
  367. if self.redrawAll: # redraw pdc and pdcVector
  368. # draw to the dc using the calculated clipping rect
  369. self.pdc.DrawToDCClipped(dc, rgn)
  370. # draw vector map layer
  371. if self.pdcVector:
  372. # decorate with GDDC (transparency)
  373. try:
  374. gcdc = wx.GCDC(dc)
  375. self.pdcVector.DrawToDCClipped(gcdc, rgn)
  376. except NotImplementedError, e:
  377. print >> sys.stderr, e
  378. self.pdcVector.DrawToDCClipped(dc, rgn)
  379. self.bufferLast = None
  380. else: # do not redraw pdc and pdcVector
  381. if self.bufferLast is None:
  382. # draw to the dc
  383. self.pdc.DrawToDC(dc)
  384. if self.pdcVector:
  385. # decorate with GDDC (transparency)
  386. try:
  387. gcdc = wx.GCDC(dc)
  388. self.pdcVector.DrawToDC(gcdc)
  389. except NotImplementedError, e:
  390. print >> sys.stderr, e
  391. self.pdcVector.DrawToDC(dc)
  392. # store buffered image
  393. # self.bufferLast = wx.BitmapFromImage(self.buffer.ConvertToImage())
  394. self.bufferLast = dc.GetAsBitmap(wx.Rect(0, 0, self.Map.width, self.Map.height))
  395. pdcLast = PseudoDC()
  396. pdcLast.DrawBitmap(self.bufferLast, 0, 0, False)
  397. pdcLast.DrawToDC(dc)
  398. # draw decorations (e.g. region box)
  399. try:
  400. gcdc = wx.GCDC(dc)
  401. self.pdcDec.DrawToDC(gcdc)
  402. except NotImplementedError, e:
  403. print >> sys.stderr, e
  404. self.pdcDec.DrawToDC(dc)
  405. # draw temporary object on the foreground
  406. ### self.pdcTmp.DrawToDCClipped(dc, rgn)
  407. self.pdcTmp.DrawToDC(dc)
  408. if switchDraw:
  409. self.redrawAll = False
  410. def OnSize(self, event):
  411. """!
  412. Scale map image so that it is
  413. the same size as the Window
  414. """
  415. Debug.msg(3, "BufferedWindow.OnSize():")
  416. # set size of the input image
  417. self.Map.ChangeMapSize(self.GetClientSize())
  418. # align extent based on center point and display resolution
  419. # this causes that image is not resized when display windows is resized
  420. # self.Map.AlignExtentFromDisplay()
  421. # Make new off screen bitmap: this bitmap will always have the
  422. # current drawing in it, so it can be used to save the image to
  423. # a file, or whatever.
  424. self.buffer = wx.EmptyBitmap(max(1, self.Map.width), max(1, self.Map.height))
  425. # get the image to be rendered
  426. self.img = self.GetImage()
  427. # update map display
  428. if self.img and self.Map.width + self.Map.height > 0: # scale image during resize
  429. self.img = self.img.Scale(self.Map.width, self.Map.height)
  430. if len(self.Map.GetListOfLayers()) > 0:
  431. self.UpdateMap()
  432. # re-render image on idle
  433. self.resize = True
  434. # reposition checkbox in statusbar
  435. self.parent.StatusbarReposition()
  436. # update statusbar
  437. self.parent.StatusbarUpdate()
  438. def OnIdle(self, event):
  439. """!
  440. Only re-render a composite map image from GRASS during
  441. idle time instead of multiple times during resizing.
  442. """
  443. if self.resize:
  444. self.UpdateMap(render=True)
  445. event.Skip()
  446. def SaveToFile(self, FileName, FileType):
  447. """!
  448. This draws the psuedo DC to a buffer that
  449. can be saved to a file.
  450. """
  451. dc = wx.BufferedPaintDC(self, self.buffer)
  452. self.pdc.DrawToDC(dc)
  453. if self.pdcVector:
  454. self.pdcVector.DrawToDC(dc)
  455. self.buffer.SaveFile(FileName, FileType)
  456. def GetOverlay(self):
  457. """!
  458. Converts rendered overlay files to wx.Image
  459. Updates self.imagedict
  460. @return list of images
  461. """
  462. imgs = []
  463. for overlay in self.Map.GetListOfLayers(l_type="overlay", l_active=True):
  464. if os.path.isfile(overlay.mapfile) and os.path.getsize(overlay.mapfile):
  465. img = wx.Image(overlay.mapfile, wx.BITMAP_TYPE_ANY)
  466. self.imagedict[img] = { 'id' : overlay.id,
  467. 'layer' : overlay }
  468. imgs.append(img)
  469. return imgs
  470. def GetImage(self):
  471. """!
  472. Converts redered map files to wx.Image
  473. Updates self.imagedict (id=99)
  474. @return wx.Image instance (map composition)
  475. """
  476. imgId = 99
  477. if self.Map.mapfile and os.path.isfile(self.Map.mapfile) and \
  478. os.path.getsize(self.Map.mapfile):
  479. img = wx.Image(self.Map.mapfile, wx.BITMAP_TYPE_ANY)
  480. else:
  481. img = None
  482. self.imagedict[img] = { 'id': imgId }
  483. return img
  484. def UpdateMap(self, render=True, renderVector=True):
  485. """!
  486. Updates the canvas anytime there is a change to the
  487. underlaying images or to the geometry of the canvas.
  488. @param render re-render map composition
  489. @param renderVector re-render vector map layer enabled for editing (used for digitizer)
  490. """
  491. start = time.clock()
  492. self.resize = False
  493. # if len(self.Map.GetListOfLayers()) == 0:
  494. # return False
  495. if self.img is None:
  496. render = True
  497. #
  498. # initialize process bar (only on 'render')
  499. #
  500. if render is True or renderVector is True:
  501. self.parent.statusbarWin['progress'].Show()
  502. if self.parent.statusbarWin['progress'].GetRange() > 0:
  503. self.parent.statusbarWin['progress'].SetValue(1)
  504. #
  505. # render background image if needed
  506. #
  507. # update layer dictionary if there has been a change in layers
  508. if self.tree and self.tree.reorder == True:
  509. self.tree.ReorderLayers()
  510. # reset flag for auto-rendering
  511. if self.tree:
  512. self.tree.rerender = False
  513. if render:
  514. # update display size
  515. self.Map.ChangeMapSize(self.GetClientSize())
  516. if self.parent.statusbarWin['resolution'].IsChecked():
  517. # use computation region resolution for rendering
  518. windres = True
  519. else:
  520. windres = False
  521. self.mapfile = self.Map.Render(force=True, mapWindow=self.parent,
  522. windres=windres)
  523. else:
  524. self.mapfile = self.Map.Render(force=False, mapWindow=self.parent)
  525. self.img = self.GetImage() # id=99
  526. #
  527. # clear pseudoDcs
  528. #
  529. for pdc in (self.pdc,
  530. self.pdcDec,
  531. self.pdcTmp):
  532. pdc.Clear()
  533. pdc.RemoveAll()
  534. #
  535. # draw background map image to PseudoDC
  536. #
  537. if not self.img:
  538. self.Draw(self.pdc, pdctype='clear')
  539. else:
  540. try:
  541. id = self.imagedict[self.img]['id']
  542. except:
  543. return False
  544. self.Draw(self.pdc, self.img, drawid=id)
  545. #
  546. # render vector map layer
  547. #
  548. digitToolbar = self.parent.toolbars['vdigit']
  549. if renderVector and digitToolbar and \
  550. digitToolbar.GetLayer():
  551. # set region
  552. self.parent.digit.driver.UpdateRegion()
  553. # re-calculate threshold for digitization tool
  554. self.parent.digit.driver.GetThreshold()
  555. # draw map
  556. self.pdcVector.Clear()
  557. self.pdcVector.RemoveAll()
  558. try:
  559. item = self.tree.FindItemByData('maplayer', digitToolbar.GetLayer())
  560. except TypeError:
  561. item = None
  562. if item and self.tree.IsItemChecked(item):
  563. self.parent.digit.driver.DrawMap()
  564. # translate tmp objects (pointer position)
  565. if digitToolbar.GetAction() == 'moveLine':
  566. if hasattr(self, "vdigitMove") and \
  567. self.vdigitMove.has_key('beginDiff'):
  568. # move line
  569. for id in self.vdigitMove['id']:
  570. self.pdcTmp.TranslateId(id,
  571. self.vdigitMove['beginDiff'][0],
  572. self.vdigitMove['beginDiff'][1])
  573. del self.vdigitMove['beginDiff']
  574. #
  575. # render overlays
  576. #
  577. for img in self.GetOverlay():
  578. # draw any active and defined overlays
  579. if self.imagedict[img]['layer'].IsActive():
  580. id = self.imagedict[img]['id']
  581. self.Draw(self.pdc, img=img, drawid=id,
  582. pdctype=self.overlays[id]['pdcType'], coords=self.overlays[id]['coords'])
  583. for id in self.textdict.keys():
  584. self.Draw(self.pdc, img=self.textdict[id], drawid=id,
  585. pdctype='text', coords=[10, 10, 10, 10])
  586. # optionally draw computational extent box
  587. self.DrawCompRegionExtent()
  588. #
  589. # redraw pdcTmp if needed
  590. #
  591. if len(self.polycoords) > 0:
  592. self.DrawLines(self.pdcTmp)
  593. if not self.parent.IsStandalone() and \
  594. self.parent.GetLayerManager().georectifying:
  595. # -> georectifier (redraw GCPs)
  596. if self.parent.toolbars['georect']:
  597. coordtype = 'gcpcoord'
  598. else:
  599. coordtype = 'mapcoord'
  600. self.parent.GetLayerManager().georectifying.DrawGCP(coordtype)
  601. #
  602. # clear measurement
  603. #
  604. if self.mouse["use"] == "measure":
  605. self.ClearLines(pdc=self.pdcTmp)
  606. self.polycoords = []
  607. self.mouse['use'] = 'pointer'
  608. self.mouse['box'] = 'point'
  609. self.mouse['end'] = [0, 0]
  610. self.SetCursor(self.parent.cursors["default"])
  611. stop = time.clock()
  612. #
  613. # hide process bar
  614. #
  615. self.parent.statusbarWin['progress'].Hide()
  616. #
  617. # update statusbar
  618. #
  619. ### self.Map.SetRegion()
  620. self.parent.StatusbarUpdate()
  621. if grass.find_file(name = 'MASK', element = 'cell')['name']:
  622. # mask found
  623. self.parent.statusbarWin['mask'].SetLabel(_('MASK'))
  624. else:
  625. self.parent.statusbarWin['mask'].SetLabel('')
  626. Debug.msg (2, "BufferedWindow.UpdateMap(): render=%s, renderVector=%s -> time=%g" % \
  627. (render, renderVector, (stop-start)))
  628. return True
  629. def DrawCompRegionExtent(self):
  630. """!
  631. Draw computational region extent in the display
  632. Display region is drawn as a blue box inside the computational region,
  633. computational region inside a display region as a red box).
  634. """
  635. if hasattr(self, "regionCoords"):
  636. compReg = self.Map.GetRegion()
  637. dispReg = self.Map.GetCurrentRegion()
  638. reg = None
  639. if self.IsInRegion(dispReg, compReg):
  640. self.polypen = wx.Pen(colour=wx.Colour(0, 0, 255, 128), width=3, style=wx.SOLID)
  641. reg = dispReg
  642. else:
  643. self.polypen = wx.Pen(colour=wx.Colour(255, 0, 0, 128),
  644. width=3, style=wx.SOLID)
  645. reg = compReg
  646. self.regionCoords = []
  647. self.regionCoords.append((reg['w'], reg['n']))
  648. self.regionCoords.append((reg['e'], reg['n']))
  649. self.regionCoords.append((reg['e'], reg['s']))
  650. self.regionCoords.append((reg['w'], reg['s']))
  651. self.regionCoords.append((reg['w'], reg['n']))
  652. # draw region extent
  653. self.DrawLines(pdc=self.pdcDec, polycoords=self.regionCoords)
  654. def IsInRegion(self, region, refRegion):
  655. """!
  656. Test if 'region' is inside of 'refRegion'
  657. @param region input region
  658. @param refRegion reference region (e.g. computational region)
  659. @return True if region is inside of refRegion
  660. @return False
  661. """
  662. if region['s'] >= refRegion['s'] and \
  663. region['n'] <= refRegion['n'] and \
  664. region['w'] >= refRegion['w'] and \
  665. region['e'] <= refRegion['e']:
  666. return True
  667. return False
  668. def EraseMap(self):
  669. """!
  670. Erase the canvas
  671. """
  672. self.Draw(self.pdc, pdctype='clear')
  673. if self.pdcVector:
  674. self.Draw(self.pdcVector, pdctype='clear')
  675. self.Draw(self.pdcDec, pdctype='clear')
  676. self.Draw(self.pdcTmp, pdctype='clear')
  677. def DragMap(self, moveto):
  678. """!
  679. Drag the entire map image for panning.
  680. """
  681. dc = wx.BufferedDC(wx.ClientDC(self))
  682. dc.SetBackground(wx.Brush("White"))
  683. dc.Clear()
  684. self.dragimg = wx.DragImage(self.buffer)
  685. self.dragimg.BeginDrag((0, 0), self)
  686. self.dragimg.GetImageRect(moveto)
  687. self.dragimg.Move(moveto)
  688. self.dragimg.DoDrawImage(dc, moveto)
  689. self.dragimg.EndDrag()
  690. return True
  691. def DragItem(self, id, event):
  692. """!
  693. Drag an overlay decoration item
  694. """
  695. if id == 99 or id == '' or id == None: return
  696. Debug.msg (5, "BufferedWindow.DragItem(): id=%d" % id)
  697. x, y = self.lastpos
  698. dx = event.GetX() - x
  699. dy = event.GetY() - y
  700. self.pdc.SetBackground(wx.Brush(self.GetBackgroundColour()))
  701. r = self.pdc.GetIdBounds(id)
  702. ### FIXME in vdigit/pseudodc.i
  703. if type(r) is list:
  704. r = wx.Rect(r[0], r[1], r[2], r[3])
  705. if id > 100: # text dragging
  706. rtop = (r[0],r[1]-r[3],r[2],r[3])
  707. r = r.Union(rtop)
  708. rleft = (r[0]-r[2],r[1],r[2],r[3])
  709. r = r.Union(rleft)
  710. self.pdc.TranslateId(id, dx, dy)
  711. r2 = self.pdc.GetIdBounds(id)
  712. if type(r2) is list:
  713. r2 = wx.Rect(r[0], r[1], r[2], r[3])
  714. if id > 100: # text
  715. self.textdict[id]['coords'] = r2
  716. r = r.Union(r2)
  717. r.Inflate(4,4)
  718. self.RefreshRect(r, False)
  719. self.lastpos = (event.GetX(), event.GetY())
  720. def MouseDraw(self, pdc=None, begin=None, end=None):
  721. """!
  722. Mouse box or line from 'begin' to 'end'
  723. If not given from self.mouse['begin'] to self.mouse['end'].
  724. """
  725. self.redrawAll = False
  726. if not pdc:
  727. return
  728. if begin is None:
  729. begin = self.mouse['begin']
  730. if end is None:
  731. end = self.mouse['end']
  732. Debug.msg (5, "BufferedWindow.MouseDraw(): use=%s, box=%s, begin=%f,%f, end=%f,%f" % \
  733. (self.mouse['use'], self.mouse['box'],
  734. begin[0], begin[1], end[0], end[1]))
  735. if self.mouse['box'] == "box":
  736. boxid = wx.ID_NEW
  737. mousecoords = [begin[0], begin[1],
  738. end[0], end[1]]
  739. r = pdc.GetIdBounds(boxid)
  740. if type(r) is list:
  741. r = wx.Rect(r[0], r[1], r[2], r[3])
  742. r.Inflate(4, 4)
  743. pdc.ClearId(boxid)
  744. self.RefreshRect(r, False)
  745. pdc.SetId(boxid)
  746. self.Draw(pdc, drawid=boxid, pdctype='box', coords=mousecoords)
  747. elif self.mouse['box'] == "line":
  748. self.lineid = wx.ID_NEW
  749. mousecoords = [begin[0], begin[1], \
  750. end[0], end[1]]
  751. x1=min(begin[0],end[0])
  752. x2=max(begin[0],end[0])
  753. y1=min(begin[1],end[1])
  754. y2=max(begin[1],end[1])
  755. r = wx.Rect(x1,y1,x2-x1,y2-y1)
  756. r.Inflate(4,4)
  757. try:
  758. pdc.ClearId(self.lineid)
  759. except:
  760. pass
  761. self.RefreshRect(r, False)
  762. pdc.SetId(self.lineid)
  763. self.Draw(pdc, drawid=self.lineid, pdctype='line', coords=mousecoords)
  764. def DrawLines(self, pdc=None, polycoords=None):
  765. """!
  766. Draw polyline in PseudoDC
  767. Set self.pline to wx.NEW_ID + 1
  768. polycoords - list of polyline vertices, geographical coordinates
  769. (if not given, self.polycoords is used)
  770. """
  771. if not pdc:
  772. pdc = self.pdcTmp
  773. if not polycoords:
  774. polycoords = self.polycoords
  775. if len(polycoords) > 0:
  776. self.plineid = wx.ID_NEW + 1
  777. # convert from EN to XY
  778. coords = []
  779. for p in polycoords:
  780. coords.append(self.Cell2Pixel(p))
  781. self.Draw(pdc, drawid=self.plineid, pdctype='polyline', coords=coords)
  782. Debug.msg (4, "BufferedWindow.DrawLines(): coords=%s, id=%s" % \
  783. (coords, self.plineid))
  784. return self.plineid
  785. return -1
  786. def DrawCross(self, pdc, coords, size, rotation=0,
  787. text=None, textAlign='lr', textOffset=(5, 5)):
  788. """!Draw cross in PseudoDC
  789. @todo implement rotation
  790. @param pdc PseudoDC
  791. @param coord center coordinates
  792. @param rotation rotate symbol
  793. @param text draw also text (text, font, color, rotation)
  794. @param textAlign alignment (default 'lower-right')
  795. @textOffset offset for text (from center point)
  796. """
  797. Debug.msg(4, "BufferedWindow.DrawCross(): pdc=%s, coords=%s, size=%d" % \
  798. (pdc, coords, size))
  799. coordsCross = ((coords[0] - size, coords[1], coords[0] + size, coords[1]),
  800. (coords[0], coords[1] - size, coords[0], coords[1] + size))
  801. self.lineid = wx.NewId()
  802. for lineCoords in coordsCross:
  803. self.Draw(pdc, drawid=self.lineid, pdctype='line', coords=lineCoords)
  804. if not text:
  805. return self.lineid
  806. if textAlign == 'ul':
  807. coord = [coords[0] - textOffset[0], coords[1] - textOffset[1], 0, 0]
  808. elif textAlign == 'ur':
  809. coord = [coords[0] + textOffset[0], coords[1] - textOffset[1], 0, 0]
  810. elif textAlign == 'lr':
  811. coord = [coords[0] + textOffset[0], coords[1] + textOffset[1], 0, 0]
  812. else:
  813. coord = [coords[0] - textOffset[0], coords[1] + textOffset[1], 0, 0]
  814. self.Draw(pdc, img=text,
  815. pdctype='text', coords=coord)
  816. return self.lineid
  817. def MouseActions(self, event):
  818. """!
  819. Mouse motion and button click notifier
  820. """
  821. if not self.processMouse:
  822. return
  823. ### if self.redrawAll is False:
  824. ### self.redrawAll = True
  825. # zoom with mouse wheel
  826. if event.GetWheelRotation() != 0:
  827. self.OnMouseWheel(event)
  828. # left mouse button pressed
  829. elif event.LeftDown():
  830. self.OnLeftDown(event)
  831. # left mouse button released
  832. elif event.LeftUp():
  833. self.OnLeftUp(event)
  834. # dragging
  835. elif event.Dragging():
  836. self.OnDragging(event)
  837. # double click
  838. elif event.ButtonDClick():
  839. self.OnButtonDClick(event)
  840. # middle mouse button pressed
  841. elif event.MiddleDown():
  842. self.OnMiddleDown(event)
  843. # right mouse button pressed
  844. elif event.RightDown():
  845. self.OnRightDown(event)
  846. # right mouse button released
  847. elif event.RightUp():
  848. self.OnRightUp(event)
  849. elif event.Moving():
  850. self.OnMouseMoving(event)
  851. # event.Skip()
  852. def OnMouseWheel(self, event):
  853. """!
  854. Mouse wheel moved
  855. """
  856. self.processMouse = False
  857. current = event.GetPositionTuple()[:]
  858. wheel = event.GetWheelRotation()
  859. Debug.msg (5, "BufferedWindow.MouseAction(): wheel=%d" % wheel)
  860. # zoom 1/2 of the screen, centered to current mouse position (TODO: settings)
  861. begin = (current[0] - self.Map.width / 4,
  862. current[1] - self.Map.height / 4)
  863. end = (current[0] + self.Map.width / 4,
  864. current[1] + self.Map.height / 4)
  865. if wheel > 0:
  866. zoomtype = 1
  867. else:
  868. zoomtype = -1
  869. # zoom
  870. self.Zoom(begin, end, zoomtype)
  871. # redraw map
  872. self.UpdateMap()
  873. ### self.OnPaint(None)
  874. # update statusbar
  875. self.parent.StatusbarUpdate()
  876. self.Refresh()
  877. self.processMouse = True
  878. # event.Skip()
  879. def OnDragging(self, event):
  880. """!
  881. Mouse dragging with left button down
  882. """
  883. Debug.msg (5, "BufferedWindow.MouseAction(): Dragging")
  884. current = event.GetPositionTuple()[:]
  885. previous = self.mouse['begin']
  886. move = (current[0] - previous[0],
  887. current[1] - previous[1])
  888. digitToolbar = self.parent.toolbars['vdigit']
  889. # dragging or drawing box with left button
  890. if self.mouse['use'] == 'pan':
  891. self.DragMap(move)
  892. # dragging decoration overlay item
  893. elif (self.mouse['use'] == 'pointer' and
  894. not digitToolbar and
  895. self.dragid != None):
  896. self.DragItem(self.dragid, event)
  897. # dragging anything else - rubber band box or line
  898. else:
  899. if (self.mouse['use'] == 'pointer' and
  900. not digitToolbar): return
  901. self.mouse['end'] = event.GetPositionTuple()[:]
  902. digitClass = self.parent.digit
  903. if (event.LeftIsDown() and
  904. not (digitToolbar and
  905. digitToolbar.GetAction() in ("moveLine",) and
  906. digitClass.driver.GetSelected() > 0)):
  907. # draw box only when left mouse button is pressed
  908. self.MouseDraw(pdc=self.pdcTmp)
  909. # event.Skip()
  910. def OnLeftDownVDigitAddLine(self, event):
  911. """!
  912. Left mouse button down - vector digitizer add new line
  913. action
  914. """
  915. digitToolbar = self.parent.toolbars['vdigit']
  916. digitClass = self.parent.digit
  917. try:
  918. mapLayer = digitToolbar.GetLayer().GetName()
  919. except:
  920. return
  921. if digitToolbar.GetAction('type') in ["point", "centroid"]:
  922. # add new point
  923. if digitToolbar.GetAction('type') == 'point':
  924. point = True
  925. else:
  926. point = False
  927. east, north = self.Pixel2Cell(self.mouse['begin'])
  928. fid = digitClass.AddPoint(mapLayer, point, east, north)
  929. if fid < 0:
  930. return
  931. self.UpdateMap(render=False) # redraw map
  932. # add new record into atribute table
  933. if UserSettings.Get(group='vdigit', key="addRecord", subkey='enabled') is True:
  934. # select attributes based on layer and category
  935. cats = { fid : {
  936. UserSettings.Get(group='vdigit', key="layer", subkey='value') :
  937. (UserSettings.Get(group='vdigit', key="category", subkey='value'), )
  938. }}
  939. posWindow = self.ClientToScreen((self.mouse['end'][0] + self.dialogOffset,
  940. self.mouse['end'][1] + self.dialogOffset))
  941. addRecordDlg = dbm_dialogs.DisplayAttributesDialog(parent=self, map=mapLayer,
  942. cats=cats,
  943. pos=posWindow,
  944. action="add")
  945. if not point:
  946. self.__geomAttrb(fid, addRecordDlg, 'area', digitClass,
  947. digitToolbar.GetLayer())
  948. self.__geomAttrb(fid, addRecordDlg, 'perimeter', digitClass,
  949. digitToolbar.GetLayer())
  950. if addRecordDlg.mapDBInfo and \
  951. addRecordDlg.ShowModal() == wx.ID_OK:
  952. sqlfile = tempfile.NamedTemporaryFile(mode="w")
  953. for sql in addRecordDlg.GetSQLString():
  954. sqlfile.file.write(sql + ";\n")
  955. sqlfile.file.flush()
  956. gcmd.RunCommand('db.execute',
  957. parent = self,
  958. quiet = True,
  959. input = sqlfile.name)
  960. if addRecordDlg.mapDBInfo:
  961. self.__updateATM()
  962. elif digitToolbar.GetAction('type') in ["line", "boundary"]:
  963. # add new point to the line
  964. self.polycoords.append(self.Pixel2Cell(event.GetPositionTuple()[:]))
  965. self.DrawLines(pdc=self.pdcTmp)
  966. def __geomAttrb(self, fid, dialog, attrb, digit, mapLayer):
  967. """!Trac geometry attributes?"""
  968. item = self.tree.FindItemByData('maplayer', mapLayer)
  969. vdigit = self.tree.GetPyData(item)[0]['vdigit']
  970. if vdigit and \
  971. vdigit.has_key('geomAttr') and \
  972. vdigit['geomAttr'].has_key(attrb):
  973. val = -1
  974. if attrb == 'length':
  975. val = digit.GetLineLength(fid)
  976. type = attrb
  977. elif attrb == 'area':
  978. val = digit.GetAreaSize(fid)
  979. type = attrb
  980. elif attrb == 'perimeter':
  981. val = digit.GetAreaPerimeter(fid)
  982. type = 'length'
  983. if val > 0:
  984. layer = int(UserSettings.Get(group='vdigit', key="layer", subkey='value'))
  985. column = vdigit['geomAttr'][attrb]['column']
  986. val = UnitsConvertValue(val, type, vdigit['geomAttr'][attrb]['units'])
  987. dialog.SetColumnValue(layer, column, val)
  988. dialog.OnReset()
  989. def __geomAttrbUpdate(self, fids):
  990. """!Update geometry atrributes of currently selected features
  991. @param fid list feature id
  992. """
  993. mapLayer = self.parent.toolbars['vdigit'].GetLayer()
  994. vectorName = mapLayer.GetName()
  995. digit = self.parent.digit
  996. item = self.tree.FindItemByData('maplayer', mapLayer)
  997. vdigit = self.tree.GetPyData(item)[0]['vdigit']
  998. if vdigit is None or not vdigit.has_key('geomAttr'):
  999. return
  1000. dbInfo = gselect.VectorDBInfo(vectorName)
  1001. sqlfile = tempfile.NamedTemporaryFile(mode="w")
  1002. for fid in fids:
  1003. for layer, cats in digit.GetLineCats(fid).iteritems():
  1004. table = dbInfo.GetTable(layer)
  1005. for attrb, item in vdigit['geomAttr'].iteritems():
  1006. val = -1
  1007. if attrb == 'length':
  1008. val = digit.GetLineLength(fid)
  1009. type = attrb
  1010. elif attrb == 'area':
  1011. val = digit.GetAreaSize(fid)
  1012. type = attrb
  1013. elif attrb == 'perimeter':
  1014. val = digit.GetAreaPerimeter(fid)
  1015. type = 'length'
  1016. if val < 0:
  1017. continue
  1018. val = UnitsConvertValue(val, type, item['units'])
  1019. for cat in cats:
  1020. sqlfile.write('UPDATE %s SET %s = %f WHERE %s = %d;\n' % \
  1021. (table, item['column'], val,
  1022. dbInfo.GetKeyColumn(layer), cat))
  1023. sqlfile.file.flush()
  1024. gcmd.RunCommand('db.execute',
  1025. parent = True,
  1026. quiet = True,
  1027. input = sqlfile.name)
  1028. def __updateATM(self):
  1029. """!Update open Attribute Table Manager
  1030. @todo: use AddDataRow() instead
  1031. """
  1032. # update ATM
  1033. digitToolbar = self.parent.toolbars['vdigit']
  1034. digitVector = digitToolbar.GetLayer().GetName()
  1035. for atm in self.lmgr.dialogs['atm']:
  1036. atmVector = atm.GetVectorName()
  1037. if atmVector == digitVector:
  1038. layer = UserSettings.Get(group='vdigit', key="layer", subkey='value')
  1039. # TODO: use AddDataRow instead
  1040. atm.LoadData(layer)
  1041. def OnLeftDownVDigitEditLine(self, event):
  1042. """!
  1043. Left mouse button down - vector digitizer edit linear feature
  1044. - add new vertex.
  1045. """
  1046. digitToolbar = self.parent.toolbars['vdigit']
  1047. digitClass = self.parent.digit
  1048. self.polycoords.append(self.Pixel2Cell(self.mouse['begin']))
  1049. self.vdigitMove['id'].append(wx.NewId())
  1050. self.DrawLines(pdc=self.pdcTmp)
  1051. def OnLeftDownVDigitMoveLine(self, event):
  1052. """!
  1053. Left mouse button down - vector digitizer move feature/vertex,
  1054. edit linear feature
  1055. """
  1056. digitToolbar = self.parent.toolbars['vdigit']
  1057. digitClass = self.parent.digit
  1058. self.vdigitMove = {}
  1059. # geographic coordinates of initial position (left-down)
  1060. self.vdigitMove['begin'] = None
  1061. # list of ids to modify
  1062. self.vdigitMove['id'] = []
  1063. # ids geographic coordinates
  1064. self.vdigitMove['coord'] = {}
  1065. if digitToolbar.GetAction() in ["moveVertex", "editLine"]:
  1066. # set pen
  1067. pcolor = UserSettings.Get(group='vdigit', key="symbol",
  1068. subkey=["highlight", "color"])
  1069. self.pen = self.polypen = wx.Pen(colour=pcolor,
  1070. width=2, style=wx.SHORT_DASH)
  1071. self.pdcTmp.SetPen(self.polypen)
  1072. def OnLeftDownVDigitDisplayCA(self, event):
  1073. """!
  1074. Left mouse button down - vector digitizer display categories
  1075. or attributes action
  1076. """
  1077. digitToolbar = self.parent.toolbars['vdigit']
  1078. digitClass = self.parent.digit
  1079. try:
  1080. mapLayer = digitToolbar.GetLayer().GetName()
  1081. except:
  1082. return
  1083. coords = self.Pixel2Cell(self.mouse['begin'])
  1084. # unselect
  1085. digitClass.driver.SetSelected([])
  1086. # select feature by point
  1087. cats = {}
  1088. if digitClass.driver.SelectLineByPoint(coords,
  1089. digitClass.GetSelectType()) is None:
  1090. return
  1091. if UserSettings.Get(group='vdigit', key='checkForDupl',
  1092. subkey='enabled'):
  1093. lines = digitClass.driver.GetSelected()
  1094. else:
  1095. lines = (digitClass.driver.GetSelected()[0],) # only first found
  1096. for line in lines:
  1097. cats[line] = digitClass.GetLineCats(line)
  1098. posWindow = self.ClientToScreen((self.mouse['end'][0] + self.dialogOffset,
  1099. self.mouse['end'][1] + self.dialogOffset))
  1100. if digitToolbar.GetAction() == "displayAttrs":
  1101. # select attributes based on coordinates (all layers)
  1102. if self.parent.dialogs['attributes'] is None:
  1103. self.parent.dialogs['attributes'] = \
  1104. dbm_dialogs.DisplayAttributesDialog(parent=self, map=mapLayer,
  1105. cats=cats,
  1106. action="update")
  1107. else:
  1108. # upgrade dialog
  1109. self.parent.dialogs['attributes'].UpdateDialog(cats=cats)
  1110. if self.parent.dialogs['attributes']:
  1111. if len(cats.keys()) > 0:
  1112. # highlight feature & re-draw map
  1113. if not self.parent.dialogs['attributes'].IsShown():
  1114. self.parent.dialogs['attributes'].Show()
  1115. else:
  1116. if self.parent.dialogs['attributes'] and \
  1117. self.parent.dialogs['attributes'].IsShown():
  1118. self.parent.dialogs['attributes'].Hide()
  1119. else: # displayCats
  1120. if self.parent.dialogs['category'] is None:
  1121. # open new dialog
  1122. dlg = VDigitCategoryDialog(parent=self,
  1123. map=mapLayer,
  1124. cats=cats,
  1125. pos=posWindow,
  1126. title=_("Update categories"))
  1127. self.parent.dialogs['category'] = dlg
  1128. else:
  1129. # update currently open dialog
  1130. self.parent.dialogs['category'].UpdateDialog(cats=cats)
  1131. if self.parent.dialogs['category']:
  1132. if len(cats.keys()) > 0:
  1133. # highlight feature & re-draw map
  1134. if not self.parent.dialogs['category'].IsShown():
  1135. self.parent.dialogs['category'].Show()
  1136. else:
  1137. if self.parent.dialogs['category'].IsShown():
  1138. self.parent.dialogs['category'].Hide()
  1139. self.UpdateMap(render=False)
  1140. def OnLeftDownVDigitCopyCA(self, event):
  1141. """!
  1142. Left mouse button down - vector digitizer copy categories
  1143. or attributes action
  1144. """
  1145. digitToolbar = self.parent.toolbars['vdigit']
  1146. digitClass = self.parent.digit
  1147. if not hasattr(self, "copyCatsList"):
  1148. self.copyCatsList = []
  1149. else:
  1150. self.copyCatsIds = []
  1151. self.mouse['box'] = 'box'
  1152. def OnLeftDownVDigitCopyLine(self, event):
  1153. """!
  1154. Left mouse button down - vector digitizer copy lines action
  1155. """
  1156. digitToolbar = self.parent.toolbars['vdigit']
  1157. digitClass = self.parent.digit
  1158. if not hasattr(self, "copyIds"):
  1159. self.copyIds = []
  1160. self.layerTmp = None
  1161. def OnLeftDownVDigitBulkLine(self, event):
  1162. """!
  1163. Left mouse button down - vector digitizer label 3d vector
  1164. lines
  1165. """
  1166. digitToolbar = self.parent.toolbars['vdigit']
  1167. digitClass = self.parent.digit
  1168. if len(self.polycoords) > 1: # start new line
  1169. self.polycoords = []
  1170. self.ClearLines(pdc=self.pdcTmp)
  1171. self.polycoords.append(self.Pixel2Cell(event.GetPositionTuple()[:]))
  1172. if len(self.polycoords) == 1:
  1173. begin = self.Pixel2Cell(self.polycoords[-1])
  1174. end = self.Pixel2Cell(self.mouse['end'])
  1175. else:
  1176. end = self.Pixel2Cell(self.polycoords[-1])
  1177. begin = self.Pixel2Cell(self.mouse['begin'])
  1178. self.DrawLines(self.pdcTmp, begin, end)
  1179. def OnLeftDown(self, event):
  1180. """!
  1181. Left mouse button pressed
  1182. """
  1183. Debug.msg (5, "BufferedWindow.OnLeftDown(): use=%s" % \
  1184. self.mouse["use"])
  1185. self.mouse['begin'] = event.GetPositionTuple()[:]
  1186. if self.mouse["use"] in ["measure", "profile"]:
  1187. # measure or profile
  1188. if len(self.polycoords) == 0:
  1189. self.mouse['end'] = self.mouse['begin']
  1190. self.polycoords.append(self.Pixel2Cell(self.mouse['begin']))
  1191. self.ClearLines(pdc=self.pdcTmp)
  1192. else:
  1193. self.mouse['begin'] = self.mouse['end']
  1194. elif self.mouse['use'] == 'zoom':
  1195. pass
  1196. #
  1197. # vector digizer
  1198. #
  1199. elif self.mouse["use"] == "pointer" and \
  1200. self.parent.toolbars['vdigit']:
  1201. digitToolbar = self.parent.toolbars['vdigit']
  1202. digitClass = self.parent.digit
  1203. try:
  1204. mapLayer = digitToolbar.GetLayer().GetName()
  1205. except:
  1206. wx.MessageBox(parent=self,
  1207. message=_("No vector map selected for editing."),
  1208. caption=_("Vector digitizer"),
  1209. style=wx.OK | wx.ICON_INFORMATION | wx.CENTRE)
  1210. event.Skip()
  1211. return
  1212. if digitToolbar.GetAction() not in ("moveVertex",
  1213. "addVertex",
  1214. "removeVertex",
  1215. "editLine"):
  1216. # set pen
  1217. self.pen = wx.Pen(colour='Red', width=2, style=wx.SHORT_DASH)
  1218. self.polypen = wx.Pen(colour='dark green', width=2, style=wx.SOLID)
  1219. if digitToolbar.GetAction() in ("addVertex",
  1220. "removeVertex",
  1221. "splitLines"):
  1222. # unselect
  1223. digitClass.driver.SetSelected([])
  1224. if digitToolbar.GetAction() == "addLine":
  1225. self.OnLeftDownVDigitAddLine(event)
  1226. elif digitToolbar.GetAction() == "editLine" and \
  1227. hasattr(self, "vdigitMove"):
  1228. self.OnLeftDownVDigitEditLine(event)
  1229. elif digitToolbar.GetAction() in ("moveLine",
  1230. "moveVertex",
  1231. "editLine") and \
  1232. not hasattr(self, "vdigitMove"):
  1233. self.OnLeftDownVDigitMoveLine(event)
  1234. elif digitToolbar.GetAction() in ("displayAttrs"
  1235. "displayCats"):
  1236. self.OnLeftDownVDigitDisplayCA(event)
  1237. elif digitToolbar.GetAction() in ("copyCats",
  1238. "copyAttrs"):
  1239. self.OnLeftDownVDigitCopyCA(event)
  1240. elif digitToolbar.GetAction() == "copyLine":
  1241. self.OnLeftDownVDigitCopyLine(event)
  1242. elif digitToolbar.GetAction() == "zbulkLine":
  1243. self.OnLeftDownVDigitBulkLine(event)
  1244. elif self.mouse['use'] == 'pointer':
  1245. # get decoration or text id
  1246. self.idlist = []
  1247. self.dragid = ''
  1248. self.lastpos = self.mouse['begin']
  1249. idlist = self.pdc.FindObjects(self.lastpos[0], self.lastpos[1],
  1250. self.hitradius)
  1251. if 99 in idlist:
  1252. idlist.remove(99)
  1253. if idlist != []:
  1254. self.dragid = idlist[0] #drag whatever is on top
  1255. else:
  1256. pass
  1257. event.Skip()
  1258. def OnLeftUpVDigitVarious(self, event):
  1259. """!
  1260. Left mouse button up - vector digitizer various actions
  1261. """
  1262. digitToolbar = self.parent.toolbars['vdigit']
  1263. digitClass = self.parent.digit
  1264. pos1 = self.Pixel2Cell(self.mouse['begin'])
  1265. pos2 = self.Pixel2Cell(self.mouse['end'])
  1266. nselected = 0
  1267. # -> delete line || move line || move vertex
  1268. if digitToolbar.GetAction() in ("moveVertex",
  1269. "editLine"):
  1270. if len(digitClass.driver.GetSelected()) == 0:
  1271. nselected = digitClass.driver.SelectLineByPoint(pos1, type=VDigit_Lines_Type)
  1272. if digitToolbar.GetAction() == "editLine":
  1273. try:
  1274. selVertex = digitClass.driver.GetSelectedVertex(pos1)[0]
  1275. except IndexError:
  1276. selVertex = None
  1277. if selVertex:
  1278. # self.UpdateMap(render=False)
  1279. ids = digitClass.driver.GetSelected(grassId=False)
  1280. # move this line to tmp layer
  1281. self.polycoords = []
  1282. for id in ids:
  1283. if id % 2: # register only vertices
  1284. e, n = self.Pixel2Cell(self.pdcVector.GetIdBounds(id)[0:2])
  1285. self.polycoords.append((e, n))
  1286. digitClass.driver.DrawSelected(False)
  1287. if selVertex < ids[-1] / 2:
  1288. # choose first or last node of line
  1289. self.vdigitMove['id'].reverse()
  1290. self.polycoords.reverse()
  1291. else:
  1292. # unselect
  1293. digitClass.driver.SetSelected([])
  1294. del self.vdigitMove
  1295. self.UpdateMap(render=False)
  1296. elif digitToolbar.GetAction() in ("copyCats",
  1297. "copyAttrs"):
  1298. if not hasattr(self, "copyCatsIds"):
  1299. # 'from' -> select by point
  1300. nselected = digitClass.driver.SelectLineByPoint(pos1, digitClass.GetSelectType())
  1301. if nselected:
  1302. self.copyCatsList = digitClass.driver.GetSelected()
  1303. else:
  1304. # -> 'to' -> select by bbox
  1305. digitClass.driver.SetSelected([])
  1306. # return number of selected features (by box/point)
  1307. nselected = digitClass.driver.SelectLinesByBox(pos1, pos2,
  1308. digitClass.GetSelectType())
  1309. if nselected == 0:
  1310. if digitClass.driver.SelectLineByPoint(pos1,
  1311. digitClass.GetSelectType()) is not None:
  1312. nselected = 1
  1313. if nselected > 0:
  1314. self.copyCatsIds = digitClass.driver.GetSelected()
  1315. elif digitToolbar.GetAction() == "queryLine":
  1316. selected = digitClass.SelectLinesByQuery(pos1, pos2)
  1317. nselected = len(selected)
  1318. if nselected > 0:
  1319. digitClass.driver.SetSelected(selected)
  1320. else:
  1321. # -> moveLine || deleteLine, etc. (select by point/box)
  1322. if digitToolbar.GetAction() == 'moveLine' and \
  1323. len(digitClass.driver.GetSelected()) > 0:
  1324. nselected = 0
  1325. else:
  1326. if digitToolbar.GetAction() == 'moveLine':
  1327. drawSeg = True
  1328. else:
  1329. drawSeg = False
  1330. nselected = digitClass.driver.SelectLinesByBox(pos1, pos2,
  1331. digitClass.GetSelectType(),
  1332. drawSeg)
  1333. if nselected == 0:
  1334. if digitClass.driver.SelectLineByPoint(pos1,
  1335. digitClass.GetSelectType()) is not None:
  1336. nselected = 1
  1337. if nselected > 0:
  1338. if digitToolbar.GetAction() in ("moveLine",
  1339. "moveVertex"):
  1340. # get pseudoDC id of objects which should be redrawn
  1341. if digitToolbar.GetAction() == "moveLine":
  1342. # -> move line
  1343. self.vdigitMove['id'] = digitClass.driver.GetSelected(grassId=False)
  1344. self.vdigitMove['coord'] = digitClass.driver.GetSelectedCoord()
  1345. else: # moveVertex
  1346. self.vdigitMove['id'] = digitClass.driver.GetSelectedVertex(pos1)
  1347. if len(self.vdigitMove['id']) == 0: # no vertex found
  1348. digitClass.driver.SetSelected([])
  1349. #
  1350. # check for duplicates
  1351. #
  1352. if UserSettings.Get(group='vdigit', key='checkForDupl', subkey='enabled') is True:
  1353. dupl = digitClass.driver.GetDuplicates()
  1354. self.UpdateMap(render=False)
  1355. if dupl:
  1356. posWindow = self.ClientToScreen((self.mouse['end'][0] + self.dialogOffset,
  1357. self.mouse['end'][1] + self.dialogOffset))
  1358. dlg = VDigitDuplicatesDialog(parent=self, data=dupl, pos=posWindow)
  1359. if dlg.ShowModal() == wx.ID_OK:
  1360. digitClass.driver.UnSelect(dlg.GetUnSelected())
  1361. # update selected
  1362. self.UpdateMap(render=False)
  1363. if digitToolbar.GetAction() != "editLine":
  1364. # -> move line || move vertex
  1365. self.UpdateMap(render=False)
  1366. else: # no vector object found
  1367. if not (digitToolbar.GetAction() in ("moveLine",
  1368. "moveVertex") and \
  1369. len(self.vdigitMove['id']) > 0):
  1370. # avoid left-click when features are already selected
  1371. self.UpdateMap(render=False, renderVector=False)
  1372. def OnLeftUpVDigitModifyLine(self, event):
  1373. """!
  1374. Left mouse button up - vector digitizer split line, add/remove
  1375. vertex action
  1376. """
  1377. digitToolbar = self.parent.toolbars['vdigit']
  1378. digitClass = self.parent.digit
  1379. pos1 = self.Pixel2Cell(self.mouse['begin'])
  1380. pointOnLine = digitClass.driver.SelectLineByPoint(pos1,
  1381. type=VDigit_Lines_Type)
  1382. if not pointOnLine:
  1383. return
  1384. if digitToolbar.GetAction() in ["splitLine", "addVertex"]:
  1385. self.UpdateMap(render=False) # highlight object
  1386. self.DrawCross(pdc=self.pdcTmp, coords=self.Cell2Pixel(pointOnLine),
  1387. size=5)
  1388. else: # removeVertex
  1389. # get only id of vertex
  1390. try:
  1391. id = digitClass.driver.GetSelectedVertex(pos1)[0]
  1392. except IndexError:
  1393. id = None
  1394. if id:
  1395. x, y = self.pdcVector.GetIdBounds(id)[0:2]
  1396. self.pdcVector.RemoveId(id)
  1397. self.UpdateMap(render=False) # highlight object
  1398. self.DrawCross(pdc=self.pdcTmp, coords=(x, y),
  1399. size=5)
  1400. else:
  1401. # unselect
  1402. digitClass.driver.SetSelected([])
  1403. self.UpdateMap(render=False)
  1404. def OnLeftUpVDigitCopyLine(self, event):
  1405. """!
  1406. Left mouse button up - vector digitizer copy feature action
  1407. """
  1408. digitToolbar = self.parent.toolbars['vdigit']
  1409. digitClass = self.parent.digit
  1410. pos1 = self.Pixel2Cell(self.mouse['begin'])
  1411. pos2 = self.Pixel2Cell(self.mouse['end'])
  1412. if UserSettings.Get(group='vdigit', key='bgmap',
  1413. subkey='value', internal=True) == '':
  1414. # no background map -> copy from current vector map layer
  1415. nselected = digitClass.driver.SelectLinesByBox(pos1, pos2,
  1416. digitClass.GetSelectType())
  1417. if nselected > 0:
  1418. # highlight selected features
  1419. self.UpdateMap(render=False)
  1420. else:
  1421. self.UpdateMap(render=False, renderVector=False)
  1422. else:
  1423. # copy features from background map
  1424. self.copyIds += digitClass.SelectLinesFromBackgroundMap(pos1, pos2)
  1425. if len(self.copyIds) > 0:
  1426. color = UserSettings.Get(group='vdigit', key='symbol',
  1427. subkey=['highlight', 'color'])
  1428. colorStr = str(color[0]) + ":" + \
  1429. str(color[1]) + ":" + \
  1430. str(color[2])
  1431. dVectTmp = ['d.vect',
  1432. 'map=%s' % UserSettings.Get(group='vdigit', key='bgmap',
  1433. subkey='value', internal=True),
  1434. 'cats=%s' % utils.ListOfCatsToRange(self.copyIds),
  1435. '-i',
  1436. 'color=%s' % colorStr,
  1437. 'fcolor=%s' % colorStr,
  1438. 'type=point,line,boundary,centroid',
  1439. 'width=2']
  1440. if not self.layerTmp:
  1441. self.layerTmp = self.Map.AddLayer(type='vector',
  1442. name=globalvar.QUERYLAYER,
  1443. command=dVectTmp)
  1444. else:
  1445. self.layerTmp.SetCmd(dVectTmp)
  1446. self.UpdateMap(render=True, renderVector=False)
  1447. else:
  1448. self.UpdateMap(render=False, renderVector=False)
  1449. self.redrawAll = None
  1450. def OnLeftUpVDigitBulkLine(self, event):
  1451. """!
  1452. Left mouse button up - vector digitizer z-bulk line action
  1453. """
  1454. digitToolbar = self.parent.toolbars['vdigit']
  1455. digitClass = self.parent.digit
  1456. # select lines to be labeled
  1457. pos1 = self.polycoords[0]
  1458. pos2 = self.polycoords[1]
  1459. nselected = digitClass.driver.SelectLinesByBox(pos1, pos2,
  1460. digitClass.GetSelectType())
  1461. if nselected > 0:
  1462. # highlight selected features
  1463. self.UpdateMap(render=False)
  1464. self.DrawLines(pdc=self.pdcTmp) # redraw temp line
  1465. else:
  1466. self.UpdateMap(render=False, renderVector=False)
  1467. def OnLeftUpVDigitConnectLine(self, event):
  1468. """!
  1469. Left mouse button up - vector digitizer connect line action
  1470. """
  1471. digitToolbar = self.parent.toolbars['vdigit']
  1472. digitClass = self.parent.digit
  1473. if len(digitClass.driver.GetSelected()) > 0:
  1474. self.UpdateMap(render=False)
  1475. def OnLeftUp(self, event):
  1476. """!
  1477. Left mouse button released
  1478. """
  1479. Debug.msg (5, "BufferedWindow.OnLeftUp(): use=%s" % \
  1480. self.mouse["use"])
  1481. self.mouse['end'] = event.GetPositionTuple()[:]
  1482. if self.mouse['use'] in ["zoom", "pan"]:
  1483. # set region in zoom or pan
  1484. begin = self.mouse['begin']
  1485. end = self.mouse['end']
  1486. if self.mouse['use'] == 'zoom':
  1487. # set region for click (zero-width box)
  1488. if begin[0] - end[0] == 0 or \
  1489. begin[1] - end[1] == 0:
  1490. # zoom 1/2 of the screen (TODO: settings)
  1491. begin = (end[0] - self.Map.width / 4,
  1492. end[1] - self.Map.height / 4)
  1493. end = (end[0] + self.Map.width / 4,
  1494. end[1] + self.Map.height / 4)
  1495. self.Zoom(begin, end, self.zoomtype)
  1496. # redraw map
  1497. self.UpdateMap(render=True)
  1498. # update statusbar
  1499. self.parent.StatusbarUpdate()
  1500. elif self.mouse["use"] == "query":
  1501. # querying
  1502. self.parent.QueryMap(self.mouse['begin'][0],self.mouse['begin'][1])
  1503. elif self.mouse["use"] == "queryVector":
  1504. # editable mode for vector map layers
  1505. self.parent.QueryVector(self.mouse['begin'][0], self.mouse['begin'][1])
  1506. # clear temp canvas
  1507. self.UpdateMap(render=False, renderVector=False)
  1508. elif self.mouse["use"] in ["measure", "profile"]:
  1509. # measure or profile
  1510. if self.mouse["use"] == "measure":
  1511. self.parent.MeasureDist(self.mouse['begin'], self.mouse['end'])
  1512. self.polycoords.append(self.Pixel2Cell(self.mouse['end']))
  1513. self.ClearLines(pdc=self.pdcTmp)
  1514. self.DrawLines(pdc=self.pdcTmp)
  1515. elif self.mouse["use"] == "pointer" and \
  1516. self.parent.GetLayerManager().georectifying:
  1517. # -> georectifying
  1518. coord = self.Pixel2Cell(self.mouse['end'])
  1519. if self.parent.toolbars['georect']:
  1520. coordtype = 'gcpcoord'
  1521. else:
  1522. coordtype = 'mapcoord'
  1523. self.parent.GetLayerManager().georectifying.SetGCPData(coordtype, coord, self)
  1524. self.UpdateMap(render=False, renderVector=False)
  1525. elif self.mouse["use"] == "pointer" and self.parent.toolbars['vdigit']:
  1526. # digitization tool active
  1527. digitToolbar = self.parent.toolbars['vdigit']
  1528. digitClass = self.parent.digit
  1529. if hasattr(self, "vdigitMove"):
  1530. if len(digitClass.driver.GetSelected()) == 0:
  1531. self.vdigitMove['begin'] = self.Pixel2Cell(self.mouse['begin']) # left down
  1532. # eliminate initial mouse moving efect
  1533. self.mouse['begin'] = self.mouse['end']
  1534. if digitToolbar.GetAction() in ("deleteLine",
  1535. "moveLine",
  1536. "moveVertex",
  1537. "copyCats",
  1538. "copyAttrs",
  1539. "editLine",
  1540. "flipLine",
  1541. "mergeLine",
  1542. "snapLine",
  1543. "queryLine",
  1544. "breakLine",
  1545. "typeConv",
  1546. "connectLine"):
  1547. self.OnLeftUpVDigitVarious(event)
  1548. elif digitToolbar.GetAction() in ("splitLine",
  1549. "addVertex",
  1550. "removeVertex"):
  1551. self.OnLeftUpVDigitModifyLine(event)
  1552. elif digitToolbar.GetAction() == "copyLine":
  1553. self.OnLeftUpVDigitCopyLine(event)
  1554. elif digitToolbar.GetAction() == "zbulkLine" and \
  1555. len(self.polycoords) == 2:
  1556. self.OnLeftUpVDigitBulkLine(event)
  1557. elif digitToolbar.GetAction() == "connectLine":
  1558. self.OnLeftUpConnectLine(event)
  1559. if len(digitClass.driver.GetSelected()) > 0:
  1560. self.redrawAll = None
  1561. elif (self.mouse['use'] == 'pointer' and
  1562. self.dragid >= 0):
  1563. # end drag of overlay decoration
  1564. if self.dragid < 99 and self.overlays.has_key(self.dragid):
  1565. self.overlays[self.dragid]['coords'] = self.pdc.GetIdBounds(self.dragid)
  1566. elif self.dragid > 100 and self.textdict.has_key(self.dragid):
  1567. self.textdict[self.dragid]['coords'] = self.pdc.GetIdBounds(self.dragid)
  1568. else:
  1569. pass
  1570. self.dragid = None
  1571. self.currtxtid = None
  1572. def OnButtonDClick(self, event):
  1573. """!
  1574. Mouse button double click
  1575. """
  1576. Debug.msg (5, "BufferedWindow.OnButtonDClick(): use=%s" % \
  1577. self.mouse["use"])
  1578. if self.mouse["use"] == "measure":
  1579. # measure
  1580. self.ClearLines(pdc=self.pdcTmp)
  1581. self.polycoords = []
  1582. self.mouse['use'] = 'pointer'
  1583. self.mouse['box'] = 'point'
  1584. self.mouse['end'] = [0, 0]
  1585. self.Refresh()
  1586. self.SetCursor(self.parent.cursors["default"])
  1587. elif self.mouse["use"] == "profile":
  1588. pass
  1589. elif self.mouse['use'] == 'pointer' and \
  1590. self.parent.toolbars['vdigit']:
  1591. # vector digitizer
  1592. pass
  1593. else:
  1594. # select overlay decoration options dialog
  1595. clickposition = event.GetPositionTuple()[:]
  1596. idlist = self.pdc.FindObjects(clickposition[0], clickposition[1], self.hitradius)
  1597. if idlist == []:
  1598. return
  1599. self.dragid = idlist[0]
  1600. # self.ovlcoords[self.dragid] = self.pdc.GetIdBounds(self.dragid)
  1601. if self.dragid > 100:
  1602. self.currtxtid = self.dragid
  1603. self.parent.OnAddText(None)
  1604. elif self.dragid == 0:
  1605. self.parent.OnAddBarscale(None)
  1606. elif self.dragid == 1:
  1607. self.parent.OnAddLegend(None)
  1608. def OnRightDown(self, event):
  1609. """!
  1610. Right mouse button pressed
  1611. """
  1612. Debug.msg (5, "BufferedWindow.OnRightDown(): use=%s" % \
  1613. self.mouse["use"])
  1614. digitToolbar = self.parent.toolbars['vdigit']
  1615. if digitToolbar:
  1616. digitClass = self.parent.digit
  1617. # digitization tool (confirm action)
  1618. if digitToolbar.GetAction() in ("moveLine",
  1619. "moveVertex") and \
  1620. hasattr(self, "vdigitMove"):
  1621. pFrom = self.vdigitMove['begin']
  1622. pTo = self.Pixel2Cell(event.GetPositionTuple())
  1623. move = (pTo[0] - pFrom[0],
  1624. pTo[1] - pFrom[1])
  1625. if digitToolbar.GetAction() == "moveLine":
  1626. # move line
  1627. if digitClass.MoveSelectedLines(move) < 0:
  1628. return
  1629. elif digitToolbar.GetAction() == "moveVertex":
  1630. # move vertex
  1631. fid = digitClass.MoveSelectedVertex(pFrom, move)
  1632. if fid < 0:
  1633. return
  1634. self.__geomAttrbUpdate([fid,])
  1635. del self.vdigitMove
  1636. event.Skip()
  1637. def OnRightUp(self, event):
  1638. """!
  1639. Right mouse button released
  1640. """
  1641. Debug.msg (5, "BufferedWindow.OnRightUp(): use=%s" % \
  1642. self.mouse["use"])
  1643. digitToolbar = self.parent.toolbars['vdigit']
  1644. if digitToolbar:
  1645. digitClass = self.parent.digit
  1646. # digitization tool (confirm action)
  1647. if digitToolbar.GetAction() == "addLine" and \
  1648. digitToolbar.GetAction('type') in ["line", "boundary"]:
  1649. # -> add new line / boundary
  1650. try:
  1651. map = digitToolbar.GetLayer().GetName()
  1652. except:
  1653. map = None
  1654. wx.MessageBox(parent=self,
  1655. message=_("No vector map selected for editing."),
  1656. caption=_("Error"), style=wx.OK | wx.ICON_ERROR | wx.CENTRE)
  1657. if map:
  1658. # mapcoords = []
  1659. # xy -> EN
  1660. # for coord in self.polycoords:
  1661. # mapcoords.append(self.Pixel2Cell(coord))
  1662. if digitToolbar.GetAction('type') == 'line':
  1663. line = True
  1664. else:
  1665. line = False
  1666. if len(self.polycoords) < 2: # ignore 'one-point' lines
  1667. return
  1668. fid = digitClass.AddLine(map, line, self.polycoords)
  1669. if fid < 0:
  1670. return
  1671. position = self.Cell2Pixel(self.polycoords[-1])
  1672. self.polycoords = []
  1673. self.UpdateMap(render=False)
  1674. self.redrawAll = True
  1675. # add new record into atribute table
  1676. if UserSettings.Get(group='vdigit', key="addRecord", subkey='enabled') and \
  1677. (line is True or \
  1678. (not line and fid > 0)):
  1679. posWindow = self.ClientToScreen((position[0] + self.dialogOffset,
  1680. position[1] + self.dialogOffset))
  1681. # select attributes based on layer and category
  1682. cats = { fid : {
  1683. UserSettings.Get(group='vdigit', key="layer", subkey='value') :
  1684. (UserSettings.Get(group='vdigit', key="category", subkey='value'), )
  1685. }}
  1686. addRecordDlg = dbm_dialogs.DisplayAttributesDialog(parent=self, map=map,
  1687. cats=cats,
  1688. pos=posWindow,
  1689. action="add")
  1690. self.__geomAttrb(fid, addRecordDlg, 'length', digitClass,
  1691. digitToolbar.GetLayer())
  1692. # auto-placing centroid
  1693. self.__geomAttrb(fid, addRecordDlg, 'area', digitClass,
  1694. digitToolbar.GetLayer())
  1695. self.__geomAttrb(fid, addRecordDlg, 'perimeter', digitClass,
  1696. digitToolbar.GetLayer())
  1697. if addRecordDlg.mapDBInfo and \
  1698. addRecordDlg.ShowModal() == wx.ID_OK:
  1699. sqlfile = tempfile.NamedTemporaryFile(mode="w")
  1700. for sql in addRecordDlg.GetSQLString():
  1701. sqlfile.file.write(sql + ";\n")
  1702. sqlfile.file.flush()
  1703. gcmd.RunCommand('db.execute',
  1704. parent = True,
  1705. quiet = True,
  1706. input = sqlfile.name)
  1707. if addRecordDlg.mapDBInfo:
  1708. self.__updateATM()
  1709. elif digitToolbar.GetAction() == "deleteLine":
  1710. # -> delete selected vector features
  1711. if digitClass.DeleteSelectedLines() < 0:
  1712. return
  1713. self.__updateATM()
  1714. elif digitToolbar.GetAction() == "splitLine":
  1715. # split line
  1716. if digitClass.SplitLine(self.Pixel2Cell(self.mouse['begin'])) < 0:
  1717. return
  1718. elif digitToolbar.GetAction() == "addVertex":
  1719. # add vertex
  1720. fid = digitClass.AddVertex(self.Pixel2Cell(self.mouse['begin']))
  1721. if fid < 0:
  1722. return
  1723. elif digitToolbar.GetAction() == "removeVertex":
  1724. # remove vertex
  1725. fid = digitClass.RemoveVertex(self.Pixel2Cell(self.mouse['begin']))
  1726. if fid < 0:
  1727. return
  1728. self.__geomAttrbUpdate([fid,])
  1729. elif digitToolbar.GetAction() in ("copyCats", "copyAttrs"):
  1730. try:
  1731. if digitToolbar.GetAction() == 'copyCats':
  1732. if digitClass.CopyCats(self.copyCatsList,
  1733. self.copyCatsIds, copyAttrb=False) < 0:
  1734. return
  1735. else:
  1736. if digitClass.CopyCats(self.copyCatsList,
  1737. self.copyCatsIds, copyAttrb=True) < 0:
  1738. return
  1739. del self.copyCatsList
  1740. del self.copyCatsIds
  1741. except AttributeError:
  1742. pass
  1743. self.__updateATM()
  1744. elif digitToolbar.GetAction() == "editLine" and \
  1745. hasattr(self, "vdigitMove"):
  1746. line = digitClass.driver.GetSelected()
  1747. if digitClass.EditLine(line, self.polycoords) < 0:
  1748. return
  1749. del self.vdigitMove
  1750. elif digitToolbar.GetAction() == "flipLine":
  1751. if digitClass.FlipLine() < 0:
  1752. return
  1753. elif digitToolbar.GetAction() == "mergeLine":
  1754. if digitClass.MergeLine() < 0:
  1755. return
  1756. elif digitToolbar.GetAction() == "breakLine":
  1757. if digitClass.BreakLine() < 0:
  1758. return
  1759. elif digitToolbar.GetAction() == "snapLine":
  1760. if digitClass.SnapLine() < 0:
  1761. return
  1762. elif digitToolbar.GetAction() == "connectLine":
  1763. if len(digitClass.driver.GetSelected()) > 1:
  1764. if digitClass.ConnectLine() < 0:
  1765. return
  1766. elif digitToolbar.GetAction() == "copyLine":
  1767. if digitClass.CopyLine(self.copyIds) < 0:
  1768. return
  1769. del self.copyIds
  1770. if self.layerTmp:
  1771. self.Map.DeleteLayer(self.layerTmp)
  1772. self.UpdateMap(render=True, renderVector=False)
  1773. del self.layerTmp
  1774. elif digitToolbar.GetAction() == "zbulkLine" and len(self.polycoords) == 2:
  1775. pos1 = self.polycoords[0]
  1776. pos2 = self.polycoords[1]
  1777. selected = digitClass.driver.GetSelected()
  1778. dlg = VDigitZBulkDialog(parent=self, title=_("Z bulk-labeling dialog"),
  1779. nselected=len(selected))
  1780. if dlg.ShowModal() == wx.ID_OK:
  1781. if digitClass.ZBulkLine(pos1, pos2, dlg.value.GetValue(),
  1782. dlg.step.GetValue()) < 0:
  1783. return
  1784. self.UpdateMap(render=False, renderVector=True)
  1785. elif digitToolbar.GetAction() == "typeConv":
  1786. # -> feature type conversion
  1787. # - point <-> centroid
  1788. # - line <-> boundary
  1789. if digitClass.TypeConvForSelectedLines() < 0:
  1790. return
  1791. if digitToolbar.GetAction() != "addLine":
  1792. # unselect and re-render
  1793. digitClass.driver.SetSelected([])
  1794. self.polycoords = []
  1795. self.UpdateMap(render=False)
  1796. self.redrawAll = True
  1797. event.Skip()
  1798. def OnMiddleDown(self, event):
  1799. """!
  1800. Middle mouse button pressed
  1801. """
  1802. digitToolbar = self.parent.toolbars['vdigit']
  1803. # digitization tool
  1804. if self.mouse["use"] == "pointer" and digitToolbar:
  1805. digitClass = self.parent.digit
  1806. if (digitToolbar.GetAction() == "addLine" and \
  1807. digitToolbar.GetAction('type') in ["line", "boundary"]) or \
  1808. digitToolbar.GetAction() == "editLine":
  1809. # add line or boundary -> remove last point from the line
  1810. try:
  1811. removed = self.polycoords.pop()
  1812. Debug.msg(4, "BufferedWindow.OnMiddleDown(): polycoords_poped=%s" % \
  1813. [removed,])
  1814. self.mouse['begin'] = self.Cell2Pixel(self.polycoords[-1])
  1815. except:
  1816. pass
  1817. if digitToolbar.GetAction() == "editLine":
  1818. # remove last vertex & line
  1819. if len(self.vdigitMove['id']) > 1:
  1820. self.vdigitMove['id'].pop()
  1821. self.UpdateMap(render=False, renderVector=False)
  1822. elif digitToolbar.GetAction() in ["deleteLine", "moveLine", "splitLine",
  1823. "addVertex", "removeVertex", "moveVertex",
  1824. "copyCats", "flipLine", "mergeLine",
  1825. "snapLine", "connectLine", "copyLine",
  1826. "queryLine", "breakLine", "typeConv"]:
  1827. # varios tools -> unselected selected features
  1828. digitClass.driver.SetSelected([])
  1829. if digitToolbar.GetAction() in ["moveLine", "moveVertex", "editLine"] and \
  1830. hasattr(self, "vdigitMove"):
  1831. del self.vdigitMove
  1832. elif digitToolbar.GetAction() == "copyCats":
  1833. try:
  1834. del self.copyCatsList
  1835. del self.copyCatsIds
  1836. except AttributeError:
  1837. pass
  1838. elif digitToolbar.GetAction() == "copyLine":
  1839. del self.copyIds
  1840. if self.layerTmp:
  1841. self.Map.DeleteLayer(self.layerTmp)
  1842. self.UpdateMap(render=True, renderVector=False)
  1843. del self.layerTmp
  1844. self.polycoords = []
  1845. self.UpdateMap(render=False) # render vector
  1846. elif digitToolbar.GetAction() == "zbulkLine":
  1847. # reset polyline
  1848. self.polycoords = []
  1849. digitClass.driver.SetSelected([])
  1850. self.UpdateMap(render=False)
  1851. self.redrawAll = True
  1852. def OnMouseMoving(self, event):
  1853. """!
  1854. Motion event and no mouse buttons were pressed
  1855. """
  1856. digitToolbar = self.parent.toolbars['vdigit']
  1857. if self.mouse["use"] == "pointer" and digitToolbar:
  1858. digitClass = self.parent.digit
  1859. self.mouse['end'] = event.GetPositionTuple()[:]
  1860. Debug.msg (5, "BufferedWindow.OnMouseMoving(): coords=%f,%f" % \
  1861. (self.mouse['end'][0], self.mouse['end'][1]))
  1862. if digitToolbar.GetAction() == "addLine" and digitToolbar.GetAction('type') in ["line", "boundary"]:
  1863. if len(self.polycoords) > 0:
  1864. self.MouseDraw(pdc=self.pdcTmp, begin=self.Cell2Pixel(self.polycoords[-1]))
  1865. elif digitToolbar.GetAction() in ["moveLine", "moveVertex", "editLine"] \
  1866. and hasattr(self, "vdigitMove"):
  1867. dx = self.mouse['end'][0] - self.mouse['begin'][0]
  1868. dy = self.mouse['end'][1] - self.mouse['begin'][1]
  1869. if len(self.vdigitMove['id']) > 0:
  1870. # draw lines on new position
  1871. if digitToolbar.GetAction() == "moveLine":
  1872. # move line
  1873. for id in self.vdigitMove['id']:
  1874. self.pdcTmp.TranslateId(id, dx, dy)
  1875. elif digitToolbar.GetAction() in ["moveVertex", "editLine"]:
  1876. # move vertex ->
  1877. # (vertex, left vertex, left line,
  1878. # right vertex, right line)
  1879. # do not draw static lines
  1880. if digitToolbar.GetAction() == "moveVertex":
  1881. self.polycoords = []
  1882. ### self.pdcTmp.TranslateId(self.vdigitMove['id'][0], dx, dy)
  1883. self.pdcTmp.RemoveId(self.vdigitMove['id'][0])
  1884. if self.vdigitMove['id'][1] > 0: # previous vertex
  1885. x, y = self.Pixel2Cell(self.pdcTmp.GetIdBounds(self.vdigitMove['id'][1])[0:2])
  1886. self.pdcTmp.RemoveId(self.vdigitMove['id'][1]+1)
  1887. self.polycoords.append((x, y))
  1888. ### x, y = self.Pixel2Cell(self.pdcTmp.GetIdBounds(self.vdigitMove['id'][0])[0:2])
  1889. self.polycoords.append(self.Pixel2Cell(self.mouse['end']))
  1890. if self.vdigitMove['id'][2] > 0: # next vertex
  1891. x, y = self.Pixel2Cell(self.pdcTmp.GetIdBounds(self.vdigitMove['id'][2])[0:2])
  1892. self.pdcTmp.RemoveId(self.vdigitMove['id'][2]-1)
  1893. self.polycoords.append((x, y))
  1894. self.ClearLines(pdc=self.pdcTmp)
  1895. self.DrawLines(pdc=self.pdcTmp)
  1896. else: # edit line
  1897. try:
  1898. if self.vdigitMove['id'][-1] > 0: # previous vertex
  1899. self.MouseDraw(pdc=self.pdcTmp,
  1900. begin=self.Cell2Pixel(self.polycoords[-1]))
  1901. except: # no line
  1902. self.vdigitMove['id'] = []
  1903. self.polycoords = []
  1904. self.Refresh() # TODO: use RefreshRect()
  1905. self.mouse['begin'] = self.mouse['end']
  1906. elif digitToolbar.GetAction() == "zbulkLine":
  1907. if len(self.polycoords) == 1:
  1908. # draw mouse moving
  1909. self.MouseDraw(self.pdcTmp)
  1910. event.Skip()
  1911. def ClearLines(self, pdc=None):
  1912. """!
  1913. Clears temporary drawn lines from PseudoDC
  1914. """
  1915. if not pdc:
  1916. pdc=self.pdcTmp
  1917. try:
  1918. pdc.ClearId(self.lineid)
  1919. pdc.RemoveId(self.lineid)
  1920. except:
  1921. pass
  1922. try:
  1923. pdc.ClearId(self.plineid)
  1924. pdc.RemoveId(self.plineid)
  1925. except:
  1926. pass
  1927. Debug.msg(4, "BufferedWindow.ClearLines(): lineid=%s, plineid=%s" %
  1928. (self.lineid, self.plineid))
  1929. ### self.Refresh()
  1930. return True
  1931. def Pixel2Cell(self, (x, y)):
  1932. """!
  1933. Convert image coordinates to real word coordinates
  1934. Input : int x, int y
  1935. Output: float x, float y
  1936. """
  1937. try:
  1938. x = int(x)
  1939. y = int(y)
  1940. except:
  1941. return None
  1942. if self.Map.region["ewres"] > self.Map.region["nsres"]:
  1943. res = self.Map.region["ewres"]
  1944. else:
  1945. res = self.Map.region["nsres"]
  1946. w = self.Map.region["center_easting"] - (self.Map.width / 2) * res
  1947. n = self.Map.region["center_northing"] + (self.Map.height / 2) * res
  1948. east = w + x * res
  1949. north = n - y * res
  1950. # extent does not correspond with whole map canvas area...
  1951. # east = self.Map.region['w'] + x * self.Map.region["ewres"]
  1952. # north = self.Map.region['n'] - y * self.Map.region["nsres"]
  1953. return (east, north)
  1954. def Cell2Pixel(self, (east, north)):
  1955. """!
  1956. Convert real word coordinates to image coordinates
  1957. """
  1958. try:
  1959. east = float(east)
  1960. north = float(north)
  1961. except:
  1962. return None
  1963. if self.Map.region["ewres"] > self.Map.region["nsres"]:
  1964. res = self.Map.region["ewres"]
  1965. else:
  1966. res = self.Map.region["nsres"]
  1967. w = self.Map.region["center_easting"] - (self.Map.width / 2) * res
  1968. n = self.Map.region["center_northing"] + (self.Map.height / 2) * res
  1969. # x = int((east - w) / res)
  1970. # y = int((n - north) / res)
  1971. x = (east - w) / res
  1972. y = (n - north) / res
  1973. return (x, y)
  1974. def Zoom(self, begin, end, zoomtype):
  1975. """!
  1976. Calculates new region while (un)zoom/pan-ing
  1977. """
  1978. x1, y1 = begin
  1979. x2, y2 = end
  1980. newreg = {}
  1981. # threshold - too small squares do not make sense
  1982. # can only zoom to windows of > 5x5 screen pixels
  1983. if abs(x2-x1) > 5 and abs(y2-y1) > 5 and zoomtype != 0:
  1984. if x1 > x2:
  1985. x1, x2 = x2, x1
  1986. if y1 > y2:
  1987. y1, y2 = y2, y1
  1988. # zoom in
  1989. if zoomtype > 0:
  1990. newreg['w'], newreg['n'] = self.Pixel2Cell((x1, y1))
  1991. newreg['e'], newreg['s'] = self.Pixel2Cell((x2, y2))
  1992. # zoom out
  1993. elif zoomtype < 0:
  1994. newreg['w'], newreg['n'] = self.Pixel2Cell((-x1 * 2, -y1 * 2))
  1995. newreg['e'], newreg['s'] = self.Pixel2Cell((self.Map.width + 2 * \
  1996. (self.Map.width - x2),
  1997. self.Map.height + 2 * \
  1998. (self.Map.height - y2)))
  1999. # pan
  2000. elif zoomtype == 0:
  2001. dx = x1 - x2
  2002. dy = y1 - y2
  2003. newreg['w'], newreg['n'] = self.Pixel2Cell((dx, dy))
  2004. newreg['e'], newreg['s'] = self.Pixel2Cell((self.Map.width + dx,
  2005. self.Map.height + dy))
  2006. # if new region has been calculated, set the values
  2007. if newreg != {}:
  2008. # LL locations
  2009. if self.parent.Map.projinfo['proj'] == 'll':
  2010. if newreg['n'] > 90.0:
  2011. newreg['n'] = 90.0
  2012. if newreg['s'] < -90.0:
  2013. newreg['s'] = -90.0
  2014. ce = newreg['w'] + (newreg['e'] - newreg['w']) / 2
  2015. cn = newreg['s'] + (newreg['n'] - newreg['s']) / 2
  2016. if hasattr(self, "vdigitMove"):
  2017. # xo = self.Cell2Pixel((self.Map.region['center_easting'], self.Map.region['center_northing']))
  2018. # xn = self.Cell2Pixel(ce, cn))
  2019. tmp = self.Pixel2Cell(self.mouse['end'])
  2020. # calculate new center point and display resolution
  2021. self.Map.region['center_easting'] = ce
  2022. self.Map.region['center_northing'] = cn
  2023. self.Map.region["ewres"] = (newreg['e'] - newreg['w']) / self.Map.width
  2024. self.Map.region["nsres"] = (newreg['n'] - newreg['s']) / self.Map.height
  2025. self.Map.AlignExtentFromDisplay()
  2026. if hasattr(self, "vdigitMove"):
  2027. tmp1 = self.mouse['end']
  2028. tmp2 = self.Cell2Pixel(self.vdigitMove['begin'])
  2029. dx = tmp1[0] - tmp2[0]
  2030. dy = tmp1[1] - tmp2[1]
  2031. self.vdigitMove['beginDiff'] = (dx, dy)
  2032. for id in self.vdigitMove['id']:
  2033. self.pdcTmp.RemoveId(id)
  2034. self.ZoomHistory(self.Map.region['n'], self.Map.region['s'],
  2035. self.Map.region['e'], self.Map.region['w'])
  2036. if self.redrawAll is False:
  2037. self.redrawAll = True
  2038. def ZoomBack(self):
  2039. """!
  2040. Zoom to previous extents in zoomhistory list
  2041. """
  2042. zoom = []
  2043. if len(self.zoomhistory) > 1:
  2044. self.zoomhistory.pop()
  2045. zoom = self.zoomhistory[len(self.zoomhistory)-1]
  2046. # (n, s, e, w)
  2047. if zoom:
  2048. # zoom to selected region
  2049. self.Map.region['center_easting'] = zoom[3] + \
  2050. (zoom[2] - zoom[3]) / 2
  2051. self.Map.region['center_northing'] = zoom[1] + \
  2052. (zoom[0] - zoom[1]) / 2
  2053. self.Map.region["ewres"] = (zoom[2] - zoom[3]) / self.Map.width
  2054. self.Map.region["nsres"] = (zoom[0] - zoom[1]) / self.Map.height
  2055. self.Map.AlignExtentFromDisplay()
  2056. # update map
  2057. self.UpdateMap()
  2058. # update statusbar
  2059. self.parent.StatusbarUpdate()
  2060. def ZoomHistory(self, n, s, e, w):
  2061. """!
  2062. Manages a list of last 10 zoom extents
  2063. Return removed history item if exists
  2064. """
  2065. removed = None
  2066. self.zoomhistory.append((n,s,e,w))
  2067. if len(self.zoomhistory) > 10:
  2068. removed = self.zoomhistory.pop(0)
  2069. if removed:
  2070. Debug.msg(4, "BufferedWindow.ZoomHistory(): hist=%s, removed=%s" %
  2071. (self.zoomhistory, removed))
  2072. else:
  2073. Debug.msg(4, "BufferedWindow.ZoomHistory(): hist=%s" %
  2074. (self.zoomhistory))
  2075. return removed
  2076. def OnZoomToMap(self, event):
  2077. """!
  2078. Set display extents to match selected raster (including NULLs)
  2079. or vector map.
  2080. """
  2081. self.ZoomToMap()
  2082. def OnZoomToRaster(self, event):
  2083. """!
  2084. Set display extents to match selected raster map (ignore NULLs)
  2085. """
  2086. self.ZoomToMap(ignoreNulls = True)
  2087. def ZoomToMap(self, layers = None, ignoreNulls = False, render = True):
  2088. """!
  2089. Set display extents to match selected raster
  2090. or vector map(s).
  2091. @param layer list of layers to be zoom to
  2092. @param ignoreNulls True to ignore null-values
  2093. @param render True to re-render display
  2094. """
  2095. zoomreg = {}
  2096. if not layers:
  2097. layers = self.GetSelectedLayer(multi = True)
  2098. if not layers:
  2099. return
  2100. rast = []
  2101. vect = []
  2102. updated = False
  2103. for l in layers:
  2104. # only raster/vector layers are currently supported
  2105. if l.type == 'raster':
  2106. rast.append(l.name)
  2107. elif l.type == 'vector':
  2108. if self.parent.digit and l.name == self.parent.digit.map and \
  2109. self.parent.digit.type == 'vdigit':
  2110. w, s, b, e, n, t = self.parent.digit.driver.GetMapBoundingBox()
  2111. self.Map.GetRegion(n=n, s=s, w=w, e=e,
  2112. update=True)
  2113. updated = True
  2114. else:
  2115. vect.append(l.name)
  2116. if not updated:
  2117. self.Map.GetRegion(rast = rast,
  2118. vect = vect,
  2119. update = True)
  2120. self.ZoomHistory(self.Map.region['n'], self.Map.region['s'],
  2121. self.Map.region['e'], self.Map.region['w'])
  2122. if render:
  2123. self.UpdateMap()
  2124. self.parent.StatusbarUpdate()
  2125. def ZoomToWind(self, event):
  2126. """!
  2127. Set display geometry to match computational
  2128. region settings (set with g.region)
  2129. """
  2130. self.Map.region = self.Map.GetRegion()
  2131. ### self.Map.SetRegion(windres=True)
  2132. self.ZoomHistory(self.Map.region['n'], self.Map.region['s'],
  2133. self.Map.region['e'], self.Map.region['w'])
  2134. self.UpdateMap()
  2135. self.parent.StatusbarUpdate()
  2136. def ZoomToDefault(self, event):
  2137. """!
  2138. Set display geometry to match default region settings
  2139. """
  2140. self.Map.region = self.Map.GetRegion(default=True)
  2141. self.Map.AdjustRegion() # aling region extent to the display
  2142. self.ZoomHistory(self.Map.region['n'], self.Map.region['s'],
  2143. self.Map.region['e'], self.Map.region['w'])
  2144. self.UpdateMap()
  2145. self.parent.StatusbarUpdate()
  2146. def DisplayToWind(self, event):
  2147. """!
  2148. Set computational region (WIND file) to
  2149. match display extents
  2150. """
  2151. tmpreg = os.getenv("GRASS_REGION")
  2152. if tmpreg:
  2153. del os.environ["GRASS_REGION"]
  2154. # We ONLY want to set extents here. Don't mess with resolution. Leave that
  2155. # for user to set explicitly with g.region
  2156. new = self.Map.AlignResolution()
  2157. gcmd.RunCommand('g.region',
  2158. parent = self,
  2159. overwrite = True,
  2160. n = new['n'],
  2161. s = new['s'],
  2162. e = new['e'],
  2163. w = new['w'],
  2164. rows = int(new['rows']),
  2165. cols = int(new['cols']))
  2166. if tmpreg:
  2167. os.environ["GRASS_REGION"] = tmpreg
  2168. def ZoomToSaved(self, event):
  2169. """!Set display geometry to match extents in
  2170. saved region file
  2171. """
  2172. dlg = gdialogs.SavedRegion(parent = self,
  2173. title = _("Zoom to saved region extents"),
  2174. loadsave='load')
  2175. if dlg.ShowModal() == wx.ID_CANCEL or not dlg.wind:
  2176. dlg.Destroy()
  2177. return
  2178. if not grass.find_file(name = dlg.wind, element = 'windows')['name']:
  2179. wx.MessageBox(parent = self,
  2180. message = _("Region <%s> not found. Operation canceled.") % dlg.wind,
  2181. caption = _("Error"), style = wx.ICON_ERROR | wx.OK | wx.CENTRE)
  2182. dlg.Destroy()
  2183. return
  2184. self.Map.GetRegion(regionName = dlg.wind,
  2185. update = True)
  2186. dlg.Destroy()
  2187. self.ZoomHistory(self.Map.region['n'],
  2188. self.Map.region['s'],
  2189. self.Map.region['e'],
  2190. self.Map.region['w'])
  2191. self.UpdateMap()
  2192. def SaveDisplayRegion(self, event):
  2193. """!
  2194. Save display extents to named region file.
  2195. """
  2196. dlg = gdialogs.SavedRegion(parent = self,
  2197. title = _("Save display extents to region file"),
  2198. loadsave='save')
  2199. if dlg.ShowModal() == wx.ID_CANCEL or not dlg.wind:
  2200. dlg.Destroy()
  2201. return
  2202. # test to see if it already exists and ask permission to overwrite
  2203. if grass.find_file(name = dlg.wind, element = 'windows')['name']:
  2204. overwrite = wx.MessageBox(parent = self,
  2205. message = _("Region file <%s> already exists. "
  2206. "Do you want to overwrite it?") % (dlg.wind),
  2207. caption = _("Warning"), style = wx.YES_NO | wx.CENTRE)
  2208. if (overwrite == wx.YES):
  2209. self.SaveRegion(dlg.wind)
  2210. else:
  2211. self.SaveRegion(dlg.wind)
  2212. dlg.Destroy()
  2213. def SaveRegion(self, wind):
  2214. """!Save region settings
  2215. @param wind region name
  2216. """
  2217. new = self.Map.GetCurrentRegion()
  2218. tmpreg = os.getenv("GRASS_REGION")
  2219. if tmpreg:
  2220. del os.environ["GRASS_REGION"]
  2221. gcmd.RunCommand('g.region',
  2222. overwrite = True,
  2223. parent = self,
  2224. flags = 'u',
  2225. n = new['n'],
  2226. s = new['s'],
  2227. e = new['e'],
  2228. w = new['w'],
  2229. rows = int(new['rows']),
  2230. cols = int(new['cols']),
  2231. save = wind)
  2232. if tmpreg:
  2233. os.environ["GRASS_REGION"] = tmpreg
  2234. def Distance(self, beginpt, endpt, screen=True):
  2235. """!Calculete distance
  2236. LL-locations not supported
  2237. @todo Use m.distance
  2238. @param beginpt first point
  2239. @param endpt second point
  2240. @param screen True for screen coordinates otherwise EN
  2241. """
  2242. x1, y1 = beginpt
  2243. x2, y2 = endpt
  2244. if screen:
  2245. dEast = (x2 - x1) * self.Map.region["ewres"]
  2246. dNorth = (y2 - y1) * self.Map.region["nsres"]
  2247. else:
  2248. dEast = (x2 - x1)
  2249. dNorth = (y2 - y1)
  2250. return (math.sqrt(math.pow((dEast),2) + math.pow((dNorth),2)), (dEast, dNorth))