Dear R-users, I am working on interpolating the station level temperature data to farm level data. I have z vector consisting of station level temperature observations and my x and y are latitude and longitude corresponding to a farm. My understanding is I can use raster combined with tps. While I am clear with the tps bit, I am not sure how I can construct the raster with teh data I have. Here is the reproducable example I made
Many thanks in advance regards, Mintewab library(fields) x <-1:20 y<- runif(20) z<- c(11, 15, 17, 2, 18, 6, 7, NA, 12, 10,21, 25, 27, 12, 28, 16, 17, NA, 12, 10) mydataset<-data.frame(z, y, z) mydataset[complete.cases(mydataset),] tpsfit <- Tps(cbind(x, y), z, scale.type="unscaled") library(raster) r <- raster(system.file("external/test.grd", package="raster")) p <- raster(r) p <- interpolate(p, tpsfit) p <- mask(p, r) plot(p) se <- interpolate(p, tpsfit, fun=predict.se) se <- mask(se, r) plot(se) ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.