<div dir="ltr"><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">Concordo Elias.<div><br></div><div>Esse é um problema que não conheço a resolução, de ter pontos em diversos fusos, por isso perguntei se alguém tinha uma luz.</div>

<div>Teria de ver qual o algoritimo que os programas de topografia fazem esse cálculo...não sei.<br><div><br></div><div>O Arthur vai ter de fazer a krigagem usando latlong mesmo, mas isso não é problema, já que ainda são dados em um plano cartesiano.</div>

</div><div><br></div><div>Hélio</div><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 27 de abril de 2013 21:35, Elias Teixeira Krainski [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4659154h42@n4.nabble.com" target="_blank">ml-node+s2285057n4659154h42@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">

        Caro Hélio, 

<br><br>Não sou especialista na área, mas acho que 

<br>não é bom separar as coordenadas de cada zona, 

<br>projetar considerando cada zona de referência e 

depois juntá-las novamente para fazer a análise.

<br>Pois, dessa forma, as distâncias entre locais 

de uma zona a locais de outra estarão erradas.

Att.

Elias.

<br><div class="im"><br>On Sat, 2013-04-27 at 21:12 -0300, Hélio Gallo Rocha wrote: 

</div><div><div class="im"><br>> Athur,

<br>> 

<br>> 

<br>> Infelizmente, desconsidere o e-mail anterior, não o CRM

<br>> 

<br>> 

<br>> Quando testei separadamenteos pontos dos fusos 22, 23 e 24, usei para

<br>> todos o valor 22, esqueci de mudar os valores dos fusos, nestas

<br>> linhas:

<br>> 

<br>> 

<br>> coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat +zone=22

<br>> +south +ellps=WGS84"));coords

<br>> coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=22 +south

<br>> +ellps=WGS84 +units=m +no_defs"));coords1

<br>> 

<br>> 

<br>> ou seja, usando o CRM postado, tem de separar os pontos de fuso 22, os

<br>> do fuso 23 e 24 e processar separadamente.

<br>> 

<br>> 

<br>> O CRM tá certo, mas tem de processar separadamente os pontos de cada

<br>> fuso

<br>> 

<br>> 

<br>> Quem dá uma luz ???

<br>> 

<br>> 

<br>> 

<br>> 

<br>> Abraço

<br>> 

<br>> 

<br>> Hélio

<br>> 

<br>> 

<br>> Em 27 de abril de 2013 19:39, Hélio Gallo Rocha

<br></div>> <<a href="http://user/SendEmail.jtp?type=node&node=4659154&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>> escreveu:

<br><div class="im">>         Arthur, fins os testes e não precisou separar os pontos por

<br>>         fuso. 

<br>>         testei da seguinte forma, fiz a conversão usando o arquivo txt

<br>>         com todos os pontos, depois separeis os pontos por fuso, 22,23

<br>>         e 24 e comparei os resultados, não houve diferença. 

<br>>         

<br>>         

<br>>         tai o CRM 

<br>>         

<br>>         

<br>>         library(geoR) 

<br>>         library(rgdal) 

<br>>         dados=read.geodata("precipita.txt

<br>>         ",coords.col=1:2,data.col=3,sep =

<br>>         "",header=T);dados;n=summary(dados)[[1]];n##o arquivo

<br>>         precipita.txt é o seu mesmo, só renomeei

<br>>         

<br>>         long= dados$coords[,1] ;long 

<br>>         lat=  dados$coords[,2] ;lat 

<br>>         coord=cbind(long,lat);coord 

<br>>         coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat

<br>>         +zone=22 +south +ellps=WGS84"));coords 

<br>>         coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=22

<br>>         +south +ellps=WGS84 +units=m +no_defs"));coords1    

<br>>         

<br>>         

<br>>         

<br>>         

<br>>         Só por curiosidade vou fazer a krigagem usando latlong e utm,

<br>>         pra ver se dá o mesmo resultado, tem de dar... 

<br>>         abraço 

<br>>         Hélio 

<br>>         

<br>>         

<br>>         Em 26 de abril de 2013 14:55, Arthur T. C. [via R-br] <ml-node

<br></div>>         +<a href="http://user/SendEmail.jtp?type=node&node=4659154&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>> escreveu: 

<br><div class="im">>         

<br>>                 A maior dúvida é o que o hélio comentou. 

<br>>                 tenho as coordenadas em três fusos distintos.  

<br>>                 Assim não dá pra projeta-los juntos. entende? 

<br>>                 

<br>>                 

<br>>                 

<br></div>>                 ______________________________________________________

<br><div class="im">>                 

<br>>                 From: [hidden email]

<br>>                 To: [hidden email]

<br>>                 Subject: RE: [R-br] Krigagem

<br>>                 

<br>>                 Date: Fri, 26 Apr 2013 14:52:22 -0300 

<br>>                 

<br>>                 

<br>>                 Desconsidere o último email. 

<br>>                 Esses dados fora retirados das normais climatológicas

<br>>                 do inmet. 

<br>>                 Perguntei a ele o sistema de coordenadas. Me disseram

>                 que é WGS84. É isso que perguntou Elias?

<br>>                 

<br>>                 

<br>>                 

<br></div>>                 ______________________________________________________

<br><div class="im">>                 

<br>>                 From: [hidden email]

<br>>                 To: [hidden email]

<br>>                 Subject: RE: [R-br] Krigagem

<br>>                 Date: Fri, 26 Apr 2013 14:47:36 -0300

<br>>                 

<br>>                 Entendi. 

>                 Sei uma maneira de fazer essa transformação no GIS.

<br>>                 Mas usando somente um fuso. Ai que esta o problema. 

<br>>                 Agradeço sua ajuda, 

<br>>                 Creio que vou ter que procurar um professor de

<br>>                 agrimensura. hehe 

<br>>                 Arthur

<br>>                 

<br>>                 

<br>>                 

<br></div>>                 ______________________________________________________

<br><div><div class="h5">>                 

<br>>                 Date: Fri, 26 Apr 2013 13:09:46 -0300

<br>>                 From: [hidden email]

<br>>                 To: [hidden email]

<br>>                 Subject: Re: [R-br] Krigagem

<br>>                 

<br>>                 Arthur 

<br>>                 Pelo que vi seus dados estão em tres fusos, 22, 23 e

<br>>                 24. 

<br>>                 Para transformar em UTM no R eu não sei de forma

<br>>                 automática todos os dados de uma vez só... 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 tirei da lista da msg de Elias T. Krainski 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 Coord <- SpatialPoints(cbind(dados$coord_x,dados

<br>>                 $coord_y), proj4string=CRS("+proj=longlat

<br>>                 +datum=WGS84"))

<br>>                 utm <- spTransform(Coord, CRS("+proj=utm +zone=22S"))

<br>>                 

<br>>                 

<br>>                 

<br>>                 e fiz o seguinte  teste 

<br>>                 d=cbind(-17,-50) # estão na zona 22S 

<br>>                  d 

<br>>                 coord<- SpatialPoints((d),

<br>>                 proj4string=CRS("+proj=longlat +datum=WGS84")) 

<br>>                 coord 

<br>>                 utm <- spTransform(coord, CRS("+proj=utm +zone=22S")) 

<br>>                 utm 

<br>>                 

<br>>                 

<br>>                 Coord <- SpatialPoints((d),

<br>>                 proj4string=CRS("+proj=longlat +zone=22S

<br>>                 +datum=WGS84")) 

<br>>                 Coord 

<br>>                 utm1 <- spTransform(Coord, CRS("+proj=utm

<br>>                 +zone=22S")) 

<br>>                 utm1 

<br>>                 

<br>>                 

<br>>                 a resposta seria W606447.747 e N8120173.351 

<br>>                 

<br>>                 

<br>>                 mas em nenhum dos testes deu certo 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                   

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 Hélio 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 Em 26 de abril de 2013 10:14, Arthur T. C. [via R-br]

<br>>                 <[hidden email]> escreveu: 

<br>>                         Verdade. Rs.  

<br>>                         Tenho os pontos em coordenadas geográficas. em

<br>>                         um datum qualquer. Por exemplo SAD 69. 

<br>>                         Necessito projetar pra saber as distâncias

<br>>                         linear entre os pontos (superficial). E não

<br>>                         cortar a esfera terrestre. 

<br>>                         É isso certo? Os dados estão em anexo. 

<br>>                         Agradeço Hélio, 

<br>>                         Arthur 

<br>>                         

<br>>                         

<br>>                         

<br></div></div>>                         ______________________________________________

<br><div><div class="h5">>                         

<br>>                         Date: Fri, 26 Apr 2013 10:00:37 -0300

<br>>                         From: [hidden email]

<br>>                         To: [hidden email]

<br>>                         Subject: Re: [R-br] Krigagem

<br>>                         

<br>>                         Arthur. 

<br>>                         Parece que está fazendo uma confusão entre

<br>>                         fuso e datum. 

<br>>                         anexe seus dados pra que eu possa dar uma

<br>>                         olhada, com o CRM que tem 

<br>>                         

<br>>                         

<br>>                         

<br>>                         

<br>>                         Hélio

<br>>                         

<br>>                         Em 26 de abril de 2013 09:24, Arthur T. C.

<br>>                         [via R-br] <[hidden email]> escreveu: 

<br>>                         

<br>>                                 Bom dia Senhores, 

<br>>                                 Estou com um pequeno problemas em

<br>>                                 relação a Krigagem e talvez os

<br>>                                 senhores tenham uma resposta. 

<br>>                                 Na verdade o problema esta em preparar

<br>>                                 os dados. Estou trabalhando com

<br>>                                 evapotranspiração no Rio Doce e possuo

<br>>                                 as coordenadas das estações

<br>>                                 meteorológicas em graus. Contudo, para

<br>>                                 a construção do semivariograma,

<br>>                                 necessito que a distância entre as

<br>>                                 mesmas seja linear. 

<br>>                                 Assim devo projeta-las. É um

<br>>                                 procedimento simples quando todos as

>                                 estações estão em um mesmo 'datum'.

<br>>                                 Contudo, as estações que estão do

<br>>                                 'datum' 23S são insuficientes para a

<br>>                                 construção do semivariograma. Assim

<br>>                                 necessito saber como faço para

<br>>                                 linearizar as distâncias entre

<br>>                                 estações meteorológicas cujos limites

>                                 extrapolam um 'datum' apenas.

<br>>                                 Entenderam? rs.  

<br>>                                 Agradeço desde já, 

<br>>                                 Arthur Calegario. 

<br>>                                 

<br>>                                 _______________________________________________ 

<br>>                                 R-br mailing list 

<br>>                                 [hidden email] 

<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

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

<br>>                                 forneça código mínimo reproduzível. 

<br>>                                 

<br>>                                 

<br></div></div>>                                 ______________________________________

<br><div class="im">>                                 If you reply to this email, your

<br>>                                 message will be added to the

<br>>                                 discussion below: 

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

>                                 To unsubscribe from R-br, click here.

<br>>                                 NAML 

<br>>                         

<br>>                         

<br>>                         

<br>>                         

<br>>                         -- 

<br>>                         Hélio Gallo Rocha

<br>>                         IFSULDEMINAS - Câmpus Muzambinho

<br>>                         

<br>>                         

<br>>                         _______________________________________________ R-br mailing list [hidden email] <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. 

<br>>                         

<br>>                         _______________________________________________ 

<br>>                         R-br mailing list 

<br>>                         [hidden email] 

<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

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

<br>>                         código mínimo reproduzível. 

<br>>                         

<br></div>>                         Dados.txt (2K) Download Attachment 

<br>>                         

<br>>                         

<br>>                         

<br>>                         ______________________________________________

<br><div class="im">>                         If you reply to this email, your message will

<br>>                         be added to the discussion below: 

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

>                         To unsubscribe from R-br, click here.

<br>>                         NAML 

<br>>                 

<br>>                 

<br>>                 

<br>>                 

<br>>                 -- 

<br>>                 Hélio Gallo Rocha

<br>>                 IFSULDEMINAS - Câmpus Muzambinho

<br>>                 

<br>>                 

<br>>                 _______________________________________________ R-br

<br>>                 mailing list [hidden email]

<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> 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>>                 _______________________________________________ 

<br>>                 R-br mailing list 

<br>>                 [hidden email] 

<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

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

<br>>                 mínimo reproduzível. 

<br>>                 

<br>>                 

<br></div>>                 ______________________________________________________

<br><div class="im">>                 If you reply to this email, your message will be added

<br>>                 to the discussion below: 

<br>>                 <a href="http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659145.html" rel="nofollow" link="external" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659145.html</a> 

>                 To unsubscribe from R-br, click here.

<br>>                 NAML 

<br>>         

<br>>         

<br>>         

<br>>         

<br>>         -- 

<br>>         Hélio Gallo Rocha

<br>>         IFSULDEMINAS - Câmpus Muzambinho

<br>>         

<br>> 

<br>> 

<br>> 

<br>> 

<br>> -- 

<br>> Hélio Gallo Rocha

<br>> IFSULDEMINAS - Câmpus Muzambinho

<br>> 

<br></div><div class="im">> _______________________________________________

<br>> R-br mailing list

<br>> <a href="http://user/SendEmail.jtp?type=node&node=4659154&i=2" 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 class="im"><br><br>_______________________________________________

<br>R-br mailing list

<br><a href="http://user/SendEmail.jtp?type=node&node=4659154&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.

        <br>

        <br>

        <hr noshade size="1" color="#cccccc">

        </div><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-Krigagem-tp4659130p4659154.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659154.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>