To compute a kriging, it is necessary the spatial covariance matrix. The spatial covariance could computed between to sets of points X and Y with different dimension and the result it is not necessarily a square matrix.
Kov(X, Y, Kern, symmetric = FALSE)
First set of spatial points.
Second set of spatial points.
Kernel function.
If result of computation will be a square matrix, the time computation can be improved setting this parameter to TRUE, the default is FALSE.
The spatial covariance matrix.
For a complete application you can check the documentation of Krig
.