<div dir="ltr">Quando trabalho com dados na escala continental eu uso a projeção Azimutal equivalente de Lambert com meridiano centrado na região de interesse, por exemplo para América do Sul em -60. Acho que pode ser uma alternativa.<div>
<br><div><div>require(sp)</div><div>require(rgdal)</div><div>p <- read.table("C:/Users/Jonatan/Downloads/Dados.txt", head = T,blank.lines.skip=T)</div><div><br></div><div>pll <- SpatialPointsDataFrame(p[,1:2], data=p[3], proj4string=CRS("+proj=longlat +datum=WGS84"))</div>
<div><br></div><div>## projecao azimutal de área equivalente de Lambert (meridiano central em -60 para America do Sul)</div><div>projlaea <- "+proj=laea +lat_0=-15 +lon_0=-60 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs"</div>
<div><br></div><div>## dados reprojetados em laea, unidade em m </div><div>plaea <- spTransform(x=pll, CRSobj=CRS(projlaea))</div><div><br></div><div>## fazer Krigagem</div><div>...</div><div><br></div><div>## reprojetar objeto com resultado para latlon se desejavel (usando ultimo comando)</div>
</div></div><div>resll <- spTransform(x=resultadoKrigagem, CRSobj=CRS("+proj=longlat +datum=WGS84"))<br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/4/26 Arthur Telles Calegario <span dir="ltr"><<a href="mailto:tcalegario@hotmail.com" target="_blank">tcalegario@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">A maior dúvida é o que o hélio comentou.<div>tenho as coordenadas em três fusos distintos. </div><div>Assim não dá pra projeta-los juntos. entende?<br><div><br><div><div class="im"><div></div><hr>From: <a href="mailto:tcalegario@hotmail.com" target="_blank">tcalegario@hotmail.com</a><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>Subject: RE: [R-br] Krigagem<br></div>Date: Fri, 26 Apr 2013 14:52:22 -0300<div><div class="h5"><br><br>
<div dir="ltr">Desconsidere o último email.<div>Esses dados fora retirados das normais climatológicas do inmet.</div><div>Perguntei a ele o sistema de coordenadas. Me disseram que é WGS84. É isso que perguntou Elias?<br>
<br><div><div></div><hr>From: <a href="mailto:tcalegario@hotmail.com" target="_blank">tcalegario@hotmail.com</a><br>To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>Subject: RE: [R-br] Krigagem<br>
Date: Fri, 26 Apr 2013 14:47:36 -0300<br><br>
<div dir="ltr">Entendi.<div>Sei uma maneira de fazer essa transformação no GIS. Mas usando somente um fuso. Ai que esta o problema.</div><div>Agradeço sua ajuda,</div><div>Creio que vou ter que procurar um professor de agrimensura. hehe</div>
<div>Arthur<br><br><div><div></div><hr>Date: Fri, 26 Apr 2013 13:09:46 -0300<br>From: <a href="mailto:heliogallorocha@gmail.com" target="_blank">heliogallorocha@gmail.com</a><br>To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Krigagem<br><br><div dir="ltr">Arthur<div>Pelo que vi seus dados estão em tres fusos, 22, 23 e 24.</div><div>Para transformar em UTM no R eu não sei de forma automática todos os dados de uma vez só...</div>
<div><br></div>
<div><br></div><div>tirei da lista da msg de <span style="font-size:16px;font-family:'times new roman','new york',times,serif">Elias T. Krainski</span></div><div><span style="font-size:16px;font-family:'times new roman','new york',times,serif"><br>
</span></div><div><span style="font-size:16px;font-family:'times new roman','new york',times,serif"><br></span></div><div><span style="font-size:16px;font-family:'times new roman','new york',times,serif">Coord <- SpatialPoints(cbind(dados$coord_x,dados$coord_y), proj4string=CRS("+proj=longlat +datum=WGS84"))<br>
<span style="font-weight:bold;background-color:rgb(255,255,153)">utm</span> <- spTransform(Coord, CRS("+proj=<span style="font-weight:bold;background-color:rgb(255,255,153)">utm</span> </span><span style="font-size:12pt;font-family:'times new roman','new york',times,serif">+zone=22S</span><span style="font-size:12pt;font-family:'times new roman','new york',times,serif">"))</span><br>
</div><div><span style="font-size:16px;font-family:'times new roman','new york',times,serif"><br></span></div><div><span style="font-size:16px;font-family:'times new roman','new york',times,serif">e fiz o seguinte teste
</span></div><div><div><font face="times new roman, new york, times, serif" size="3">d=cbind(-17,-50) # estão na zona 22S</font></div><div><font face="times new roman, new york, times, serif" size="3"> d</font></div>
<div><font face="times new roman, new york, times, serif" size="3">coord<- SpatialPoints((d), proj4string=CRS("+proj=longlat +datum=WGS84"))</font></div><div><font face="times new roman, new york, times, serif" size="3">coord</font></div>
<div><font face="times new roman, new york, times, serif" size="3">utm <- spTransform(coord, CRS("+proj=utm +zone=22S"))</font></div><div><font face="times new roman, new york, times, serif" size="3">utm</font></div>
<div><font face="times new roman, new york, times, serif" size="3"><br></font></div><div><font face="times new roman, new york, times, serif" size="3">Coord <- SpatialPoints((d), proj4string=CRS("+proj=longlat +zone=22S +datum=WGS84"))</font></div>
<div><font face="times new roman, new york, times, serif" size="3">Coord</font></div><div><font face="times new roman, new york, times, serif" size="3">utm1 <- spTransform(Coord, CRS("+proj=utm +zone=22S"))</font></div>
<div><font face="times new roman, new york, times, serif" size="3">utm1</font></div><div><font face="times new roman, new york, times, serif" size="3"><br></font></div><div><font face="times new roman, new york, times, serif" size="3">a resposta seria W606447.747 e N8120173.351</font></div>
<div><font face="times new roman, new york, times, serif" size="3"><br></font></div><div><font face="times new roman, new york, times, serif" size="3">mas em nenhum dos testes deu certo</font></div><table border="0" cellpadding="0" cellspacing="0" width="161" style="border-collapse:collapse;width:121pt">
<tbody><tr height="17" style="min-height:12.75pt"><td height="17" align="right" width="84" style="min-height:12.75pt;width:63pt"><br></td>
<td align="right" width="77" style="width:58pt"><br></td></tr></tbody></table> </div><div><br></div><div><br></div><div>Hélio</div><div><br></div><div><br>
</div><div><br></div><div><br></div></div><div><br><br><div>Em 26 de abril de 2013 10:14, Arthur T. C. [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4659134h22@n4.nabble.com" target="_blank">ml-node+s2285057n4659134h22@n4.nabble.com</a>></span> escreveu:<br>
<blockquote style="border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Verdade. Rs. <div>Tenho os pontos em coordenadas geográficas. em um datum qualquer. Por exemplo SAD 69.</div><div>Necessito projetar pra saber as distâncias linear entre os pontos (superficial). E não cortar a esfera terrestre.</div>
<div>É isso certo? Os dados estão em anexo.</div><div>Agradeço Hélio,</div><div>Arthur<br><div><div><br><div><div></div><hr>Date: Fri, 26 Apr 2013 10:00:37 -0300<br>From: <a href="http://user/SendEmail.jtp?type=node&node=4659134&i=0" rel="nofollow" target="_blank">[hidden email]</a><br>
To: <a href="http://user/SendEmail.jtp?type=node&node=4659134&i=1" rel="nofollow" target="_blank">[hidden email]</a><br>Subject: Re: [R-br] Krigagem<br><br><div dir="ltr"><div>Arthur.<div>
Parece que está fazendo uma confusão entre fuso e datum.</div><div>anexe seus dados pra que eu possa dar uma olhada, com o CRM que tem</div><div><br></div><div><br></div></div><div>
Hélio<br><br><div>Em 26 de abril de 2013 09:24, Arthur T. C. [via R-br] <span dir="ltr"><<a href="http://user/SendEmail.jtp?type=node&node=4659134&i=2" rel="nofollow" target="_blank">[hidden email]</a>></span> escreveu:<div>
<br>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em">
<div dir="ltr">Bom dia Senhores,<div>Estou com um pequeno problemas em relação a Krigagem e talvez os senhores tenham uma resposta.</div><div>Na verdade o problema esta em preparar os dados. Estou trabalhando com evapotranspiração no Rio Doce e possuo as coordenadas das estações meteorológicas em graus. Contudo, para a construção do semivariograma, necessito que a distância entre as mesmas seja linear.</div>
<div>Assim devo projeta-las. É um procedimento simples quando todos as estações estão em um mesmo 'datum'. Contudo, as estações que estão do 'datum' 23S são insuficientes para a construção do semivariograma. Assim necessito saber como faço para linearizar as distâncias entre estações meteorológicas cujos limites extrapolam um 'datum' apenas. Entenderam? rs. </div>
<div>Agradeço desde já,</div><div>Arthur Calegario.</div> </div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4659130&i=0" rel="nofollow" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
<a href="http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130.html" rel="nofollow" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130.html</a>
</div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;line-height:1.5em">
To unsubscribe from R-br, <a rel="nofollow">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble:email.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble:email.naml-instant_emails%21nabble:email.naml-send_instant_email%21nabble:email.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></blockquote></div></div><div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div></div></div>
<br>_______________________________________________
R-br mailing list
<a href="http://user/SendEmail.jtp?type=node&node=4659134&i=3" rel="nofollow" target="_blank">[hidden email]</a>
<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>
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.</div></div></div></div> </div><div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4659134&i=4" rel="nofollow" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div>
<div><br><img> <strong>Dados.txt</strong> (2K) <a href="http://r-br.2285057.n4.nabble.com/attachment/4659134/0/Dados.txt" rel="nofollow" target="_blank">Download Attachment</a></div>
<br>
<br>
<hr noshade size="1" color="#cccccc">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div>
<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-Krigagem-tp4659130p4659134.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659134.html</a>
</div><div><div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;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:email.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble:email.naml-instant_emails%21nabble:email.naml-send_instant_email%21nabble:email.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>
<br>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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.</div></div> </div></div></div> </div></div>
</div></div></div></div> </div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">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>
<div>###############################################################<br>## Jônatan Dupont Tatsch<br>## *Pesquisador Associado Recém Doutor </div><div>## do Programa de Pós Graduação em Meteorologia<br>## Universidade Federal de Santa Maria (UFSM)<br>
## Santa Maria/RS - Brasil</div><div>## Telefone: (55) 3301 2084</div>## <a href="http://www.ufsm.br/meteorologia/pos" target="_blank">www.ufsm.br/meteorologia/pos</a><div>## * Pesquisador colaborador do:<br>## Laboratório de Clima e Biosfera LCB/IAG/USP <br>
## Núcleo de Apoio a Pesquisas em Mudanças Climáticas (INCLINE)<br>## <a href="http://www.incline.iag.usp.br" target="_blank">www.incline.iag.usp.br</a></div>###############################################################
</div>