frame.py 54 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398
  1. """!
  2. @package mapdisp.frame
  3. @brief Map display with toolbar for various display management
  4. functions, and additional toolbars (vector digitizer, 3d view).
  5. Can be used either from Layer Manager or as d.mon backend.
  6. Classes:
  7. - mapdisp::MapFrame
  8. (C) 2006-2012 by the GRASS Development Team
  9. This program is free software under the GNU General Public License
  10. (>=v2). Read the file COPYING that comes with GRASS for details.
  11. @author Michael Barton
  12. @author Jachym Cepicky
  13. @author Martin Landa <landa.martin gmail.com>
  14. @author Vaclav Petras <wenzeslaus gmail.com> (SingleMapFrame, handlers support)
  15. @author Anna Kratochvilova <kratochanna gmail.com> (SingleMapFrame)
  16. @author Stepan Turek <stepan.turek seznam.cz> (handlers support)
  17. """
  18. import os
  19. import sys
  20. import math
  21. import copy
  22. from core import globalvar
  23. import wx
  24. import wx.aui
  25. sys.path.append(os.path.join(globalvar.ETCWXDIR, "icons"))
  26. sys.path.append(os.path.join(globalvar.ETCDIR, "python"))
  27. from core import globalvar
  28. from core.render import EVT_UPDATE_PRGBAR
  29. from vdigit.toolbars import VDigitToolbar
  30. from mapdisp.toolbars import MapToolbar, NvizIcons
  31. from mapdisp.gprint import PrintOptions
  32. from core.gcmd import GError, GMessage, RunCommand
  33. from dbmgr.dialogs import DisplayAttributesDialog
  34. from core.utils import ListOfCatsToRange, GetLayerNameFromCmd
  35. from gui_core.dialogs import GetImageHandlers, ImageSizeDialog, DecorationDialog, TextLayerDialog
  36. from core.debug import Debug
  37. from core.settings import UserSettings
  38. from gui_core.mapdisp import SingleMapFrame
  39. from mapdisp.mapwindow import BufferedWindow
  40. from modules.histogram import HistogramFrame
  41. from wxplot.histogram import Histogram2Frame
  42. from wxplot.profile import ProfileFrame
  43. from wxplot.scatter import ScatterFrame
  44. from mapdisp import statusbar as sb
  45. import grass.script as grass
  46. haveCtypes = False
  47. class MapFrame(SingleMapFrame):
  48. """!Main frame for map display window. Drawing takes place in
  49. child double buffered drawing window.
  50. """
  51. def __init__(self, parent = None, title = _("GRASS GIS - Map display"),
  52. toolbars = ["map"], tree = None, notebook = None, lmgr = None,
  53. page = None, Map = None, auimgr = None, name = 'MapWindow', **kwargs):
  54. """!Main map display window with toolbars, statusbar and
  55. BufferedWindow (map canvas)
  56. @param toolbars array of activated toolbars, e.g. ['map', 'digit']
  57. @param tree reference to layer tree
  58. @param notebook control book ID in Layer Manager
  59. @param lmgr Layer Manager
  60. @param page notebook page with layer tree
  61. @param Map instance of render.Map
  62. @param auimgs AUI manager
  63. @param name frame name
  64. @param kwargs wx.Frame attributes
  65. """
  66. SingleMapFrame.__init__(self, parent = parent, title = title,
  67. Map = Map, auimgr = auimgr, name = name, **kwargs)
  68. self._layerManager = lmgr # Layer Manager object
  69. self.tree = tree # Layer Manager layer tree object
  70. self.page = page # Notebook page holding the layer tree
  71. self.layerbook = notebook # Layer Manager layer tree notebook
  72. #
  73. # Add toolbars
  74. #
  75. for toolb in toolbars:
  76. self.AddToolbar(toolb)
  77. #
  78. # Add statusbar
  79. #
  80. # items for choice
  81. self.statusbarItems = [sb.SbCoordinates,
  82. sb.SbRegionExtent,
  83. sb.SbCompRegionExtent,
  84. sb.SbShowRegion,
  85. sb.SbAlignExtent,
  86. sb.SbResolution,
  87. sb.SbDisplayGeometry,
  88. sb.SbMapScale,
  89. sb.SbGoTo,
  90. sb.SbProjection]
  91. self.statusbarItemsHiddenInNviz = (sb.SbAlignExtent,
  92. sb.SbDisplayGeometry,
  93. sb.SbShowRegion,
  94. sb.SbResolution,
  95. sb.SbMapScale)
  96. # create statusbar and its manager
  97. statusbar = self.CreateStatusBar(number = 4, style = 0)
  98. statusbar.SetStatusWidths([-5, -2, -1, -1])
  99. self.statusbarManager = sb.SbManager(mapframe = self, statusbar = statusbar)
  100. # fill statusbar manager
  101. self.statusbarManager.AddStatusbarItemsByClass(self.statusbarItems, mapframe = self, statusbar = statusbar)
  102. self.statusbarManager.AddStatusbarItem(sb.SbMask(self, statusbar = statusbar, position = 2))
  103. self.statusbarManager.AddStatusbarItem(sb.SbRender(self, statusbar = statusbar, position = 3))
  104. self.statusbarManager.Update()
  105. #
  106. # Init map display (buffered DC & set default cursor)
  107. #
  108. self.MapWindow2D = BufferedWindow(self, id = wx.ID_ANY,
  109. Map = self.Map, frame = self, tree = self.tree, lmgr = self._layerManager)
  110. # default is 2D display mode
  111. self.MapWindow = self.MapWindow2D
  112. self.MapWindow.SetCursor(self.cursors["default"])
  113. # used by vector digitizer
  114. self.MapWindowVDigit = None
  115. # used by Nviz (3D display mode)
  116. self.MapWindow3D = None
  117. #
  118. # initialize region values
  119. #
  120. self._initMap(Map = self.Map)
  121. #
  122. # Bind various events
  123. #
  124. self.Bind(wx.EVT_ACTIVATE, self.OnFocus)
  125. self.Bind(wx.EVT_CLOSE, self.OnCloseWindow)
  126. self.Bind(EVT_UPDATE_PRGBAR, self.OnUpdateProgress)
  127. #
  128. # Update fancy gui style
  129. #
  130. self._mgr.AddPane(self.MapWindow, wx.aui.AuiPaneInfo().CentrePane().
  131. Dockable(False).BestSize((-1,-1)).Name('2d').
  132. CloseButton(False).DestroyOnClose(True).
  133. Layer(0))
  134. self._mgr.Update()
  135. #
  136. # Init print module and classes
  137. #
  138. self.printopt = PrintOptions(self, self.MapWindow)
  139. #
  140. # Init zoom history
  141. #
  142. self.MapWindow.ZoomHistory(self.Map.region['n'],
  143. self.Map.region['s'],
  144. self.Map.region['e'],
  145. self.Map.region['w'])
  146. #
  147. # Re-use dialogs
  148. #
  149. self.dialogs = {}
  150. self.dialogs['attributes'] = None
  151. self.dialogs['category'] = None
  152. self.dialogs['barscale'] = None
  153. self.dialogs['legend'] = None
  154. self.dialogs['vnet'] = None
  155. self.decorationDialog = None # decoration/overlays
  156. def GetMapWindow(self):
  157. return self.MapWindow
  158. def _addToolbarVDigit(self):
  159. """!Add vector digitizer toolbar
  160. """
  161. from vdigit.main import haveVDigit, VDigit
  162. if not haveVDigit:
  163. from vdigit import errorMsg
  164. msg = _("Unable to start wxGUI vector digitizer.\nDo you want to start "
  165. "TCL/TK digitizer (v.digit) instead?\n\n"
  166. "Details: %s" % errorMsg)
  167. self.toolbars['map'].combo.SetValue(_("2D view"))
  168. dlg = wx.MessageDialog(parent = self,
  169. message = msg,
  170. caption=_("Vector digitizer failed"),
  171. style = wx.YES_NO | wx.CENTRE)
  172. if dlg.ShowModal() == wx.ID_YES:
  173. mapName = self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].GetName()
  174. self._layerManager.goutput.RunCmd(['v.digit', 'map=%s' % mapName],
  175. switchPage = False)
  176. dlg.Destroy()
  177. self.toolbars['map'].combo.SetValue(_("2D view"))
  178. return
  179. if self._layerManager:
  180. log = self._layerManager.goutput
  181. else:
  182. log = None
  183. if not self.MapWindowVDigit:
  184. from vdigit.mapwindow import VDigitWindow
  185. self.MapWindowVDigit = VDigitWindow(self, id = wx.ID_ANY, frame = self,
  186. Map = self.Map, tree = self.tree,
  187. lmgr = self._layerManager)
  188. self.MapWindowVDigit.Show()
  189. self._mgr.AddPane(self.MapWindowVDigit, wx.aui.AuiPaneInfo().CentrePane().
  190. Dockable(False).BestSize((-1,-1)).Name('vdigit').
  191. CloseButton(False).DestroyOnClose(True).
  192. Layer(0))
  193. self.MapWindow = self.MapWindowVDigit
  194. if self._mgr.GetPane('2d').IsShown():
  195. self._mgr.GetPane('2d').Hide()
  196. elif self._mgr.GetPane('3d').IsShown():
  197. self._mgr.GetPane('3d').Hide()
  198. self._mgr.GetPane('vdigit').Show()
  199. self.toolbars['vdigit'] = VDigitToolbar(parent = self, MapWindow = self.MapWindow,
  200. digitClass = VDigit, layerTree = self.tree,
  201. log = log)
  202. self.MapWindowVDigit.SetToolbar(self.toolbars['vdigit'])
  203. self._mgr.AddPane(self.toolbars['vdigit'],
  204. wx.aui.AuiPaneInfo().
  205. Name("vdigittoolbar").Caption(_("Vector Digitizer Toolbar")).
  206. ToolbarPane().Top().Row(1).
  207. LeftDockable(False).RightDockable(False).
  208. BottomDockable(False).TopDockable(True).
  209. CloseButton(False).Layer(2).
  210. BestSize((self.toolbars['vdigit'].GetBestSize())))
  211. # change mouse to draw digitized line
  212. self.MapWindow.mouse['box'] = "point"
  213. self.MapWindow.zoomtype = 0
  214. self.MapWindow.pen = wx.Pen(colour = 'red', width = 2, style = wx.SOLID)
  215. self.MapWindow.polypen = wx.Pen(colour = 'green', width = 2, style = wx.SOLID)
  216. def AddNviz(self):
  217. """!Add 3D view mode window
  218. """
  219. from nviz.main import haveNviz, GLWindow
  220. # check for GLCanvas and OpenGL
  221. if not haveNviz:
  222. self.toolbars['map'].combo.SetValue(_("2D view"))
  223. GError(parent = self,
  224. message = _("Unable to switch to 3D display mode.\nThe Nviz python extension "
  225. "was not found or loaded properly.\n"
  226. "Switching back to 2D display mode.\n\nDetails: %s" % nviz.errorMsg))
  227. return
  228. # disable 3D mode for other displays
  229. for page in range(0, self._layerManager.GetLayerNotebook().GetPageCount()):
  230. mapdisp = self._layerManager.GetLayerNotebook().GetPage(page).maptree.GetMapDisplay()
  231. if self._layerManager.GetLayerNotebook().GetPage(page) != self._layerManager.currentPage:
  232. if '3D' in mapdisp.toolbars['map'].combo.GetString(1):
  233. mapdisp.toolbars['map'].combo.Delete(1)
  234. self.toolbars['map'].Enable2D(False)
  235. # add rotate tool to map toolbar
  236. self.toolbars['map'].InsertTool((('rotate', NvizIcons['rotate'],
  237. self.OnRotate, wx.ITEM_CHECK, 7),)) # 7 is position
  238. self.toolbars['map'].InsertTool((('flyThrough', NvizIcons['flyThrough'],
  239. self.OnFlyThrough, wx.ITEM_CHECK, 8),))
  240. self.toolbars['map'].ChangeToolsDesc(mode2d = False)
  241. # update status bar
  242. self.statusbarManager.HideStatusbarChoiceItemsByClass(self.statusbarItemsHiddenInNviz)
  243. self.statusbarManager.SetMode(0)
  244. # erase map window
  245. self.MapWindow.EraseMap()
  246. self._layerManager.goutput.WriteCmdLog(_("Starting 3D view mode..."),
  247. switchPage = False)
  248. self.SetStatusText(_("Please wait, loading data..."), 0)
  249. # create GL window
  250. if not self.MapWindow3D:
  251. self.MapWindow3D = GLWindow(self, id = wx.ID_ANY, frame = self,
  252. Map = self.Map, tree = self.tree, lmgr = self._layerManager)
  253. self.MapWindow = self.MapWindow3D
  254. self.MapWindow.SetCursor(self.cursors["default"])
  255. # add Nviz notebookpage
  256. self._layerManager.AddNvizTools()
  257. # switch from MapWindow to MapWindowGL
  258. self._mgr.GetPane('2d').Hide()
  259. self._mgr.AddPane(self.MapWindow3D, wx.aui.AuiPaneInfo().CentrePane().
  260. Dockable(False).BestSize((-1,-1)).Name('3d').
  261. CloseButton(False).DestroyOnClose(True).
  262. Layer(0))
  263. self.MapWindow3D.Show()
  264. self.MapWindow3D.ResetViewHistory()
  265. self.MapWindow3D.UpdateView(None)
  266. else:
  267. self.MapWindow = self.MapWindow3D
  268. os.environ['GRASS_REGION'] = self.Map.SetRegion(windres = True, windres3 = True)
  269. self.MapWindow3D.GetDisplay().Init()
  270. del os.environ['GRASS_REGION']
  271. # switch from MapWindow to MapWindowGL
  272. self._mgr.GetPane('2d').Hide()
  273. self._mgr.GetPane('3d').Show()
  274. # add Nviz notebookpage
  275. self._layerManager.AddNvizTools()
  276. self.MapWindow3D.ResetViewHistory()
  277. for page in ('view', 'light', 'fringe', 'constant', 'cplane', 'animation'):
  278. self._layerManager.nviz.UpdatePage(page)
  279. self.MapWindow3D.overlays = self.MapWindow2D.overlays
  280. self.MapWindow3D.textdict = self.MapWindow2D.textdict
  281. # update overlays needs to be called after because getClientSize
  282. # is called during update and it must give reasonable values
  283. wx.CallAfter(self.MapWindow3D.UpdateOverlays)
  284. self.SetStatusText("", 0)
  285. self._mgr.Update()
  286. def RemoveNviz(self):
  287. """!Restore 2D view"""
  288. try:
  289. self.toolbars['map'].RemoveTool(self.toolbars['map'].rotate)
  290. self.toolbars['map'].RemoveTool(self.toolbars['map'].flyThrough)
  291. except AttributeError:
  292. pass
  293. # update status bar
  294. self.statusbarManager.ShowStatusbarChoiceItemsByClass(self.statusbarItemsHiddenInNviz)
  295. self.statusbarManager.SetMode(UserSettings.Get(group = 'display',
  296. key = 'statusbarMode',
  297. subkey = 'selection'))
  298. self.SetStatusText(_("Please wait, unloading data..."), 0)
  299. self._layerManager.goutput.WriteCmdLog(_("Switching back to 2D view mode..."),
  300. switchPage = False)
  301. if self.MapWindow3D:
  302. self.MapWindow3D.OnClose(event = None)
  303. # switch from MapWindowGL to MapWindow
  304. self._mgr.GetPane('2d').Show()
  305. self._mgr.GetPane('3d').Hide()
  306. self.MapWindow = self.MapWindow2D
  307. # remove nviz notebook page
  308. self._layerManager.RemoveNvizTools()
  309. try:
  310. self.MapWindow2D.overlays = self.MapWindow3D.overlays
  311. self.MapWindow2D.textdict = self.MapWindow3D.textdict
  312. except AttributeError:
  313. pass
  314. self.MapWindow.UpdateMap()
  315. self._mgr.Update()
  316. def AddToolbar(self, name):
  317. """!Add defined toolbar to the window
  318. Currently known toolbars are:
  319. - 'map' - basic map toolbar
  320. - 'vdigit' - vector digitizer
  321. - 'gcpdisp' - GCP Manager Display
  322. """
  323. # default toolbar
  324. if name == "map":
  325. self.toolbars['map'] = MapToolbar(self, self.Map)
  326. self._mgr.AddPane(self.toolbars['map'],
  327. wx.aui.AuiPaneInfo().
  328. Name("maptoolbar").Caption(_("Map Toolbar")).
  329. ToolbarPane().Top().Name('mapToolbar').
  330. LeftDockable(False).RightDockable(False).
  331. BottomDockable(False).TopDockable(True).
  332. CloseButton(False).Layer(2).
  333. BestSize((self.toolbars['map'].GetBestSize())))
  334. # vector digitizer
  335. elif name == "vdigit":
  336. self._addToolbarVDigit()
  337. self._mgr.Update()
  338. def RemoveToolbar (self, name):
  339. """!Removes defined toolbar from the window
  340. @todo Only hide, activate by calling AddToolbar()
  341. """
  342. # cannot hide main toolbar
  343. if name == "map":
  344. return
  345. self._mgr.DetachPane(self.toolbars[name])
  346. self.toolbars[name].Destroy()
  347. self.toolbars.pop(name)
  348. if name == 'vdigit':
  349. self._mgr.GetPane('vdigit').Hide()
  350. self._mgr.GetPane('2d').Show()
  351. self.MapWindow = self.MapWindow2D
  352. self.toolbars['map'].combo.SetValue(_("2D view"))
  353. self.toolbars['map'].Enable2D(True)
  354. self._mgr.Update()
  355. def IsPaneShown(self, name):
  356. """!Check if pane (toolbar, mapWindow ...) of given name is currently shown"""
  357. if self._mgr.GetPane(name).IsOk():
  358. return self._mgr.GetPane(name).IsShown()
  359. return False
  360. def OnUpdateProgress(self, event):
  361. """!Update progress bar info
  362. """
  363. self.GetProgressBar().SetValue(event.value)
  364. event.Skip()
  365. def OnFocus(self, event):
  366. """!Change choicebook page to match display.
  367. """
  368. # change bookcontrol page to page associated with display
  369. if self.page:
  370. pgnum = self.layerbook.GetPageIndex(self.page)
  371. if pgnum > -1:
  372. self.layerbook.SetSelection(pgnum)
  373. self._layerManager.currentPage = self.layerbook.GetCurrentPage()
  374. event.Skip()
  375. def RemoveQueryLayer(self):
  376. """!Removes temporary map layers (queries)"""
  377. qlayer = self.GetMap().GetListOfLayers(l_name = globalvar.QUERYLAYER)
  378. for layer in qlayer:
  379. self.GetMap().DeleteLayer(layer)
  380. def OnRender(self, event):
  381. """!Re-render map composition (each map layer)
  382. """
  383. self.RemoveQueryLayer()
  384. # delete tmp lines
  385. if self.MapWindow.mouse["use"] in ("measure",
  386. "profile"):
  387. self.MapWindow.polycoords = []
  388. self.MapWindow.ClearLines()
  389. # deselect features in vdigit
  390. if self.GetToolbar('vdigit'):
  391. if self.MapWindow.digit:
  392. self.MapWindow.digit.GetDisplay().SetSelected([])
  393. self.MapWindow.UpdateMap(render = True, renderVector = True)
  394. else:
  395. self.MapWindow.UpdateMap(render = True)
  396. # update statusbar
  397. self.StatusbarUpdate()
  398. def OnPointer(self, event):
  399. """!Pointer button clicked
  400. """
  401. if self.GetMapToolbar():
  402. if event:
  403. self.SwitchTool(self.toolbars['map'], event)
  404. self.toolbars['map'].action['desc'] = ''
  405. self.MapWindow.mouse['use'] = "pointer"
  406. self.MapWindow.mouse['box'] = "point"
  407. # change the cursor
  408. if self.GetToolbar('vdigit'):
  409. # digitization tool activated
  410. self.MapWindow.SetCursor(self.cursors["cross"])
  411. # reset mouse['box'] if needed
  412. if self.toolbars['vdigit'].GetAction() in ['addLine']:
  413. if self.toolbars['vdigit'].GetAction('type') in ['point', 'centroid']:
  414. self.MapWindow.mouse['box'] = 'point'
  415. else: # line, boundary
  416. self.MapWindow.mouse['box'] = 'line'
  417. elif self.toolbars['vdigit'].GetAction() in ['addVertex', 'removeVertex', 'splitLine',
  418. 'editLine', 'displayCats', 'queryMap',
  419. 'copyCats']:
  420. self.MapWindow.mouse['box'] = 'point'
  421. else: # moveLine, deleteLine
  422. self.MapWindow.mouse['box'] = 'box'
  423. else:
  424. self.MapWindow.SetCursor(self.cursors["default"])
  425. def OnRotate(self, event):
  426. """!Rotate 3D view
  427. """
  428. if self.GetMapToolbar():
  429. self.SwitchTool(self.toolbars['map'], event)
  430. self.toolbars['map'].action['desc'] = ''
  431. self.MapWindow.mouse['use'] = "rotate"
  432. # change the cursor
  433. self.MapWindow.SetCursor(self.cursors["hand"])
  434. def OnFlyThrough(self, event):
  435. """!Fly-through mode
  436. """
  437. if self.GetMapToolbar():
  438. self.SwitchTool(self.toolbars['map'], event)
  439. self.toolbars['map'].action['desc'] = ''
  440. self.MapWindow.mouse['use'] = "fly"
  441. # change the cursor
  442. self.MapWindow.SetCursor(self.cursors["hand"])
  443. self.MapWindow.SetFocus()
  444. def OnZoomRegion(self, event):
  445. """!Zoom to region
  446. """
  447. self.Map.getRegion()
  448. self.Map.getResolution()
  449. self.UpdateMap()
  450. # event.Skip()
  451. def OnAlignRegion(self, event):
  452. """!Align region
  453. """
  454. if not self.Map.alignRegion:
  455. self.Map.alignRegion = True
  456. else:
  457. self.Map.alignRegion = False
  458. # event.Skip()
  459. def SaveToFile(self, event):
  460. """!Save map to image
  461. """
  462. if self.IsPaneShown('3d'):
  463. filetype = "TIF file (*.tif)|*.tif|PPM file (*.ppm)|*.ppm"
  464. ltype = [{ 'ext' : 'tif', 'type' : 'tif' },
  465. { 'ext' : 'ppm', 'type' : 'ppm' }]
  466. else:
  467. img = self.MapWindow.img
  468. if not img:
  469. GMessage(parent = self,
  470. message = _("Nothing to render (empty map). Operation canceled."))
  471. return
  472. filetype, ltype = GetImageHandlers(img)
  473. # get size
  474. dlg = ImageSizeDialog(self)
  475. dlg.CentreOnParent()
  476. if dlg.ShowModal() != wx.ID_OK:
  477. dlg.Destroy()
  478. return
  479. width, height = dlg.GetValues()
  480. dlg.Destroy()
  481. # get filename
  482. dlg = wx.FileDialog(parent = self,
  483. message = _("Choose a file name to save the image "
  484. "(no need to add extension)"),
  485. wildcard = filetype,
  486. style = wx.SAVE | wx.FD_OVERWRITE_PROMPT)
  487. if dlg.ShowModal() == wx.ID_OK:
  488. path = dlg.GetPath()
  489. if not path:
  490. dlg.Destroy()
  491. return
  492. base, ext = os.path.splitext(path)
  493. fileType = ltype[dlg.GetFilterIndex()]['type']
  494. extType = ltype[dlg.GetFilterIndex()]['ext']
  495. if ext != extType:
  496. path = base + '.' + extType
  497. self.MapWindow.SaveToFile(path, fileType,
  498. width, height)
  499. dlg.Destroy()
  500. def PrintMenu(self, event):
  501. """
  502. Print options and output menu for map display
  503. """
  504. point = wx.GetMousePosition()
  505. printmenu = wx.Menu()
  506. # Add items to the menu
  507. setup = wx.MenuItem(printmenu, wx.ID_ANY, _('Page setup'))
  508. printmenu.AppendItem(setup)
  509. self.Bind(wx.EVT_MENU, self.printopt.OnPageSetup, setup)
  510. preview = wx.MenuItem(printmenu, wx.ID_ANY, _('Print preview'))
  511. printmenu.AppendItem(preview)
  512. self.Bind(wx.EVT_MENU, self.printopt.OnPrintPreview, preview)
  513. doprint = wx.MenuItem(printmenu, wx.ID_ANY, _('Print display'))
  514. printmenu.AppendItem(doprint)
  515. self.Bind(wx.EVT_MENU, self.printopt.OnDoPrint, doprint)
  516. # Popup the menu. If an item is selected then its handler
  517. # will be called before PopupMenu returns.
  518. self.PopupMenu(printmenu)
  519. printmenu.Destroy()
  520. def OnCloseWindow(self, event):
  521. """!Window closed.
  522. Also close associated layer tree page
  523. """
  524. pgnum = None
  525. self.Map.Clean()
  526. # close edited map and 3D tools properly
  527. if self.GetToolbar('vdigit'):
  528. maplayer = self.toolbars['vdigit'].GetLayer()
  529. if maplayer:
  530. self.toolbars['vdigit'].OnExit()
  531. if self.IsPaneShown('3d'):
  532. self.RemoveNviz()
  533. if not self._layerManager:
  534. self.Destroy()
  535. elif self.page:
  536. pgnum = self.layerbook.GetPageIndex(self.page)
  537. if pgnum > -1:
  538. self.layerbook.DeletePage(pgnum)
  539. def Query(self, x, y, layers):
  540. """!Query selected layers.
  541. Calls QueryMap in case of raster or more vectors,
  542. or QueryVector in case of one vector with db connection.
  543. @param x,y coordinates
  544. @param layers selected tree item layers
  545. """
  546. num = 0
  547. filteredLayers = []
  548. for layer in layers:
  549. ltype = self.tree.GetPyData(layer)[0]['maplayer'].GetType()
  550. if ltype in ('raster', 'rgb', 'his',
  551. 'vector', 'thememap', 'themechart'):
  552. filteredLayers.append(layer)
  553. if not filteredLayers:
  554. GMessage(parent = self,
  555. message = _('No raster or vector map layer selected for querying.'))
  556. return
  557. layers = filteredLayers
  558. # set query snap distance for v.what at map unit equivalent of 10 pixels
  559. qdist = 10.0 * ((self.Map.region['e'] - self.Map.region['w']) / self.Map.width)
  560. east, north = self.MapWindow.Pixel2Cell((x, y))
  561. posWindow = self.ClientToScreen((x + self.MapWindow.dialogOffset,
  562. y + self.MapWindow.dialogOffset))
  563. isRaster = False
  564. nVectors = 0
  565. isDbConnection = False
  566. for l in layers:
  567. maplayer = self.tree.GetPyData(l)[0]['maplayer']
  568. if maplayer.GetType() == 'raster':
  569. isRaster = True
  570. break
  571. if maplayer.GetType() == 'vector':
  572. nVectors += 1
  573. isDbConnection = grass.vector_db(maplayer.GetName())
  574. if not self.IsPaneShown('3d'):
  575. if isRaster or nVectors > 1 or not isDbConnection:
  576. self.QueryMap(east, north, qdist, layers)
  577. else:
  578. self.QueryVector(east, north, qdist, posWindow, layers[0])
  579. else:
  580. if isRaster:
  581. self.MapWindow.QuerySurface(x, y)
  582. if nVectors > 1 or not isDbConnection:
  583. self.QueryMap(east, north, qdist, layers)
  584. elif nVectors == 1:
  585. self.QueryVector(east, north, qdist, posWindow, layers[0])
  586. def QueryMap(self, east, north, qdist, layers):
  587. """!Query raster or vector map layers by r/v.what
  588. @param east,north coordinates
  589. @param qdist query distance
  590. @param layers selected tree items
  591. """
  592. rast = list()
  593. vect = list()
  594. rcmd = ['r.what', '--v']
  595. vcmd = ['v.what', '--v']
  596. for layer in layers:
  597. ltype = self.tree.GetPyData(layer)[0]['maplayer'].GetType()
  598. dcmd = self.tree.GetPyData(layer)[0]['cmd']
  599. name, found = GetLayerNameFromCmd(dcmd)
  600. if not found:
  601. continue
  602. if ltype == 'raster':
  603. rast.append(name)
  604. elif ltype in ('rgb', 'his'):
  605. for iname in name.split('\n'):
  606. rast.append(iname)
  607. elif ltype in ('vector', 'thememap', 'themechart'):
  608. vect.append(name)
  609. # use display region settings instead of computation region settings
  610. self.tmpreg = os.getenv("GRASS_REGION")
  611. os.environ["GRASS_REGION"] = self.Map.SetRegion(windres = False)
  612. # build query commands for any selected rasters and vectors
  613. if rast:
  614. rcmd.append('-f')
  615. rcmd.append('-n')
  616. rcmd.append('map=%s' % ','.join(rast))
  617. rcmd.append('coordinates=%f,%f' % (float(east), float(north)))
  618. if vect:
  619. # check for vector maps open to be edited
  620. digitToolbar = self.GetToolbar('vdigit')
  621. if digitToolbar:
  622. lmap = digitToolbar.GetLayer().GetName()
  623. for name in vect:
  624. if lmap == name:
  625. self._layerManager.goutput.WriteWarning(_("Vector map <%s> "
  626. "opened for editing - skipped.") % map)
  627. vect.remove(name)
  628. if len(vect) < 1:
  629. self._layerManager.goutput.WriteCmdLog(_("Nothing to query."))
  630. return
  631. vcmd.append('-a')
  632. vcmd.append('map=%s' % ','.join(vect))
  633. vcmd.append('layer=%s' % ','.join(['-1'] * len(vect)))
  634. vcmd.append('coordinates=%f,%f' % (float(east), float(north)))
  635. vcmd.append('distance=%f' % float(qdist))
  636. Debug.msg(1, "QueryMap(): raster=%s vector=%s" % (','.join(rast),
  637. ','.join(vect)))
  638. # parse query command(s)
  639. if rast and not self.IsPaneShown('3d'):
  640. self._layerManager.goutput.RunCmd(rcmd,
  641. compReg = False,
  642. onDone = self._QueryMapDone)
  643. if vect:
  644. self._layerManager.goutput.RunCmd(vcmd,
  645. onDone = self._QueryMapDone)
  646. def _QueryMapDone(self, cmd, returncode):
  647. """!Restore settings after querying (restore GRASS_REGION)
  648. @param returncode command return code
  649. """
  650. if hasattr(self, "tmpreg"):
  651. if self.tmpreg:
  652. os.environ["GRASS_REGION"] = self.tmpreg
  653. elif 'GRASS_REGION' in os.environ:
  654. del os.environ["GRASS_REGION"]
  655. elif 'GRASS_REGION' in os.environ:
  656. del os.environ["GRASS_REGION"]
  657. if hasattr(self, "tmpreg"):
  658. del self.tmpreg
  659. def QueryVector(self, east, north, qdist, posWindow, layer):
  660. """!Query vector map layer features
  661. Attribute data of selected vector object are displayed in GUI dialog.
  662. Data can be modified (On Submit)
  663. """
  664. mapName = self.tree.GetPyData(layer)[0]['maplayer'].name
  665. if self.tree.GetPyData(layer)[0]['maplayer'].GetMapset() != \
  666. grass.gisenv()['MAPSET']:
  667. mode = 'display'
  668. else:
  669. mode = 'update'
  670. if self.dialogs['attributes'] is None:
  671. dlg = DisplayAttributesDialog(parent = self.MapWindow,
  672. map = mapName,
  673. query = ((east, north), qdist),
  674. pos = posWindow,
  675. action = mode)
  676. self.dialogs['attributes'] = dlg
  677. else:
  678. # selection changed?
  679. if not self.dialogs['attributes'].mapDBInfo or \
  680. self.dialogs['attributes'].mapDBInfo.map != mapName:
  681. self.dialogs['attributes'].UpdateDialog(map = mapName, query = ((east, north), qdist),
  682. action = mode)
  683. else:
  684. self.dialogs['attributes'].UpdateDialog(query = ((east, north), qdist),
  685. action = mode)
  686. if not self.dialogs['attributes'].IsFound():
  687. self._layerManager.goutput.WriteLog(_('Nothing found.'))
  688. cats = self.dialogs['attributes'].GetCats()
  689. qlayer = None
  690. if not self.IsPaneShown('3d') and self.IsAutoRendered():
  691. try:
  692. qlayer = self.Map.GetListOfLayers(l_name = globalvar.QUERYLAYER)[0]
  693. except IndexError:
  694. pass
  695. if self.dialogs['attributes'].mapDBInfo and cats:
  696. if not self.IsPaneShown('3d') and self.IsAutoRendered():
  697. # highlight feature & re-draw map
  698. if qlayer:
  699. qlayer.SetCmd(self.AddTmpVectorMapLayer(mapName, cats,
  700. useId = False,
  701. addLayer = False))
  702. else:
  703. qlayer = self.AddTmpVectorMapLayer(mapName, cats, useId = False)
  704. # set opacity based on queried layer
  705. opacity = self.tree.GetPyData(layer)[0]['maplayer'].GetOpacity(float = True)
  706. qlayer.SetOpacity(opacity)
  707. self.MapWindow.UpdateMap(render = False, renderVector = False)
  708. if not self.dialogs['attributes'].IsShown():
  709. self.dialogs['attributes'].Show()
  710. else:
  711. if qlayer:
  712. self.Map.DeleteLayer(qlayer)
  713. self.MapWindow.UpdateMap(render = False, renderVector = False)
  714. if self.dialogs['attributes'].IsShown():
  715. self.dialogs['attributes'].Hide()
  716. def OnQuery(self, event):
  717. """!Query tools menu"""
  718. if self.GetMapToolbar():
  719. self.SwitchTool(self.toolbars['map'], event)
  720. action = self.toolbars['map'].GetAction()
  721. self.toolbars['map'].action['desc'] = 'queryMap'
  722. self.MapWindow.mouse['use'] = "query"
  723. if not self.IsStandalone():
  724. # switch to output console to show query results
  725. self._layerManager.notebook.SetSelectionByName('output')
  726. self.MapWindow.mouse['box'] = "point"
  727. self.MapWindow.zoomtype = 0
  728. # change the cursor
  729. self.MapWindow.SetCursor(self.cursors["cross"])
  730. def AddTmpVectorMapLayer(self, name, cats, useId = False, addLayer = True):
  731. """!Add temporal vector map layer to map composition
  732. @param name name of map layer
  733. @param useId use feature id instead of category
  734. """
  735. # color settings from ATM
  736. color = UserSettings.Get(group = 'atm', key = 'highlight', subkey = 'color')
  737. colorStr = str(color[0]) + ":" + \
  738. str(color[1]) + ":" + \
  739. str(color[2])
  740. # icon used in vector display and its size
  741. icon = ''
  742. size = 0
  743. vparam = self.tree.GetPyData(self.tree.layer_selected)[0]['cmd']
  744. for p in vparam:
  745. if '=' in p:
  746. parg,pval = p.split('=', 1)
  747. if parg == 'icon': icon = pval
  748. elif parg == 'size': size = float(pval)
  749. pattern = ["d.vect",
  750. "map=%s" % name,
  751. "color=%s" % colorStr,
  752. "fcolor=%s" % colorStr,
  753. "width=%d" % UserSettings.Get(group = 'atm', key = 'highlight', subkey = 'width')]
  754. if icon != '':
  755. pattern.append('icon=%s' % icon)
  756. if size > 0:
  757. pattern.append('size=%i' % size)
  758. if useId:
  759. cmd = pattern
  760. cmd.append('-i')
  761. cmd.append('cats=%s' % str(cats))
  762. else:
  763. cmd = []
  764. for layer in cats.keys():
  765. cmd.append(copy.copy(pattern))
  766. lcats = cats[layer]
  767. cmd[-1].append("layer=%d" % layer)
  768. cmd[-1].append("cats=%s" % ListOfCatsToRange(lcats))
  769. if addLayer:
  770. if useId:
  771. return self.Map.AddLayer(type = 'vector', name = globalvar.QUERYLAYER, command = cmd,
  772. l_active = True, l_hidden = True, l_opacity = 1.0)
  773. else:
  774. return self.Map.AddLayer(type = 'command', name = globalvar.QUERYLAYER, command = cmd,
  775. l_active = True, l_hidden = True, l_opacity = 1.0)
  776. else:
  777. return cmd
  778. def OnMeasure(self, event):
  779. """!Init measurement routine that calculates map distance
  780. along transect drawn on map display
  781. """
  782. self.totaldist = 0.0 # total measured distance
  783. self.SwitchTool(self.toolbars['map'], event)
  784. # switch Layer Manager to output console to show measure results
  785. self._layerManager.notebook.SetSelectionByName('output')
  786. # change mouse to draw line for measurement
  787. self.MapWindow.mouse['use'] = "measure"
  788. self.MapWindow.mouse['box'] = "line"
  789. self.MapWindow.zoomtype = 0
  790. self.MapWindow.pen = wx.Pen(colour = 'red', width = 2, style = wx.SHORT_DASH)
  791. self.MapWindow.polypen = wx.Pen(colour = 'green', width = 2, style = wx.SHORT_DASH)
  792. # change the cursor
  793. self.MapWindow.SetCursor(self.cursors["pencil"])
  794. # initiating output
  795. style = self._layerManager.goutput.cmdOutput.StyleWarning
  796. self._layerManager.goutput.WriteLog(_('Click and drag with left mouse button '
  797. 'to measure.%s'
  798. 'Double click with left button to clear.') % \
  799. (os.linesep), style)
  800. if self.Map.projinfo['proj'] != 'xy':
  801. units = self.Map.projinfo['units']
  802. self._layerManager.goutput.WriteCmdLog(_('Measuring distance') + ' ('
  803. + units + '):')
  804. else:
  805. self._layerManager.goutput.WriteCmdLog(_('Measuring distance:'))
  806. if self.Map.projinfo['proj'] == 'll':
  807. try:
  808. import grass.lib.gis as gislib
  809. global haveCtypes
  810. haveCtypes = True
  811. gislib.G_begin_distance_calculations()
  812. except ImportError, e:
  813. self._layerManager.goutput.WriteWarning(_('Geodesic distance is not yet '
  814. 'supported by this tool.\n'
  815. 'Reason: %s' % e))
  816. def MeasureDist(self, beginpt, endpt):
  817. """!Calculate map distance from screen distance
  818. and print to output window
  819. """
  820. self._layerManager.notebook.SetSelectionByName('output')
  821. dist, (north, east) = self.MapWindow.Distance(beginpt, endpt)
  822. dist = round(dist, 3)
  823. d, dunits = self.FormatDist(dist)
  824. self.totaldist += dist
  825. td, tdunits = self.FormatDist(self.totaldist)
  826. strdist = str(d)
  827. strtotdist = str(td)
  828. if self.Map.projinfo['proj'] == 'xy' or 'degree' not in self.Map.projinfo['unit']:
  829. angle = int(math.degrees(math.atan2(north,east)) + 0.5)
  830. angle = 180 - angle
  831. if angle < 0:
  832. angle = 360 + angle
  833. mstring = '%s = %s %s\n%s = %s %s\n%s = %d %s\n%s' \
  834. % (_('segment'), strdist, dunits,
  835. _('total distance'), strtotdist, tdunits,
  836. _('bearing'), angle, _('deg'),
  837. '-' * 60)
  838. else:
  839. mstring = '%s = %s %s\n%s = %s %s\n%s' \
  840. % (_('segment'), strdist, dunits,
  841. _('total distance'), strtotdist, tdunits,
  842. '-' * 60)
  843. self._layerManager.goutput.WriteLog(mstring)
  844. return dist
  845. def OnProfile(self, event):
  846. """!Launch profile tool
  847. """
  848. raster = []
  849. if self.tree.layer_selected and \
  850. self.tree.GetPyData(self.tree.layer_selected)[0]['type'] == 'raster':
  851. raster.append(self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].name)
  852. win = ProfileFrame(parent = self, rasterList = raster)
  853. win.CentreOnParent()
  854. win.Show()
  855. # Open raster select dialog to make sure that a raster (and
  856. # the desired raster) is selected to be profiled
  857. win.OnSelectRaster(None)
  858. def FormatDist(self, dist):
  859. """!Format length numbers and units in a nice way,
  860. as a function of length. From code by Hamish Bowman
  861. Grass Development Team 2006"""
  862. mapunits = self.Map.projinfo['units']
  863. if mapunits == 'metres':
  864. mapunits = 'meters'
  865. outunits = mapunits
  866. dist = float(dist)
  867. divisor = 1.0
  868. # figure out which units to use
  869. if mapunits == 'meters':
  870. if dist > 2500.0:
  871. outunits = 'km'
  872. divisor = 1000.0
  873. else: outunits = 'm'
  874. elif mapunits == 'feet':
  875. # nano-bug: we match any "feet", but US Survey feet is really
  876. # 5279.9894 per statute mile, or 10.6' per 1000 miles. As >1000
  877. # miles the tick markers are rounded to the nearest 10th of a
  878. # mile (528'), the difference in foot flavours is ignored.
  879. if dist > 5280.0:
  880. outunits = 'miles'
  881. divisor = 5280.0
  882. else:
  883. outunits = 'ft'
  884. elif 'degree' in mapunits and \
  885. not haveCtypes:
  886. if dist < 1:
  887. outunits = 'min'
  888. divisor = (1/60.0)
  889. else:
  890. outunits = 'deg'
  891. else:
  892. outunits = 'meters'
  893. # format numbers in a nice way
  894. if (dist/divisor) >= 2500.0:
  895. outdist = round(dist/divisor)
  896. elif (dist/divisor) >= 1000.0:
  897. outdist = round(dist/divisor,1)
  898. elif (dist/divisor) > 0.0:
  899. outdist = round(dist/divisor,int(math.ceil(3-math.log10(dist/divisor))))
  900. else:
  901. outdist = float(dist/divisor)
  902. return (outdist, outunits)
  903. def OnHistogramPyPlot(self, event):
  904. """!Init PyPlot histogram display canvas and tools
  905. """
  906. raster = []
  907. for layer in self.tree.GetSelections():
  908. if self.tree.GetPyData(layer)[0]['maplayer'].GetType() != 'raster':
  909. continue
  910. raster.append(self.tree.GetPyData(layer)[0]['maplayer'].GetName())
  911. win = Histogram2Frame(parent = self, rasterList = raster)
  912. win.CentreOnParent
  913. win.Show()
  914. # Open raster select dialog to make sure that a raster (and the desired raster)
  915. # is selected to be histogrammed
  916. win.OnSelectRaster(None)
  917. def OnScatterplot(self, event):
  918. """!Init PyPlot scatterplot display canvas and tools
  919. """
  920. raster = []
  921. for layer in self.tree.GetSelections():
  922. if self.tree.GetPyData(layer)[0]['maplayer'].GetType() != 'raster':
  923. continue
  924. raster.append(self.tree.GetPyData(layer)[0]['maplayer'].GetName())
  925. win = ScatterFrame(parent = self, rasterList = raster)
  926. win.CentreOnParent()
  927. win.Show()
  928. # Open raster select dialog to make sure that at least 2 rasters (and the desired rasters)
  929. # are selected to be plotted
  930. win.OnSelectRaster(None)
  931. def OnHistogram(self, event):
  932. """!Init histogram display canvas and tools
  933. """
  934. win = HistogramFrame(self)
  935. win.CentreOnParent()
  936. win.Show()
  937. win.Refresh()
  938. win.Update()
  939. def OnAddBarscale(self, event):
  940. """!Handler for scale/arrow map decoration menu selection.
  941. """
  942. if self.dialogs['barscale']:
  943. return
  944. self.SwitchTool(self.toolbars['map'], event)
  945. id = 0 # unique index for overlay layer
  946. # If location is latlon, only display north arrow (scale won't work)
  947. # proj = self.Map.projinfo['proj']
  948. # if proj == 'll':
  949. # barcmd = 'd.barscale -n'
  950. # else:
  951. # barcmd = 'd.barscale'
  952. # decoration overlay control dialog
  953. self.dialogs['barscale'] = \
  954. DecorationDialog(parent = self, title = _('Scale and North arrow'),
  955. size = (350, 200),
  956. style = wx.DEFAULT_DIALOG_STYLE | wx.CENTRE,
  957. cmd = ['d.barscale', 'at=0,95'],
  958. ovlId = id,
  959. name = 'barscale',
  960. checktxt = _("Show/hide scale and North arrow"),
  961. ctrltxt = _("scale object"))
  962. self.dialogs['barscale'].CentreOnParent()
  963. ### dialog cannot be show as modal - in the result d.barscale is not selectable
  964. ### self.dialogs['barscale'].ShowModal()
  965. self.dialogs['barscale'].Show()
  966. self.MapWindow.mouse['use'] = 'pointer'
  967. def OnAddLegend(self, event):
  968. """!Handler for legend map decoration menu selection.
  969. """
  970. if self.dialogs['legend']:
  971. return
  972. id = 1 # index for overlay layer in render
  973. cmd = ['d.legend', 'at=5,50,2,5']
  974. if self.tree and self.tree.layer_selected and \
  975. self.tree.GetPyData(self.tree.layer_selected)[0]['type'] == 'raster':
  976. cmd.append('map=%s' % self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].name)
  977. # Decoration overlay control dialog
  978. self.dialogs['legend'] = \
  979. DecorationDialog(parent = self, title = ('Legend'),
  980. size = (350, 200),
  981. style = wx.DEFAULT_DIALOG_STYLE | wx.CENTRE,
  982. cmd = cmd,
  983. ovlId = id,
  984. name = 'legend',
  985. checktxt = _("Show/hide legend"),
  986. ctrltxt = _("legend object"))
  987. self.dialogs['legend'].CentreOnParent()
  988. ### dialog cannot be show as modal - in the result d.legend is not selectable
  989. ### self.dialogs['legend'].ShowModal()
  990. self.dialogs['legend'].Show()
  991. self.MapWindow.mouse['use'] = 'pointer'
  992. def OnAddText(self, event):
  993. """!Handler for text decoration menu selection.
  994. """
  995. self.SwitchTool(self.toolbars['map'], event)
  996. if self.MapWindow.dragid > -1:
  997. id = self.MapWindow.dragid
  998. self.MapWindow.dragid = -1
  999. else:
  1000. # index for overlay layer in render
  1001. if len(self.MapWindow.textdict.keys()) > 0:
  1002. id = max(self.MapWindow.textdict.keys()) + 1
  1003. else:
  1004. id = 101
  1005. self.dialogs['text'] = TextLayerDialog(parent = self, ovlId = id,
  1006. title = _('Add text layer'),
  1007. size = (400, 200))
  1008. self.dialogs['text'].CenterOnParent()
  1009. # If OK button pressed in decoration control dialog
  1010. if self.dialogs['text'].ShowModal() == wx.ID_OK:
  1011. text = self.dialogs['text'].GetValues()['text']
  1012. active = self.dialogs['text'].GetValues()['active']
  1013. # delete object if it has no text or is not active
  1014. if text == '' or active == False:
  1015. try:
  1016. self.MapWindow2D.pdc.ClearId(id)
  1017. self.MapWindow2D.pdc.RemoveId(id)
  1018. del self.MapWindow.textdict[id]
  1019. if self.IsPaneShown('3d'):
  1020. self.MapWindow3D.UpdateOverlays()
  1021. self.MapWindow.UpdateMap()
  1022. else:
  1023. self.MapWindow2D.UpdateMap(render = False, renderVector = False)
  1024. except:
  1025. pass
  1026. return
  1027. self.MapWindow.textdict[id] = self.dialogs['text'].GetValues()
  1028. if self.IsPaneShown('3d'):
  1029. self.MapWindow3D.UpdateOverlays()
  1030. self.MapWindow3D.UpdateMap()
  1031. else:
  1032. self.MapWindow2D.pdc.ClearId(id)
  1033. self.MapWindow2D.pdc.SetId(id)
  1034. self.MapWindow2D.UpdateMap(render = False, renderVector = False)
  1035. self.MapWindow.mouse['use'] = 'pointer'
  1036. def OnAddArrow(self, event):
  1037. """!Handler for north arrow menu selection.
  1038. Opens Appearance page of nviz notebook.
  1039. """
  1040. self._layerManager.nviz.SetPage('decoration')
  1041. self.MapWindow3D.SetDrawArrow((70, 70))
  1042. def GetOptData(self, dcmd, type, params, propwin):
  1043. """!Callback method for decoration overlay command generated by
  1044. dialog created in menuform.py
  1045. """
  1046. # Reset comand and rendering options in render.Map. Always render decoration.
  1047. # Showing/hiding handled by PseudoDC
  1048. self.Map.ChangeOverlay(ovltype = type, type = 'overlay', name = '', command = dcmd,
  1049. l_active = True, l_render = False)
  1050. self.params[type] = params
  1051. self.propwin[type] = propwin
  1052. def OnZoomToMap(self, event):
  1053. """!Set display extents to match selected raster (including
  1054. NULLs) or vector map.
  1055. """
  1056. layers = None
  1057. if self.IsStandalone():
  1058. layers = self.MapWindow.GetMap().GetListOfLayers(l_active = False)
  1059. self.MapWindow.ZoomToMap(layers = layers)
  1060. def OnZoomToRaster(self, event):
  1061. """!Set display extents to match selected raster map (ignore NULLs)
  1062. """
  1063. self.MapWindow.ZoomToMap(ignoreNulls = True)
  1064. def OnZoomToSaved(self, event):
  1065. """!Set display geometry to match extents in
  1066. saved region file
  1067. """
  1068. self.MapWindow.ZoomToSaved()
  1069. def OnDisplayToWind(self, event):
  1070. """!Set computational region (WIND file) to match display
  1071. extents
  1072. """
  1073. self.MapWindow.DisplayToWind()
  1074. def SaveDisplayRegion(self, event):
  1075. """!Save display extents to named region file.
  1076. """
  1077. self.MapWindow.SaveDisplayRegion()
  1078. def OnZoomMenu(self, event):
  1079. """!Popup Zoom menu
  1080. """
  1081. point = wx.GetMousePosition()
  1082. zoommenu = wx.Menu()
  1083. # Add items to the menu
  1084. zoomwind = wx.MenuItem(zoommenu, wx.ID_ANY, _('Zoom to computational region'))
  1085. zoommenu.AppendItem(zoomwind)
  1086. self.Bind(wx.EVT_MENU, self.OnZoomToWind, zoomwind)
  1087. zoomdefault = wx.MenuItem(zoommenu, wx.ID_ANY, _('Zoom to default region'))
  1088. zoommenu.AppendItem(zoomdefault)
  1089. self.Bind(wx.EVT_MENU, self.OnZoomToDefault, zoomdefault)
  1090. zoomsaved = wx.MenuItem(zoommenu, wx.ID_ANY, _('Zoom to saved region'))
  1091. zoommenu.AppendItem(zoomsaved)
  1092. self.Bind(wx.EVT_MENU, self.OnZoomToSaved, zoomsaved)
  1093. savewind = wx.MenuItem(zoommenu, wx.ID_ANY, _('Set computational region from display extent'))
  1094. zoommenu.AppendItem(savewind)
  1095. self.Bind(wx.EVT_MENU, self.OnDisplayToWind, savewind)
  1096. savezoom = wx.MenuItem(zoommenu, wx.ID_ANY, _('Save display geometry to named region'))
  1097. zoommenu.AppendItem(savezoom)
  1098. self.Bind(wx.EVT_MENU, self.SaveDisplayRegion, savezoom)
  1099. # Popup the menu. If an item is selected then its handler
  1100. # will be called before PopupMenu returns.
  1101. self.PopupMenu(zoommenu)
  1102. zoommenu.Destroy()
  1103. def SetProperties(self, render = False, mode = 0, showCompExtent = False,
  1104. constrainRes = False, projection = False, alignExtent = True):
  1105. """!Set properies of map display window"""
  1106. self.SetProperty('render', render)
  1107. self.statusbarManager.SetMode(mode)
  1108. self.StatusbarUpdate()
  1109. self.SetProperty('region', showCompExtent)
  1110. self.SetProperty('alignExtent', alignExtent)
  1111. self.SetProperty('resolution', constrainRes)
  1112. self.SetProperty('projection', projection)
  1113. def IsStandalone(self):
  1114. """!Check if Map display is standalone"""
  1115. if self._layerManager:
  1116. return False
  1117. return True
  1118. def GetLayerManager(self):
  1119. """!Get reference to Layer Manager
  1120. @return window reference
  1121. @return None (if standalone)
  1122. """
  1123. return self._layerManager
  1124. def GetMapToolbar(self):
  1125. """!Returns toolbar with zooming tools"""
  1126. return self.toolbars['map']
  1127. def OnVnet(self, event):
  1128. """!Dialog for v.net* modules
  1129. """
  1130. if self.dialogs['vnet']:
  1131. return
  1132. else:
  1133. try:
  1134. from vnet.dialogs import VNETDialog
  1135. except ImportError:
  1136. GError(_("Extension <%s> not available, run '%s' to install it.") % \
  1137. ('wx.vnet', 'g.extension -s extension=wx.vnet'),
  1138. parent = self, showTraceback = False)
  1139. return
  1140. self.dialogs['vnet'] = VNETDialog(parent = self)
  1141. self.dialogs['vnet'].CenterOnScreen()
  1142. self.dialogs['vnet'].Show()
  1143. def SwitchTool(self, toolbar, event):
  1144. """!Calls UpdateTools to manage connected toolbars"""
  1145. self.UpdateTools(event)
  1146. SingleMapFrame.SwitchTool(self, toolbar, event)
  1147. def UpdateTools(self, event):
  1148. """!Method deals with relations of toolbars and other
  1149. elements"""
  1150. # untoggles button in other toolbars
  1151. for toolbar in self.toolbars.itervalues():
  1152. if hasattr(event, 'GetEventObject') == True:
  1153. if event.GetEventObject() == toolbar:
  1154. continue
  1155. toolbar.ToggleTool(toolbar.action['id'], False)
  1156. toolbar.action['id'] = -1
  1157. toolbar.OnTool(None)
  1158. # mouse settings
  1159. self.MapWindow.mouse['box'] = 'point'
  1160. self.MapWindow.mouse['use'] = 'pointer'
  1161. # untoggles button in add legend dialog
  1162. if self.dialogs['legend']:
  1163. if hasattr(event ,'GetEventObject'):
  1164. if event.GetEventObject().GetId() == \
  1165. self.dialogs['legend'].FindWindowByName('resize').GetId():
  1166. return
  1167. self.dialogs['legend'].FindWindowByName('resize').SetValue(0)