dtm_kriging: Spatial Interpolation Algorithm

In lidR: Airborne LiDAR Data Manipulation and Visualization for Forestry Applications

Spatial Interpolation Algorithm

Description

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.

Usage

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.

See Also

Other dtm algorithms: dtm_idw, dtm_tin

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)

lidR documentation built on June 8, 2025, 9:39 p.m.