[R-br] Salvar como Geotiff objeto krige.conv do geoR

[Éder Comunello]

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 em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]


Em 30 de dezembro de 2013 19:42, Éder Comunello
<comunello.eder em 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.html#a4660867
>
> --
> Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131231/91f0b263/attachment.html>