# NOT RUN {
LASfile <- system.file("extdata", "Megaplot.laz", package="lidR")
las <- readLAS(LASfile)
plot(las)
# Outliers in intensity break the color range. Use the trim parameter.
col <- heat.colors(50)
plot(las, color = "Intensity", colorPalette = col)
plot(las, color = "Intensity", colorPalette = col, trim = 70)
# single file catalog using data provided in lidR
ctg = readLAScatalog(LASfile)
plot(ctg)
# }
Run the code above in your browser using DataLab