Note: This document is for an older version of GRASS GIS that will be discontinued soon. You should upgrade, and read the current manual page.
Note: This addon document is for an older version of GRASS GIS that will be discontinued soon. You should upgrade your GRASS GIS installation, and read the current addon manual page.
The computational region of a LiDAR point file can be determined by scanning the file using the -s flag.
Optionally r.in.pdal creates an estimated vector footprint area map of the LAS file when using the footprint parameter (the footprint is generated by PDAL).
Since a new raster map is created during the binning, the binning of points depends on the computational region settings (extent and resolution) which is by default set to the extent of the LiDAR input file (see more about binning below). The resulting raster resolution can be specified with the parameter resolution.
r.in.pdal is designed for processing massive point cloud datasets, for example raw LiDAR or sidescan sonar swath data. It has been tested with large datasets.
For details concerning raster binning see the manual page of r.in.lidar.
# check metadata pdal info --summary lidar_raleigh_nc_spm_height_feet.las # scan extent and exit r.in.pdal input=lidar_raleigh_nc_spm_height_feet.las output=lidar_raleigh -s # scan extent (g.region style) and exit r.in.pdal input=lidar_raleigh_nc_spm_height_feet.las output=lidar_raleigh -s -g # n=228500 s=215000.01 w=633370.82 e=645000 t=558.87 b=88.5 # set computation region to this extent g.region n=228500 s=215000.01 w=633370.82 e=645000 -p # import while aligning pixel geometry to existing "elevation" 10m res. raster map # specifying EPSG manually because SRS information is missing in this LAS file r.in.pdal input=lidar_raleigh_nc_spm_height_feet.las raster_reference=elevation \ resolution=10 output=lidar_raleigh method=mean # optionally: footprint=lidar_raleigh_footprint # visualize d.mon wx0 g.list vector d.rast lidar_raleigh d.vect streets_wake # analyse differences between DEM and rasterized point cloud # LAS files come with height in US feet units r.mapcalc "diff = elevation - lidar_raleigh * 0.3048006096012192" r.univar -e diff
# pulling official PDAL docker image docker pull pdal/pdal # using PDAL docker as command with mounted data volume # (caution: the LAZ file has to be stored in the mounted folder! (here: $(pwd))) export pdal_docker="docker run --rm -v $(pwd):/data -t pdal/pdal pdal" # scan extent and exit r.in.pdal input=/data/simple.laz output=lidar -s pdal_cmd="$pdal_docker" # scan extent (g.region style) and exit r.in.pdal input=/data/simple.laz output=lidar -s -g pdal_cmd="$pdal_docker" # n=853535.43 s=848899.7 w=635619.85 e=638982.55 t=586.38 b=406.59 # set computation region to this extent g.region n=853535.43 s=848899.7 w=635619.85 e=638982.55 -p # import data r.in.pdal input=/data/simple.laz output=lidar_perc95 method=percentile pth=95 pdal_cmd="$pdal_docker" r.univar lidar_perc95
Documentation: Markus Neteler, mundialis GmbH & Co. KG
Available at: r.in.pdal source code (history)
Latest change: Monday Jun 28 07:54:09 2021 in commit: 1cfc0af029a35a5d6c7dae5ca7204d0eb85dbc55
Note: This document is for an older version of GRASS GIS that will be discontinued soon. You should upgrade, and read the current manual page.
Note: This addon document is for an older version of GRASS GIS that will be discontinued soon. You should upgrade your GRASS GIS installation, and read the current addon manual page.
Main index | Raster index | Topics index | Keywords index | Graphical index | Full index
© 2003-2023 GRASS Development Team, GRASS GIS 7.8.9dev Reference Manual