lidR (version 4.1.0)

plot_3d: Add a spatial object to a point cloud scene

Description

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.

Usage

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

if (FALSE) {
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)
}