<div dir="ltr">Athur,<div><br></div><div>Infelizmente, desconsidere o e-mail anterior, não o CRM</div><div><br></div><div style>Quando testei separadamenteos pontos dos fusos 22, 23 e 24, usei para todos o valor 22, esqueci de mudar os valores dos fusos, nestas linhas:</div>
<div style><br></div><div style><div>coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat +zone=<font color="#ff00ff" size="4"><b>22</b> </font>+south +ellps=WGS84"));coords</div><div>coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=<font color="#ff00ff" size="4"><b>22</b> </font>+south +ellps=WGS84 +units=m +no_defs"));coords1</div>
<div><br></div><div style>ou seja, usando o CRM postado, tem de separar os pontos de fuso 22, os do fuso 23 e 24 e processar separadamente.</div><div style><br></div><div style>O CRM tá certo, mas tem de processar separadamente os pontos de cada fuso</div>
<div style><br></div><div style>Quem dá uma luz ???</div><div style><br></div><div style><br></div><div style>Abraço</div><div style><br></div><div style>Hélio</div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
Em 27 de abril de 2013 19:39, Hélio Gallo Rocha <span dir="ltr"><<a href="mailto:heliogallorocha@gmail.com" target="_blank">heliogallorocha@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Arthur, fins os testes e não precisou separar os pontos por fuso.<div>testei da seguinte forma, fiz a conversão usando o arquivo txt com todos os pontos, depois separeis os pontos por fuso, 22,23 e 24 e comparei os resultados, não houve diferença.</div>

<div><br></div><div>tai o CRM</div><div><br></div><div><div>library(geoR)</div><div>library(rgdal)</div><div>dados=read.geodata("precipita.txt ",coords.col=1:2,data.col=3,sep = "",header=T);dados;n=summary(dados)[[1]];n##o arquivo precipita.txt é o seu mesmo, só renomeei<br>

</div><div>long= dados$coords[,1] ;long</div><div>lat=  dados$coords[,2] ;lat</div><div>coord=cbind(long,lat);coord</div><div>coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat +zone=22 +south +ellps=WGS84"));coords</div>

<div>coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=22 +south +ellps=WGS84 +units=m +no_defs"));coords1   </div><div><br></div><div><br></div><div>Só por curiosidade vou fazer a krigagem usando latlong e utm, pra ver se dá o mesmo resultado, tem de dar...</div>

<div>abraço</div><div>Hélio</div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 26 de abril de 2013 14:55, Arthur T. C. [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4659145h63@n4.nabble.com" target="_blank">ml-node+s2285057n4659145h63@n4.nabble.com</a>></span> escreveu:<div>
<div class="h5"><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

        


<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><div></div><hr>From: <a href="http://user/SendEmail.jtp?type=node&node=4659145&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>

To: <a href="http://user/SendEmail.jtp?type=node&node=4659145&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>Subject: RE: [R-br] Krigagem<br></div>Date: Fri, 26 Apr 2013 14:52:22 -0300<div>

<div><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="http://user/SendEmail.jtp?type=node&node=4659145&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>To: <a href="http://user/SendEmail.jtp?type=node&node=4659145&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</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="http://user/SendEmail.jtp?type=node&node=4659145&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>

To: <a href="http://user/SendEmail.jtp?type=node&node=4659145&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</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="http://user/SendEmail.jtp?type=node&node=4659145&i=6" rel="nofollow" link="external" target="_blank">[hidden email]</a>></span> escreveu:<br>


<blockquote style="border-left:2px solid #cccccc;padding:0 1em">

        


<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" link="external" target="_blank">[hidden email]</a><br>


To: <a href="http://user/SendEmail.jtp?type=node&node=4659134&i=1" rel="nofollow" link="external" 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" link="external" 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" link="external" target="_blank">[hidden email]</a>
<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>
        <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" link="external" 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" link="external">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" link="external" 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" link="external" target="_blank">[hidden email]</a>
<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></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" link="external" target="_blank">[hidden email]</a>
<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.</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" link="external" 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" rel="nofollow" link="external" 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 rel="nofollow" link="external">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" link="external" 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="http://user/SendEmail.jtp?type=node&node=4659145&i=7" rel="nofollow" link="external" 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></div>

</div></div></div></div>                                          </div><div><div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4659145&i=8" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<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>
        <br>
        <hr noshade size="1" color="#cccccc">
        </div></div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div><div>
                <div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
                </div></div><a href="http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659145.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659145.html</a>
        </div><div><div>
        <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></div></div><div><div class="h5"><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>