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 <[email protected]>[email protected]> Dourados, MS - [22 16.5'S, 54 49'W] Em 30 de dezembro de 2013 19:42, Éder Comunello <[email protected]>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.html#a4660867 > > -- > Éder Comunello <c <[email protected]>[email protected]> > Dourados, MS - [22 16.5'S, 54 49'W] > >
_______________________________________________ R-br mailing list [email protected] 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.
