<div dir="ltr">Alexandre<div><br></div><div><br></div><div style>Não sei se já tinha este resultado com valores de kc$predict, mas consegui colocar numa tabela as coordenadas(var1 e var2) e o kc$predict.</div><div style><br>
</div><div style>Não tive muito tempo de trabalhar nisso, mas o ideal era exportar para shape, mas não consegui.</div><div style><br></div><div style>Quem puder ajudar, se não teremos de lançar mão em soft pago para gerar um shape.</div>
<div style><br></div><div style>ai vai o CRM:</div><div style><br></div><div style><div>require(geoR)</div><div>library(raster)</div><div>ml <- likfit(s100, ini = c(1,0.5), fix.nugget = T)</div><div><br></div><div># defining the grid</div>
<div>pred.grid <- expand.grid(seq(0,1, l=100), seq(0,1, l=100))# ;pred.grid</div><div>#pred.grid</div><div># kriging calculations</div><div>kc<- krige.conv(s100, loc = pred.grid, krige = krige.control(obj.m = ml))</div>
<div><br></div><div>#recorte=matrix(c(.2,.8,.6,0,.1,.2,.7,.4),ncol=2) #c(0.2, 0.6, 0.2, 1.1, 0.2, 0.3, 1.0, 1.1, ncol=2)</div><div>#recorte</div><div>#plot(recorte)</div><div>ima=image(kc, loc = pred.grid,col=gray(seq(1,0.1,l=10)))</div>
<div>##tutorial do paulo para incrementar o mapa sde krigagem</div><div>contour(kc)</div><div>image(kc)</div><div>contour(kc, add=TRUE, nlev=21)</div><div>persp(kc, theta=20, phi=20)</div><div>contour(kc, filled=TRUE)</div>
<div>contour(kc, filled=TRUE, color=terrain.colors)</div><div>contour(kc, filled=TRUE, col=gray(seq(1,0,l=21)))</div><div>coordinates(pred.grid) <- ~ Var1+Var2;plot(pred.grid)## Transformando as coordenadas do grid</div>
<div>#fim</div><div>kc$predict<-as.data.frame(kc$predict)### Valores dos pixels preditos</div><div>kc_output<-SpatialPixelsDataFrame(points=pred.grid, data=kc$predict)### Criando o objeto SpatialPixelsDataFrame</div>
<div>#kc_output</div><div>#plot(x,y,xlab="predição da krigagem",ylab="dados")</div><div>edit(data.frame(kc_output))</div><div><br></div><div><br></div><div style>Seria bom conseguir apenas os valores de uma área com borda, também não consegui</div>
<div style><br></div><div style><br></div></div><div style>Abraço</div><div style>Hélio</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 11 de abril de 2013 10:37, ASANTOS [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4658942h70@n4.nabble.com" target="_blank">ml-node+s2285057n4658942h70@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hélio,<br>
<br>
Postei o problema no R-sig-geo internacional (<a href="https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/018003.html" style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:left;font-style:normal;font-weight:normal;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:'Times New Roman';word-spacing:0px" rel="nofollow" link="external" target="_blank"><span> </span>Extract a pixel values in
regular grid on a krigging object<span> </span></a><a name="13df950f7f42972d_18003" style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:left;font-style:normal;font-weight:normal;line-height:normal;text-transform:none;font-size:medium;white-space:normal;font-family:'Times New Roman';word-spacing:0px" rel="nofollow" link="external"> </a>), acompanhe os
desdobramentos em:<br>
<br>
<a href="https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/subject.html" rel="nofollow" link="external" target="_blank">https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/subject.html</a><br>
<br>
Abraço,<br>
<br>
<br>
Alexandre<br>
<br>
<br>
<div>Em 11/04/2013 07:12, Hélio Gallo Rocha
escreveu:<br>
</div>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" type="cite">
<div dir="ltr"><div class="im">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><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="http://user/SendEmail.jtp?type=node&node=4658942&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span>
escreveu:<br>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" class="gmail_quote"><div class="im">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" rel="nofollow" link="external" target="_blank">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br>
</div>
Phone: <a moz-do-not-send="true"
href="tel:%28612%29%20323%209898" value="<a href="tel:%2B16123239898" value="+16123239898" target="_blank">+16123239898</a>"
target="_blank"><a href="tel:%28612%29%20323%209898" value="+16123239898" target="_blank">(612) 323 9898</a><br>
<br>
<br>
________________________________<div class="im"><br>
From: Hélio Gallo Rocha <<a href="http://user/SendEmail.jtp?type=node&node=4658942&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
To: r-br <<a href="http://user/SendEmail.jtp?type=node&node=4658942&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br></div><div class="im">
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>
</div><a href="http://user/SendEmail.jtp?type=node&node=4658942&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a><div class="im"><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>)
e forneça código mínimo reproduzível.<br>
</div></blockquote>
</div><div class="im">
<br>
<br clear="all">
<div><br>
</div>
-- <br>
Hélio Gallo Rocha<br>
IFSULDEMINAS - Câmpus Muzambinho<br>
</div></div>
</div>
<br>
<fieldset></fieldset>
<br>
<pre><div class="im">_______________________________________________
R-br mailing list
<a href="http://user/SendEmail.jtp?type=node&node=4658942&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a>
</div><div class="im"><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div></pre>
</blockquote>
<br>
<pre cols="72">--
======================================================================
Alexandre dos Santos
Proteção Florestal
Coordenador do curso Técnico em Florestas
Vice Coordenador do curso de Engenharia 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: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
<a href="http://user/SendEmail.jtp?type=node&node=4658942&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<a href="http://user/SendEmail.jtp?type=node&node=4658942&i=6" rel="nofollow" link="external" target="_blank">[hidden email]</a>
======================================================================
</pre><div class="im">
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4658942&i=7" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br></div><div class="im"><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" 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" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
</div><hr noshade size="1" color="#cccccc">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div class="im">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div><a href="http://r-br.2285057.n4.nabble.com/R-br-raster-e-shape-file-tp4658937p4658942.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-raster-e-shape-file-tp4658937p4658942.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>