library(sp)
set.seed(1)
d <- data.frame(lon = runif(n=500,min = -179, max = 179),
lat = runif(n=500,min = -90, max = 90),
z = rnorm(500))
coordinates(d) <- ~lon + lat
slot(d, "proj4string") = CRS("+proj=longlat")
### Now create basis functions on sphere
G <- auto_basis(manifold = sphere(),data=d,
nres = 2,prune=15,
type = "bisquare",
subsamp = 20000)
nres(G)
Run the code above in your browser using DataLab