<div dir="ltr">Thiago, <div>Obrigado pelo retorno.</div><div><br></div><div>Consegui plotar, mas não consegui extrair os dados kc$predict com coordenadas para fazer uma recomendação de calagem</div><div>
Mas dá mais uma luz.</div><div><br></div><div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br></div><div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 10 de abril de 2013 22:40, Thiago V. dos Santos <span dir="ltr"><<a href="mailto:thi_veloso@yahoo.com.br" target="_blank">thi_veloso@yahoo.com.br</a>></span> escreveu:<br>


<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hélio,<br>
<br>
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.<br>
<br>
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?<br>



<br>
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:<br>



<br>
#--------------------------------------<br>
<br>
# First, add kriging values to the coordinates list<br>
pred.grid$krig<-kc$predict<br>
<br>
# Then, create a raster<br>
r <- rasterFromXYZ(pred.grid)<br>
<br>
# Last, check the result<br>
plot(r)<br>
<br>
# Converting to shapefile is straightforward<br>
s <- rasterToPolygons(r)<br>
<br>
# Plot shapefile (BEWARE - pretty slow)<br>
plot(s)<br>
#--------------------------------------<br>
<br>
 <br>
<br>
Saudações,<br>
--<br>
Thiago V. dos Santos<br>
PhD student<br>
Land and Atmospheric Science<br>
University of Minnesota<br>
<a href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm" target="_blank">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br>
Phone: <a href="tel:%28612%29%20323%209898" value="+16123239898" target="_blank">(612) 323 9898</a><br>
<br>
<br>
________________________________<br>
 From: Hélio Gallo Rocha <<a href="mailto:heliogallorocha@gmail.com" target="_blank">heliogallorocha@gmail.com</a>><br>
To: r-br <<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>><br>
Sent: Wednesday, April 10, 2013 2:27 PM<br>
Subject: [R-br] raster e shape file<br>
<div><div><br>
<br>
<br>
Caros,<br>
<br>
Também estou as voltas com resultados de krigagen.<br>
Estou acompanhando a msg do Alexandre Santos.<br>
<br>
Já pesquisei bastante, mas não achei solução para transformar o resultado de uma krigagem em um raster ou num shape file.<br>
<br>
<br>
<br>
Agradeço desde já<br>
<br>
--<br>
Hélio Gallo Rocha<br>
IFSULDEMINAS - Câmpus Muzambinho<br>
<br>
</div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div></div>