<div dir="ltr"><div>Alexandre, </div><div><br></div><div>Uma forma de fazer:</div><div><br></div><div><font face="courier new, monospace">### <code r></font></div><div><font face="courier new, monospace">require(geoR); require(rgdal); require(raster)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"># Predição por krigagem ordinária </font></div><div><font face="courier new, monospace">kloci   <- expand.grid(seq(0,1,l=21), seq(0,1,l=21))</font></div>
<div><font face="courier new, monospace">kcontr <- krige.control(cov.pars=c(1, .25))</font></div><div><font face="courier new, monospace">kc <- krige.conv(s100, loc=kloci, krige=kcontr) </font></div><div><font face="courier new, monospace">image(kc, asp=1)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">pred <- cbind(kloci,kc$pred)     ### data.frame com dados da predição</font></div><div><font face="courier new, monospace">coordinates(pred) <- ~Var1+Var2  ### transforma em SPointsDF</font></div>
<div><font face="courier new, monospace">gridded(pred) = TRUE             ### transforma em SPixelsDF</font></div><div><font face="courier new, monospace">spplot(pred)</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">writeGDAL(pred, 'teste.tif', drivername = 'GTiff')</font></div><div><font face="courier new, monospace">### </code></font></div></div><div class="gmail_extra">
<br clear="all"><div><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br>
</div></div>
<br><br><div class="gmail_quote">Em 30 de dezembro de 2013 19:42, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Alexandre,<br><br>Acredito que seja necessário transformar o objeto da predição em SPDF ou raster primeiro. O exemplo que postei no link que segue pode ajudar...<br><br><a href="http://r-br.2285057.n4.nabble.com/R-br-de-cada-area-na-krigagem-td4660853.html#a4660867" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-de-cada-area-na-krigagem-td4660853.html#a4660867</a><span class="HOEnZb"><font color="#888888"><br>

<br>-- <br><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br>

</div><br>
</font></span></blockquote></div><br></div>