# NOT RUN {
LASfile <- system.file("extdata", "Topography.laz", package="lidR")
lidar = readLAS(LASfile)
plot(lidar)
# --- First option: compute the DTM with grid_terrain -----
dtm = grid_terrain(lidar, method = "delaunay")
lidar_norm = lasnormalize(lidar, dtm)
plot(dtm)
plot(lidar_norm)
# }
# NOT RUN {
# --- Second option: read the DTM from a file -----
dtm = raster::raster(terrain.tiff)
lidar_norm = lidar - dtm
plot(lidar_norm)
# }
Run the code above in your browser using DataLab