"""! @package mapdisp.py @brief GIS map display canvas, with toolbar for various display management functions, and additional toolbars (vector digitizer, 3d view). Can be used either from Layer Manager or as p.mon backend. Classes: - MapFrame - MapApp Usage: python mapdisp.py monitor-identifier /path/to/command/file (C) 2006-2011 by the GRASS Development Team This program is free software under the GNU General Public License (>=v2). Read the file COPYING that comes with GRASS for details. @author Michael Barton @author Jachym Cepicky @author Martin Landa """ import os import sys import glob import math import tempfile import copy import globalvar import wx import wx.aui try: import subprocess except: CompatPath = os.path.join(globalvar.ETCWXDIR) sys.path.append(CompatPath) from compat import subprocess gmpath = os.path.join(globalvar.ETCWXDIR, "icons") sys.path.append(gmpath) grassPath = os.path.join(globalvar.ETCDIR, "python") sys.path.append(grassPath) import render import toolbars import menuform import gselect import disp_print import gcmd import dbm import dbm_dialogs import histogram import profile import globalvar import utils import gdialogs from grass.script import core as grass from debug import Debug from icon import Icons from preferences import globalSettings as UserSettings from mapdisp_command import Command from mapdisp_window import BufferedWindow # for standalone app cmdfilename = None haveCtypes = False class MapFrame(wx.Frame): """!Main frame for map display window. Drawing takes place in child double buffered drawing window. """ def __init__(self, parent = None, id = wx.ID_ANY, title = _("GRASS GIS - Map display"), style = wx.DEFAULT_FRAME_STYLE, toolbars = ["map"], tree = None, notebook = None, lmgr = None, page = None, Map = None, auimgr = None, **kwargs): """!Main map display window with toolbars, statusbar and DrawWindow @param toolbars array of activated toolbars, e.g. ['map', 'digit'] @param tree reference to layer tree @param notebook control book ID in Layer Manager @param lmgr Layer Manager @param page notebook page with layer tree @param Map instance of render.Map @param auimgs AUI manager @param kwargs wx.Frame attribures """ self._layerManager = lmgr # Layer Manager object self.Map = Map # instance of render.Map self.tree = tree # Layer Manager layer tree object self.page = page # Notebook page holding the layer tree self.layerbook = notebook # Layer Manager layer tree notebook self.parent = parent if 'name' not in kwargs: kwargs['name'] = 'MapWindow' wx.Frame.__init__(self, parent, id, title, style = style, **kwargs) # available cursors self.cursors = { # default: cross # "default" : wx.StockCursor(wx.CURSOR_DEFAULT), "default" : wx.StockCursor(wx.CURSOR_ARROW), "cross" : wx.StockCursor(wx.CURSOR_CROSS), "hand" : wx.StockCursor(wx.CURSOR_HAND), "pencil" : wx.StockCursor(wx.CURSOR_PENCIL), "sizenwse": wx.StockCursor(wx.CURSOR_SIZENWSE) } # # set the size & system icon # self.SetClientSize(self.GetSize()) self.iconsize = (16, 16) self.SetIcon(wx.Icon(os.path.join(globalvar.ETCICONDIR, 'grass_map.ico'), wx.BITMAP_TYPE_ICO)) # # Fancy gui # self._mgr = wx.aui.AuiManager(self) # # Add toolbars # self.toolbars = { 'map' : None, 'vdigit' : None, 'georect' : None, 'gcpdisp' : None, 'nviz' : None } for toolb in toolbars: self.AddToolbar(toolb) # # Add statusbar # self.statusbar = self.CreateStatusBar(number = 4, style = 0) self.statusbar.SetStatusWidths([-5, -2, -1, -1]) self.statusbarWin = dict() self.statusbarWin['toggle'] = wx.Choice(self.statusbar, wx.ID_ANY, choices = globalvar.MAP_DISPLAY_STATUSBAR_MODE) self.statusbarWin['toggle'].SetSelection(UserSettings.Get(group = 'display', key = 'statusbarMode', subkey = 'selection')) self.statusbar.Bind(wx.EVT_CHOICE, self.OnToggleStatus, self.statusbarWin['toggle']) # auto-rendering checkbox self.statusbarWin['render'] = wx.CheckBox(parent = self.statusbar, id = wx.ID_ANY, label = _("Render")) self.statusbar.Bind(wx.EVT_CHECKBOX, self.OnToggleRender, self.statusbarWin['render']) self.statusbarWin['render'].SetValue(UserSettings.Get(group = 'display', key = 'autoRendering', subkey = 'enabled')) self.statusbarWin['render'].SetToolTip(wx.ToolTip (_("Enable/disable auto-rendering"))) # show region self.statusbarWin['region'] = wx.CheckBox(parent = self.statusbar, id = wx.ID_ANY, label = _("Show computational extent")) self.statusbar.Bind(wx.EVT_CHECKBOX, self.OnToggleShowRegion, self.statusbarWin['region']) self.statusbarWin['region'].SetValue(False) self.statusbarWin['region'].Hide() self.statusbarWin['region'].SetToolTip(wx.ToolTip (_("Show/hide computational " "region extent (set with g.region). " "Display region drawn as a blue box inside the " "computational region, " "computational region inside a display region " "as a red box)."))) # set resolution self.statusbarWin['resolution'] = wx.CheckBox(parent = self.statusbar, id = wx.ID_ANY, label = _("Constrain display resolution to computational settings")) self.statusbar.Bind(wx.EVT_CHECKBOX, self.OnToggleResolution, self.statusbarWin['resolution']) self.statusbarWin['resolution'].SetValue(UserSettings.Get(group = 'display', key = 'compResolution', subkey = 'enabled')) self.statusbarWin['resolution'].Hide() self.statusbarWin['resolution'].SetToolTip(wx.ToolTip (_("Constrain display resolution " "to computational region settings. " "Default value for new map displays can " "be set up in 'User GUI settings' dialog."))) # map scale self.statusbarWin['mapscale'] = wx.ComboBox(parent = self.statusbar, id = wx.ID_ANY, style = wx.TE_PROCESS_ENTER, size = (150, -1)) self.statusbarWin['mapscale'].SetItems(['1:1000', '1:5000', '1:10000', '1:25000', '1:50000', '1:100000', '1:1000000']) self.statusbarWin['mapscale'].Hide() self.statusbar.Bind(wx.EVT_TEXT_ENTER, self.OnChangeMapScale, self.statusbarWin['mapscale']) self.statusbar.Bind(wx.EVT_COMBOBOX, self.OnChangeMapScale, self.statusbarWin['mapscale']) # go to self.statusbarWin['goto'] = wx.TextCtrl(parent = self.statusbar, id = wx.ID_ANY, value = "", style = wx.TE_PROCESS_ENTER, size = (300, -1)) self.statusbarWin['goto'].Hide() self.statusbar.Bind(wx.EVT_TEXT_ENTER, self.OnGoTo, self.statusbarWin['goto']) # projection self.statusbarWin['projection'] = wx.CheckBox(parent = self.statusbar, id = wx.ID_ANY, label = _("Use defined projection")) self.statusbarWin['projection'].SetValue(False) size = self.statusbarWin['projection'].GetSize() self.statusbarWin['projection'].SetMinSize((size[0] + 150, size[1])) self.statusbarWin['projection'].SetToolTip(wx.ToolTip (_("Reproject coordinates displayed " "in the statusbar. Projection can be " "defined in GUI preferences dialog " "(tab 'Display')"))) self.statusbarWin['projection'].Hide() # mask self.statusbarWin['mask'] = wx.StaticText(parent = self.statusbar, id = wx.ID_ANY, label = '') self.statusbarWin['mask'].SetForegroundColour(wx.Colour(255, 0, 0)) # on-render gauge self.statusbarWin['progress'] = wx.Gauge(parent = self.statusbar, id = wx.ID_ANY, range = 0, style = wx.GA_HORIZONTAL) self.statusbarWin['progress'].Hide() self.StatusbarReposition() # reposition statusbar # # Init map display (buffered DC & set default cursor) # self.MapWindow2D = BufferedWindow(self, id = wx.ID_ANY, Map = self.Map, tree = self.tree, lmgr = self._layerManager) # default is 2D display mode self.MapWindow = self.MapWindow2D self.MapWindow.SetCursor(self.cursors["default"]) # used by vector digitizer self.MapWindowVDigit = None # used by Nviz (3D display mode) self.MapWindow3D = None # # initialize region values # self._initDisplay() # # Bind various events # self.Bind(wx.EVT_ACTIVATE, self.OnFocus) self.Bind(wx.EVT_CLOSE, self.OnCloseWindow) self.Bind(render.EVT_UPDATE_PRGBAR, self.OnUpdateProgress) # # Update fancy gui style # self._mgr.AddPane(self.MapWindow, wx.aui.AuiPaneInfo().CentrePane(). Dockable(False).BestSize((-1,-1)). CloseButton(False).DestroyOnClose(True). Layer(0)) self._mgr.Update() # # Init print module and classes # self.printopt = disp_print.PrintOptions(self, self.MapWindow) # # Init zoom history # self.MapWindow.ZoomHistory(self.Map.region['n'], self.Map.region['s'], self.Map.region['e'], self.Map.region['w']) # # Re-use dialogs # self.dialogs = {} self.dialogs['attributes'] = None self.dialogs['category'] = None self.dialogs['barscale'] = None self.dialogs['legend'] = None self.decorationDialog = None # decoration/overlays def _addToolbarVDigit(self): """!Add vector digitizer toolbar """ from vdigit import haveVDigit if not haveVDigit: from vdigit import errorMsg msg = _("Unable to start wxGUI vector digitizer.\nDo you want to start " "TCL/TK digitizer (v.digit) instead?\n\n" "Details: %s" % errorMsg) self.toolbars['map'].combo.SetValue(_("2D view")) dlg = wx.MessageDialog(parent = self, message = msg, caption=_("Vector digitizer failed"), style = wx.YES_NO | wx.CENTRE) if dlg.ShowModal() == wx.ID_YES: mapName = self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].GetName() self._layerManager.goutput.RunCmd(['v.digit', 'map=%s' % mapName], switchPage = False) dlg.Destroy() self.toolbars['map'].combo.SetValue(_("2D view")) return if self._layerManager: log = self._layerManager.goutput else: log = None if not self.MapWindowVDigit: from mapdisp_vdigit import VDigitWindow self.MapWindowVDigit = VDigitWindow(self, id = wx.ID_ANY, Map = self.Map, tree = self.tree, lmgr = self._layerManager) self.MapWindowVDigit.Show() self.MapWindow = self.MapWindowVDigit self._mgr.DetachPane(self.MapWindow2D) self.MapWindow2D.Hide() self.toolbars['vdigit'] = toolbars.VDigitToolbar(parent = self, mapcontent = self.Map, layerTree = self.tree, log = log) self.MapWindowVDigit.SetToolbar(self.toolbars['vdigit']) self._mgr.AddPane(self.MapWindowVDigit, wx.aui.AuiPaneInfo().CentrePane(). Dockable(False).BestSize((-1,-1)). CloseButton(False).DestroyOnClose(True). Layer(0)) self._mgr.AddPane(self.toolbars['vdigit'], wx.aui.AuiPaneInfo(). Name("vdigittoolbar").Caption(_("Vector Digitizer Toolbar")). ToolbarPane().Top().Row(1). LeftDockable(False).RightDockable(False). BottomDockable(False).TopDockable(True). CloseButton(False).Layer(2). BestSize((self.toolbars['vdigit'].GetBestSize()))) # change mouse to draw digitized line self.MapWindow.mouse['box'] = "point" self.MapWindow.zoomtype = 0 self.MapWindow.pen = wx.Pen(colour = 'red', width = 2, style = wx.SOLID) self.MapWindow.polypen = wx.Pen(colour = 'green', width = 2, style = wx.SOLID) def _addToolbarNviz(self): """!Add 3D view mode toolbar """ import nviz # check for GLCanvas and OpenGL if not nviz.haveNviz: self.toolbars['map'].combo.SetValue (_("2D view")) gcmd.GError(parent = self, message = _("Unable to switch to 3D display mode.\nThe Nviz python extension " "was not found or loaded properly.\n" "Switching back to 2D display mode.\n\nDetails: %s" % nviz.errorMsg)) return # add Nviz toolbar and disable 2D display mode tools self.toolbars['nviz'] = toolbars.NvizToolbar(self, self.Map) self.toolbars['map'].Enable2D(False) # update status bar self.statusbarWin['toggle'].Enable(False) # erase map window self.MapWindow.EraseMap() self._layerManager.goutput.WriteCmdLog(_("Starting 3D view mode...")) self.statusbar.SetStatusText(_("Please wait, loading data..."), 0) # create GL window & NVIZ toolbar if not self.MapWindow3D: self.MapWindow3D = nviz.GLWindow(self, id = wx.ID_ANY, Map = self.Map, tree = self.tree, lmgr = self._layerManager) self.MapWindow = self.MapWindow3D self.MapWindow.SetCursor(self.cursors["default"]) # add Nviz notebookpage self._layerManager.AddNviz() self.MapWindow3D.OnPaint(None) # -> LoadData self.MapWindow3D.Show() self.MapWindow3D.UpdateView(None) else: self.MapWindow = self.MapWindow3D # add Nviz notebookpage self._layerManager.AddNviz() self._layerManager.nviz.UpdatePage('view') self._layerManager.nviz.UpdatePage('light') # switch from MapWindow to MapWindowGL # add nviz toolbar self._mgr.DetachPane(self.MapWindow2D) self.MapWindow2D.Hide() self._mgr.AddPane(self.MapWindow3D, wx.aui.AuiPaneInfo().CentrePane(). Dockable(False).BestSize((-1,-1)). CloseButton(False).DestroyOnClose(True). Layer(0)) self._mgr.AddPane(self.toolbars['nviz'], wx.aui.AuiPaneInfo(). Name("nviztoolbar").Caption(_("3D View Toolbar")). ToolbarPane().Top().Row(1). LeftDockable(False).RightDockable(False). BottomDockable(False).TopDockable(True). CloseButton(False).Layer(2). BestSize((self.toolbars['nviz'].GetBestSize()))) self.SetStatusText("", 0) def AddToolbar(self, name): """!Add defined toolbar to the window Currently known toolbars are: - 'map' - basic map toolbar - 'vdigit' - vector digitizer - 'gcpdisp' - GCP Manager Display - 'georect' - georectifier - 'nviz' - 3D view mode """ # default toolbar if name == "map": self.toolbars['map'] = toolbars.MapToolbar(self, self.Map) self._mgr.AddPane(self.toolbars['map'], wx.aui.AuiPaneInfo(). Name("maptoolbar").Caption(_("Map Toolbar")). ToolbarPane().Top(). LeftDockable(False).RightDockable(False). BottomDockable(False).TopDockable(True). CloseButton(False).Layer(2). BestSize((self.toolbars['map'].GetBestSize()))) # vector digitizer elif name == "vdigit": self._addToolbarVDigit() # georectifier elif name == "georect": self.toolbars['georect'] = toolbars.GRToolbar(self, self.Map) self._mgr.AddPane(self.toolbars['georect'], wx.aui.AuiPaneInfo(). Name("georecttoolbar").Caption(_("Georectification Toolbar")). ToolbarPane().Top(). LeftDockable(False).RightDockable(False). BottomDockable(False).TopDockable(True). CloseButton(False).Layer(2). BestSize((self.toolbars['georect'].GetBestSize()))) # nviz elif name == "nviz": self._addToolbarNviz() self._mgr.Update() def RemoveToolbar (self, name): """!Removes defined toolbar from the window @todo Only hide, activate by calling AddToolbar() """ # cannot hide main toolbar if name == "map": return self._mgr.DetachPane(self.toolbars[name]) self.toolbars[name].Destroy() self.toolbars[name] = None if name == 'vdigit': self._mgr.DetachPane(self.MapWindowVDigit) self.MapWindowVDigit.Hide() self.MapWindow2D.Show() self._mgr.AddPane(self.MapWindow2D, wx.aui.AuiPaneInfo().CentrePane(). Dockable(False).BestSize((-1,-1)). CloseButton(False).DestroyOnClose(True). Layer(0)) self.MapWindow = self.MapWindow2D elif name == 'nviz': # unload data # self.MapWindow3D.Reset() # switch from MapWindowGL to MapWindow self._mgr.DetachPane(self.MapWindow3D) self.MapWindow3D.Hide() self.MapWindow2D.Show() self._mgr.AddPane(self.MapWindow2D, wx.aui.AuiPaneInfo().CentrePane(). Dockable(False).BestSize((-1,-1)). CloseButton(False).DestroyOnClose(True). Layer(0)) self.MapWindow = self.MapWindow2D # remove nviz notebook page self._layerManager.RemoveNviz() self.MapWindow.UpdateMap() self.toolbars['map'].combo.SetValue(_("2D view")) self.toolbars['map'].Enable2D(True) self.statusbarWin['toggle'].Enable(True) self._mgr.Update() def _initDisplay(self): """!Initialize map display, set dimensions and map region """ if not grass.find_program('g.region', ['--help']): sys.exit(_("GRASS module '%s' not found. Unable to start map " "display window.") % 'g.region') self.width, self.height = self.GetClientSize() Debug.msg(2, "MapFrame._initDisplay():") self.Map.ChangeMapSize(self.GetClientSize()) self.Map.region = self.Map.GetRegion() # g.region -upgc # self.Map.SetRegion() # adjust region to match display window def OnUpdateProgress(self, event): """!Update progress bar info """ self.statusbarWin['progress'].SetValue(event.value) event.Skip() def OnFocus(self, event): """ Change choicebook page to match display. Or set display for georectifying """ if self._layerManager and \ self._layerManager.georectifying: # in georectifying session; display used to get geographic # coordinates for GCPs self.OnPointer(event) else: # change bookcontrol page to page associated with display if self.page: pgnum = self.layerbook.GetPageIndex(self.page) if pgnum > -1: self.layerbook.SetSelection(pgnum) event.Skip() def OnDraw(self, event): """!Re-display current map composition """ self.MapWindow.UpdateMap(render = False) def OnRender(self, event): """!Re-render map composition (each map layer) """ # delete tmp map layers (queries) qlayer = self.Map.GetListOfLayers(l_name = globalvar.QUERYLAYER) for layer in qlayer: self.Map.DeleteLayer(layer) # delete tmp lines if self.MapWindow.mouse["use"] in ("measure", "profile"): self.MapWindow.polycoords = [] self.MapWindow.ClearLines() # deselect features in vdigit if self.toolbars['vdigit']: if self.MapWindow.digit: self.MapWindow.digit.GetDisplay().SetSelected([]) self.MapWindow.UpdateMap(render = True, renderVector = True) else: self.MapWindow.UpdateMap(render = True) # update statusbar self.StatusbarUpdate() def OnPointer(self, event): """!Pointer button clicked """ if self.toolbars['map']: if event: self.toolbars['map'].OnTool(event) self.toolbars['map'].action['desc'] = '' self.MapWindow.mouse['use'] = "pointer" self.MapWindow.mouse['box'] = "point" # change the cursor if self.toolbars['vdigit']: # digitization tool activated self.MapWindow.SetCursor(self.cursors["cross"]) # reset mouse['box'] if needed if self.toolbars['vdigit'].GetAction() in ['addLine']: if self.toolbars['vdigit'].GetAction('type') in ['point', 'centroid']: self.MapWindow.mouse['box'] = 'point' else: # line, boundary self.MapWindow.mouse['box'] = 'line' elif self.toolbars['vdigit'].GetAction() in ['addVertex', 'removeVertex', 'splitLine', 'editLine', 'displayCats', 'queryMap', 'copyCats']: self.MapWindow.mouse['box'] = 'point' else: # moveLine, deleteLine self.MapWindow.mouse['box'] = 'box' elif self._layerManager and self._layerManager.georectifying: self.MapWindow.SetCursor(self.cursors["cross"]) else: self.MapWindow.SetCursor(self.cursors["default"]) def OnZoomIn(self, event): """ Zoom in the map. Set mouse cursor, zoombox attributes, and zoom direction """ if self.toolbars['map']: self.toolbars['map'].OnTool(event) self.toolbars['map'].action['desc'] = '' self.MapWindow.mouse['use'] = "zoom" self.MapWindow.mouse['box'] = "box" self.MapWindow.zoomtype = 1 self.MapWindow.pen = wx.Pen(colour = 'Red', width = 2, style = wx.SHORT_DASH) # change the cursor self.MapWindow.SetCursor(self.cursors["cross"]) def OnZoomOut(self, event): """ Zoom out the map. Set mouse cursor, zoombox attributes, and zoom direction """ if self.toolbars['map']: self.toolbars['map'].OnTool(event) self.toolbars['map'].action['desc'] = '' self.MapWindow.mouse['use'] = "zoom" self.MapWindow.mouse['box'] = "box" self.MapWindow.zoomtype = -1 self.MapWindow.pen = wx.Pen(colour = 'Red', width = 2, style = wx.SHORT_DASH) # change the cursor self.MapWindow.SetCursor(self.cursors["cross"]) def OnZoomBack(self, event): """ Zoom last (previously stored position) """ self.MapWindow.ZoomBack() def OnPan(self, event): """ Panning, set mouse to drag """ if self.toolbars['map']: self.toolbars['map'].OnTool(event) self.toolbars['map'].action['desc'] = '' self.MapWindow.mouse['use'] = "pan" self.MapWindow.mouse['box'] = "pan" self.MapWindow.zoomtype = 0 # change the cursor self.MapWindow.SetCursor(self.cursors["hand"]) def OnErase(self, event): """ Erase the canvas """ self.MapWindow.EraseMap() def OnZoomRegion(self, event): """ Zoom to region """ self.Map.getRegion() self.Map.getResolution() self.UpdateMap() # event.Skip() def OnAlignRegion(self, event): """ Align region """ if not self.Map.alignRegion: self.Map.alignRegion = True else: self.Map.alignRegion = False # event.Skip() def OnToggleRender(self, event): """!Enable/disable auto-rendering """ if self.statusbarWin['render'].GetValue(): self.OnRender(None) def IsAutoRendered(self): """!Check if auto-rendering is enabled""" return self.statusbarWin['render'].IsChecked() def OnToggleShowRegion(self, event): """!Show/Hide extent in map canvas """ if self.statusbarWin['region'].GetValue(): # show extent self.MapWindow.regionCoords = [] else: del self.MapWindow.regionCoords # redraw map if auto-rendering is enabled if self.statusbarWin['render'].GetValue(): self.OnRender(None) def OnToggleResolution(self, event): """ Use resolution of computation region settings for redering image instead of display resolution """ # redraw map if auto-rendering is enabled if self.statusbarWin['render'].GetValue(): self.OnRender(None) def OnToggleStatus(self, event): """ Toggle status text """ self.StatusbarUpdate() def OnChangeMapScale(self, event): """ Map scale changed by user """ scale = event.GetString() try: if scale[:2] != '1:': raise ValueError value = int(scale[2:]) except ValueError: self.statusbarWin['mapscale'].SetValue('1:%ld' % int(self.mapScaleValue)) return dEW = value * (self.Map.region['cols'] / self.ppm[0]) dNS = value * (self.Map.region['rows'] / self.ppm[1]) self.Map.region['n'] = self.Map.region['center_northing'] + dNS / 2. self.Map.region['s'] = self.Map.region['center_northing'] - dNS / 2. self.Map.region['w'] = self.Map.region['center_easting'] - dEW / 2. self.Map.region['e'] = self.Map.region['center_easting'] + dEW / 2. # add to zoom history self.MapWindow.ZoomHistory(self.Map.region['n'], self.Map.region['s'], self.Map.region['e'], self.Map.region['w']) # redraw a map self.MapWindow.UpdateMap() self.statusbarWin['mapscale'].SetFocus() def OnGoTo(self, event): """ Go to position """ try: if self.statusbarWin['projection'].IsChecked(): if not UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'proj4'): self.statusbar.SetStatusText(_("Projection not defined (check the settings)"), 0) else: # reproject values projIn = UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'proj4') projOut = gcmd.RunCommand('g.proj', flags = 'jf', read = True) proj = projIn.split(' ')[0].split('=')[1] if proj in ('ll', 'latlong', 'longlat'): e, n = self.statusbarWin['goto'].GetValue().split(';') e, n = utils.DMS2Deg(e, n) proj, coord1 = utils.ReprojectCoordinates(coord = (e, n), projIn = projIn, projOut = projOut, flags = 'd') e, n = coord1 else: e, n = map(float, self.statusbarWin['goto'].GetValue().split(';')) proj, coord1 = utils.ReprojectCoordinates(coord = (e, n), projIn = projIn, projOut = projOut, flags = 'd') e, n = coord1 else: if self.Map.projinfo['proj'] == 'll': e, n = self.statusbarWin['goto'].GetValue().split(';') else: e, n = map(float, self.statusbarWin['goto'].GetValue().split(';')) region = self.Map.GetCurrentRegion() if self.statusbarWin['projection'].IsChecked(): if not UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'proj4'): self.statusbar.SetStatusText(_("Projection not defined (check the settings)"), 0) else: region['center_easting'], region['center_northing'] = e, n else: if self.Map.projinfo['proj'] == 'll': region['center_easting'], region['center_northing'] = utils.DMS2Deg(e, n) else: region['center_easting'], region['center_northing'] = e, n except ValueError: region = self.Map.GetCurrentRegion() precision = int(UserSettings.Get(group = 'projection', key = 'format', subkey = 'precision')) format = UserSettings.Get(group = 'projection', key = 'format', subkey = 'll') if self.Map.projinfo['proj'] == 'll' and format == 'DMS': self.statusbarWin['goto'].SetValue("%s" % utils.Deg2DMS(region['center_easting'], region['center_northing'], precision = precision)) else: self.statusbarWin['goto'].SetValue("%.*f; %.*f" % \ (precision, region['center_easting'], precision, region['center_northing'])) return dn = (region['nsres'] * region['rows']) / 2. region['n'] = region['center_northing'] + dn region['s'] = region['center_northing'] - dn de = (region['ewres'] * region['cols']) / 2. region['e'] = region['center_easting'] + de region['w'] = region['center_easting'] - de self.Map.AdjustRegion() # add to zoom history self.MapWindow.ZoomHistory(region['n'], region['s'], region['e'], region['w']) # redraw a map self.MapWindow.UpdateMap() self.statusbarWin['goto'].SetFocus() def StatusbarUpdate(self): """!Update statusbar content""" self.statusbarWin['region'].Hide() self.statusbarWin['resolution'].Hide() self.statusbarWin['mapscale'].Hide() self.statusbarWin['goto'].Hide() self.statusbarWin['projection'].Hide() self.mapScaleValue = self.ppm = None if self.statusbarWin['toggle'].GetSelection() == 0: # Coordinates self.statusbar.SetStatusText("", 0) # enable long help self.StatusbarEnableLongHelp() elif self.statusbarWin['toggle'].GetSelection() in (1, 2): # Extent sel = self.statusbarWin['toggle'].GetSelection() if sel == 1: region = self.Map.region else: region = self.Map.GetRegion() # computation region precision = int(UserSettings.Get(group = 'projection', key = 'format', subkey = 'precision')) format = UserSettings.Get(group = 'projection', key = 'format', subkey = 'll') if self.statusbarWin['projection'].IsChecked(): if not UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'proj4'): self.statusbar.SetStatusText(_("Projection not defined (check the settings)"), 0) else: projOut = UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'proj4') proj, coord1 = utils.ReprojectCoordinates(coord = (region["w"], region["s"]), projOut = projOut, flags = 'd') proj, coord2 = utils.ReprojectCoordinates(coord = (region["e"], region["n"]), projOut = projOut, flags = 'd') if sel == 2: proj, coord3 = utils.ReprojectCoordinates(coord = (0.0, 0.0), projOut = projOut, flags = 'd') proj, coord4 = utils.ReprojectCoordinates(coord = (region["ewres"], region["nsres"]), projOut = projOut, flags = 'd') if coord1 and coord2: if proj in ('ll', 'latlong', 'longlat') and format == 'DMS': w, s = utils.Deg2DMS(coord1[0], coord1[1], string = False, precision = precision) e, n = utils.Deg2DMS(coord2[0], coord2[1], string = False, precision = precision) if sel == 1: self.statusbar.SetStatusText("%s - %s, %s - %s" % (w, e, s, n), 0) else: ewres, nsres = utils.Deg2DMS(abs(coord3[0]) - abs(coord4[0]), abs(coord3[1]) - abs(coord4[1]), string = False, hemisphere = False, precision = precision) self.statusbar.SetStatusText("%s - %s, %s - %s (%s, %s)" % (w, e, s, n, ewres, nsres), 0) else: w, s = coord1 e, n = coord2 if sel == 1: self.statusbar.SetStatusText("%.*f - %.*f, %.*f - %.*f" % (precision, w, precision, e, precision, s, precision, n), 0) else: ewres, nsres = coord3 self.statusbar.SetStatusText("%.*f - %.*f, %.*f - %.*f (%.*f, %.*f)" % (precision, w, precision, e, precision, s, precision, n, precision, ewres, precision, nsres), 0) else: self.statusbar.SetStatusText(_("Error in projection (check the settings)"), 0) else: if self.Map.projinfo['proj'] == 'll' and format == 'DMS': w, s = utils.Deg2DMS(region["w"], region["s"], string = False, precision = precision) e, n = utils.Deg2DMS(region["e"], region["n"], string = False, precision = precision) if sel == 1: self.statusbar.SetStatusText("%s - %s, %s - %s" % (w, e, s, n), 0) else: ewres, nsres = utils.Deg2DMS(region['ewres'], region['nsres'], string = False, precision = precision) self.statusbar.SetStatusText("%s - %s, %s - %s (%s, %s)" % (w, e, s, n, ewres, nsres), 0) else: w, s = region["w"], region["s"] e, n = region["e"], region["n"] if sel == 1: self.statusbar.SetStatusText("%.*f - %.*f, %.*f - %.*f" % (precision, w, precision, e, precision, s, precision, n), 0) else: ewres, nsres = region['ewres'], region['nsres'] self.statusbar.SetStatusText("%.*f - %.*f, %.*f - %.*f (%.*f, %.*f)" % (precision, w, precision, e, precision, s, precision, n, precision, ewres, precision, nsres), 0) # enable long help self.StatusbarEnableLongHelp() elif self.statusbarWin['toggle'].GetSelection() == 3: # Show comp. extent self.statusbar.SetStatusText("", 0) self.statusbarWin['region'].Show() # disable long help self.StatusbarEnableLongHelp(False) elif self.statusbarWin['toggle'].GetSelection() == 4: # Display mode self.statusbar.SetStatusText("", 0) self.statusbarWin['resolution'].Show() # disable long help self.StatusbarEnableLongHelp(False) elif self.statusbarWin['toggle'].GetSelection() == 5: # Display geometry self.statusbar.SetStatusText("rows=%d; cols=%d; nsres=%.2f; ewres=%.2f" % (self.Map.region["rows"], self.Map.region["cols"], self.Map.region["nsres"], self.Map.region["ewres"]), 0) # enable long help self.StatusbarEnableLongHelp() elif self.statusbarWin['toggle'].GetSelection() == 6: # Map scale # TODO: need to be fixed... ### screen X region problem ### user should specify ppm dc = wx.ScreenDC() dpSizePx = wx.DisplaySize() # display size in pixels dpSizeMM = wx.DisplaySizeMM() # display size in mm (system) dpSizeIn = (dpSizeMM[0] / 25.4, dpSizeMM[1] / 25.4) # inches sysPpi = dc.GetPPI() comPpi = (dpSizePx[0] / dpSizeIn[0], dpSizePx[1] / dpSizeIn[1]) ppi = comPpi # pixel per inch self.ppm = ((ppi[0] / 2.54) * 100, # pixel per meter (ppi[1] / 2.54) * 100) Debug.msg(4, "MapFrame.StatusbarUpdate(mapscale): size: px=%d,%d mm=%f,%f " "in=%f,%f ppi: sys=%d,%d com=%d,%d; ppm=%f,%f" % \ (dpSizePx[0], dpSizePx[1], dpSizeMM[0], dpSizeMM[1], dpSizeIn[0], dpSizeIn[1], sysPpi[0], sysPpi[1], comPpi[0], comPpi[1], self.ppm[0], self.ppm[1])) region = self.Map.region heightCm = region['rows'] / self.ppm[1] * 100 widthCm = region['cols'] / self.ppm[0] * 100 Debug.msg(4, "MapFrame.StatusbarUpdate(mapscale): width_cm=%f, height_cm=%f" % (widthCm, heightCm)) xscale = (region['e'] - region['w']) / (region['cols'] / self.ppm[0]) yscale = (region['n'] - region['s']) / (region['rows'] / self.ppm[1]) scale = (xscale + yscale) / 2. Debug.msg(3, "MapFrame.StatusbarUpdate(mapscale): xscale=%f, yscale=%f -> scale=%f" % \ (xscale, yscale, scale)) self.statusbar.SetStatusText("") try: self.statusbarWin['mapscale'].SetValue("1:%ld" % (scale + 0.5)) except TypeError: pass self.mapScaleValue = scale self.statusbarWin['mapscale'].Show() # disable long help self.StatusbarEnableLongHelp(False) elif self.statusbarWin['toggle'].GetSelection() == 7: # go to self.statusbar.SetStatusText("") region = self.Map.GetCurrentRegion() precision = int(UserSettings.Get(group = 'projection', key = 'format', subkey = 'precision')) format = UserSettings.Get(group = 'projection', key = 'format', subkey = 'll') if self.statusbarWin['projection'].IsChecked(): if not UserSettings.Get(group='projection', key='statusbar', subkey='proj4'): self.statusbar.SetStatusText(_("Projection not defined (check the settings)"), 0) else: proj, coord = utils.ReprojectCoordinates(coord = (region['center_easting'], region['center_northing']), projOut = UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'proj4'), flags = 'd') if coord: if proj in ('ll', 'latlong', 'longlat') and format == 'DMS': self.statusbarWin['goto'].SetValue("%s" % utils.Deg2DMS(coord[0], coord[1], precision = precision)) else: self.statusbarWin['goto'].SetValue("%.*f; %.*f" % (precision, coord[0], precision, coord[1])) else: self.statusbar.SetStatusText(_("Error in projection (check the settings)"), 0) else: if self.Map.projinfo['proj'] == 'll' and format == 'DMS': self.statusbarWin['goto'].SetValue("%s" % utils.Deg2DMS(region['center_easting'], region['center_northing'], precision = precision)) else: self.statusbarWin['goto'].SetValue("%.*f; %.*f" % (precision, region['center_easting'], precision, region['center_northing'])) self.statusbarWin['goto'].Show() # disable long help self.StatusbarEnableLongHelp(False) elif self.statusbarWin['toggle'].GetSelection() == 8: # projection self.statusbar.SetStatusText("") epsg = UserSettings.Get(group = 'projection', key = 'statusbar', subkey = 'epsg') if epsg: label = '%s (EPSG: %s)' % (_("Use defined projection"), epsg) self.statusbarWin['projection'].SetLabel(label) else: self.statusbarWin['projection'].SetLabel(_("Use defined projection")) self.statusbarWin['projection'].Show() # disable long help self.StatusbarEnableLongHelp(False) else: self.statusbar.SetStatusText("", 1) def StatusbarEnableLongHelp(self, enable = True): """!Enable/disable toolbars long help""" for toolbar in self.toolbars.itervalues(): if toolbar: toolbar.EnableLongHelp(enable) def StatusbarReposition(self): """!Reposition checkbox in statusbar""" # reposition checkbox widgets = [(0, self.statusbarWin['region']), (0, self.statusbarWin['resolution']), (0, self.statusbarWin['mapscale']), (0, self.statusbarWin['progress']), (0, self.statusbarWin['projection']), (1, self.statusbarWin['toggle']), (2, self.statusbarWin['mask']), (3, self.statusbarWin['render'])] for idx, win in widgets: rect = self.statusbar.GetFieldRect(idx) if idx == 0: # show region / mapscale / process bar # -> size wWin, hWin = win.GetBestSize() if win == self.statusbarWin['progress']: wWin = rect.width - 6 # -> position # if win == self.statusbarWin['region']: # x, y = rect.x + rect.width - wWin, rect.y - 1 # align left # else: x, y = rect.x + 3, rect.y - 1 w, h = wWin, rect.height + 2 else: # choice || auto-rendering x, y = rect.x, rect.y - 1 w, h = rect.width, rect.height + 2 if idx == 2: # mask x += 5 y += 4 elif idx == 3: # render x += 5 win.SetPosition((x, y)) win.SetSize((w, h)) def SaveToFile(self, event): """!Save map to image """ if self.toolbars['nviz']: filetype = "PPM file (*.ppm)|*.ppm|TIF file (*.tif)|*.tif" ltype = [{ 'ext' : 'ppm', 'type' : -1 }, { 'ext' : 'tif', 'type' : wx.BITMAP_TYPE_TIF }] else: img = self.MapWindow.img if not img: gcmd.GMessage(parent = self, message = _("Nothing to render (empty map). Operation canceled.")) return filetype, ltype = gdialogs.GetImageHandlers(img) # get size dlg = gdialogs.ImageSizeDialog(self) dlg.CentreOnParent() if dlg.ShowModal() != wx.ID_OK: dlg.Destroy() return width, height = dlg.GetValues() dlg.Destroy() # get filename dlg = wx.FileDialog(parent = self, message = _("Choose a file name to save the image " "(no need to add extension)"), wildcard = filetype, style = wx.SAVE | wx.FD_OVERWRITE_PROMPT) if dlg.ShowModal() == wx.ID_OK: path = dlg.GetPath() if not path: dlg.Destroy() return base, ext = os.path.splitext(path) fileType = ltype[dlg.GetFilterIndex()]['type'] extType = ltype[dlg.GetFilterIndex()]['ext'] if ext != extType: path = base + '.' + extType self.MapWindow.SaveToFile(path, fileType, width, height) dlg.Destroy() def PrintMenu(self, event): """ Print options and output menu for map display """ point = wx.GetMousePosition() printmenu = wx.Menu() # Add items to the menu setup = wx.MenuItem(printmenu, wx.ID_ANY, _('Page setup')) printmenu.AppendItem(setup) self.Bind(wx.EVT_MENU, self.printopt.OnPageSetup, setup) preview = wx.MenuItem(printmenu, wx.ID_ANY, _('Print preview')) printmenu.AppendItem(preview) self.Bind(wx.EVT_MENU, self.printopt.OnPrintPreview, preview) doprint = wx.MenuItem(printmenu, wx.ID_ANY, _('Print display')) printmenu.AppendItem(doprint) self.Bind(wx.EVT_MENU, self.printopt.OnDoPrint, doprint) # Popup the menu. If an item is selected then its handler # will be called before PopupMenu returns. self.PopupMenu(printmenu) printmenu.Destroy() def OnCloseWindow(self, event): """!Window closed. Also close associated layer tree page """ pgnum = None self.Map.Clean() # close edited map and 3D tools properly if self.toolbars['vdigit']: maplayer = self.toolbars['vdigit'].GetLayer() if maplayer: self.toolbars['vdigit'].OnExit() if self.toolbars['nviz']: self.toolbars['nviz'].OnExit() if not self._layerManager: self.Destroy() elif self.page: pgnum = self.layerbook.GetPageIndex(self.page) if pgnum > -1: self.layerbook.DeletePage(pgnum) def GetRender(self): """!Returns current instance of render.Map() """ return self.Map def GetWindow(self): """!Get map window""" return self.MapWindow def OnNvizQuerySurface(self, event): """!Query current surface in 3D view mode""" if self.toolbars['map'].GetAction() == 'nvizQuerySurface': self.toolbars['map'].SelectDefault(event) return self.toolbars['map'].action['desc'] = 'nvizQuerySurface' self.MapWindow.mouse['use'] = "nvizQuerySurface" self._OnQuery() def OnNvizQueryVector(self, event): """!Query current vector in 3D view mode""" if self.toolbars['map'].GetAction() == 'nvizQueryVector': self.toolbars['map'].SelectDefault(event) return self.toolbars['map'].action['desc'] = 'nvizQueryVector' self.MapWindow.mouse['use'] = "nvizQueryVector" self._OnQuery() def QueryMap(self, x, y): """!Query raster or vector map layers by r/v.what @param x,y coordinates """ # set query snap distance for v.what at map unit equivalent of 10 pixels qdist = 10.0 * ((self.Map.region['e'] - self.Map.region['w']) / self.Map.width) east, north = self.MapWindow.Pixel2Cell((x, y)) if not self.IsStandalone(): num = 0 for layer in self.tree.GetSelections(): ltype = self.tree.GetPyData(layer)[0]['maplayer'].GetType() if ltype in ('raster', 'rgb', 'his', 'vector', 'thememap', 'themechart'): num += 1 if num < 1: gcmd.GMessage(parent = self, message = _('No raster or vector map layer selected for querying.')) return rast = list() vect = list() rcmd = ['r.what', '--v'] vcmd = ['v.what', '--v'] if self.IsStandalone(): pass else: for layer in self.tree.GetSelections(): ltype = self.tree.GetPyData(layer)[0]['maplayer'].GetType() dcmd = self.tree.GetPyData(layer)[0]['cmd'] name, found = utils.GetLayerNameFromCmd(dcmd) if not found: continue if ltype == 'raster': rast.append(name) elif ltype in ('rgb', 'his'): for iname in name.split('\n'): rast.append(iname) elif ltype in ('vector', 'thememap', 'themechart'): vect.append(name) # use display region settings instead of computation region settings self.tmpreg = os.getenv("GRASS_REGION") os.environ["GRASS_REGION"] = self.Map.SetRegion(windres = False) # build query commands for any selected rasters and vectors if rast: rcmd.append('-f') rcmd.append('-n') rcmd.append('input=%s' % ','.join(rast)) rcmd.append('east_north=%f,%f' % (float(east), float(north))) if vect: # check for vector maps open to be edited digitToolbar = self.toolbars['vdigit'] if digitToolbar: lmap = digitToolbar.GetLayer().GetName() for name in vect: if lmap == name: self._layerManager.goutput.WriteWarning(_("Vector map <%s> " "opened for editing - skipped.") % map) vect.remove(name) if len(vect) < 1: self._layerManager.goutput.WriteCmdLog(_("Nothing to query.")) return vcmd.append('-a') vcmd.append('map=%s' % ','.join(vect)) vcmd.append('east_north=%f,%f' % (float(east), float(north))) vcmd.append('distance=%f' % float(qdist)) Debug.msg(1, "QueryMap(): raster=%s vector=%s" % (','.join(rast), ','.join(vect))) # parse query command(s) if not self.IsStandalone(): if rast: self._layerManager.goutput.RunCmd(rcmd, compReg = False, onDone = self._QueryMapDone) if vect: self._layerManager.goutput.RunCmd(vcmd, onDone = self._QueryMapDone) else: if rast: gcmd.RunCommand(rcmd) if vect: gcmd.RunCommand(vcmd) def _QueryMapDone(self, cmd, returncode): """!Restore settings after querying (restore GRASS_REGION) @param returncode command return code """ if hasattr(self, "tmpreg"): if self.tmpreg: os.environ["GRASS_REGION"] = self.tmpreg elif 'GRASS_REGION' in os.environ: del os.environ["GRASS_REGION"] elif 'GRASS_REGION' in os.environ: del os.environ["GRASS_REGION"] if hasattr(self, "tmpreg"): del self.tmpreg def QueryVector(self, x, y): """!Query vector map layer features Attribute data of selected vector object are displayed in GUI dialog. Data can be modified (On Submit) """ if not self.tree.layer_selected or \ self.tree.GetPyData(self.tree.layer_selected)[0]['type'] != 'vector': gcmd.GMessage(parent = self, message = _("No map layer selected for querying.")) return posWindow = self.ClientToScreen((x + self.MapWindow.dialogOffset, y + self.MapWindow.dialogOffset)) qdist = 10.0 * ((self.Map.region['e'] - self.Map.region['w']) / self.Map.width) east, north = self.MapWindow.Pixel2Cell((x, y)) mapName = self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].name if self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].GetMapset() != \ grass.gisenv()['MAPSET']: mode = 'display' else: mode = 'update' if self.dialogs['attributes'] is None: dlg = dbm_dialogs.DisplayAttributesDialog(parent = self.MapWindow, map = mapName, query = ((east, north), qdist), pos = posWindow, action = mode) self.dialogs['attributes'] = dlg else: # selection changed? if not self.dialogs['attributes'].mapDBInfo or \ self.dialogs['attributes'].mapDBInfo.map != mapName: self.dialogs['attributes'].UpdateDialog(map = mapName, query = ((east, north), qdist), action = mode) else: self.dialogs['attributes'].UpdateDialog(query = ((east, north), qdist), action = mode) if not self.dialogs['attributes'].IsFound(): self._layerManager.goutput.WriteLog(_('Nothing found.')) cats = self.dialogs['attributes'].GetCats() try: qlayer = self.Map.GetListOfLayers(l_name = globalvar.QUERYLAYER)[0] except IndexError: qlayer = None if self.dialogs['attributes'].mapDBInfo and cats: # highlight feature & re-draw map if qlayer: qlayer.SetCmd(self.AddTmpVectorMapLayer(mapName, cats, useId = False, addLayer = False)) else: qlayer = self.AddTmpVectorMapLayer(mapName, cats, useId = False) # set opacity based on queried layer opacity = self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].GetOpacity(float = True) qlayer.SetOpacity(opacity) self.MapWindow.UpdateMap(render = False, renderVector = False) if not self.dialogs['attributes'].IsShown(): self.dialogs['attributes'].Show() else: if qlayer: self.Map.DeleteLayer(qlayer) self.MapWindow.UpdateMap(render = False, renderVector = False) if self.dialogs['attributes'].IsShown(): self.dialogs['attributes'].Hide() def OnQuery(self, event): """!Query tools menu""" if self.toolbars['map']: self.toolbars['map'].OnTool(event) action = self.toolbars['map'].GetAction() if self.toolbars['nviz']: toolsmenu = wx.Menu() raster = wx.MenuItem(parentMenu = toolsmenu, id = wx.ID_ANY, text = _("Query surface (raster map)"), kind = wx.ITEM_CHECK) toolsmenu.AppendItem(raster) self.Bind(wx.EVT_MENU, self.OnNvizQuerySurface, raster) if action == "nvizQuerySurface": raster.Check(True) vector = wx.MenuItem(parentMenu = toolsmenu, id = wx.ID_ANY, text = _("Query vector map"), kind = wx.ITEM_CHECK) toolsmenu.AppendItem(vector) self.Bind(wx.EVT_MENU, self.OnNvizQueryVector, vector) if action == "nvizQueryVector": vector.Check(True) self.PopupMenu(toolsmenu) toolsmenu.Destroy() else: self.toolbars['map'].action['desc'] = 'queryMap' self.MapWindow.mouse['use'] = "query" if not self.IsStandalone(): # switch to output console to show query results self._layerManager.notebook.SetSelectionByName('output') self.MapWindow.mouse['box'] = "point" self.MapWindow.zoomtype = 0 # change the cursor self.MapWindow.SetCursor(self.cursors["cross"]) def AddTmpVectorMapLayer(self, name, cats, useId = False, addLayer = True): """!Add temporal vector map layer to map composition @param name name of map layer @param useId use feature id instead of category """ # color settings from ATM color = UserSettings.Get(group = 'atm', key = 'highlight', subkey = 'color') colorStr = str(color[0]) + ":" + \ str(color[1]) + ":" + \ str(color[2]) # icon used in vector display and its size icon = '' size = 0 vparam = self.tree.GetPyData(self.tree.layer_selected)[0]['cmd'] for p in vparam: if '=' in p: parg,pval = p.split('=') if parg == 'icon': icon = pval elif parg == 'size': size = int(pval) pattern = ["d.vect", "map=%s" % name, "color=%s" % colorStr, "fcolor=%s" % colorStr, "width=%d" % UserSettings.Get(group = 'atm', key = 'highlight', subkey = 'width')] if icon != '': pattern.append('icon=%s' % icon) if size > 0: pattern.append('size=%i' % size) if useId: cmd = pattern cmd.append('-i') cmd.append('cats=%s' % str(cats)) else: cmd = [] for layer in cats.keys(): cmd.append(copy.copy(pattern)) lcats = cats[layer] cmd[-1].append("layer=%d" % layer) cmd[-1].append("cats=%s" % utils.ListOfCatsToRange(lcats)) if addLayer: if useId: return self.Map.AddLayer(type = 'vector', name = globalvar.QUERYLAYER, command = cmd, l_active = True, l_hidden = True, l_opacity = 1.0) else: return self.Map.AddLayer(type = 'command', name = globalvar.QUERYLAYER, command = cmd, l_active = True, l_hidden = True, l_opacity = 1.0) else: return cmd def OnAnalyze(self, event): """!Analysis tools menu """ point = wx.GetMousePosition() toolsmenu = wx.Menu() icons = Icons['displayWindow'] # Add items to the menu measure = wx.MenuItem(toolsmenu, wx.ID_ANY, icons["measure"].GetLabel()) measure.SetBitmap(icons["measure"].GetBitmap(self.iconsize)) toolsmenu.AppendItem(measure) self.Bind(wx.EVT_MENU, self.OnMeasure, measure) profile = wx.MenuItem(toolsmenu, wx.ID_ANY, icons["profile"].GetLabel()) profile.SetBitmap(icons["profile"].GetBitmap(self.iconsize)) toolsmenu.AppendItem(profile) self.Bind(wx.EVT_MENU, self.Profile, profile) histogram = wx.MenuItem(toolsmenu, wx.ID_ANY, icons["histogram"].GetLabel()) histogram.SetBitmap(icons["histogram"].GetBitmap(self.iconsize)) toolsmenu.AppendItem(histogram) self.Bind(wx.EVT_MENU, self.Histogram, histogram) # Popup the menu. If an item is selected then its handler # will be called before PopupMenu returns. self.PopupMenu(toolsmenu) toolsmenu.Destroy() def OnMeasure(self, event): """!Init measurement routine that calculates map distance along transect drawn on map display """ self.totaldist = 0.0 # total measured distance # switch Layer Manager to output console to show measure results self._layerManager.notebook.SetSelectionByName('output') # change mouse to draw line for measurement self.MapWindow.mouse['use'] = "measure" self.MapWindow.mouse['box'] = "line" self.MapWindow.zoomtype = 0 self.MapWindow.pen = wx.Pen(colour = 'red', width = 2, style = wx.SHORT_DASH) self.MapWindow.polypen = wx.Pen(colour = 'green', width = 2, style = wx.SHORT_DASH) # change the cursor self.MapWindow.SetCursor(self.cursors["pencil"]) # initiating output style = self._layerManager.goutput.cmd_output.StyleWarning self._layerManager.goutput.WriteLog(_('Click and drag with left mouse button ' 'to measure.%s' 'Double click with left button to clear.') % \ (os.linesep), style) if self.Map.projinfo['proj'] != 'xy': units = self.Map.projinfo['units'] self._layerManager.goutput.WriteCmdLog(_('Measuring distance') + ' (' + units + '):') else: self._layerManager.goutput.WriteCmdLog(_('Measuring distance:')) if self.Map.projinfo['proj'] == 'll': try: import grass.lib.gis as gislib global haveCtypes haveCtypes = True gislib.G_begin_distance_calculations() except ImportError, e: self._layerManager.goutput.WriteWarning(_('Geodesic distance is not yet ' 'supported by this tool.\n' 'Reason: %s' % e)) def MeasureDist(self, beginpt, endpt): """!Calculate map distance from screen distance and print to output window """ self._layerManager.notebook.SetSelectionByName('output') dist, (north, east) = self.MapWindow.Distance(beginpt, endpt) dist = round(dist, 3) d, dunits = self.FormatDist(dist) self.totaldist += dist td, tdunits = self.FormatDist(self.totaldist) strdist = str(d) strtotdist = str(td) if self.Map.projinfo['proj'] == 'xy' or 'degree' not in self.Map.projinfo['unit']: angle = int(math.degrees(math.atan2(north,east)) + 0.5) angle = 180 - angle if angle < 0: angle = 360 + angle mstring = '%s = %s %s\n%s = %s %s\n%s = %d %s\n%s' \ % (_('segment'), strdist, dunits, _('total distance'), strtotdist, tdunits, _('bearing'), angle, _('deg'), '-' * 60) else: mstring = '%s = %s %s\n%s = %s %s\n%s' \ % (_('segment'), strdist, dunits, _('total distance'), strtotdist, tdunits, '-' * 60) self._layerManager.goutput.WriteLog(mstring) return dist def Profile(self, event): """!Init profile canvas and tools """ raster = [] if self.tree.layer_selected and \ self.tree.GetPyData(self.tree.layer_selected)[0]['type'] == 'raster': raster.append(self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].name) self.profile = profile.ProfileFrame(self, id = wx.ID_ANY, pos = wx.DefaultPosition, size = (700,300), style = wx.DEFAULT_FRAME_STYLE, rasterList = raster) self.profile.Show() # Open raster select dialog to make sure that a raster (and the desired raster) # is selected to be profiled self.profile.OnSelectRaster(None) def FormatDist(self, dist): """!Format length numbers and units in a nice way, as a function of length. From code by Hamish Bowman Grass Development Team 2006""" mapunits = self.Map.projinfo['units'] if mapunits == 'metres': mapunits = 'meters' outunits = mapunits dist = float(dist) divisor = 1.0 # figure out which units to use if mapunits == 'meters': if dist > 2500.0: outunits = 'km' divisor = 1000.0 else: outunits = 'm' elif mapunits == 'feet': # nano-bug: we match any "feet", but US Survey feet is really # 5279.9894 per statute mile, or 10.6' per 1000 miles. As >1000 # miles the tick markers are rounded to the nearest 10th of a # mile (528'), the difference in foot flavours is ignored. if dist > 5280.0: outunits = 'miles' divisor = 5280.0 else: outunits = 'ft' elif 'degree' in mapunits and \ not haveCtypes: if dist < 1: outunits = 'min' divisor = (1/60.0) else: outunits = 'deg' else: outunits = 'meters' # format numbers in a nice way if (dist/divisor) >= 2500.0: outdist = round(dist/divisor) elif (dist/divisor) >= 1000.0: outdist = round(dist/divisor,1) elif (dist/divisor) > 0.0: outdist = round(dist/divisor,int(math.ceil(3-math.log10(dist/divisor)))) else: outdist = float(dist/divisor) return (outdist, outunits) def Histogram(self, event): """!Init histogram display canvas and tools """ self.histogram = histogram.HistFrame(self, id = wx.ID_ANY, size = globalvar.HIST_WINDOW_SIZE, style = wx.DEFAULT_FRAME_STYLE) #show new display self.histogram.Show() self.histogram.Refresh() self.histogram.Update() def OnDecoration(self, event): """!Decorations overlay menu """ point = wx.GetMousePosition() decmenu = wx.Menu() icons = Icons['displayWindow'] # Add items to the menu AddScale = wx.MenuItem(decmenu, wx.ID_ANY, icons["addBarscale"].GetLabel()) AddScale.SetBitmap(icons["addBarscale"].GetBitmap(self.iconsize)) decmenu.AppendItem(AddScale) self.Bind(wx.EVT_MENU, self.OnAddBarscale, AddScale) AddLegend = wx.MenuItem(decmenu, wx.ID_ANY, icons["addLegend"].GetLabel()) AddLegend.SetBitmap(icons["addLegend"].GetBitmap(self.iconsize)) decmenu.AppendItem(AddLegend) self.Bind(wx.EVT_MENU, self.OnAddLegend, AddLegend) AddText = wx.MenuItem(decmenu, wx.ID_ANY, icons["addText"].GetLabel()) AddText.SetBitmap(icons["addText"].GetBitmap(self.iconsize)) decmenu.AppendItem(AddText) self.Bind(wx.EVT_MENU, self.OnAddText, AddText) # Popup the menu. If an item is selected then its handler # will be called before PopupMenu returns. self.PopupMenu(decmenu) decmenu.Destroy() def OnAddBarscale(self, event): """!Handler for scale/arrow map decoration menu selection. """ if self.dialogs['barscale']: return id = 0 # unique index for overlay layer # If location is latlon, only display north arrow (scale won't work) # proj = self.Map.projinfo['proj'] # if proj == 'll': # barcmd = 'd.barscale -n' # else: # barcmd = 'd.barscale' # decoration overlay control dialog self.dialogs['barscale'] = \ gdialogs.DecorationDialog(parent = self, title = _('Scale and North arrow'), size = (350, 200), style = wx.DEFAULT_DIALOG_STYLE | wx.CENTRE, cmd = ['d.barscale', 'at=0,5'], ovlId = id, name = 'barscale', checktxt = _("Show/hide scale and North arrow"), ctrltxt = _("scale object")) self.dialogs['barscale'].CentreOnParent() ### dialog cannot be show as modal - in the result d.barscale is not selectable ### self.dialogs['barscale'].ShowModal() self.dialogs['barscale'].Show() self.MapWindow.mouse['use'] = 'pointer' def OnAddLegend(self, event): """!Handler for legend map decoration menu selection. """ if self.dialogs['legend']: return id = 1 # index for overlay layer in render cmd = ['d.legend', 'at=5,50,2,5'] if self.tree.layer_selected and \ self.tree.GetPyData(self.tree.layer_selected)[0]['type'] == 'raster': cmd.append('map=%s' % self.tree.GetPyData(self.tree.layer_selected)[0]['maplayer'].name) # Decoration overlay control dialog self.dialogs['legend'] = \ gdialogs.DecorationDialog(parent = self, title = ('Legend'), size = (350, 200), style = wx.DEFAULT_DIALOG_STYLE | wx.CENTRE, cmd = cmd, ovlId = id, name = 'legend', checktxt = _("Show/hide legend"), ctrltxt = _("legend object")) self.dialogs['legend'].CentreOnParent() ### dialog cannot be show as modal - in the result d.legend is not selectable ### self.dialogs['legend'].ShowModal() self.dialogs['legend'].Show() self.MapWindow.mouse['use'] = 'pointer' def OnAddText(self, event): """!Handler for text decoration menu selection. """ if self.MapWindow.dragid > -1: id = self.MapWindow.dragid else: # index for overlay layer in render if len(self.MapWindow.textdict.keys()) > 0: id = self.MapWindow.textdict.keys()[-1] + 1 else: id = 101 self.dialogs['text'] = gdialogs.TextLayerDialog(parent = self, ovlId = id, title = _('Add text layer'), size = (400, 200)) self.dialogs['text'].CenterOnParent() # If OK button pressed in decoration control dialog if self.dialogs['text'].ShowModal() == wx.ID_OK: text = self.dialogs['text'].GetValues()['text'] active = self.dialogs['text'].GetValues()['active'] coords, w, h = self.MapWindow.TextBounds(self.dialogs['text'].GetValues()) # delete object if it has no text or is not active if text == '' or active == False: try: self.MapWindow.pdc.ClearId(id) self.MapWindow.pdc.RemoveId(id) del self.MapWindow.textdict[id] except: pass return self.MapWindow.pdc.ClearId(id) self.MapWindow.pdc.SetId(id) self.MapWindow.textdict[id] = self.dialogs['text'].GetValues() self.MapWindow.Draw(self.MapWindow.pdcDec, img = self.MapWindow.textdict[id], drawid = id, pdctype = 'text', coords = coords) self.MapWindow.UpdateMap(render = False, renderVector = False) self.MapWindow.mouse['use'] = 'pointer' def GetOptData(self, dcmd, type, params, propwin): """!Callback method for decoration overlay command generated by dialog created in menuform.py """ # Reset comand and rendering options in render.Map. Always render decoration. # Showing/hiding handled by PseudoDC self.Map.ChangeOverlay(ovltype = type, type = 'overlay', name = '', command = dcmd, l_active = True, l_render = False) self.params[type] = params self.propwin[type] = propwin def OnZoomToMap(self, event): """!Set display extents to match selected raster (including NULLs) or vector map. """ self.MapWindow.ZoomToMap() def OnZoomToRaster(self, event): """!Set display extents to match selected raster map (ignore NULLs) """ self.MapWindow.ZoomToMap(ignoreNulls = True) def OnZoomToWind(self, event): """!Set display geometry to match computational region settings (set with g.region) """ self.MapWindow.ZoomToWind() def OnZoomToDefault(self, event): """!Set display geometry to match default region settings """ self.MapWindow.ZoomToDefault() def OnZoomToSaved(self, event): """!Set display geometry to match extents in saved region file """ self.MapWindow.ZoomToSaved() def OnDisplayToWind(self, event): """!Set computational region (WIND file) to match display extents """ self.MapWindow.DisplayToWind() def SaveDisplayRegion(self, event): """!Save display extents to named region file. """ self.MapWindow.SaveDisplayRegion() def OnZoomMenu(self, event): """!Popup Zoom menu """ point = wx.GetMousePosition() zoommenu = wx.Menu() # Add items to the menu zoomwind = wx.MenuItem(zoommenu, wx.ID_ANY, _('Zoom to computational region (set with g.region)')) zoommenu.AppendItem(zoomwind) self.Bind(wx.EVT_MENU, self.OnZoomToWind, zoomwind) zoomdefault = wx.MenuItem(zoommenu, wx.ID_ANY, _('Zoom to default region')) zoommenu.AppendItem(zoomdefault) self.Bind(wx.EVT_MENU, self.OnZoomToDefault, zoomdefault) zoomsaved = wx.MenuItem(zoommenu, wx.ID_ANY, _('Zoom to saved region')) zoommenu.AppendItem(zoomsaved) self.Bind(wx.EVT_MENU, self.OnZoomToSaved, zoomsaved) savewind = wx.MenuItem(zoommenu, wx.ID_ANY, _('Set computational region from display extent')) zoommenu.AppendItem(savewind) self.Bind(wx.EVT_MENU, self.OnDisplayToWind, savewind) savezoom = wx.MenuItem(zoommenu, wx.ID_ANY, _('Save display geometry to named region')) zoommenu.AppendItem(savezoom) self.Bind(wx.EVT_MENU, self.SaveDisplayRegion, savezoom) # Popup the menu. If an item is selected then its handler # will be called before PopupMenu returns. self.PopupMenu(zoommenu) zoommenu.Destroy() def SetProperties(self, render = False, mode = 0, showCompExtent = False, constrainRes = False, projection = False): """!Set properies of map display window""" self.statusbarWin['render'].SetValue(render) self.statusbarWin['toggle'].SetSelection(mode) self.StatusbarUpdate() self.statusbarWin['region'].SetValue(showCompExtent) self.statusbarWin['resolution'].SetValue(constrainRes) self.statusbarWin['projection'].SetValue(projection) if showCompExtent: self.MapWindow.regionCoords = [] def IsStandalone(self): """!Check if Map display is standalone""" if self._layerManager: return False return True def GetLayerManager(self): """!Get reference to Layer Manager @return window reference @return None (if standalone) """ return self._layerManager # end of class MapFrame class MapApp(wx.App): def OnInit(self): wx.InitAllImageHandlers() if __name__ == "__main__": Map = render.Map() # instance of Map class to render GRASS display output to PPM file else: Map = None self.mapFrm = MapFrame(parent = None, id = wx.ID_ANY, Map = Map, size = globalvar.MAP_WINDOW_SIZE) #self.SetTopWindow(Map) self.mapFrm.Show() if __name__ == "__main__": # redraw map, if new command appears self.redraw = False status = Command(self, Map, cmdfilename) status.start() self.timer = wx.PyTimer(self.watcher) # check each 0.1s self.timer.Start(100) return 1 def OnExit(self): if __name__ == "__main__": # stop the timer self.timer.Stop() # terminate thread (a bit ugly) os.system("""!echo "quit" >> %s""" % (cmdfilename)) def watcher(self): """!Redraw, if new layer appears""" if self.redraw: self.mapFrm.OnDraw(None) self.redraw = False return # end of class MapApp if __name__ == "__main__": ###### SET command variable if len(sys.argv) != 3: print __doc__ sys.exit() title = sys.argv[1] cmdfilename = sys.argv[2] import gettext gettext.install('grasswxpy', os.path.join(os.getenv("GISBASE"), 'locale'), unicode = True) print >> sys.stderr, "\nStarting monitor <%s>...\n" % (title) gm_map = MapApp(0) # set title gm_map.mapFrm.SetTitle(_("GRASS GIS Map Display: " + title + " - Location: " + grass.gisenv()["LOCATION_NAME"])) gm_map.MainLoop() if grass.gisenv().has_key("MONITOR"): os.system("d.mon sel=%s" % grass.gisenv()["MONITOR"]) os.remove(cmdfilename) os.system("""!g.gisenv set="GRASS_PYCMDFILE" """) print >> sys.stderr, "\nStoping monitor <%s>...\n" % (title) sys.exit(0)