<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>To keep the voxels between the upper DEM and lower bedrock surface, I am trying to this by step using mapcalc.  as the DEM (GTADEM) is a 2d raster, and the 3d raster (jan27551ten100) is '3D', it seems that r3.mapcalc can only work with 3d raster as it does not recognize GTADEM.  the inverse happens when using r.mapcalc.  if mapcalc cannot work with a mix of 2d and 3d rasters, how can i go around this?  </div><div><br></div><div>r3.mapcalc LessThanDEM=jan27551tens100 <= GTADEM                                <br></div><div dir="ltr"><div dir="ltr"><div dir="ltr">Invalid map <GTADEM></div><div dir="ltr">Parse error</div><div dir="ltr">ERROR: parse error</div><div dir="ltr"><br></div><div dir="ltr"><div dir="ltr">r.mapcalc LessThanDEM=jan27551tens100 <= GTADEM                                 </div><div dir="ltr">Invalid map <jan27551tens100></div><div dir="ltr">Parse error</div><div dir="ltr">ERROR: parse error</div><div dir="ltr">(Fri Feb 01 22:31</div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mer. 30 janv. 2019 à 21:23, Francois Chartier <<a href="mailto:fra.chartier@gmail.com">fra.chartier@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I have a 3D raster showing an interpolation of soil texture using v.vol.rst.  The resulting raster is rectangular (with a resolution of 5x5x1 m), and i would like to remove cells above the DEM and below top of bedrock.  I looked at r3mapcalc, as well as tried r3.out.vtk selecting top and bottom surface (on first tab) with and without selecting -s.  </div><div>I am not sure if this is possible.  </div><div><br></div><div>Thanks,</div><div>Francois</div></div>
</blockquote></div>