kriging

0th

Percentile

Spatial Interpolation Algorithm

This function is made to be used in grid_terrain or lasground. 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. If NULL it performs an ordinary or weighted least squares prediction.

k

numeric. Number of k-nearest neighbours. Default 10.

See Also

Other spatial interpolation algorithms: knnidw, tin

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

# plot(las)

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

plot(dtm, col = terrain.colors(50))
plot_dtm3d(dtm)
# }
Documentation reproduced from package lidR, version 2.0.0, License: GPL-3

Community examples

Looks like there are no examples yet.