# kriging

From lidR v2.0.0
by Jean-Romain Roussel

##### 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

##### 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.