[R-br] IDW

[Elias T. Krainski]

Hélio,

Se nao ha' dependencia espacial nao faz sentido interpolar. O IDW pode 
ser resumido numa conta matricial. Exemplo:

n = 30;
coo = cbind(runif(n)*2, runif(n))
y = rnorm(n)

x0 = seq(0, 1, length=101)
y0 = seq(0, 2, length=201)
grid = as.matrix(expand.grid(x0, y0))

p = 5
w = 1/(outer(grid[,1], coo[,1], '-')^2 +
            outer(grid[,2], coo[,2], '-')^2)^(p/2)
w = w/rowSums(w)

result = matrix(w%*%y, length(x0))

require(fields)
image.plot(list(x=x0, y=y0, z=result), asp=1)
points(coo, cex=0.1+(y-min(y))/diff(range(y)))

Elias.

On 28/05/14 21:21, Hélio Gallo Rocha wrote:
> Caros Listeiros.
>
> Estou utilizando o pacote geoR para análise geoestatística.
> Numa das áreas em estudo não apresentou dependência espacial.
> Sei que no pacote geosat contempla o método IDW.
> Alguém teria um CRM com este método?
>
> Agradeço desde já.
>
>
> -- 
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140529/a32eda8c/attachment.html>