plot_3d function

Add a spatial object to a point cloud scene

Add a raster (raster, stars terra) object that represents a digital terrain model or a SpatialPointsDataFrame or sf that represents tree tops to a point cloud scene. To add elements to a scene with a point cloud plotted with the function plot from lidR, the functions add_* take as first argument the output of the plot function (see examples), because the plot function does not plot the actual coordinates of the point cloud, but offset values. See function plot and its argument clear_artifacts for more details. It works only with rgl i.e. backend = "rgl" which is the default.

plot_dtm3d(dtm, bg = "black", clear_artifacts = TRUE, ...)

add_dtm3d(x, dtm, ...)

add_treetops3d(x, ttops, z = "Z", ...)

add_flightlines3d(x, flightlines, z = "Z", ...)

Arguments

• dtm: An object of the class RasterLayer or stars or SpatRaster
• bg: The color for the background. Default is black.
• clear_artifacts: logical. It is a known and documented issue that 3D visualisation with rgl displays artifacts. The points and lines are inaccurately positioned in the space and thus the rendering may look false or weird. This is because rgl computes with single precision float. To fix this, the objects are shifted to (0,0) to reduce the number of digits needed to represent their coordinates. The drawback is that the objects are not plotted at their actual coordinates.
• ...: Supplementary parameters for surface3d or spheres3d .
• x: The output of the function plot used with a LAS object.
• ttops: A SpatialPointsDataFrame or sf/sfc that contains tree tops coordinates.
• z: character. The name of the attribute that contains the height of the tree tops or of the flightlines. Only for XY geometries Ignored if the input have XYZ geometries
• flightlines: A SpatialPointsDataFrame or sf that contains flightlines coordinates.

Examples

## Not run:

LASfile <- system.file("extdata", "Topography.laz", package="lidR")
las <- readLAS(LASfile)

dtm <- rasterize_terrain(las, algorithm = tin())
ttops <- locate_trees(las, lmf(ws = 5))

plot_dtm3d(dtm)

x <- plot(las)
add_dtm3d(x, dtm)
add_treetops3d(x, ttops)

plot(las) |> add_dtm3d(dtm) |> add_treetops3d(ttops)
## End(Not run)