'\" -*- coding: UTF-8 -*-
.if \n(.g .ds T< \\FC
.if \n(.g .ds T> \\F[\n[.fam]]
.de URL
\\$2 \(la\\$1\(ra\\$3
..
.if \n(.g .mso www.tmac
.TH pklas2img 1 "11 June 2020" "" ""
.SH NAME
pklas2img \- Rasterize LAS/LAZ point clouds with filtering/compositing options
.SH SYNOPSIS
'nh
.fi
.ad l
\fBpklas2img\fR \kx
.if (\nx>(\n(.l/2)) .nr x (\n(.l/5)
'in \n(.iu+\nxu
\fB-i\fR \fIlasfile\fR \fB-o\fR \fIoutput\fR [\fIoptions\fR] [\fIadvanced options\fR]
'in \n(.iu-\nxu
.ad b
'hy
.SH DESCRIPTION
\fBpklas2img\fR converts a las/laz point cloud into a
gridded raster dataset.
The implementation is based on
.URL http://www.liblas.org/ liblas
API.
You can define the bounding box, grid cell size and spatial reference
set.
The composite rule for multiple returns within a single grid cell can be
set with the option \*(T<\fB\-comp\fR\*(T>.
The default attribute is z (height), but can also be intensity
(if available), the scan angle rank (\*(T<\fB\-n\fR\*(T> angle), the
return number (\*(T<\fB\-n\fR\*(T> return) or the total number of returns
in that grid cell (\*(T<\fB\-n\fR\*(T> nreturn).
To select specific returns only, set the option \*(T<\fB\-fir\fR\*(T>
(first, last, single, multiple, or all).
.SH OPTIONS
.TP 
\*(T<\fB\-i\fR\*(T> \fIfilename\fR, \*(T<\fB\-\-input\fR\*(T> \fIfilename\fR
Input las file
.TP 
\*(T<\fB\-o\fR\*(T> \fIfilename\fR, \*(T<\fB\-\-output\fR\*(T> \fIfilename\fR
Output image file
.TP 
\*(T<\fB\-n\fR\*(T> \fIattribute\fR, \*(T<\fB\-\-name\fR\*(T> \fIattribute\fR
names of the attribute to select: intensity, angle, return, nreturn, spacing, z
.TP 
\*(T<\fB\-comp\fR\*(T> \fIcomposite\fR, \*(T<\fB\-\-comp\fR\*(T> \fIcomposite\fR
composite for multiple points in cell (min, max, absmin, absmax,
median, mean, sum, first, last, profile (percentile height values),
number (point density)).
Last: overwrite cells with latest point
.TP 
\*(T<\fB\-fir\fR\*(T> \fItype\fR, \*(T<\fB\-\-filter\fR\*(T> \fItype\fR
filter las points (first,last,single,multiple,all).
.TP 
\*(T<\fB\-angle_min\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-angle_min\fR\*(T> \fIvalue\fR
minimum scan angle to read points.
.TP 
\*(T<\fB\-angle_max\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-angle_max\fR\*(T> \fIvalue\fR
maximum scan angle to read points.
.TP 
\*(T<\fB\-a_srs\fR\*(T> \fIEPSG:number\fR, \*(T<\fB\-\-a_srs\fR\*(T> \fIEPSG:number\fR
assign the projection for the output file in epsg code,
e.g., epsg:3035 for European LAEA projection
.TP 
\*(T<\fB\-ulx\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-ulx\fR\*(T> \fIvalue\fR
Upper left x value bounding box
(in geocoordinates if georef is true).
0 is read from input file
.TP 
\*(T<\fB\-uly\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-uly\fR\*(T> \fIvalue\fR
Upper left y value bounding box
(in geocoordinates if georef is true).
0 is read from input file
.TP 
\*(T<\fB\-lrx\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-lrx\fR\*(T> \fIvalue\fR
Lower right x value bounding box
(in geocoordinates if georef is true).
0 is read from input file
.TP 
\*(T<\fB\-lry\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-lry\fR\*(T> \fIvalue\fR
Lower right y value bounding box
(in geocoordinates if georef is true).
0 is read from input file
.TP 
\*(T<\fB\-dx\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-dx\fR\*(T> \fIvalue\fR
Output resolution in x (in meter)
.TP 
\*(T<\fB\-dy\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-dy\fR\*(T> \fIvalue\fR
Output resolution in y (in meter)
.TP 
\*(T<\fB\-ot\fR\*(T> \fItype\fR, \*(T<\fB\-\-otype\fR\*(T> \fItype\fR
Data type for output image
({Byte / Int16 / UInt16 / UInt32 / Int32 / Float32 / Float64 / CInt16 / CInt32 / CFloat32 / CFloat64}).
Empty string: inherit type from input image
.TP 
\*(T<\fB\-of\fR\*(T> \fIout_format\fR, \*(T<\fB\-\-oformat\fR\*(T> \fIout_format\fR
Output image format (see also
\fBgdal_translate\fR(1)).
.TP 
\*(T<\fB\-ret\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-ret\fR\*(T> \fIvalue\fR
number(s) of returns to include
.TP 
\*(T<\fB\-class\fR\*(T> \fInumber\fR, \*(T<\fB\-\-class\fR\*(T> \fInumber\fR
classes to keep:
0 (created, never classified),
1 (unclassified),
2 (ground),
3 (low vegetation),
4 (medium vegetation),
5 (high vegetation),
6 (building),
7 (low point, noise),
8 (model key-point),
9 (water),
10 (reserved),
11 (reserved),
12 (overlap)
.TP 
\*(T<\fB\-v\fR\*(T> \fIlevel\fR, \*(T<\fB\-\-verbose\fR\*(T> \fIlevel\fR
verbose mode
.PP
Advanced options
.TP 
\*(T<\fB\-nbin\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-nbin\fR\*(T> \fIvalue\fR
Number of percentile bins for calculating percentile height value
profile (=number of output bands)
.TP 
\*(T<\fB\-nodata\fR\*(T> \fIvalue\fR, \*(T<\fB\-\-nodata\fR\*(T> \fIvalue\fR
nodata value to put in image if not valid
.TP 
\*(T<\fB\-co\fR\*(T> \fIoption\fR, \*(T<\fB\-\-co\fR\*(T> \fIoption\fR
Creation option for output file.
Multiple options can be specified.
.TP 
\*(T<\fB\-ct\fR\*(T> \fIfilename\fR, \*(T<\fB\-\-ct\fR\*(T> \fIfilename\fR
colour table in ASCII format having 5 columns:
id R G B ALFA (0: transparent, 255: solid)