[R-br] IDW

[Elias T. Krainski]

Corrigindo os limites x-y:

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

x0 = seq(0, 2, length=201)
y0 = seq(0, 1, length=101)
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)))

On 29/05/14 08:21, Elias T. Krainski wrote:
> 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.
>
>
>
> _______________________________________________
> 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/ce15c478/attachment.html>