Spatial Interpolation Algorithm
This function is made to be used in grid_terrain or normalize_height. It implements an algorithm for spatial interpolation. Interpolation is done using a k-nearest neighbour (KNN) approach with an inverse-distance weighting (IDW).
knnidw(k = 10, p = 2, rmax = 50)
k |
integer. Number of k-nearest neighbours. Default 10. |
p |
numeric. Power for inverse-distance weighting. Default 2. |
rmax |
numeric. Maximum radius where to search for knn. Default 50. |
LASfile <- system.file("extdata", "Topography.laz", package="lidR") las = readLAS(LASfile) #plot(las) dtm = grid_terrain(las, algorithm = knnidw(k = 6L, p = 2)) #plot(dtm, col = terrain.colors(50)) #plot_dtm3d(dtm)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.