.TH r3.in.lidar 1grass "" "GRASS 7.2.0" "Grass User's Manual"
.SH NAME
\fI\fBr3.in.lidar\fR\fR  \- Creates a 3D raster map from LAS LiDAR points
.SH KEYWORDS
3D raster, import, LIDAR
.SH SYNOPSIS
\fBr3.in.lidar\fR
.br
\fBr3.in.lidar \-\-help\fR
.br
\fBr3.in.lidar\fR [\-\fBdv\fR] \fBinput\fR=\fIname\fR \fBn\fR=\fIname\fR \fBsum\fR=\fIname\fR \fBmean\fR=\fIname\fR \fBproportional_n\fR=\fIname\fR \fBproportional_sum\fR=\fIname\fR  [\fBreturn_filter\fR=\fIstring\fR]   [\fBclass_filter\fR=\fIinteger\fR[,\fIinteger\fR,...]]   [\fBbase_raster\fR=\fIname\fR]   [\-\-\fBoverwrite\fR]  [\-\-\fBhelp\fR]  [\-\-\fBverbose\fR]  [\-\-\fBquiet\fR]  [\-\-\fBui\fR]
.SS Flags:
.IP "\fB\-d\fR" 4m
.br
Use base raster actual resolution instead of computational region
.IP "\fB\-v\fR" 4m
.br
Use only valid points
.br
Points invalid according to APSRS LAS specification will be filtered out
.IP "\fB\-\-overwrite\fR" 4m
.br
Allow output files to overwrite existing files
.IP "\fB\-\-help\fR" 4m
.br
Print usage summary
.IP "\fB\-\-verbose\fR" 4m
.br
Verbose module output
.IP "\fB\-\-quiet\fR" 4m
.br
Quiet module output
.IP "\fB\-\-ui\fR" 4m
.br
Force launching GUI dialog
.SS Parameters:
.IP "\fBinput\fR=\fIname\fR \fB[required]\fR" 4m
.br
LAS input file
.br
LiDAR input file in LAS format (*.las or *.laz)
.IP "\fBn\fR=\fIname\fR \fB[required]\fR" 4m
.br
Count of points per cell
.br
Name for output 3D raster map
.IP "\fBsum\fR=\fIname\fR \fB[required]\fR" 4m
.br
Sum of values of point intensities per cell
.br
Name for output 3D raster map
.IP "\fBmean\fR=\fIname\fR \fB[required]\fR" 4m
.br
Mean of point intensities per cell
.br
Name for output 3D raster map
.IP "\fBproportional_n\fR=\fIname\fR \fB[required]\fR" 4m
.br
3D raster map of proportional point count
.br
Point count per 3D cell divided by point count per vertical column
.IP "\fBproportional_sum\fR=\fIname\fR \fB[required]\fR" 4m
.br
3D raster map of proportional sum of values
.br
Sum of values per 3D cell divided by sum of values per vertical column
.IP "\fBreturn_filter\fR=\fIstring\fR" 4m
.br
Only import points of selected return type
.br
If not specified, all points are imported
.br
Options: \fIfirst, last, mid\fR
.IP "\fBclass_filter\fR=\fIinteger[,\fIinteger\fR,...]\fR" 4m
.br
Only import points of selected class(es)
.br
Input is comma separated integers. If not specified, all points are imported.
.IP "\fBbase_raster\fR=\fIname\fR" 4m
.br
Subtract raster values from the z coordinates
.br
The scale for z is applied beforehand, the filter afterwards
.SH DESCRIPTION
.PP
.PP
\fI
Figure: Proportional count of points per 3D cell. When 50% of all
points in a vertical column fall into a given 3D cell, the value
is 0.5. Here, the green color was assigned to 0.5, red to 1 and
yellow to 0. The figure shows vertical slices and green color
indicates high vegetation while red color indicates bare ground.
\fR
.SH NOTES
.RS 4n
.IP \(bu 4n
This module is highly experimental. Don\(cqt rely on its
functionality or interface. Please report issues on the mailing
list or in the bug tracker.
.IP \(bu 4n
No projection check or reprojection is performed, make sure you
are using data in the right projection for your GRASS Location.
.IP \(bu 4n
Selection of points according to return or class is not yet
supported.
.IP \(bu 4n
All outputs are currently mandatory.
.IP \(bu 4n
Some temporary maps are created but not cleaned up.
.IP \(bu 4n
Expects points to have intensity (random result for related
outputs when they don\(cqt).
.RE
.SH EXAMPLES
Set the region according to a 2D raster and adding 3D minimum
(bottom), maximum (top) and vertical (top\-bottom) resolution.
.br
.nf
\fC
g.region rast=secref b=80 t=160 tbres=5 \-p3
\fR
.fi
Now, \fIr3.in.lidar\fR will create the 3D raster of the size
given by the computation region:
.br
.nf
\fC
r3.in.lidar input=points.las n=points_n sum=points_sum \(rs
    mean=points_mean proportional_n=points_n_prop \(rs
    proportional_sum=points_sum_prop
\fR
.fi
.SS Point density vertical structure reduced to the terrain
Create ground raster:
.br
.nf
\fC
r.in.lidar input=points.las output=ground method=mean class_filter=2
\fR
.fi
Set vertical extent of computational region to (relative) coordinates
above ground:
.br
.nf
\fC
g.region rast=secref b=0 t=47 \-p3
\fR
.fi
Compute point density:
.br
.nf
\fC
r3.in.lidar input=points.las n=points_n sum=points_sum \(rs
    mean=points_mean proportional_n=points_n_prop \(rs
    proportional_sum=points_sum_prop \(rs
    base_raster=ground
\fR
.fi
.SH SEE ALSO
\fI
r3.in.xyz,
r.in.lidar,
v.in.lidar,
r.to.rast3,
r3.mapcalc,
g.region
\fR
.SH REFERENCES
.RS 4n
.IP \(bu 4n
V. Petras, A. Petrasova, J. Jeziorska, H. Mitasova (2016):
\fIProcessing UAV and lidar point clouds in GRASS GIS\fR.
XXIII ISPRS Congress 2016
[ISPRS Archives,
ResearchGate]
.IP \(bu 4n
ASPRS LAS format
.IP \(bu 4n
LAS library
.IP \(bu 4n
LAS library C API documentation
.RE
.SH AUTHOR
Vaclav Petras, NCSU GeoForAll Lab
.PP
\fILast changed: $Date: 2016\-08\-22 04:16:33 +0200 (Mon, 22 Aug 2016) $\fR
.SH SOURCE CODE
.PP
Available at: r3.in.lidar source code (history)
.PP
Main index |
3D raster index |
Topics index |
Keywords index |
Graphical index |
Full index
.PP
© 2003\-2016
GRASS Development Team,
GRASS GIS 7.2.0 Reference Manual