[R-br] raster e shape file

[Thiago V. dos Santos]

Hélio,

Algumas sugestões foram postadas por mim e pelo Paulo, que envolvem converter o resultado da krigagem em um raster usando as funções do sp ou do raster.

Certas perguntas podem ajudar na solução do problema: qual o nome do objeto produzido pela krigagem em que você está interessado? qual a classe desse objeto? onde estão as coordenadas dele? como converter um objeto dessa classe para um objeto da classe raster?

Seguindo o código que eu Alexandre postou recentemente, eu faria uma matriz ou dataframe a partir de lat lon e krig (baseado no output da krigagem) e montaria um raster usando a função rasterFromXYZ. Uma vez em raster, a conversão para polígono (shapefile) é direta. Veja código abaixo, que pode ser acoplado ao CMR da última mensagem do Alexandre:

#--------------------------------------

# First, add kriging values to the coordinates list
pred.grid$krig<-kc$predict

# Then, create a raster
r <- rasterFromXYZ(pred.grid)

# Last, check the result
plot(r)

# Converting to shapefile is straightforward
s <- rasterToPolygons(r)

# Plot shapefile (BEWARE - pretty slow)
plot(s)
#--------------------------------------

 

Saudações,
--
Thiago V. dos Santos
PhD student
Land and Atmospheric Science
University of Minnesota
http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm
Phone: (612) 323 9898


________________________________
 From: Hélio Gallo Rocha <heliogallorocha em gmail.com>
To: r-br <r-br em listas.c3sl.ufpr.br> 
Sent: Wednesday, April 10, 2013 2:27 PM
Subject: [R-br] raster e shape file
 


Caros,

Também estou as voltas com resultados de krigagen.
Estou acompanhando a msg do Alexandre Santos.

Já pesquisei bastante, mas não achei solução para transformar o resultado de uma krigagem em um raster ou num shape file.



Agradeço desde já

-- 
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho

_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
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.