1r3.in.xyz(1) GRASS GIS User's Manual r3.in.xyz(1)
2
3
4
6 r3.in.xyz - Create a 3D raster map from an assemblage of many coordi‐
7 nates using univariate statistics
8
10 raster3d, import, voxel, LIDAR, statistics, conversion, aggregation,
11 binning
12
14 r3.in.xyz
15 r3.in.xyz --help
16 r3.in.xyz [-sgi] input=name output=name [method=string]
17 [type=string] [separator=character] [x=integer] [y=integer]
18 [z=integer] [value_column=integer] [vrange=min,max]
19 [vscale=float] [percent=integer] [pth=integer] [trim=float]
20 [workers=integer] [--overwrite] [--help] [--verbose] [--quiet]
21 [--ui]
22
23 Flags:
24 -s
25 Scan data file for extent then exit
26
27 -g
28 In scan mode, print using shell script style
29
30 -i
31 Ignore broken lines
32
33 --overwrite
34 Allow output files to overwrite existing files
35
36 --help
37 Print usage summary
38
39 --verbose
40 Verbose module output
41
42 --quiet
43 Quiet module output
44
45 --ui
46 Force launching GUI dialog
47
48 Parameters:
49 input=name [required]
50 ASCII file containing input data
51
52 output=name [required]
53 Name for output raster map
54
55 method=string
56 Statistic to use for raster values
57 Options: n, min, max, range, sum, mean, stddev, variance,
58 coeff_var, median, percentile, skewness, trimmean
59 Default: mean
60
61 type=string
62 Storage type for resultant raster map
63 Options: float, double
64 Default: float
65
66 separator=character
67 Field separator
68 Special characters: pipe, comma, space, tab, newline
69 Default: pipe
70
71 x=integer
72 Column number of x coordinates in input file (first column is 1)
73 Default: 1
74
75 y=integer
76 Column number of y coordinates in input file
77 Default: 2
78
79 z=integer
80 Column number of z coordinates in input file
81 Default: 3
82
83 value_column=integer
84 Column number of data values in input file
85 If not given or set to 0, the data points’ z-values are used
86 Default: 0
87
88 vrange=min,max
89 Filter range for value column data (min,max)
90
91 vscale=float
92 Scaling factor to apply to value column data
93 Default: 1.0
94
95 percent=integer
96 Percent of map to keep in memory
97 Options: 1-100
98 Default: 100
99
100 pth=integer
101 pth percentile of the values
102 Options: 1-100
103
104 trim=float
105 Discard <trim> percent of the smallest and <trim> percent of the
106 largest observations
107 Options: 0-50
108
109 workers=integer
110 Number of parallel processes to launch
111 Options: 1-256
112 Default: 1
113
115 r3.in.xyz imports sparse XYZ data from an ASCII file into a 3D raster
116 map (voxels). It does this by running the r.in.xyz module multiple
117 times for different z-ranges and then assembling the slices with
118 r.to.rast3.
119
120 See the r.in.xyz help page for general parameter usage and tips.
121
122 The map is created using the rows, columns, and depths set by current
123 region settings. Be sure to check and adjust these with the g.region
124 module before performing the import.
125
126 You may either use the z-value as the data value for the voxel (e.g.
127 with the ’n’ statistic), or alternately scan another column for the
128 data values to bin into the voxels. This alternate data column can be
129 both filtered by range and have a scaling factor applied to it.
130
132 The 2D and 3D horizontal region resolutions must match. See the EXAM‐
133 PLES section below.
134
135 Unlike r.in.xyz, reading from stdin and z-scaling are not possible.
136 Filtering by z-range is accomplished by setting the 3D region.
137
138 To enable parallel processing support, set the workers= option to match
139 the number of CPUs or CPU-cores available on your system. Alterna‐
140 tively, the WORKERS environment variable can be set to the number of
141 concurrent processes desired.
142
143 Points falling exactly on a vertical bound will belong to the depth
144 band below them, except for points exactly on the top bound, which will
145 belong to the top-most slice.
146
147 The script is expected to be nearly as efficient as if it was fully
148 written in C.
149
151 Using the Serpent Mound dataset. (see the GRASS LiDAR wiki page)
152 #scan dataset for extent:
153 r3.in.xyz -s in=Serpent_Mound_Model_LAS_Data.txt out=dummy \
154 x=1 y=2 z=3 separator=space
155 # set the 2D and 3D regions:
156 g.region n=4323641.57 s=4320942.61 w=289020.90 e=290106.02 res=1 -a
157 g.region b=166 t=216 res3=1 tbres=5 -3 -p
158 r3.in.xyz in=Serpent_Mound_Model_LAS_Data.txt out=serpent3D \
159 method=mean x=1 y=2 z=3 separator=space type=float
160 The same, but aggregate and store backscatter strength from column 5
161 into voxels in instead of the z-value:
162 r3.in.xyz in=Serpent_Mound_Model_LAS_Data.txt out=serpent3D.bakscat \
163 method=mean x=1 y=2 z=3 val=5 separator=space type=float
164
166 r.to.rast3 always creates a double output map regardless of input.
167
169 g.region, r.in.xyz, r.to.rast3
170
172 Hamish Bowman
173 Dunedin, New Zealand
174
176 Available at: r3.in.xyz source code (history)
177
178 Main index | 3D raster index | Topics index | Keywords index | Graphi‐
179 cal index | Full index
180
181 © 2003-2020 GRASS Development Team, GRASS GIS 7.8.5 Reference Manual
182
183
184
185GRASS 7.8.5 r3.in.xyz(1)