# load lidar elevation data
ctg <- lidR::readLAScatalog("/path/to/your/laz_data/")
las <- lidR::readLAS(ctg, select = "xyz")
dem <- lidR::grid_canopy(las, res = 1, p2r())
# load lidar elevation data
ctg <- lidR::readLAScatalog("/path/to/your/laz_data/")
las <- lidR::readLAS(ctg, select = "xyz")
dem <- lidR::grid_canopy(las, res = 1, p2r())