kriging: Spatial Interpolation Algorithm

Description Usage Arguments See Also Examples

View source: R/algorithm-spi.r

Description

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

1
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

Examples

1
2
3
4
5
6
7
8
9
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)

Jean-Romain/lidR documentation built on Nov. 28, 2019, 2:15 a.m.