g.region raster=elevation.10m -p g.region res=200 res3=200 t=2000 b=0 tbres=20 -p # Write the values of raster map soils based on the elevation of elevation.10m # to the 3D map volev r.to.rast3elev in=soils elev=elevation.10m out=volev # Write the values of map soils based on the elevation of elevation.10m # to the 3D map volev_l and fill the lower cells with the soils map values r.to.rast3elev in=soils elev=elevation.10m out=volev_l -l # Write the values of map soils based on the elevation of elevation.10m # to the 3D map volev_u and fill the upper cells with the soils map values r.to.rast3elev in=soils elev=elevation.10m out=volev_u -u # Example with multiple elevation maps. ## first we need three support maps r.mapcalc expression="one = 1" r.mapcalc expression="two = 2" r.mapcalc expression="three = 3" ## Now we generate the new evelation maps r.mapcalc expression="elev_mid = elevation.10m - 500" r.mapcalc expression="elev_bottom = elevation.10m - 1000" ## Now fill the lower cells below the elevation maps with the values one, two and three r.to.rast3elev -l input=one,two,three elevation=elevation.10m,elev_mid,elev_bottom output=threelayer ## Export the map for visualization with paraview (http://www.paraview.org) ## By default the null value is -9999.99, we adjust it to 0.0 for ## better visualization r3.out.vtk null=0.0 input=threelayer output=/tmp/threelayer.vtk # Start paraview paraview --data=/tmp/threelayer.vtk # Note: First you need to choose the surface representation style and # then color by "threelayer" in paraview.
Available at: r.to.rast3elev source code (history)
Latest change: Monday Feb 27 13:45:25 2023 in commit: a82501dc85294207e8e58437c3ec7f1056465ed4
Main index | Raster index | Topics index | Keywords index | Graphical index | Full index
© 2003-2024 GRASS Development Team, GRASS GIS 8.5.0dev Reference Manual