dtm_kriging function

Spatial Interpolation Algorithm

This function is made to be used in rasterize_terrain or normalize_height . It implements an algorithm for spatial interpolation. Spatial interpolation is based on universal kriging using the krige() function from gstat. This method combines the KNN approach with the kriging approach. For each point of interest it kriges the terrain using the k-nearest neighbour ground points. This method is more difficult to manipulate but it is also the most advanced method for interpolating spatial data.

kriging(model = gstat::vgm(0.59, "Sph", 874), k = 10L)

Arguments

• model: A variogram model computed with vgm() from package gstat. If NULL it performs an ordinary or weighted least squares prediction.
• k: numeric. Number of k-nearest neighbours. Default 10.

Examples

## Not run:

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

plot(las)

dtm = rasterize_terrain(las, algorithm = kriging())

plot(dtm)
plot_dtm3d(dtm)
## End(Not run)

See Also

Other dtm algorithms: dtm_idw, dtm_tin