Rasterize LAS/LAZ point clouds with filtering/compositing options

SYNOPSIS

Usage: pklas2img -i input.las -o output.tif

Options: [-n attribute] [-comp method] [-fir type] [-a_srs] [-ulx value -uly value -lrx value -lry value] [-dx value -dy value] [-ot type] [-of format] [-ret value]* [-class number]*

Advanced options: [-nbin value] [-nodata value] [-co option]* [-ct colortable]

Description

The utility pklas2img converts a las/laz point cloud into a gridded raster dataset. The implementation is based on 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 -comp. The default attribute is z (heiht), but can also be intensity (if available), scan angle rank (-n angle), the return number (-n return) or the total number of returns in that grid cell (-n nreturn). To select specific returns only, set the option -fir (first, last, single, multiple, or all).

Options

• use either -short or --long options (both --long=value and --long value are supported)
• short option -h shows basic options only, long option --help shows all options

  short	long	type	default	description
  i	input	std::string		Input las file
  n	name	std::string	z	names of the attribute to select: intensity, angle, return, nreturn, spacing, z
  ret	ret	unsigned short		number(s) of returns to include
  class	class	unsigned short		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)
  comp	comp	std::string	last	composite for multiple points in cell (min, max, absmin, absmax, median, mean, sum, first, last, profile (percentile height values), percentile, number (point density)). Last: overwrite cells with latest point
  fir	filter	std::string	all	filter las points (first,last,single,multiple,all).
  angle_min	angle_min	unsigned short		minimum scan angle to read points.
  angle_max	angle_max	unsigned short		maximum scan angle to read points.
  o	output	std::string		Output image file
  a_srs	a_srs	std::string		assign the projection for the output file in epsg code, e.g., epsg:3035 for European LAEA projection
  ulx	ulx	double	0	Upper left x value bounding box (in geocoordinates if georef is true). 0 is read from input file
  uly	uly	double	0	Upper left y value bounding box (in geocoordinates if georef is true). 0 is read from input file
  lrx	lrx	double	0	Lower right x value bounding box (in geocoordinates if georef is true). 0 is read from input file
  lry	lry	double	0	Lower right y value bounding box (in geocoordinates if georef is true). 0 is read from input file
  ot	otype	std::string	Byte	Data type for output image ({Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/CInt16/CInt32/CFloat32/CFloat64}). Empty string: inherit type from input image
  of	oformat	std::string	GTiff	Output image format (see also gdal_translate).
  dx	dx	double	1	Output resolution in x (in meter)
  dy	dy	double	1	Output resolution in y (in meter)
  nbin	nbin	short	10	Number of percentile bins for calculating percentile height value profile (=number of output bands)
  perc	perc	double	95	Percentile value used for rule percentile
  nodata	nodata	short	0	nodata value to put in image
  co	co	std::string		Creation option for output file. Multiple options can be specified.
  ct	ct	std::string		color table (file with 5 columns: id R G B ALFA (0: transparent, 255: solid)

  pklas2img -i lasfile -o output