← Back to branch summary

~ubuntu-branches/ubuntu/vivid/grass/vivid-proposed

« back to all changes in this revision

Viewing changes to scripts/r3.in.xyz/r3.in.xyz.html

  • Committer: Package Import Robot
  • Author(s): Bas Couwenberg
  • Date: 2015-02-20 23:12:08 UTC
  • mfrom: (8.2.6 experimental)
  • Revision ID: package-import@ubuntu.com-20150220231208-1u6qvqm84v430b10
Tags: 7.0.0-1~exp1
* New upstream release.
* Update python-ctypes-ternary.patch to use if/else instead of and/or.
* Drop check4dev patch, rely on upstream check.
* Add build dependency on libpq-dev to grass-dev for libpq-fe.h.
* Drop patches applied upstream, refresh remaining patches.
* Update symlinks for images switched from jpg to png.

Show diffs side-by-side

added added

removed removed

Lines of Context:
scripts/r3.in.xyz/r3.in.xyz.html
 
1
<h2>DESCRIPTION</h2>
 
2
 
 
3
<em>r3.in.xyz</em> imports sparse XYZ data from an ASCII file into
 
4
a 3D raster map (voxels). It does this by running the <em>r.in.xyz</em>
 
5
module multiple times for different z-ranges and then assembling the
 
6
slices with <em>r.to.rast3</em>.
 
7
<p>
 
8
See the <a href="r.in.xyz.html">r.in.xyz</a> help page for general
 
9
parameter usage and tips.
 
10
<p>
 
11
The map is created using the rows, columns, and depths set by
 
12
current region settings. Be sure to check and adjust these with
 
13
the <em>g.region</em> module before performing the import.
 
14
<p>
 
15
You may either use the z-value as the data value for the voxel
 
16
(e.g. with the 'n' statistic), or alternately scan another
 
17
column for the data values to bin into the voxels. This alternate
 
18
data column can be both filtered by range and have a scaling
 
19
factor applied to it.
 
20
 
 
21
 
 
22
<h2>NOTES</h2>
 
23
 
 
24
The 2D and 3D horizontal region resolutions must match. See the
 
25
EXAMPLES section below.
 
26
<p>
 
27
Unlike <em>r.in.xyz</em>, reading from stdin and z-scaling are not
 
28
possible. Filtering by z-range is accomplished by setting the 3D region.
 
29
<p>
 
30
To enable parallel processing support, set the <b>workers=</b> option
 
31
to match the number of CPUs or CPU-cores available on your system.
 
32
Alternatively, the <tt>WORKERS</tt> environment variable can be set
 
33
to the number of concurrent processes desired.
 
34
<p>
 
35
Points falling exactly on a vertical bound will belong to the depth
 
36
band below them, except for points exactly on the top bound, which will
 
37
belong to the top-most slice.
 
38
<p>
 
39
The script is expected to be nearly as efficient as if it was fully
 
40
written in C.
 
41
 
 
42
 
 
43
<h2>EXAMPLE</h2>
 
44
 
 
45
Using the Serpent Mound dataset. (see the
 
46
 <a href="http://grasswiki.osgeo.org/wiki/LIDAR">GRASS LiDAR wiki page</a>)
 
47
 
 
48
<div class="code"><pre>
 
49
  #scan dataset for extent:
 
50
  r3.in.xyz -s in=Serpent_Mound_Model_LAS_Data.txt out=dummy \
 
51
     x=1 y=2 z=3 separator=space
 
52
 
 
53
  # set the 2D and 3D regions:
 
54
  g.region n=4323641.57 s=4320942.61 w=289020.90 e=290106.02 res=1 -a
 
55
  g.region b=166 t=216 res3=1 tbres=5 -3 -p
 
56
 
 
57
  r3.in.xyz in=Serpent_Mound_Model_LAS_Data.txt out=serpent3D \
 
58
     method=mean x=1 y=2 z=3 separator=space type=float
 
59
</pre></div>
 
60
 
 
61
The same, but aggregate and store backscatter strength from column 5
 
62
into voxels in instead of the z-value:
 
63
<div class="code"><pre>
 
64
  r3.in.xyz in=Serpent_Mound_Model_LAS_Data.txt out=serpent3D.bakscat \
 
65
     method=mean x=1 y=2 z=3 val=5 separator=space type=float
 
66
</pre></div>
 
67
 
 
68
 
 
69
<h2>BUGS</h2>
 
70
 
 
71
<em>r.to.rast3</em> always creates a <tt>double</tt> output map
 
72
regardless of input.
 
73
 
 
74
 
 
75
<h2>SEE ALSO</h2>
 
76
 
 
77
<em>
 
78
<a href="g.region.html">g.region</a>,
 
79
<a href="r.in.xyz.html">r.in.xyz</a>,
 
80
<a href="r.to.rast3.html">r.to.rast3</a>
 
81
</em>
 
82
 
 
83
 
 
84
<h2>AUTHOR</h2>
 
85
 
 
86
Hamish Bowman<br>
 
87
<i>Dunedin, New Zealand</i>
 
88
 
 
89
<p>
 
90
<i>Last changed: $Date: 2015-01-25 18:56:33 +0100 (Sun, 25 Jan 2015) $</i>