nviz_mapdisp.py 62 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398139914001401140214031404140514061407140814091410141114121413141414151416141714181419142014211422142314241425142614271428142914301431143214331434143514361437143814391440144114421443144414451446144714481449145014511452145314541455145614571458145914601461146214631464146514661467146814691470147114721473147414751476147714781479148014811482148314841485148614871488148914901491149214931494149514961497149814991500150115021503150415051506150715081509151015111512151315141515151615171518151915201521152215231524152515261527152815291530153115321533153415351536153715381539154015411542154315441545154615471548154915501551155215531554155515561557155815591560156115621563156415651566156715681569157015711572157315741575157615771578157915801581158215831584158515861587
  1. """!
  2. @package nviz_mapdisp.py
  3. @brief wxGUI 3D view mode (map canvas)
  4. This module implements 3D visualization mode for map display.
  5. List of classes:
  6. - NvizThread
  7. - GLWindow
  8. (C) 2008-2011 by the GRASS Development Team
  9. This program is free software under the GNU General Public
  10. License (>=v2). Read the file COPYING that comes with GRASS
  11. for details.
  12. @author Martin Landa <landa.martin gmail.com> (Google SoC 2008/2010)
  13. @author Anna Kratochvilova <KratochAnna seznam.cz> (Google SoC 2011)
  14. """
  15. import os
  16. import sys
  17. import time
  18. import copy
  19. import math
  20. from threading import Thread
  21. import wx
  22. import wx.lib.scrolledpanel as scrolled
  23. from wx.lib.newevent import NewEvent
  24. from wx import glcanvas
  25. import gcmd
  26. import globalvar
  27. from debug import Debug
  28. from mapdisp_window import MapWindow
  29. from goutput import wxCmdOutput
  30. from preferences import globalSettings as UserSettings
  31. from workspace import Nviz as NvizDefault
  32. import wxnviz
  33. wxUpdateProperties, EVT_UPDATE_PROP = NewEvent()
  34. wxUpdateView, EVT_UPDATE_VIEW = NewEvent()
  35. wxUpdateLight, EVT_UPDATE_LIGHT = NewEvent()
  36. wxUpdateCPlane, EVT_UPDATE_CPLANE = NewEvent()
  37. class NvizThread(Thread):
  38. def __init__(self, log, progressbar, window):
  39. Thread.__init__(self)
  40. Debug.msg(5, "NvizThread.__init__():")
  41. self.log = log
  42. self.progressbar = progressbar
  43. self.window = window
  44. self._display = None
  45. self.setDaemon(True)
  46. def run(self):
  47. self._display = wxnviz.Nviz(self.log, self.progressbar)
  48. def GetDisplay(self):
  49. """!Get display instance"""
  50. return self._display
  51. class GLWindow(MapWindow, glcanvas.GLCanvas):
  52. """!OpenGL canvas for Map Display Window"""
  53. def __init__(self, parent, id = wx.ID_ANY,
  54. Map = None, tree = None, lmgr = None):
  55. self.parent = parent # MapFrame
  56. Debug.msg(5, "GLCanvas.__init__(): begin")
  57. glcanvas.GLCanvas.__init__(self, parent, id)
  58. MapWindow.__init__(self, parent, id,
  59. Map, tree, lmgr)
  60. self.Hide()
  61. self.init = False
  62. self.initView = False
  63. # render mode
  64. self.render = { 'quick' : False,
  65. # do not render vector lines in quick mode
  66. 'vlines' : False,
  67. 'vpoints' : False }
  68. self.mouse = {
  69. 'use': 'default'
  70. }
  71. self.cursors = {
  72. 'default' : wx.StockCursor(wx.CURSOR_ARROW),
  73. 'cross' : wx.StockCursor(wx.CURSOR_CROSS),
  74. }
  75. # list of loaded map layers (layer tree items)
  76. self.layers = list()
  77. # list of constant surfaces
  78. self.constants = list()
  79. # list of cutting planes
  80. self.cplanes = list()
  81. # list of query points
  82. self.qpoints = list()
  83. #
  84. # use display region instead of computational
  85. #
  86. os.environ['GRASS_REGION'] = self.Map.SetRegion()
  87. #
  88. # create nviz instance
  89. #
  90. if self.lmgr:
  91. self.log = self.lmgr.goutput
  92. logerr = self.lmgr.goutput.cmd_stderr
  93. logmsg = self.lmgr.goutput.cmd_output
  94. else:
  95. self.log = logmsg = sys.stdout
  96. logerr = sys.stderr
  97. self.nvizThread = NvizThread(logerr,
  98. self.parent.statusbarWin['progress'],
  99. logmsg)
  100. self.nvizThread.start()
  101. time.sleep(.1)
  102. self._display = self.nvizThread.GetDisplay()
  103. # GRASS_REGION needed only for initialization
  104. del os.environ['GRASS_REGION']
  105. self.img = wx.Image(self.Map.mapfile, wx.BITMAP_TYPE_ANY)
  106. # size of MapWindow, to avoid resizing if size is the same
  107. self.size = (0,0)
  108. #
  109. # default values
  110. #
  111. self.view = copy.deepcopy(UserSettings.Get(group = 'nviz', key = 'view')) # copy
  112. self.iview = UserSettings.Get(group = 'nviz', key = 'view', internal = True)
  113. self.nvizDefault = NvizDefault()
  114. self.light = copy.deepcopy(UserSettings.Get(group = 'nviz', key = 'light')) # copy
  115. self.Bind(wx.EVT_ERASE_BACKGROUND, self.OnEraseBackground)
  116. self.Bind(wx.EVT_SIZE, self.OnSize)
  117. self.Bind(wx.EVT_PAINT, self.OnPaint)
  118. self.Bind(wx.EVT_LEFT_UP, self.OnLeftUp)
  119. self.Bind(wx.EVT_MOUSE_EVENTS, self.OnMouseAction)
  120. self.Bind(wx.EVT_MOTION, self.OnMotion)
  121. self.Bind(EVT_UPDATE_PROP, self.UpdateMapObjProperties)
  122. self.Bind(EVT_UPDATE_VIEW, self.UpdateView)
  123. self.Bind(EVT_UPDATE_LIGHT, self.UpdateLight)
  124. self.Bind(EVT_UPDATE_CPLANE, self.UpdateCPlane)
  125. self.Bind(wx.EVT_CLOSE, self.OnClose)
  126. Debug.msg(5, "GLCanvas.__init__(): end")
  127. #cplanes cannot be initialized now
  128. wx.CallAfter(self.InitCPlanes)
  129. def InitCPlanes(self):
  130. """!Initialize cutting planes list"""
  131. for i in range(self._display.GetCPlanesCount()):
  132. cplane = copy.deepcopy(UserSettings.Get(group = 'nviz', key = 'cplane'))
  133. self.cplanes.append(cplane)
  134. def OnClose(self, event):
  135. # cleanup when window actually closes (on quit) and not just is hidden
  136. self.Reset()
  137. def OnEraseBackground(self, event):
  138. pass # do nothing, to avoid flashing on MSW
  139. def OnSize(self, event):
  140. size = self.GetClientSize()
  141. if self.size != size \
  142. and self.GetContext():
  143. Debug.msg(3, "GLCanvas.OnSize(): w = %d, h = %d" % \
  144. (size.width, size.height))
  145. self.SetCurrent()
  146. self._display.ResizeWindow(size.width,
  147. size.height)
  148. self.size = size
  149. event.Skip()
  150. def OnPaint(self, event):
  151. Debug.msg(1, "GLCanvas.OnPaint()")
  152. dc = wx.PaintDC(self)
  153. self.DoPaint()
  154. def DoPaint(self):
  155. self.SetCurrent()
  156. if not self.initView:
  157. self._display.InitView()
  158. self.initView = True
  159. self.LoadDataLayers()
  160. self.UnloadDataLayers()
  161. if not self.init:
  162. self.ResetView()
  163. if hasattr(self.lmgr, "nviz"):
  164. self.lmgr.nviz.UpdatePage('view')
  165. self.lmgr.nviz.UpdatePage('light')
  166. self.lmgr.nviz.UpdatePage('cplane')
  167. layer = self.GetSelectedLayer()
  168. if layer:
  169. if layer.type == 'raster':
  170. self.lmgr.nviz.UpdatePage('surface')
  171. self.lmgr.nviz.UpdatePage('fringe')
  172. elif layer.type == 'vector':
  173. self.lmgr.nviz.UpdatePage('vector')
  174. self.lmgr.nviz.UpdateSettings()
  175. # update widgets
  176. win = self.lmgr.nviz.FindWindowById( \
  177. self.lmgr.nviz.win['vector']['lines']['surface'])
  178. win.SetItems(self.GetLayerNames('raster'))
  179. self.init = True
  180. self.UpdateMap()
  181. def OnMouseAction(self, event):
  182. """!Handle mouse events"""
  183. # zoom with mouse wheel
  184. if event.GetWheelRotation() != 0:
  185. self.OnMouseWheel(event)
  186. # left mouse button pressed
  187. elif event.LeftDown():
  188. self.OnLeftDown(event)
  189. # left mouse button released
  190. elif event.LeftUp():
  191. self.OnLeftUp(event)
  192. # dragging
  193. elif event.Dragging():
  194. self.OnDragging(event)
  195. event.Skip()
  196. def OnMouseWheel(self, event):
  197. """!Change perspective"""
  198. wheel = event.GetWheelRotation()
  199. current = event.GetPositionTuple()[:]
  200. Debug.msg (5, "GLWindow.OnMouseMotion(): wheel = %d" % wheel)
  201. prev_value = self.view['persp']['value']
  202. if wheel > 0:
  203. value = -1 * self.view['persp']['step']
  204. else:
  205. value = self.view['persp']['step']
  206. self.view['persp']['value'] += value
  207. if self.view['persp']['value'] < 1:
  208. self.view['persp']['value'] = 1
  209. elif self.view['persp']['value'] > 100:
  210. self.view['persp']['value'] = 100
  211. if prev_value != self.view['persp']['value']:
  212. if hasattr(self.lmgr, "nviz"):
  213. self.lmgr.nviz.UpdateSettings()
  214. self._display.SetView(self.view['position']['x'], self.view['position']['y'],
  215. self.iview['height']['value'],
  216. self.view['persp']['value'],
  217. self.view['twist']['value'])
  218. # redraw map
  219. self.DoPaint()
  220. # update statusbar
  221. ### self.parent.StatusbarUpdate()
  222. def OnLeftDown(self, event):
  223. """!On left mouse down"""
  224. if self.mouse['use'] == "lookHere":
  225. pos = event.GetPosition()
  226. size = self.GetClientSize()
  227. self._display.LookHere(pos[0], size[1] - pos[1])
  228. self.Refresh(False)
  229. focus = self._display.GetFocus()
  230. for i, coord in enumerate(('x', 'y', 'z')):
  231. self.iview['focus'][coord] = focus[i]
  232. toggle = self.lmgr.nviz.FindWindowByName('here')
  233. toggle.SetValue(False)
  234. self.mouse['use'] = 'default'
  235. self.SetCursor(self.cursors['default'])
  236. if self.mouse['use'] == "cplane":
  237. pos = event.GetPosition()
  238. size = self.GetClientSize()
  239. x, y, z = self._display.SetCPlaneInteractively(pos[0], size[1] - pos[1])
  240. if x is not None:
  241. idx = self._display.GetCPlaneCurrent()
  242. self.cplanes[idx]['position']['x'] = x
  243. self.cplanes[idx]['position']['y'] = y
  244. self.render['quick'] = True
  245. event.Skip()
  246. def OnDragging(self, event):
  247. if self.mouse['use'] == "cplane":
  248. idx = self._display.GetCPlaneCurrent()
  249. pos = event.GetPosition()
  250. size = self.GetClientSize()
  251. x, y, z = self._display.SetCPlaneInteractively(pos[0], size[1] - pos[1])
  252. if x is not None:
  253. self.cplanes[idx]['position']['x'] = x
  254. self.cplanes[idx]['position']['y'] = y
  255. event.Skip()
  256. def Pixel2Cell(self, (x, y)):
  257. """!Convert image coordinates to real word coordinates
  258. @param x, y image coordinates
  259. @return easting, northing
  260. @return None on error
  261. """
  262. size = self.GetClientSize()
  263. # UL -> LL
  264. sid, x, y, z = self._display.GetPointOnSurface(x, size[1] - y)
  265. if not sid:
  266. return None
  267. return (x, y)
  268. def OnLeftUp(self, event):
  269. self.ReleaseMouse()
  270. if self.mouse["use"] == "nvizQuerySurface":
  271. self.OnQuerySurface(event)
  272. elif self.mouse["use"] == "nvizQueryVector":
  273. self.OnQueryVector(event)
  274. elif self.mouse["use"] == 'cplane':
  275. self.lmgr.nviz.OnCPlaneChangeDone(None)
  276. idx = self._display.GetCPlaneCurrent()
  277. self.lmgr.nviz.UpdateCPlanePage(idx)
  278. self.lmgr.nviz.FindWindowByName('cplaneHere').SetValue(False)
  279. self.mouse['use'] = 'default'
  280. self.SetCursor(self.cursors['default'])
  281. def OnQuerySurface(self, event):
  282. """!Query surface on given position"""
  283. result = self._display.QueryMap(event.GetX(), event.GetY())
  284. if result:
  285. self.qpoints.append((result['x'], result['y'], result['z']))
  286. self.log.WriteLog("%-30s: %.3f" % (_("Easting"), result['x']))
  287. self.log.WriteLog("%-30s: %.3f" % (_("Northing"), result['y']))
  288. self.log.WriteLog("%-30s: %.3f" % (_("Elevation"), result['z']))
  289. self.log.WriteLog("%-30s: %s" % (_("Surface map elevation"), result['elevation']))
  290. self.log.WriteLog("%-30s: %s" % (_("Surface map color"), result['color']))
  291. if len(self.qpoints) > 1:
  292. prev = self.qpoints[-2]
  293. curr = self.qpoints[-1]
  294. dxy = math.sqrt(pow(prev[0]-curr[0], 2) +
  295. pow(prev[1]-curr[1], 2))
  296. dxyz = math.sqrt(pow(prev[0]-curr[0], 2) +
  297. pow(prev[1]-curr[1], 2) +
  298. pow(prev[2]-curr[2], 2))
  299. self.log.WriteLog("%-30s: %.3f" % (_("XY distance from previous"), dxy))
  300. self.log.WriteLog("%-30s: %.3f" % (_("XYZ distance from previous"), dxyz))
  301. self.log.WriteLog("%-30s: %.3f" % (_("Distance along surface"),
  302. self._display.GetDistanceAlongSurface(result['id'],
  303. (curr[0], curr[1]),
  304. (prev[0], prev[1]),
  305. useExag = False)))
  306. self.log.WriteLog("%-30s: %.3f" % (_("Distance along exag. surface"),
  307. self._display.GetDistanceAlongSurface(result['id'],
  308. (curr[0], curr[1]),
  309. (prev[0], prev[1]),
  310. useExag = True)))
  311. self.log.WriteCmdLog('-' * 80)
  312. else:
  313. self.log.WriteLog(_("No point on surface"))
  314. self.log.WriteCmdLog('-' * 80)
  315. def OnQueryVector(self, event):
  316. """!Query vector on given position"""
  317. self.log.WriteWarning(_("Function not implemented yet"))
  318. self.log.WriteCmdLog('-' * 80)
  319. def UpdateView(self, event):
  320. """!Change view settings"""
  321. data = self.view
  322. self._display.SetView(data['position']['x'], data['position']['y'],
  323. self.iview['height']['value'],
  324. data['persp']['value'],
  325. data['twist']['value'])
  326. if event and event.zExag and 'value' in data['z-exag']:
  327. self._display.SetZExag(data['z-exag']['value'])
  328. if self.iview['focus']['x'] != -1:
  329. self._display.SetFocus(self.iview['focus']['x'], self.iview['focus']['y'],
  330. self.iview['focus']['z'])
  331. if event:
  332. event.Skip()
  333. def UpdateLight(self, event):
  334. """!Change light settings"""
  335. data = self.light
  336. self._display.SetLight(x = data['position']['x'], y = data['position']['y'],
  337. z = data['position']['z'] / 100., color = data['color'],
  338. bright = data['bright'] / 100.,
  339. ambient = data['ambient'] / 100.)
  340. self._display.DrawLightingModel()
  341. if hasattr(event, 'refresh'):
  342. self.Refresh(False)
  343. def UpdateMap(self, render = True):
  344. """!Updates the canvas anytime there is a change to the
  345. underlaying images or to the geometry of the canvas.
  346. @param render re-render map composition
  347. """
  348. start = time.clock()
  349. self.resize = False
  350. if self.render['quick'] is False:
  351. self.parent.statusbarWin['progress'].Show()
  352. self.parent.statusbarWin['progress'].SetRange(2)
  353. self.parent.statusbarWin['progress'].SetValue(0)
  354. if self.render['quick'] is False:
  355. self.parent.statusbarWin['progress'].SetValue(1)
  356. self._display.Draw(False, -1)
  357. elif self.render['quick'] is True:
  358. # quick
  359. mode = wxnviz.DRAW_QUICK_SURFACE | wxnviz.DRAW_QUICK_VOLUME
  360. if self.render['vlines']:
  361. mode |= wxnviz.DRAW_QUICK_VLINES
  362. if self.render['vpoints']:
  363. mode |= wxnviz.DRAW_QUICK_VPOINTS
  364. self._display.Draw(True, mode)
  365. else: # None -> reuse last rendered image
  366. pass # TODO
  367. self.SwapBuffers()
  368. # draw fringe after SwapBuffers, otherwise it don't have to be visible
  369. # on some computers
  370. if self.render['quick'] is False:
  371. self._display.DrawFringe()
  372. stop = time.clock()
  373. if self.render['quick'] is False:
  374. self.parent.statusbarWin['progress'].SetValue(2)
  375. # hide process bar
  376. self.parent.statusbarWin['progress'].Hide()
  377. Debug.msg(3, "GLWindow.UpdateMap(): quick = %d, -> time = %g" % \
  378. (self.render['quick'], (stop-start)))
  379. def EraseMap(self):
  380. """!Erase the canvas
  381. """
  382. self._display.EraseMap()
  383. self.SwapBuffers()
  384. def IsLoaded(self, item):
  385. """!Check if layer (item) is already loaded
  386. @param item layer item
  387. """
  388. layer = self.tree.GetPyData(item)[0]['maplayer']
  389. data = self.tree.GetPyData(item)[0]['nviz']
  390. if not data:
  391. return 0
  392. if layer.type == 'raster':
  393. if 'object' not in data['surface']:
  394. return 0
  395. elif layer.type == 'vector':
  396. if 'object' not in data['vlines'] and \
  397. 'object' not in data['points']:
  398. return 0
  399. return 1
  400. def _GetDataLayers(self, item, litems):
  401. """!Return get list of enabled map layers"""
  402. # load raster & vector maps
  403. while item and item.IsOk():
  404. type = self.tree.GetPyData(item)[0]['type']
  405. if type == 'group':
  406. subItem = self.tree.GetFirstChild(item)[0]
  407. self._GetDataLayers(subItem, litems)
  408. item = self.tree.GetNextSibling(item)
  409. if not item.IsChecked() or \
  410. type not in ('raster', 'vector', '3d-raster'):
  411. item = self.tree.GetNextSibling(item)
  412. continue
  413. litems.append(item)
  414. item = self.tree.GetNextSibling(item)
  415. def LoadDataLayers(self):
  416. """!Load raster/vector from current layer tree
  417. @todo volumes
  418. """
  419. if not self.tree:
  420. return
  421. listOfItems = []
  422. item = self.tree.GetFirstChild(self.tree.root)[0]
  423. self._GetDataLayers(item, listOfItems)
  424. start = time.time()
  425. while(len(listOfItems) > 0):
  426. item = listOfItems.pop()
  427. type = self.tree.GetPyData(item)[0]['type']
  428. if item in self.layers:
  429. continue
  430. # "raster (double click to set properties)" - tries to load this
  431. # layer - no idea how to fix it
  432. if ' ' in self.tree.GetPyData(item)[0]['maplayer'].name:
  433. return
  434. try:
  435. if type == 'raster':
  436. self.LoadRaster(item)
  437. elif type == '3d-raster':
  438. self.LoadRaster3d(item)
  439. elif type == 'vector':
  440. # data = self.tree.GetPyData(item)[0]['nviz']
  441. # vecType = []
  442. # if data and 'vector' in data:
  443. # for v in ('lines', 'points'):
  444. # if data['vector'][v]:
  445. # vecType.append(v)
  446. layer = self.tree.GetPyData(item)[0]['maplayer']
  447. npoints, nlines, nfeatures, mapIs3D = self.lmgr.nviz.VectorInfo(layer)
  448. if npoints > 0:
  449. self.LoadVector(item, points = True)
  450. if nlines > 0:
  451. self.LoadVector(item, points = False)
  452. except gcmd.GException, e:
  453. GError(parent = self,
  454. message = e.value)
  455. self.init = False
  456. stop = time.time()
  457. Debug.msg(3, "GLWindow.LoadDataLayers(): time = %f" % (stop-start))
  458. def UnloadDataLayers(self):
  459. """!Unload any layers that have been deleted from layer tree"""
  460. if not self.tree:
  461. return
  462. listOfItems = []
  463. item = self.tree.GetFirstChild(self.tree.root)[0]
  464. self._GetDataLayers(item, listOfItems)
  465. start = time.time()
  466. for layer in self.layers:
  467. if layer not in listOfItems:
  468. ltype = self.tree.GetPyData(layer)[0]['type']
  469. try:
  470. if ltype == 'raster':
  471. self.UnloadRaster(layer)
  472. elif ltype == '3d-raster':
  473. self.UnloadRaster3d(layer)
  474. elif ltype == 'vector':
  475. self.UnloadVector(layer, True)
  476. self.UnloadVector(layer, False)
  477. self.UpdateView(None)
  478. except gcmd.GException, e:
  479. gcmd.GError(parent = self,
  480. message = e.value)
  481. self.lmgr.nviz.UpdateSettings()
  482. stop = time.time()
  483. Debug.msg(3, "GLWindow.UnloadDataLayers(): time = %f" % (stop-start))
  484. def SetVectorSurface(self, data):
  485. """!Set reference surfaces of vector"""
  486. data['mode']['surface'] = {}
  487. data['mode']['surface']['value'] = list()
  488. data['mode']['surface']['show'] = list()
  489. for name in self.GetLayerNames('raster'):
  490. data['mode']['surface']['value'].append(name)
  491. data['mode']['surface']['show'].append(True)
  492. def SetVectorFromCmd(self, item, data):
  493. """!Set 3D view properties from cmd (d.vect)
  494. @param item Layer Tree item
  495. @param nviz data
  496. """
  497. cmd = self.tree.GetPyData(item)[0]['cmd']
  498. if cmd[0] != 'd.vect':
  499. return
  500. for opt in cmd[1:]:
  501. try:
  502. key, value = opt.split('=')
  503. except ValueError:
  504. continue
  505. if key == 'color':
  506. data['lines']['color']['value'] = value
  507. data['points']['color']['value'] = value
  508. def SetMapObjProperties(self, item, id, nvizType):
  509. """!Set map object properties
  510. Properties must be afterwards updated by
  511. UpdateMapObjProperties().
  512. @param item layer item
  513. @param id nviz layer id (or -1)
  514. @param nvizType nviz data type (surface, points, vector)
  515. """
  516. if nvizType != 'constant':
  517. type = self.tree.GetPyData(item)[0]['maplayer'].type
  518. # reference to original layer properties (can be None)
  519. data = self.tree.GetPyData(item)[0]['nviz']
  520. else:
  521. type = nvizType
  522. data = self.constants[item]
  523. if not data:
  524. # init data structure
  525. if nvizType != 'constant':
  526. self.tree.GetPyData(item)[0]['nviz'] = {}
  527. data = self.tree.GetPyData(item)[0]['nviz']
  528. if type == 'raster':
  529. # reset to default properties
  530. data[nvizType] = self.nvizDefault.SetSurfaceDefaultProp()
  531. elif type == 'vector':
  532. # reset to default properties (lines/points)
  533. data['vector'] = self.nvizDefault.SetVectorDefaultProp()
  534. self.SetVectorFromCmd(item, data['vector'])
  535. self.SetVectorSurface(data['vector']['points'])
  536. self.SetVectorSurface(data['vector']['lines'])
  537. elif type == '3d-raster':
  538. # reset to default properties
  539. data[nvizType] = self.nvizDefault.SetVolumeDefaultProp()
  540. elif type == 'constant':
  541. data['constant'] = self.nvizDefault.SetConstantDefaultProp()
  542. else:
  543. # complete data (use default values), not sure if this is necessary
  544. if type == 'raster':
  545. if not data['surface']:
  546. data['surface'] = self.nvizDefault.SetSurfaceDefaultProp()
  547. if type == 'vector':
  548. if not data['vector']['lines']:
  549. self.nvizDefault.SetVectorLinesDefaultProp(data['vector']['lines'])
  550. if not data['vector']['points']:
  551. self.nvizDefault.SetVectorPointsDefaultProp(data['vector']['points'])
  552. # set updates
  553. for sec in data.keys():
  554. for sec1 in data[sec].keys():
  555. if sec1 == 'position':
  556. data[sec][sec1]['update'] = None
  557. continue
  558. for sec2 in data[sec][sec1].keys():
  559. if sec2 != 'all':
  560. data[sec][sec1][sec2]['update'] = None
  561. event = wxUpdateProperties(data = data)
  562. wx.PostEvent(self, event)
  563. # set id
  564. if id > 0:
  565. if type in ('raster', '3d-raster'):
  566. data[nvizType]['object'] = { 'id' : id,
  567. 'init' : False }
  568. elif type == 'vector':
  569. data['vector'][nvizType]['object'] = { 'id' : id,
  570. 'init' : False }
  571. elif type == 'constant':
  572. data[nvizType]['object'] = { 'id' : id,
  573. 'init' : False }
  574. return data
  575. def LoadRaster(self, item):
  576. """!Load 2d raster map and set surface attributes
  577. @param layer item
  578. """
  579. return self._loadRaster(item)
  580. def LoadRaster3d(self, item):
  581. """!Load 3d raster map and set surface attributes
  582. @param layer item
  583. """
  584. return self._loadRaster(item)
  585. def _loadRaster(self, item):
  586. """!Load 2d/3d raster map and set its attributes
  587. @param layer item
  588. """
  589. layer = self.tree.GetPyData(item)[0]['maplayer']
  590. if layer.type not in ('raster', '3d-raster'):
  591. return
  592. if layer.type == 'raster':
  593. id = self._display.LoadSurface(str(layer.name), None, None)
  594. nvizType = 'surface'
  595. errorMsg = _("Loading raster map")
  596. elif layer.type == '3d-raster':
  597. id = self._display.LoadVolume(str(layer.name), None, None)
  598. nvizType = 'volume'
  599. errorMsg = _("Loading 3d raster map")
  600. else:
  601. id = -1
  602. if id < 0:
  603. if layer.type in ('raster', '3d-raster'):
  604. self.log.WriteError("%s <%s> %s" % (errorMsg, layer.name, _("failed")))
  605. else:
  606. self.log.WriteError(_("Unsupported layer type '%s'") % layer.type)
  607. self.layers.append(item)
  608. # set default/workspace layer properties
  609. data = self.SetMapObjProperties(item, id, nvizType)
  610. # update properties
  611. event = wxUpdateProperties(data = data)
  612. wx.PostEvent(self, event)
  613. # update tools window
  614. if hasattr(self.lmgr, "nviz") and \
  615. item == self.GetSelectedLayer(type = 'item'):
  616. toolWin = self.lmgr.nviz
  617. if layer.type == 'raster':
  618. win = toolWin.FindWindowById( \
  619. toolWin.win['vector']['lines']['surface'])
  620. win.SetItems(self.GetLayerNames(layer.type))
  621. #toolWin.UpdatePage(nvizType)
  622. #toolWin.SetPage(nvizType)
  623. return id
  624. def NewConstant(self):
  625. """!Create new constant"""
  626. index = len(self.constants)
  627. try:
  628. name = self.constants[-1]['constant']['object']['name'] + 1
  629. except IndexError:
  630. name = 1
  631. data = dict()
  632. self.constants.append(data)
  633. data = self.SetMapObjProperties(item = index, id = -1, nvizType = 'constant')
  634. self.AddConstant(data, name)
  635. return name
  636. def AddConstant(self, data, name):
  637. """!Add new constant"""
  638. id = self._display.AddConstant(value = data['constant']['value'], color = data['constant']['color'])
  639. self._display.SetSurfaceRes(id, data['constant']['resolution'], data['constant']['resolution'])
  640. data['constant']['object'] = { 'id' : id,
  641. 'name': name,
  642. 'init' : False }
  643. def DeleteConstant(self, index):
  644. """!Delete constant layer"""
  645. id = self.constants[index]['constant']['object']['id']
  646. self._display.UnloadSurface(id)
  647. del self.constants[index]
  648. def SelectCPlane(self, index):
  649. """!Select cutting plane"""
  650. for plane in range (self._display.GetCPlanesCount()):
  651. if plane == index:
  652. self._display.SelectCPlane(plane)
  653. else:
  654. self._display.UnselectCPlane(plane)
  655. def UpdateCPlane(self, event):
  656. """!Change cutting plane settings"""
  657. current = event.current
  658. for each in event.update:
  659. if each == 'rotation':
  660. self._display.SetCPlaneRotation(0, self.cplanes[current]['rotation']['tilt'],
  661. self.cplanes[current]['rotation']['rot'])
  662. if each == 'position':
  663. self._display.SetCPlaneTranslation(self.cplanes[current]['position']['x'],
  664. self.cplanes[current]['position']['y'],
  665. self.cplanes[current]['position']['z'])
  666. if each == 'shading':
  667. self._display.SetFenceColor(self.cplanes[current]['shading'])
  668. def UnloadRaster(self, item):
  669. """!Unload 2d raster map
  670. @param layer item
  671. """
  672. return self._unloadRaster(item)
  673. def UnloadRaster3d(self, item):
  674. """!Unload 3d raster map
  675. @param layer item
  676. """
  677. return self._unloadRaster(item)
  678. def _unloadRaster(self, item):
  679. """!Unload 2d/3d raster map
  680. @param item layer item
  681. """
  682. layer = self.tree.GetPyData(item)[0]['maplayer']
  683. if layer.type not in ('raster', '3d-raster'):
  684. return
  685. data = self.tree.GetPyData(item)[0]['nviz']
  686. if layer.type == 'raster':
  687. nvizType = 'surface'
  688. unloadFn = self._display.UnloadSurface
  689. errorMsg = _("Unable to unload raster map")
  690. successMsg = _("Raster map")
  691. else:
  692. nvizType = 'volume'
  693. unloadFn = self._display.UnloadVolume
  694. errorMsg = _("Unable to unload 3d raster map")
  695. successMsg = _("3d raster map")
  696. id = data[nvizType]['object']['id']
  697. if unloadFn(id) == 0:
  698. self.log.WriteError("%s <%s>" % (errorMsg, layer.name))
  699. else:
  700. self.log.WriteLog("%s <%s> %s" % (successMsg, layer.name, _("unloaded successfully")))
  701. data[nvizType].pop('object')
  702. self.layers.remove(item)
  703. # update tools window
  704. if hasattr(self.lmgr, "nviz") and \
  705. layer.type == 'raster':
  706. toolWin = self.lmgr.nviz
  707. win = toolWin.FindWindowById( \
  708. toolWin.win['vector']['lines']['surface'])
  709. win.SetItems(self.GetLayerNames(layer.type))
  710. def LoadVector(self, item, points = None):
  711. """!Load 2D or 3D vector map overlay
  712. @param item layer item
  713. @param points True to load points, False to load lines, None
  714. to load both
  715. """
  716. layer = self.tree.GetPyData(item)[0]['maplayer']
  717. if layer.type != 'vector':
  718. return
  719. # set default properties
  720. if points is None:
  721. self.SetMapObjProperties(item, -1, 'lines')
  722. self.SetMapObjProperties(item, -1, 'points')
  723. vecTypes = ('points', 'lines')
  724. elif points:
  725. self.SetMapObjProperties(item, -1, 'points')
  726. vecTypes = ('points', )
  727. else:
  728. self.SetMapObjProperties(item, -1, 'lines')
  729. vecTypes = ('lines', )
  730. id = -1
  731. for vecType in vecTypes:
  732. if vecType == 'lines':
  733. id = self._display.LoadVector(str(layer.GetName()), False)
  734. else:
  735. id = self._display.LoadVector(str(layer.GetName()), True)
  736. if id < 0:
  737. self.log.WriteError(_("Loading vector map <%(name)s> (%(type)s) failed") % \
  738. { 'name' : layer.name, 'type' : vecType })
  739. # update layer properties
  740. self.SetMapObjProperties(item, id, vecType)
  741. self.layers.append(item)
  742. # update properties
  743. data = self.tree.GetPyData(item)[0]['nviz']
  744. event = wxUpdateProperties(data = data)
  745. wx.PostEvent(self, event)
  746. # update tools window
  747. if hasattr(self.lmgr, "nviz") and \
  748. item == self.GetSelectedLayer(type = 'item'):
  749. toolWin = self.lmgr.nviz
  750. toolWin.UpdatePage('vector')
  751. ### toolWin.SetPage('vector')
  752. return id
  753. def UnloadVector(self, item, points = None):
  754. """!Unload vector map overlay
  755. @param item layer item
  756. @param points,lines True to unload given feature type
  757. """
  758. layer = self.tree.GetPyData(item)[0]['maplayer']
  759. data = self.tree.GetPyData(item)[0]['nviz']['vector']
  760. # if vecType is None:
  761. # vecType = []
  762. # for v in ('lines', 'points'):
  763. # if UserSettings.Get(group = 'nviz', key = 'vector',
  764. # subkey = [v, 'show']):
  765. # vecType.append(v)
  766. if points is None:
  767. vecTypes = ('points', 'lines')
  768. elif points:
  769. vecTypes = ('points', )
  770. else:
  771. vecTypes = ('lines', )
  772. for vecType in vecTypes:
  773. if 'object' not in data[vecType]:
  774. continue
  775. id = data[vecType]['object']['id']
  776. if vecType == 'lines':
  777. ret = self._display.UnloadVector(id, False)
  778. else:
  779. ret = self._display.UnloadVector(id, True)
  780. if ret == 0:
  781. self.log.WriteError(_("Unable to unload vector map <%(name)s> (%(type)s)") % \
  782. { 'name': layer.name, 'type' : vecType })
  783. else:
  784. self.log.WriteLog(_("Vector map <%(name)s> (%(type)s) unloaded successfully") % \
  785. { 'name' : layer.name, 'type' : vecType })
  786. data[vecType].pop('object')
  787. ### self.layers.remove(id)
  788. def Reset(self):
  789. """!Reset (unload data)"""
  790. for item in self.layers:
  791. type = self.tree.GetPyData(item)[0]['maplayer'].type
  792. if type == 'raster':
  793. self.UnloadRaster(item)
  794. elif type == '3d-raster':
  795. self.UnloadRaster3d(item)
  796. elif type == 'vector':
  797. self.UnloadVector(item)
  798. self.init = False
  799. def OnZoomToMap(self, event):
  800. """!Set display extents to match selected raster or vector
  801. map or volume.
  802. @todo vector, volume
  803. """
  804. layer = self.GetSelectedLayer()
  805. if layer is None:
  806. return
  807. Debug.msg (3, "GLWindow.OnZoomToMap(): layer = %s, type = %s" % \
  808. (layer.name, layer.type))
  809. self._display.SetViewportDefault()
  810. def ResetView(self):
  811. """!Reset to default view"""
  812. self.view['z-exag']['value'], \
  813. self.iview['height']['value'], \
  814. self.iview['height']['min'], \
  815. self.iview['height']['max'] = self._display.SetViewDefault()
  816. self.view['z-exag']['min'] = 0
  817. self.view['z-exag']['max'] = self.view['z-exag']['value'] * 10
  818. self.view['position']['x'] = UserSettings.Get(group = 'nviz', key = 'view',
  819. subkey = ('position', 'x'))
  820. self.view['position']['y'] = UserSettings.Get(group = 'nviz', key = 'view',
  821. subkey = ('position', 'y'))
  822. self.view['persp']['value'] = UserSettings.Get(group = 'nviz', key = 'view',
  823. subkey = ('persp', 'value'))
  824. self.view['twist']['value'] = UserSettings.Get(group = 'nviz', key = 'view',
  825. subkey = ('twist', 'value'))
  826. self._display.LookAtCenter()
  827. focus = self.iview['focus']
  828. focus['x'], focus['y'], focus['z'] = self._display.GetFocus()
  829. event = wxUpdateView(zExag = False)
  830. wx.PostEvent(self, event)
  831. def UpdateMapObjProperties(self, event):
  832. """!Generic method to update data layer properties"""
  833. data = event.data
  834. if 'surface' in data:
  835. try:
  836. id = data['surface']['object']['id']
  837. except KeyError:
  838. pass
  839. self.UpdateSurfaceProperties(id, data['surface'])
  840. # -> initialized
  841. data['surface']['object']['init'] = True
  842. elif 'constant' in data:
  843. id = data['constant']['object']['id']
  844. self.UpdateConstantProperties(id, data['constant'])
  845. # -> initialized
  846. data['constant']['object']['init'] = True
  847. elif 'volume' in data:
  848. id = data['volume']['object']['id']
  849. self.UpdateVolumeProperties(id, data['volume'])
  850. # -> initialized
  851. data['volume']['object']['init'] = True
  852. elif 'vector' in data:
  853. for type in ('lines', 'points'):
  854. if 'object' in data['vector'][type]:
  855. id = data['vector'][type]['object']['id']
  856. self.UpdateVectorProperties(id, data['vector'], type)
  857. # -> initialized
  858. data['vector'][type]['object']['init'] = True
  859. def UpdateConstantProperties(self, id, data):
  860. """!Update surface map object properties"""
  861. self._display.SetSurfaceColor(id = id, map = False, value = data['color'])
  862. self._display.SetSurfaceTopo(id = id, map = False, value = data['value'])
  863. self._display.SetSurfaceRes(id, data['resolution'], data['resolution'])
  864. def UpdateSurfaceProperties(self, id, data):
  865. """!Update surface map object properties"""
  866. # surface attributes
  867. for attrb in ('color', 'mask',
  868. 'transp', 'shine'):
  869. if attrb not in data['attribute'] or \
  870. 'update' not in data['attribute'][attrb]:
  871. continue
  872. map = data['attribute'][attrb]['map']
  873. value = data['attribute'][attrb]['value']
  874. if map is None: # unset
  875. # only optional attributes
  876. if attrb == 'mask':
  877. # TODO: invert mask
  878. # TODO: broken in NVIZ
  879. self._display.UnsetSurfaceMask(id)
  880. elif attrb == 'transp':
  881. self._display.UnsetSurfaceTransp(id)
  882. else:
  883. if type(value) == type('') and \
  884. len(value) <= 0: # ignore empty values (TODO: warning)
  885. continue
  886. if attrb == 'color':
  887. self._display.SetSurfaceColor(id, map, str(value))
  888. elif attrb == 'mask':
  889. # TODO: invert mask
  890. # TODO: broken in NVIZ
  891. self._display.SetSurfaceMask(id, False, str(value))
  892. elif attrb == 'transp':
  893. self._display.SetSurfaceTransp(id, map, str(value))
  894. elif attrb == 'shine':
  895. self._display.SetSurfaceShine(id, map, str(value))
  896. data['attribute'][attrb].pop('update')
  897. # draw res
  898. if 'update' in data['draw']['resolution']:
  899. coarse = data['draw']['resolution']['coarse']
  900. fine = data['draw']['resolution']['fine']
  901. if data['draw']['all']:
  902. self._display.SetSurfaceRes(-1, fine, coarse)
  903. else:
  904. self._display.SetSurfaceRes(id, fine, coarse)
  905. data['draw']['resolution'].pop('update')
  906. # draw style
  907. if 'update' in data['draw']['mode']:
  908. if data['draw']['mode']['value'] < 0: # need to calculate
  909. data['draw']['mode']['value'] = \
  910. self.nvizDefault.GetDrawMode(mode = data['draw']['mode']['desc']['mode'],
  911. style = data['draw']['mode']['desc']['style'],
  912. shade = data['draw']['mode']['desc']['shading'],
  913. string = True)
  914. style = data['draw']['mode']['value']
  915. if data['draw']['all']:
  916. self._display.SetSurfaceStyle(-1, style)
  917. else:
  918. self._display.SetSurfaceStyle(id, style)
  919. data['draw']['mode'].pop('update')
  920. # wire color
  921. if 'update' in data['draw']['wire-color']:
  922. color = data['draw']['wire-color']['value']
  923. if data['draw']['all']:
  924. self._display.SetWireColor(-1, str(color))
  925. else:
  926. self._display.SetWireColor(id, str(color))
  927. data['draw']['wire-color'].pop('update')
  928. # position
  929. if 'update' in data['position']:
  930. x = data['position']['x']
  931. y = data['position']['y']
  932. z = data['position']['z']
  933. self._display.SetSurfacePosition(id, x, y, z)
  934. data['position'].pop('update')
  935. data['draw']['all'] = False
  936. def UpdateVolumeProperties(self, id, data, isosurfId = None):
  937. """!Update volume (isosurface/slice) map object properties"""
  938. if 'update' in data['draw']['resolution']:
  939. self._display.SetIsosurfaceRes(id, data['draw']['resolution']['value'])
  940. data['draw']['resolution'].pop('update')
  941. if 'update' in data['draw']['shading']:
  942. if data['draw']['shading']['value'] < 0: # need to calculate
  943. data['draw']['shading']['value'] = \
  944. self.nvizDefault.GetDrawMode(shade = data['draw']['shading'],
  945. string = False)
  946. data['draw']['shading'].pop('update')
  947. #
  948. # isosurface attributes
  949. #
  950. isosurfId = 0
  951. for isosurf in data['isosurface']:
  952. for attrb in ('color', 'mask',
  953. 'transp', 'shine', 'emit'):
  954. if attrb not in isosurf or \
  955. 'update' not in isosurf[attrb]:
  956. continue
  957. map = isosurf[attrb]['map']
  958. value = isosurf[attrb]['value']
  959. if map is None: # unset
  960. # only optional attributes
  961. if attrb == 'mask':
  962. # TODO: invert mask
  963. # TODO: broken in NVIZ
  964. self._display.UnsetIsosurfaceMask(id, isosurfId)
  965. elif attrb == 'transp':
  966. self._display.UnsetIsosurfaceTransp(id, isosurfId)
  967. elif attrb == 'emit':
  968. self._display.UnsetIsosurfaceEmit(id, isosurfId)
  969. else:
  970. if type(value) == type('') and \
  971. len(value) <= 0: # ignore empty values (TODO: warning)
  972. continue
  973. elif attrb == 'color':
  974. self._display.SetIsosurfaceColor(id, isosurfId, map, str(value))
  975. elif attrb == 'mask':
  976. # TODO: invert mask
  977. # TODO: broken in NVIZ
  978. self._display.SetIsosurfaceMask(id, isosurfId, False, str(value))
  979. elif attrb == 'transp':
  980. self._display.SetIsosurfaceTransp(id, isosurfId, map, str(value))
  981. elif attrb == 'shine':
  982. self._display.SetIsosurfaceShine(id, isosurfId, map, str(value))
  983. elif attrb == 'emit':
  984. self._display.SetIsosurfaceEmit(id, isosurfId, map, str(value))
  985. isosurf[attrb].pop('update')
  986. isosurfId += 1
  987. def UpdateVectorProperties(self, id, data, type):
  988. """!Update vector layer properties
  989. @param id layer id
  990. @param data properties
  991. @param type lines/points
  992. """
  993. if type == 'points':
  994. self.UpdateVectorPointsProperties(id, data[type])
  995. else:
  996. self.UpdateVectorLinesProperties(id, data[type])
  997. def UpdateVectorLinesProperties(self, id, data):
  998. """!Update vector line map object properties"""
  999. # mode
  1000. if 'update' in data['color'] or \
  1001. 'update' in data['width'] or \
  1002. 'update' in data['mode']:
  1003. width = data['width']['value']
  1004. color = data['color']['value']
  1005. if data['mode']['type'] == 'flat':
  1006. flat = True
  1007. if 'surface' in data['mode']:
  1008. data['mode'].pop('surface')
  1009. else:
  1010. flat = False
  1011. self._display.SetVectorLineMode(id, color,
  1012. width, flat)
  1013. if 'update' in data['color']:
  1014. data['color'].pop('update')
  1015. if 'update' in data['width']:
  1016. data['width'].pop('update')
  1017. # height
  1018. if 'update' in data['height']:
  1019. self._display.SetVectorLineHeight(id,
  1020. data['height']['value'])
  1021. data['height'].pop('update')
  1022. # surface
  1023. if 'surface' in data['mode'] and 'update' in data['mode']:
  1024. for item in range(len(data['mode']['surface']['value'])):
  1025. for type in ('raster', 'constant'):
  1026. sid = self.GetLayerId(type = type,
  1027. name = data['mode']['surface']['value'][item])
  1028. if sid > -1:
  1029. if data['mode']['surface']['show'][item]:
  1030. self._display.SetVectorLineSurface(id, sid)
  1031. else:
  1032. self._display.UnsetVectorLineSurface(id, sid)
  1033. break
  1034. if 'update' in data['mode']:
  1035. data['mode'].pop('update')
  1036. def UpdateVectorPointsProperties(self, id, data):
  1037. """!Update vector point map object properties"""
  1038. if 'update' in data['size'] or \
  1039. 'update' in data['width'] or \
  1040. 'update' in data['marker'] or \
  1041. 'update' in data['color']:
  1042. ret = self._display.SetVectorPointMode(id, data['color']['value'],
  1043. data['width']['value'], float(data['size']['value']),
  1044. data['marker']['value'] + 1)
  1045. error = None
  1046. if ret == -1:
  1047. error = _("Vector point layer not found (id = %d)") % id
  1048. elif ret == -2:
  1049. error = _("Unable to set data layer properties (id = %d)") % id
  1050. if error:
  1051. raise gcmd.GException(_("Setting data layer properties failed.\n\n%s") % error)
  1052. for prop in ('size', 'width', 'marker', 'color'):
  1053. if 'update' in data[prop]:
  1054. data[prop].pop('update')
  1055. # height
  1056. if 'update' in data['height']:
  1057. self._display.SetVectorPointHeight(id,
  1058. data['height']['value'])
  1059. data['height'].pop('update')
  1060. # surface
  1061. if 'update' in data['mode']:
  1062. for item in range(len(data['mode']['surface']['value'])):
  1063. for type in ('raster', 'constant'):
  1064. sid = self.GetLayerId(type = type,
  1065. name = data['mode']['surface']['value'][item])
  1066. if sid > -1:
  1067. if data['mode']['surface']['show'][item]:
  1068. self._display.SetVectorPointSurface(id, sid)
  1069. else:
  1070. self._display.UnsetVectorPointSurface(id, sid)
  1071. break
  1072. data['mode'].pop('update')
  1073. def GetLayerNames(self, type):
  1074. """!Return list of map layer names of given type"""
  1075. layerName = []
  1076. if type == 'constant':
  1077. for item in self.constants:
  1078. layerName.append(_("constant#") + str(item['constant']['object']['name']))
  1079. else:
  1080. for item in self.layers:
  1081. mapLayer = self.tree.GetPyData(item)[0]['maplayer']
  1082. if type != mapLayer.GetType():
  1083. continue
  1084. layerName.append(mapLayer.GetName())
  1085. return layerName
  1086. def GetLayerId(self, type, name, vsubtyp = None):
  1087. """!Get layer object id or -1"""
  1088. if len(name) < 1:
  1089. return -1
  1090. if type == 'constant':
  1091. for item in self.constants:
  1092. if _("constant#") + str(item['constant']['object']['name']) == name:
  1093. return item['constant']['object']['id']
  1094. return self.constants
  1095. for item in self.layers:
  1096. mapLayer = self.tree.GetPyData(item)[0]['maplayer']
  1097. if type != mapLayer.GetType() or \
  1098. name != mapLayer.GetName():
  1099. continue
  1100. data = self.tree.GetPyData(item)[0]['nviz']
  1101. if type == 'raster':
  1102. return data['surface']['object']['id']
  1103. elif type == 'vector':
  1104. if vsubtyp == 'vpoint':
  1105. return data['vector']['points']['object']['id']
  1106. elif vsubtyp == 'vline':
  1107. return data['vector']['lines']['object']['id']
  1108. elif type == '3d-raster':
  1109. return data['volume']['object']['id']
  1110. return -1
  1111. def NvizCmdCommand(self):
  1112. """!Generate command for nviz_cmd according to current state"""
  1113. cmd = 'nviz_cmd '
  1114. rasters = []
  1115. vectors = []
  1116. for item in self.layers:
  1117. if self.tree.GetPyData(item)[0]['type'] == 'raster':
  1118. rasters.append(item)
  1119. elif self.tree.GetPyData(item)[0]['type'] == 'vector':
  1120. vectors.append(item)
  1121. if not rasters and not self.constants:
  1122. return _("At least one raster map required")
  1123. # elevation_map/elevation_value
  1124. if self.constants:
  1125. subcmd = "elevation_value="
  1126. for constant in self.constants:
  1127. subcmd += "%d," % constant['constant']['value']
  1128. subcmd = subcmd.strip(', ') + ' '
  1129. cmd += subcmd
  1130. if rasters:
  1131. subcmd = "elevation_map="
  1132. for item in rasters:
  1133. subcmd += "%s," % self.tree.GetPyData(item)[0]['maplayer'].GetName()
  1134. subcmd = subcmd.strip(', ') + ' '
  1135. cmd += subcmd
  1136. #
  1137. # draw mode
  1138. #
  1139. cmdMode = "mode="
  1140. cmdFine = "resolution_fine="
  1141. cmdCoarse = "resolution_coarse="
  1142. cmdShading = "shading="
  1143. cmdStyle = "style="
  1144. cmdWire = "wire_color="
  1145. # test -a flag
  1146. flag_a = "-a "
  1147. nvizDataFirst = self.tree.GetPyData(rasters[0])[0]['nviz']['surface']['draw']
  1148. for item in rasters:
  1149. nvizData = self.tree.GetPyData(item)[0]['nviz']['surface']['draw']
  1150. if nvizDataFirst != nvizData:
  1151. flag_a = ""
  1152. cmd += flag_a
  1153. for item in rasters:
  1154. nvizData = self.tree.GetPyData(item)[0]['nviz']['surface']['draw']
  1155. cmdMode += "%s," % nvizData['mode']['desc']['mode']
  1156. cmdFine += "%s," % nvizData['resolution']['fine']
  1157. cmdCoarse += "%s," % nvizData['resolution']['coarse']
  1158. cmdShading += "%s," % nvizData['mode']['desc']['shading']
  1159. cmdStyle += "%s," % nvizData['mode']['desc']['style']
  1160. cmdWire += "%s," % nvizData['wire-color']['value']
  1161. for item in self.constants:
  1162. cmdMode += "fine,"
  1163. cmdFine += "%s," % item['constant']['resolution']
  1164. cmdCoarse += "%s," % item['constant']['resolution']
  1165. cmdShading += "gouraud,"
  1166. cmdStyle += "surface,"
  1167. cmdWire += "0:0:0,"
  1168. mode = []
  1169. for subcmd in (cmdMode, cmdFine, cmdCoarse, cmdShading, cmdStyle, cmdWire):
  1170. if flag_a:
  1171. mode.append(subcmd.split(',')[0] + ' ')
  1172. else:
  1173. subcmd = subcmd.strip(', ') + ' '
  1174. cmd += subcmd
  1175. if flag_a:# write only meaningful possibilities
  1176. cmd += mode[0]
  1177. if 'fine' in mode[0]:
  1178. cmd += mode[1]
  1179. elif 'coarse' in mode[0]:
  1180. cmd += mode[2]
  1181. elif 'both' in mode[0]:
  1182. cmd += mode[2]
  1183. cmd += mode[1]
  1184. if 'flat' in mode[3]:
  1185. cmd += mode[3]
  1186. if 'wire' in mode[4]:
  1187. cmd += mode[4]
  1188. if 'coarse' in mode[0] or 'both' in mode[0] and 'wire' in mode[3]:
  1189. cmd += mode[5]
  1190. #
  1191. # attributes
  1192. #
  1193. cmdColorMap = "color_map="
  1194. cmdColorVal = "color="
  1195. for item in rasters:
  1196. nvizData = self.tree.GetPyData(item)[0]['nviz']['surface']['attribute']
  1197. if 'color' not in nvizData:
  1198. cmdColorMap += "%s," % self.tree.GetPyData(item)[0]['maplayer'].GetName()
  1199. else:
  1200. if nvizData['color']['map']:
  1201. cmdColorMap += "%s," % nvizData['color']['value']
  1202. else:
  1203. cmdColorVal += "%s," % nvizData['color']['value']
  1204. #TODO
  1205. # transparency, shine, mask
  1206. for item in self.constants:
  1207. cmdColorVal += "%s," % item['constant']['color']
  1208. if cmdColorMap.split("=")[1]:
  1209. cmd += cmdColorMap.strip(', ') + ' '
  1210. if cmdColorVal.split("=")[1]:
  1211. cmd += cmdColorVal.strip(', ') + ' '
  1212. cmd += "\\\n"
  1213. #
  1214. # vlines
  1215. #
  1216. if vectors:
  1217. cmdLines = cmdLWidth = cmdLHeight = cmdLColor = cmdLMode = cmdLPos = \
  1218. cmdPoints = cmdPWidth = cmdPSize = cmdPColor = cmdPMarker = cmdPPos = cmdPLayer = ""
  1219. markers = ['x', 'box', 'sphere', 'cube', 'diamond',
  1220. 'dec_tree', 'con_tree', 'aster', 'gyro', 'histogram']
  1221. for vector in vectors:
  1222. npoints, nlines, nfeatures, mapIs3D = self.lmgr.nviz.VectorInfo(
  1223. self.tree.GetPyData(vector)[0]['maplayer'])
  1224. nvizData = self.tree.GetPyData(vector)[0]['nviz']['vector']
  1225. if nlines > 0:
  1226. cmdLines += "%s," % self.tree.GetPyData(vector)[0]['maplayer'].GetName()
  1227. cmdLWidth += "%d," % nvizData['lines']['width']['value']
  1228. cmdLHeight += "%d," % nvizData['lines']['height']['value']
  1229. cmdLColor += "%s," % nvizData['lines']['color']['value']
  1230. cmdLMode += "%s," % nvizData['lines']['mode']['type']
  1231. cmdLPos += "0,0,%d," % nvizData['lines']['height']['value']
  1232. if npoints > 0:
  1233. cmdPoints += "%s," % self.tree.GetPyData(vector)[0]['maplayer'].GetName()
  1234. cmdPWidth += "%d," % nvizData['points']['width']['value']
  1235. cmdPSize += "%d," % nvizData['points']['size']['value']
  1236. cmdPColor += "%s," % nvizData['points']['color']['value']
  1237. cmdPMarker += "%s," % markers[nvizData['points']['marker']['value']]
  1238. cmdPPos += "0,0,%d," % nvizData['points']['height']['value']
  1239. cmdPLayer += "1,1,"
  1240. if cmdLines:
  1241. cmd += "vline=" + cmdLines.strip(',') + ' '
  1242. cmd += "vline_width=" + cmdLWidth.strip(',') + ' '
  1243. cmd += "vline_color=" + cmdLColor.strip(',') + ' '
  1244. cmd += "vline_height=" + cmdLHeight.strip(',') + ' '
  1245. cmd += "vline_mode=" + cmdLMode.strip(',') + ' '
  1246. cmd += "vline_position=" + cmdLPos.strip(',') + ' '
  1247. if cmdPoints:
  1248. cmd += "vpoint=" + cmdPoints.strip(',') + ' '
  1249. cmd += "vpoint_width=" + cmdPWidth.strip(',') + ' '
  1250. cmd += "vpoint_color=" + cmdPColor.strip(',') + ' '
  1251. cmd += "vpoint_size=" + cmdPSize.strip(',') + ' '
  1252. cmd += "vpoint_marker=" + cmdPMarker.strip(',') + ' '
  1253. cmd += "vpoint_position=" + cmdPPos.strip(',') + ' '
  1254. cmd += "vpoint_layer=" + cmdPLayer.strip(',') + ' '
  1255. cmd += "\\\n"
  1256. #
  1257. # cutting planes
  1258. #
  1259. cplane = self.lmgr.nviz.FindWindowById(self.lmgr.nviz.win['cplane']['planes']).GetStringSelection()
  1260. try:
  1261. planeIndex = int(cplane.split()[1])
  1262. except IndexError:
  1263. planeIndex = None
  1264. if planeIndex is not None:
  1265. shading = ['clear', 'top', 'bottom', 'blend', 'shaded']
  1266. cmd += "cplane=%d " % planeIndex
  1267. cmd += "cplane_rotation=%d " % self.cplanes[planeIndex]['rotation']['rot']
  1268. cmd += "cplane_tilt=%d " % self.cplanes[planeIndex]['rotation']['tilt']
  1269. cmd += "cplane_position=%d,%d,%d " % (self.cplanes[planeIndex]['position']['x'],
  1270. self.cplanes[planeIndex]['position']['y'],
  1271. self.cplanes[planeIndex]['position']['z'])
  1272. cmd += "cplane_shading=%s " % shading[self.cplanes[planeIndex]['shading']]
  1273. cmd += "\\\n"
  1274. #
  1275. # viewpoint
  1276. #
  1277. subcmd = "position=%.2f,%.2f " % (self.view['position']['x'], self.view['position']['y'])
  1278. subcmd += "height=%d " % (self.iview['height']['value'])
  1279. subcmd += "perspective=%d " % (self.view['persp']['value'])
  1280. subcmd += "twist=%d " % (self.view['twist']['value'])
  1281. subcmd += "zexag=%d " % (self.view['z-exag']['value'])
  1282. subcmd += "focus=%d,%d,%d " % (self.iview['focus']['x'],self.iview['focus']['y'],self.iview['focus']['z'])
  1283. cmd += subcmd
  1284. # background
  1285. subcmd = "bgcolor=%d:%d:%d " % (self.view['background']['color'])
  1286. if self.view['background']['color'] != (255, 255, 255):
  1287. cmd += subcmd
  1288. cmd += "\\\n"
  1289. # light
  1290. subcmd = "light_position=%.2f,%.2f,%.2f " % (self.light['position']['x'],
  1291. self.light['position']['y'],
  1292. self.light['position']['z']/100.)
  1293. subcmd += "light_brightness=%d " % (self.light['bright'])
  1294. subcmd += "light_ambient=%d " % (self.light['ambient'])
  1295. subcmd += "light_color=%d:%d:%d " % (self.light['color'])
  1296. cmd += subcmd
  1297. cmd += "\\\n"
  1298. # fringe
  1299. toolWindow = self.lmgr.nviz
  1300. direction = ''
  1301. for dir in ('nw', 'ne', 'sw', 'se'):
  1302. if toolWindow.FindWindowById(toolWindow.win['fringe'][dir]).IsChecked():
  1303. direction += "%s," % dir
  1304. if direction:
  1305. subcmd = "fringe=%s " % (direction.strip(','))
  1306. color = toolWindow.FindWindowById(toolWindow.win['fringe']['color']).GetValue()
  1307. subcmd += "fringe_color=%d:%d:%d " % (color[0], color[1], color[2])
  1308. subcmd += "fringe_elevation=%d " % (toolWindow.FindWindowById(toolWindow.win['fringe']['elev']).GetValue())
  1309. cmd += subcmd
  1310. cmd += "\\\n"
  1311. # output
  1312. subcmd = 'output=nviz_output '
  1313. subcmd += 'format=ppm '
  1314. subcmd += 'size=%d,%d ' % self.GetClientSizeTuple()
  1315. cmd += subcmd
  1316. return cmd
  1317. def OnNvizCmd(self):
  1318. """!Generate and write command to command output"""
  1319. self.log.WriteLog(self.NvizCmdCommand(), switchPage = True)
  1320. def SaveToFile(self, FileName, FileType, width, height):
  1321. """!This draws the DC to a buffer that can be saved to a file.
  1322. @todo fix BufferedPaintDC
  1323. @param FileName file name
  1324. @param FileType type of bitmap
  1325. @param width image width
  1326. @param height image height
  1327. """
  1328. self._display.SaveToFile(FileName, width, height)
  1329. # pbuffer = wx.EmptyBitmap(max(1, self.Map.width), max(1, self.Map.height))
  1330. # dc = wx.BufferedPaintDC(self, pbuffer)
  1331. # dc.Clear()
  1332. # self.SetCurrent()
  1333. # self._display.Draw(False, -1)
  1334. # pbuffer.SaveFile(FileName, FileType)
  1335. # self.SwapBuffers()
  1336. def GetDisplay(self):
  1337. """!Get display instance"""
  1338. return self._display
  1339. def ZoomToMap(self):
  1340. """!Reset view
  1341. """
  1342. self.lmgr.nviz.OnResetView(None)