1r3.cross.rast(1) Grass User's Manual r3.cross.rast(1)
2
3
4
6 r3.cross.rast - Creates cross section 2D raster map from 3D raster map
7 based on 2D elevation map
8
10 raster3d, profile, raster, voxel
11
13 r3.cross.rast
14 r3.cross.rast --help
15 r3.cross.rast [-m] input=string elevation=string output=string
16 [--overwrite] [--help] [--verbose] [--quiet] [--ui]
17
18 Flags:
19 -m
20 Use 3D raster mask (if exists) with input map
21
22 --overwrite
23 Allow output files to overwrite existing files
24
25 --help
26 Print usage summary
27
28 --verbose
29 Verbose module output
30
31 --quiet
32 Quiet module output
33
34 --ui
35 Force launching GUI dialog
36
37 Parameters:
38 input=string [required]
39 Input 3D raster map for cross section
40
41 elevation=string [required]
42 2D elevation map used to create the cross section map
43
44 output=string [required]
45 Resulting cross section 2D raster map
46
48 This module creates a cross section 2D map from one 3D raster volume
49 map based on a 2D elevation map. It checks if the value of the eleva‐
50 tion map is located in the z-coordinate space of the 3D map. If so, the
51 3D voxel value for this position is transferred to the related cross
52 section output map cell, otherwise the NULL value is set.
53
54 How r3.cross.rast works
55
56
57 If the 2D and 3D region settings are different, the 2D resolution will
58 be adjust to the 3D resolution.
59
61 To create a cut plane elevation map use r.mapcalc. Some examples:
62
63 · To create a cut plane elevation map in x direction type
64 r.mapcalc "cutplane = col()*x",
65 x be the value for the elevation. If the range of col() is 1
66 ... 10, the elevation map has the range 1 ... 10 if x == 1 and
67 if x == 10 the range 10 ... 100
68
69 · To create a cut plane elevation map in y direction type
70 r.mapcalc "cutplane = row()*x",
71 x be the value for the elevation. If the range of col() is 1
72 ... 10, the elevation map has the range 1 ... 10 if x == 1 and
73 if x == 10 the range 10 ... 100
74
75 · The user can also make a cut in y and x direction with r.map‐
76 calc by using
77 r.mapcalc "cutplane = (row()+col())*x"
78
80 Simple Spearfish example
81 g.region -d
82 g.region res=150 res3=150 t=1000 b=0 tbres=100
83 # synthetic data, could be geological structures:
84 r3.mapcalc "map3d = sin(row())+sin(col())+sin(depth()*depth())"
85 #create a cutplane map
86 r.mapcalc "cutplane = col()*10"
87 #create the cross section map
88 r3.cross.rast input=map3d elevation=cutplane output=crosssection
89
91 g.region, r.mapcalc, r3.mapcalc, r3.to.rast
92
94 Sören Gebbert
95
96 Last changed: $Date: 2013-03-26 21:28:19 +0100 (Tue, 26 Mar 2013) $
97
99 Available at: r3.cross.rast source code (history)
100
101 Main index | 3D raster index | Topics index | Keywords index | Graphi‐
102 cal index | Full index
103
104 © 2003-2019 GRASS Development Team, GRASS GIS 7.4.4 Reference Manual
105
106
107
108GRASS 7.4.4 r3.cross.rast(1)