Salvar como Geotiff objeto krige.conv do geoR

Boa tarde Pessoal, Gostaria de pegar o resultado de uma krigagem realizada com a função krige.conv() do pacote geoR e salvar como uma imagem Geotiff, abaixo fiz um CRM com a função writeGDAL(), sem sucesso, alguém teria alguma ideia para resolver isto, obrigado, require(geoR) loci <- expand.grid(seq(0,1,l=21), seq(0,1,l=21)) # Predição por krigagem ordinária kc <- krige.conv(s100, loc=loci, krige=krige.control(cov.pars=c(1, .25))) require(rgdal) geoTiff <- tempfile() writeGDAL(kc["predict"], geoTiff , drivername="GTiff") GDALinfo(tf) -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

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]

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]
participantes (2)
-
ASANTOS
-
Éder Comunello