Spatial Interpolation Algorithm
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.
tin(..., extrapolate = knnidw(3, 1, 50))
... |
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. |
LASfile <- system.file("extdata", "Topography.laz", package="lidR")
las = readLAS(LASfile, filter = "-inside 273450 5274350 273550 5274450")
#plot(las)
dtm = grid_terrain(las, algorithm = tin())
plot(dtm, col = terrain.colors(50))
#plot_dtm3d(dtm)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.