Thiago, Obrigado pelo retorno. Consegui plotar, mas não consegui extrair os dados kc$predict com coordenadas para fazer uma recomendação de calagem Mas dá mais uma luz. -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho Em 10 de abril de 2013 22:40, Thiago V. dos Santos <thi_veloso@yahoo.com.br>escreveu:
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/inde... Phone: (612) 323 9898
________________________________ From: Hélio Gallo Rocha <heliogallorocha@gmail.com> To: r-br <r-br@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@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.
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho