<div dir="auto">Thanks Moritz<div dir="auto"><br></div><div dir="auto"><div dir="auto" style="font-family:sans-serif">My goal is to summarize the steps for 3d interpolation at the end of the process. it took me a year to get here and still have a few more steps to do which i will contact the User list for.  </div><div dir="auto" style="font-family:sans-serif">I am actually doing this as part of a 2nd master thesis with le Cnam in Paris as continuing education, so i will formalize the process.</div><div dir="auto" style="font-family:sans-serif"><br></div></div><div dir="auto">The voxels above the dem are not Null values in this interpolation as the vvolrst interpolates the data within the full 'footprint' of the g.region; is there a possibility to constrain the interpolation to in between the DEM and lower surface, by setting a more customized Region than a Block?</div><div dir="auto"><br></div><div dir="auto"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Feb 6, 2019 04:45, "Moritz Lennert" <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi François,<br>
<br>
Thanks for sharing the solution you found. From your first mail it was <br>
clear to me what these voxels above the DEM contained. NULL values ?<br>
<br>
It would be great if you could summarize all your experience working on <br>
3D data in GRASS GIS in the wiki, at least by amending <br>
<a href="https://grasswiki.osgeo.org/wiki/Help_with_3D" rel="noreferrer noreferrer" target="_blank">https://grasswiki.osgeo.org/wiki/Help_with_3D</a> with your additions or <br>
corrections ? This would be very helpful many others.<br>
<br>
Moritz<div class="quoted-text"><br>
<br>
On 6/02/19 03:46, Francois Chartier wrote:<br>
> See attached screenshot, I was able to remove voxels above the DEM and <br>
> below bedrock surface using r3.out.vtk and selecting export as <br></div>
> *pointdata (p)* not voxels.  If exported as a Voxel from r3.out.vtk, the <br><div class="quoted-text">
> voxel file did not open properly in Paraview could be due to size of <br>
> file. I changed option to export as points and it is working.<br>
> grid resolution is 5 x 5 x 1 m over 4km x 4 km x 65 m thickness approx. <br>
> (actual thickness varies between DEM and top of rock).<br>
> <br>
> r3.out.vtk -p -s input=jan27551tens100@Toronto output=C:\Users\Francois <br>
> Chartier\Documents\GeoModel\WorkingFile\Feb6_VTK_Pts top=GTADEM@Toronto <br>
> bottom=rock@PERMANENT<br>
> <br>
> Le mar. 5 févr. 2019 à 21:20, Francois Chartier <<a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">fra.chartier@gmail.com</a> <br></div>
> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">fra.chartier@gmail.com</a>>> a écrit :<div class="quoted-text"><br>
> <br>
>     See attached screenshot, I was able to remove voxels above the DEM<br>
>     and below bedrock surface using r3.out.vtk and selecting export as<br></div>
>     *pointdata (p)* not voxels.  If exported as a Voxel from r3.out.vtk,<div class="quoted-text"><br>
>     the voxel file did not open properly in Paraview could be due to<br>
>     size of file. I changed option to export as points and it is working.<br>
>     grid resolution is 5 x 5 x 1 m over 4km x 4 km x 65 m thickness<br>
>     approx. (actual thickness varies between DEM and top of rock).<br>
> <br>
>     r3.out.vtk -p -s input=jan27551tens100@Toronto<br>
>     output=C:\Users\Francois<br>
>     Chartier\Documents\GeoModel\WorkingFile\Feb6_VTK_Pts<br>
>     top=GTADEM@Toronto bottom=rock@PERMANENT<br>
> <br>
> <br>
>     Le dim. 3 févr. 2019 à 16:42, Francois Chartier<br></div>
>     <<a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">fra.chartier@gmail.com</a> <mailto:<a href="mailto:fra.chartier@gmail.com" target="_blank" rel="noreferrer">fra.chartier@gmail.com</a>>> a écrit :<div class="quoted-text"><br>
> <br>
>         Hi,<br>
> <br>
>         I would like to remove voxels located at a higher elevation<br>
>         than  DEM.<br>
>         I have been trying r3.mapcalc without sucess.<br>
>         Does the XY resolutions of both raster have to the same?<br>
>         any ideas how to do this?<br>
>         Thanks<br>
>         Francois<br>
> <br>
> <br></div>
> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org" target="_blank" rel="noreferrer">grass-user@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
> <br>
<br>
<br>
</blockquote></div><br></div>