[GRASS-user] 3D raster remove voxels above DEM

[Francois Chartier]

See attached screenshot, I was able to remove voxels above the DEM and
below bedrock surface using r3.out.vtk and selecting export as *pointdata
(p)* not voxels.  If exported as a Voxel from r3.out.vtk, the voxel file
did not open properly in Paraview could be due to size of file. I changed
option to export as points and it is working.
grid resolution is 5 x 5 x 1 m over 4km x 4 km x 65 m thickness approx.
(actual thickness varies between DEM and top of rock).

r3.out.vtk -p -s input=jan27551tens100 at Toronto output=C:\Users\Francois
Chartier\Documents\GeoModel\WorkingFile\Feb6_VTK_Pts top=GTADEM at Toronto
bottom=rock at PERMANENT

Le mar. 5 févr. 2019 à 21:20, Francois Chartier <fra.chartier at gmail.com> a
écrit :

> See attached screenshot, I was able to remove voxels above the DEM and
> below bedrock surface using r3.out.vtk and selecting export as *pointdata
> (p)* not voxels.  If exported as a Voxel from r3.out.vtk, the voxel file
> did not open properly in Paraview could be due to size of file. I changed
> option to export as points and it is working.
> grid resolution is 5 x 5 x 1 m over 4km x 4 km x 65 m thickness approx.
> (actual thickness varies between DEM and top of rock).
>
> r3.out.vtk -p -s input=jan27551tens100 at Toronto output=C:\Users\Francois
> Chartier\Documents\GeoModel\WorkingFile\Feb6_VTK_Pts top=GTADEM at Toronto
> bottom=rock at PERMANENT
>
>
> Le dim. 3 févr. 2019 à 16:42, Francois Chartier <fra.chartier at gmail.com>
> a écrit :
>
>> Hi,
>>
>> I would like to remove voxels located at a higher elevation than  DEM.
>> I have been trying r3.mapcalc without sucess.
>> Does the XY resolutions of both raster have to the same?
>> any ideas how to do this?
>> Thanks
>> Francois
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190205/09f5dd3b/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: overall3dview.jpg
Type: image/jpeg
Size: 166858 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190205/09f5dd3b/attachment-0001.jpg>