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@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.
_______________________________________________ R-br mailing list R-br@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.