Alexandre, Uma forma de fazer: ### <code r> require(geoR); require(rgdal); require(raster) # Predição por krigagem ordinária kloci <- expand.grid(seq(0,1,l=21), seq(0,1,l=21)) kcontr <- krige.control(cov.pars=c(1, .25)) kc <- krige.conv(s100, loc=kloci, krige=kcontr) image(kc, asp=1) pred <- cbind(kloci,kc$pred) ### data.frame com dados da predição coordinates(pred) <- ~Var1+Var2 ### transforma em SPointsDF gridded(pred) = TRUE ### transforma em SPixelsDF spplot(pred) writeGDAL(pred, 'teste.tif', drivername = 'GTiff') ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 30 de dezembro de 2013 19:42, Éder Comunello <comunello.eder@gmail.com>escreveu:
Alexandre,
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...
http://r-br.2285057.n4.nabble.com/R-br-de-cada-area-na-krigagem-td4660853.ht...
-- Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]