<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Olá pessoal..<div>Tenho outra dúvida..</div><div>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...</div><div><div>Minha dúvida é: como fazer a predição espacial para os pontos não amostrados?</div><div>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)</div><div><br></div><div>Segue rotina.. </div><div>require (geoR)</div><div>library(sp)</div><div>require(gstat)</div><div>require(splancs)</div><div>require(classInt)</div><div><br></div><div>#lendo os dados</div><div>h <- read.table("phca.txt ",head=T)</div><div>names(h)</div><div>class(h)</div><div>coordinates(h) = ~x+y</div><div><br></div><div>g = gstat(NULL, "ph", ph~1, h)</div><div>g = gstat(g, "ca", ca~1, h)</div><div>plot(variogram(g))</div><div>vm = variogram(g)</div><div>variogram(g)</div><div><br></div><div>##Ajuste do modelo esferico</div><div>sph <- fit.lmc(vm,g, vgm(0.3, "Sph", 200, 0.2))</div><div>sph</div><div><br></div><div>#modelo ajustado </div><div>plot(vm,sph)</div><div><br></div><div>##Predição para pontos não amostrados</div><div>#construção de uma malha fina</div><div>borda <- read.table("borda.txt") </div><div>borda </div><div>plot(borda)</div><div>polygon(borda)</div><div>apply(borda,2,range) </div><div>gr1 <- expand.grid(x=seq(239182.8, 241027.9, by=5),y=seq(7236744, 7238593, by=5))</div><div>gi <- polygrid(gr1, bor=borda)</div><div>plot(gi)</div><div><br></div><div>Att, Jacqueline</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div> </div></body>
</html>