This code generates high resolution surface models and other data from a LiDAR point cloud making use of FUSION LDV (available only for Windows) and the OsGeo/QGIS python environment.
Input data required:
- Lidar point cloud (currently NH - Nationell höjddatabas)
- Building vector foot print (polygon) in shapefile format
- Vector polygon representing domain (shapefile) [optional]
- Water, downloaded from Open Street Map - (Quick OSM plugin required) or as a polygon vector layer [optional]
- Elevated rail structures and bridges as polygon layers [optional]
Output data generated:
- DSM: Digital surface model (ground and building heights - masl)
- CDSM: Canopy Digital Surface Model (tall vegetation - magl)
- DEM: Digital Elevation Model (ground heights only - masl)
- Landcover (7-classes) [optional]
- LAI (Leaf Area Index) [optional]
Main file: makedsmfromlidar_loop.py
Remarks
- If you have .laz-files you need to make sure that FUSION can read .laz-files. Go to the FUSION-manual to see how.
- If you have .laz-files, you need to set the parameter updateLAS to 'yes' the first time you run the code. This will unpack your .laz-files to .las.