Olá pessoal..

Tenho outra dúvida..

Construí o semivariograma cruzado, ajustei um modelo e quero fazer a predição para pontos que não foram amostrados para isso construí uma malha fina...

Minha dúvida é: como fazer a predição espacial para os pontos não amostrados?

Lembrando que são 2 variáveis, sendo que para predizer os valores de uma variável será utilizado os valores da outra variável (cokrigagem)

Segue rotina..

require (geoR)

library(sp)

require(gstat)

require(splancs)

require(classInt)

#lendo os dados

h <- read.table("phca.txt ",head=T)

names(h)

class(h)

coordinates(h) = ~x+y

g = gstat(NULL, "ph", ph~1, h)

g = gstat(g, "ca", ca~1, h)

plot(variogram(g))

vm = variogram(g)

variogram(g)

##Ajuste do modelo esferico

sph <- fit.lmc(vm,g, vgm(0.3, "Sph", 200, 0.2))

sph

#modelo ajustado

plot(vm,sph)

##Predição para pontos não amostrados

#construção de uma malha fina

borda <- read.table("borda.txt")

borda

plot(borda)

polygon(borda)

apply(borda,2,range)

gr1 <- expand.grid(x=seq(239182.8, 241027.9, by=5),y=seq(7236744, 7238593, by=5))

gi <- polygrid(gr1, bor=borda)

plot(gi)

Att, Jacqueline