<div class="gmail_quote">Hi Benjamin,<br><br><div class="gmail_quote">2009/8/28 Benjamin Ducke <span dir="ltr"><<a href="mailto:benjamin.ducke@oxfordarch.co.uk" target="_blank">benjamin.ducke@oxfordarch.co.uk</a>></span><div class="im">
<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
This is great! Did you consider writing a little plug-on for one<br>
of the VTK viewers, like ParaView or VisIt, to allow direct access<br>
to GRASS data from there?</blockquote></div><div><br>No.I do not consider to implement plug-ins for ParaView or Visit for now.<br>This may change in the far future.<br>First i need to validate the library concept and then i will implement full raster,<br>
vector, volume and grass sql database support .. and this will take a while.<br><br>Then i need to understand the vtkGeo and projection concept in VTK 5.4+ . <br><br>So the plugin development does not have a high priority on my todo list. :)<br>
<br>Cheers<br><font color="#888888">Soeren<br><br></font></div><div><div></div><div class="h5"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
<br>
Cheers,<br>
<br>
Ben<br>
<br>
Soeren Gebbert wrote:<br>
> Hi,<br>
> JFYI i have implemented a small demo library which connects parts of the<br>
> grass gis library<br>
> with the Visual-Tool-Kit VTK.<br>
> It is called vtkGRASSBridge.<br>
><br>
> Additionally the VTK Python and Java wrapping mechanism have been<br>
> enabled. So the user can use the<br>
> library from C++, Python and Java.<br>
><br>
> Features:<br>
> * The library can read and write raster maps row by row, as well as<br>
> regions and the raster history.<br>
> * It can list files of the grass database.<br>
> * It can convert grass raster maps into the vtkImageData format and the<br>
> other way round.<br>
> Missing Features:<br>
> * Full raster support (categories, colors, ..)<br>
> * Imagery support (groups, ...)<br>
> * Full vector support ....<br>
> * ....<br>
><br>
> You can use every vtkDataSet processing algorithm available in VTK (more<br>
> than 50 i think)<br>
> to process the converted raster map and write it back to the grass<br>
> raster database.<br>
><br>
> You can use the VTK render functionality to display the raster map.<br>
><br>
> There are examples and tests in the library showing this functionality<br>
> in Python.<br>
><br>
> The documentation of the library can be found here:<br>
> <a href="http://www-pool.math.tu-berlin.de/%7Esoeren/grass/files/vtkGRASSBridge/doc/html/index.html" target="_blank">http://www-pool.math.tu-berlin.de/~soeren/grass/files/vtkGRASSBridge/doc/html/index.html</a><br>
><br>
> The GPL source code is located here:<br>
> <a href="http://www-pool.math.tu-berlin.de/%7Esoeren/grass/files/vtkGRASSBridge/vtkGRASSBridge_0.1.tar.gz" target="_blank">http://www-pool.math.tu-berlin.de/~soeren/grass/files/vtkGRASSBridge/vtkGRASSBridge_0.1.tar.gz</a><br>
><br>
> Here is a small Python code example to read and write grass raster maps<br>
> in Python (within a grass session):<br>
><br>
> ##############################################<br>
> #include the VTK and vtkGRASSBridge Python libraries<br>
> from libvtkCommonPython import *<br>
> from libvtkFilteringPython import *<br>
> from libvtkGraphicsPython import *<br>
> from libvtkRenderingPython import *<br>
> from libvtkIOPython import *<br>
> from libvtkImagingPython import *<br>
> from libvtkGRASSBridgeIOPython import *<br>
> from libvtkGRASSBridgeCommonPython import *<br>
><br>
> #Initiate grass<br>
> init = vtkGRASSInit()<br>
><br>
> #The region settings<br>
> region = vtkGRASSRegion()<br>
> region.ReadCurrentRegion()<br>
> region.SetCols(10)<br>
> region.SetRows(10)<br>
> region.AdjustRegion()<br>
><br>
> #Create a new map and write some data into it<br>
> writer = vtkGRASSRasterMapWriter()<br>
> writer.SetMapTypeToCELL()<br>
> writer.SetRegion(region)<br>
> writer.UseUserDefinedRegion()<br>
> writer.OpenMap("test_cell")<br>
><br>
> #this is the row buffer<br>
> data = vtkIntArray()<br>
> data.SetNumberOfTuples(writer.GetNumberOfCols())<br>
><br>
> #iterate over each raster cell<br>
> for i in range(writer.GetNumberOfRows()):<br>
> for j in range(writer.GetNumberOfCols()):<br>
> data.SetTuple1(j, i + j + 100)<br>
> writer.PutNextRow(data)<br>
><br>
> #close the map<br>
> writer.CloseMap()<br>
><br>
> #now reopen the map for reading<br>
> reader = vtkGRASSRasterMapReader()<br>
> reader.OpenMap("test_cell")<br>
><br>
> #print the cells to stdout<br>
> print " "<br>
> for i in range(reader.GetNumberOfRows()):<br>
> data = reader.GetRow(i)<br>
> for j in range(reader.GetNumberOfCols()):<br>
> print data.GetTuple1(j),<br>
> print " "<br>
><br>
> #print the range of the raster map<br>
> val = [0,0]<br>
> reader.GetRange(val);<br>
> print "Range ", val<br>
><br>
> #print the history<br>
> print reader.GetHistory()<br>
> reader.CloseMap()<br>
><br>
> #############################################<br>
><br>
><br>
> Have fun<br>
> Soeren<br>
><br>
><br>
> ------------------------------------------------------------------------<br>
><br>
> _______________________________________________<br>
> grass-dev mailing list<br>
> <a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-dev</a><br>
<br>
<br>
--<br>
Benjamin Ducke<br>
Senior Applications Support and Development Officer<br>
<br>
Oxford Archaeological Unit Limited<br>
Janus House<br>
Osney Mead<br>
OX2 0ES<br>
Oxford, U.K.<br>
<br>
Tel: +44 (0)1865 263 800 (switchboard)<br>
Tel: +44 (0)1865 980 758 (direct)<br>
Fax :+44 (0)1865 793 496<br>
<a href="mailto:benjamin.ducke@oxfordarch.co.uk" target="_blank">benjamin.ducke@oxfordarch.co.uk</a><br>
<br>
<br>
<br>
<br>
------<br>
Files attached to this email may be in ISO 26300 format (OASIS Open Document Format). If you have difficulty opening them, please visit <a href="http://iso26300.info" target="_blank">http://iso26300.info</a> for more information.<br>
<br>
</blockquote></div></div></div><br>
</div><br>