lidR (version 3.0.3)

tin: Spatial Interpolation Algorithm

Description

This function is made to be used in grid_terrain or normalize_height. It implements an algorithm for spatial interpolation. Spatial interpolation is based on a Delaunay triangulation, which performs a linear interpolation within each triangle. There are usually a few points outside the convex hull, determined by the ground points at the very edge of the dataset, that cannot be interpolated with a triangulation. Extrapolation is done using the nearest neighbour approach.

Usage

tin(..., extrapolate = knnidw(1, 1, 50))

Arguments

...

unused

extrapolate

There are usually a few points outside the convex hull, determined by the ground points at the very edge of the dataset, that cannot be interpolated with a triangulation. Extrapolation is done using the nearest neighbour approach by default using knnidw.

See Also

Other spatial interpolation algorithms: knnidw(), kriging()

Examples

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

# plot(las)

dtm = grid_terrain(las, algorithm = tin())

plot(dtm, col = terrain.colors(50))
plot_dtm3d(dtm)
# }