https://github.com/cran/fields
Tip revision: 6c8b30169bba182a68765ee3cb9b4e2ef7d38332 authored by Doug Nychka on 16 November 2011, 00:00:00 UTC
version 6.6.3
version 6.6.3
Tip revision: 6c8b301
fastTps.R
# fields, Tools for spatial data
# Copyright 2004-2011, Institute for Mathematics Applied Geosciences
# University Corporation for Atmospheric Research
# Licensed under the GPL -- www.gpl.org/licenses/gpl.html
"fastTps" <- function(x, Y, m = NULL, p = NULL, theta,
lon.lat = FALSE, ...) {
x <- as.matrix(x)
d <- ncol(x)
if (is.null(p)) {
if (is.null(m)) {
m <- max(c(2, ceiling(d/2 + 0.1)))
}
p <- (2 * m - d)
if (p <= 0) {
stop(" m is too small you must have 2*m -d >0")
}
}
# special arguments to send to the wendland covariance/taper function.
# see nearest.dist for some explanation of 'method'
cov.args <- list(k = p, Dist.args = list(method = ifelse(!lon.lat,
"euclidean", "greatcircle")))
Tpscall <- match.call()
Tpscall$cov.function <- "Compactly supported Wendland function "
mKrig(x, Y, cov.function = "wendland.cov", m = m, cov.args = cov.args,
theta = theta, find.trA = TRUE, ...)
}