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 <[email protected]>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/index.htm > Phone: (612) 323 9898 > > > ________________________________ > From: Hélio Gallo Rocha <[email protected]> > To: r-br <[email protected]> > 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 > [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. > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ 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.
