Hello,
It's not really clear from your message whether y1, y2, y3 are the y
part of x/y coordinates, and what is the exact structure of your data
frame. Without knowing this, the code below is to be adapted to your needs
##load the packages
library(gstat)
library(raster)
#Load the data that are needed
setwd("f:/Grids")
r2 <- raster("idgd_wgs.asc")
y.list <- list(y1<-runif(10), y2 <- runif(10), y3 <- runif(10))
da1<-data.frame(x=1:10, z=1:10)
for (i in 1:length(y.list) {
da.temp<-data.frame(da1, y=y.list[[i]])
IDW.temp <- gstat(id = "z", formula = z~1, locations = ~x+y,
data=da.temp, set=list(idp =2))
IDWRaster.temp <- interpolate(r2, IDW.temp)
#add here one or two command lines to save the objects you want to
save, or check
rm(da.temp, IDW.temp, IDWRaster.temp) #remove the temporary elements
(always safer)
}
HTH
Alex
Dear r users,
I have the following script that I would like to run repeatedly for
different values of y ( the list of which I created below).
y.list <- list(y1=runif(10), y2 <- runif(10), y3 <- runif(10))
and here is the script I want to run for different values of y.
############
x <- 1:10
y <- 1:10
z<- 1:10
da1<-data.frame(x,y,z)
***** this file is an aready existing ascii file that I use to form the
raster
library(raster)
setwd("f:/Grids")
r2 <- raster("idgd_wgs.asc")
library(raster)
library(gstat)
IDW <- gstat(id = "z", formula = z~1, locations = ~x+y, data=da1,
set=list(idp =2))
IDWRaster <- interpolate(r2, IDW)
Any suggestion appreciated.
Thanks.
minti
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
