[R-sig-Geo] how to get weights for inverse distance weighting

Edzer Pebesma edzer.pebesma at uni-muenster.de
Sat Apr 8 00:35:08 CEST 2017


Scott,

you may want to have a look at gstat::idw0 , which expresses the idw
problem in a hand full of expressions, essentially

idw0 = function (data, newdata, y, idp = 2) {
    s = coordinates(data)
    s0 = coordinates(newdata)
    D = 1/(spDists(s0, s)^idp)
    sumD = apply(D, 1, sum)
    D %*% y/sumD
}


On 08/04/17 00:01, Waichler, Scott R wrote:
> Hello, I am trying to get inverse distance weights to estimate values on a regular grid from a set of data points, over a sequence of times.  The locations of the data points don't vary with time, but their values will with each instance, and in general some are NA.  I want to determine the number of unique weight matrices needed to do the IDW estimation across time.  I'm looking at gstat and spdep but am still having trouble with some of the basics . . .I should be more fluent with sp, I know.  The weight matrix has grid points for rows and data points for columns.
> 
> x.grid <- y.grid <- seq(0, 100, by=10)  # regular grid, locations I want to estimate z at
> N <- 10  # number of data points
> x.data <- runif(0, 100, n=N)  # set locations of the data points
> y.data <- runif(0, 100, n=N)
> 
> w <- list()  # 
> # Loop thru 20 instances (a sequence of months in my real problem)
> for(i in 1:20) {
>   # set the z.data values to something between 0 and 1
>   z.data <- runif(0, 1, n=N)
>   # set some of the z.data values to NA
>   ind.na <- sample.int(N, size=round(runif(1, 0, N)))
>   z.data[ind.na] <- NA
> 
>   # For this time, I want to get the matrix of weights (>= 0) that would be used in an inverse
>   # distance weighting interpolation.  I would like to specify parameters like nmax, nmin, omax, maxdist, and force
>   # as used in gstat().   Rows are grid points, columns are data points.
>   #w[[i]] <- matrix(data= ????, ncol=N, byrow=T)
> }
> 
> Finally, how can I determine the number of unique w[[i]] matrices?  I know for my toy problem it will
> in general be 20, but in my real problem it will be less than the number of times.
> 
> Thank you,
> Scott Waichler
> Pacific Northwest National Laboratory
> Richland, WA    USA
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma
Institute for Geoinformatics  (ifgi),  University of Münster
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software:   http://www.jstatsoft.org/
Computers & Geosciences:   http://elsevier.com/locate/cageo/

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 473 bytes
Desc: OpenPGP digital signature
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20170408/38bb272e/attachment.sig>


More information about the R-sig-Geo mailing list