1r.to.rast3elev(1) GRASS GIS User's Manual r.to.rast3elev(1)
2
3
4
6 r.to.rast3elev - Creates a 3D volume map based on 2D elevation and
7 value raster maps.
8
10 raster, conversion, raster3d, voxel
11
13 r.to.rast3elev
14 r.to.rast3elev --help
15 r.to.rast3elev [-ulm] input=name[,name,...] elevation=name[,name,...]
16 output=name [upper=float] [lower=float] [tilesize=integer]
17 [--overwrite] [--help] [--verbose] [--quiet] [--ui]
18
19 Flags:
20 -u
21 Use the input map values to fill the upper cells
22
23 -l
24 Use the input map values to fill the lower cells
25
26 -m
27 Use 3D raster mask (if exists) with input map
28
29 --overwrite
30 Allow output files to overwrite existing files
31
32 --help
33 Print usage summary
34
35 --verbose
36 Verbose module output
37
38 --quiet
39 Quiet module output
40
41 --ui
42 Force launching GUI dialog
43
44 Parameters:
45 input=name[,name,...]Â [required]
46 Name of input raster map(s)
47
48 elevation=name[,name,...]Â [required]
49 Name of input elevation raster map(s)
50
51 output=name [required]
52 Name for output 3D raster map
53
54 upper=float
55 The value to fill the upper cells, default is null
56
57 lower=float
58 The value to fill the lower cells, default is null
59
60 tilesize=integer
61 The maximum tile size in kilo bytes. Default is 32KB.
62 Default: 32
63
65 Creates a 3D volume map based on 2D elevation and value raster maps.
66 If the 2d and 3d region settings are different, the 2d resolution will
67 be adjust to the 3d resolution.
68
69 How r.to.rast3elev works
70
71
73 The hight of the 2D elevation maps will be used to verify the position
74 within the 3D region. If the cell value of the elevation raster maps is
75 located within the 3D region, the cell value of the appropriate 2D in‐
76 put raster maps will be written to the associated 3D cell. There are
77 flags and options to fill the upper and lower 3D cells with a specific
78 value, or the input raster maps values.
79
81 Simple Spearfish sample data set example
82 g.region raster=elevation.10m -p
83 g.region res=200 res3=200 t=2000 b=0 tbres=20 -p
84 # Write the values of raster map soils based on the elevation of elevation.10m
85 # to the 3D map volev
86 r.to.rast3elev in=soils elev=elevation.10m out=volev
87 # Write the values of map soils based on the elevation of elevation.10m
88 # to the 3D map volev_l and fill the lower cells with the soils map values
89 r.to.rast3elev in=soils elev=elevation.10m out=volev_l -l
90 # Write the values of map soils based on the elevation of elevation.10m
91 # to the 3D map volev_u and fill the upper cells with the soils map values
92 r.to.rast3elev in=soils elev=elevation.10m out=volev_u -u
93 # Example with multiple elevation maps.
94 ## first we need three support maps
95 r.mapcalc expression="one = 1"
96 r.mapcalc expression="two = 2"
97 r.mapcalc expression="three = 3"
98 ## Now we generate the new evelation maps
99 r.mapcalc expression="elev_mid = elevation.10m - 500"
100 r.mapcalc expression="elev_bottom = elevation.10m - 1000"
101 ## Now fill the lower cells below the elevation maps with the values one, two and three
102 r.to.rast3elev -l input=one,two,three elevation=elevation.10m,elev_mid,elev_bottom output=threelayer
103 ## Export the map for visualization with paraview (http://www.paraview.org)
104 ## By default the null value is -9999.99, we adjust it to 0.0 for
105 ## better visualization
106 r3.out.vtk null=0.0 input=threelayer output=/tmp/threelayer.vtk
107 # Start paraview
108 paraview --data=/tmp/threelayer.vtk
109 # Note: First you need to choose the surface representation style and
110 # then color by "threelayer" in paraview.
111
113 r.to.rast3, r3.cross.rast, g.region
114
116 Soeren Gebbert
117
119 Available at: r.to.rast3elev source code (history)
120
121 Accessed: Saturday Jan 21 20:39:13 2023
122
123 Main index | Raster index | Topics index | Keywords index | Graphical
124 index | Full index
125
126 © 2003-2023 GRASS Development Team, GRASS GIS 8.2.1 Reference Manual
127
128
129
130GRASS 8.2.1 r.to.rast3elev(1)