<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-2">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7233.28">
<TITLE>ODP: [GRASS-user] voxel from dems</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->

<P><FONT SIZE=2>hi<BR>
thanks Soren for your sugestions i will try r.to.rast3<BR>
best<BR>
maciej<BR>
<BR>
Maciej Tomaszczyk<BR>
Polish Geological Institute<BR>
Geological Mapping Department<BR>
<BR>
<BR>
<BR>
-----Wiadomość oryginalna-----<BR>
Od: Soeren Gebbert [<A HREF="mailto:soerengebbert@gmx.de">mailto:soerengebbert@gmx.de</A>]<BR>
Wysłano: Śr 2006-07-26 00:47<BR>
Do: grassuser@grass.itc.it<BR>
DW: Maciej Tomaszczyk<BR>
Temat: Re: [GRASS-user] voxel from dems<BR>
<BR>
Dear Maciej,<BR>
<BR>
On Tuesday 25 July 2006 21:47, Maciej Tomaszczyk wrote:<BR>
&gt; hello all<BR>
&gt;<BR>
&gt; i try to make a grid3d geological model from a drillhole data<BR>
&gt; first i generated a couple of dems which represent a geological bounduary,<BR>
&gt; (using Z values of geological bounduary from drillholes),<BR>
&gt; after that i used r3.out.vtk and everythink looked ok in paraview, but<BR>
&gt; i got a big problem with some geological structures wich are not more less planar<BR>
&gt; as salt diapir or overturned folds, its mean that in some drillhole profills more<BR>
&gt; than one Z value exist (which represent the same bounduary)&nbsp;<BR>
&gt; and i am not able to make a DEM using both Z values.<BR>
<BR>
&gt;There is no way to create such volume surfaces with r3.out.vtk + dems,<BR>
&gt;because of the limits of dems to have only one Z value per cell.<BR>
<BR>
&gt;But there are maybe some other ways:<BR>
<BR>
&gt;One way would be to modell the geological structures as 3d faces in<BR>
&gt;grass and export this with v.out.vtk. But this is very .... hard.<BR>
<BR>
&gt;The second way would be to modell the geological structures as raster<BR>
&gt;map slices at a specific hight and put them together with r.to.rast3.<BR>
&gt;take a look at:<BR>
&gt;<A HREF="http://www-pool.math.tu-berlin.de/~soeren/grass/modules/screenshots/r.to.rast3.png">http://www-pool.math.tu-berlin.de/~soeren/grass/modules/screenshots/r.to.rast3.png</A><BR>
<BR>
&gt;Now you can export this data with r3.out.vtk and use the threshold or isosurface<BR>
&gt;filter of Paraview to extract the geological structure. To use these Filter you will need<BR>
&gt;VTK point data (export with r3.otu.vtk or the Cell-Data-To-Point-Data filter of Paraview).<BR>
<BR>
&gt;For now there is no way to import the generated VTK isosufaces into grass.<BR>
&gt;And it would be nice if GRASS would have such a threshold or isosurface modules like VTK.<BR>
&gt;But this is just a wish... .<BR>
&gt;Well .... i have some ideas in my mind how to implement this. ;)<BR>
<BR>
&gt;Best regards<BR>
&gt;Soeren<BR>
<BR>
&gt; maybe some of you got the same problem<BR>
&gt; can you help me solve it<BR>
&gt; best<BR>
&gt; maciej&nbsp;<BR>
&gt;<BR>
&gt; Maciej Tomaszczyk<BR>
&gt; Polish Geological Institute<BR>
&gt; Geological Mapping Department<BR>
&gt;<BR>
&gt;<BR>
<BR>
</FONT>
</P>

</BODY>
</HTML>