[R-br] Krigagem

Hélio Gallo Rocha heliogallorocha em gmail.com
Sábado Abril 27 22:22:24 BRT 2013


Concordo Elias.

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.
Teria de ver qual o algoritimo que os programas de topografia fazem esse
cálculo...não sei.

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.

Hélio



Em 27 de abril de 2013 21:35, Elias Teixeira Krainski [via R-br] <
ml-node+s2285057n4659154h42 em n4.nabble.com> escreveu:

> Caro Hélio,
>
> Não sou especialista na área, mas acho que
> não é bom separar as coordenadas de cada zona,
> projetar considerando cada zona de referência e
> depois juntá-las novamente para fazer a análise.
> Pois, dessa forma, as distâncias entre locais
> de uma zona a locais de outra estarão erradas.
>
> Att.
> Elias.
>
> On Sat, 2013-04-27 at 21:12 -0300, Hélio Gallo Rocha wrote:
>
> > Athur,
> >
> >
> > Infelizmente, desconsidere o e-mail anterior, não o CRM
> >
> >
> > 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:
> >
> >
> > coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat +zone=22
> > +south +ellps=WGS84"));coords
> > coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=22 +south
> > +ellps=WGS84 +units=m +no_defs"));coords1
> >
> >
> > ou seja, usando o CRM postado, tem de separar os pontos de fuso 22, os
> > do fuso 23 e 24 e processar separadamente.
> >
> >
> > O CRM tá certo, mas tem de processar separadamente os pontos de cada
> > fuso
> >
> >
> > Quem dá uma luz ???
> >
> >
> >
> >
> > Abraço
> >
> >
> > Hélio
> >
> >
> > Em 27 de abril de 2013 19:39, Hélio Gallo Rocha
> > <[hidden email] <http://user/SendEmail.jtp?type=node&node=4659154&i=0>>
> escreveu:
> >         Arthur, fins os testes e não precisou separar os pontos por
> >         fuso.
> >         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.
> >
> >
> >         tai o CRM
> >
> >
> >         library(geoR)
> >         library(rgdal)
> >         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
> >
> >         long= dados$coords[,1] ;long
> >         lat=  dados$coords[,2] ;lat
> >         coord=cbind(long,lat);coord
> >         coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat
> >         +zone=22 +south +ellps=WGS84"));coords
> >         coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=22
> >         +south +ellps=WGS84 +units=m +no_defs"));coords1
> >
> >
> >
> >
> >         Só por curiosidade vou fazer a krigagem usando latlong e utm,
> >         pra ver se dá o mesmo resultado, tem de dar...
> >         abraço
> >         Hélio
> >
> >
> >         Em 26 de abril de 2013 14:55, Arthur T. C. [via R-br] <ml-node
> >         +[hidden email]<http://user/SendEmail.jtp?type=node&node=4659154&i=1>>
> escreveu:
> >
> >                 A maior dúvida é o que o hélio comentou.
> >                 tenho as coordenadas em três fusos distintos.
> >                 Assim não dá pra projeta-los juntos. entende?
> >
> >
> >
> >                 ______________________________________________________
> >
> >                 From: [hidden email]
> >                 To: [hidden email]
> >                 Subject: RE: [R-br] Krigagem
> >
> >                 Date: Fri, 26 Apr 2013 14:52:22 -0300
> >
> >
> >                 Desconsidere o último email.
> >                 Esses dados fora retirados das normais climatológicas
> >                 do inmet.
> >                 Perguntei a ele o sistema de coordenadas. Me disseram
> >                 que é WGS84. É isso que perguntou Elias?
> >
> >
> >
> >                 ______________________________________________________
> >
> >                 From: [hidden email]
> >                 To: [hidden email]
> >                 Subject: RE: [R-br] Krigagem
> >                 Date: Fri, 26 Apr 2013 14:47:36 -0300
> >
> >                 Entendi.
> >                 Sei uma maneira de fazer essa transformação no GIS.
> >                 Mas usando somente um fuso. Ai que esta o problema.
> >                 Agradeço sua ajuda,
> >                 Creio que vou ter que procurar um professor de
> >                 agrimensura. hehe
> >                 Arthur
> >
> >
> >
> >                 ______________________________________________________
> >
> >                 Date: Fri, 26 Apr 2013 13:09:46 -0300
> >                 From: [hidden email]
> >                 To: [hidden email]
> >                 Subject: Re: [R-br] Krigagem
> >
> >                 Arthur
> >                 Pelo que vi seus dados estão em tres fusos, 22, 23 e
> >                 24.
> >                 Para transformar em UTM no R eu não sei de forma
> >                 automática todos os dados de uma vez só...
> >
> >
> >
> >
> >                 tirei da lista da msg de Elias T. Krainski
> >
> >
> >
> >
> >                 Coord <- SpatialPoints(cbind(dados$coord_x,dados
> >                 $coord_y), proj4string=CRS("+proj=longlat
> >                 +datum=WGS84"))
> >                 utm <- spTransform(Coord, CRS("+proj=utm +zone=22S"))
> >
> >
> >
> >                 e fiz o seguinte  teste
> >                 d=cbind(-17,-50) # estão na zona 22S
> >                  d
> >                 coord<- SpatialPoints((d),
> >                 proj4string=CRS("+proj=longlat +datum=WGS84"))
> >                 coord
> >                 utm <- spTransform(coord, CRS("+proj=utm +zone=22S"))
> >                 utm
> >
> >
> >                 Coord <- SpatialPoints((d),
> >                 proj4string=CRS("+proj=longlat +zone=22S
> >                 +datum=WGS84"))
> >                 Coord
> >                 utm1 <- spTransform(Coord, CRS("+proj=utm
> >                 +zone=22S"))
> >                 utm1
> >
> >
> >                 a resposta seria W606447.747 e N8120173.351
> >
> >
> >                 mas em nenhum dos testes deu certo
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >                 Hélio
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >
> >                 Em 26 de abril de 2013 10:14, Arthur T. C. [via R-br]
> >                 <[hidden email]> escreveu:
> >                         Verdade. Rs.
> >                         Tenho os pontos em coordenadas geográficas. em
> >                         um datum qualquer. Por exemplo SAD 69.
> >                         Necessito projetar pra saber as distâncias
> >                         linear entre os pontos (superficial). E não
> >                         cortar a esfera terrestre.
> >                         É isso certo? Os dados estão em anexo.
> >                         Agradeço Hélio,
> >                         Arthur
> >
> >
> >
> >                         ______________________________________________
> >
> >                         Date: Fri, 26 Apr 2013 10:00:37 -0300
> >                         From: [hidden email]
> >                         To: [hidden email]
> >                         Subject: Re: [R-br] Krigagem
> >
> >                         Arthur.
> >                         Parece que está fazendo uma confusão entre
> >                         fuso e datum.
> >                         anexe seus dados pra que eu possa dar uma
> >                         olhada, com o CRM que tem
> >
> >
> >
> >
> >                         Hélio
> >
> >                         Em 26 de abril de 2013 09:24, Arthur T. C.
> >                         [via R-br] <[hidden email]> escreveu:
> >
> >                                 Bom dia Senhores,
> >                                 Estou com um pequeno problemas em
> >                                 relação a Krigagem e talvez os
> >                                 senhores tenham uma resposta.
> >                                 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.
> >                                 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.
> >                                 Agradeço desde já,
> >                                 Arthur Calegario.
> >
> >
> _______________________________________________
> >                                 R-br mailing list
> >                                 [hidden email]
> >
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >                                 Leia o guia de postagem
> >                                 (http://www.leg.ufpr.br/r-br-guia) e
> >                                 forneça código mínimo reproduzível.
> >
> >
> >                                 ______________________________________
> >                                 If you reply to this email, your
> >                                 message will be added to the
> >                                 discussion below:
> >
> http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130.html
> >                                 To unsubscribe from R-br, click here.
> >                                 NAML
> >
> >
> >
> >
> >                         --
> >                         Hélio Gallo Rocha
> >                         IFSULDEMINAS - Câmpus Muzambinho
> >
> >
> >                         _______________________________________________
> R-br mailing list [hidden email]
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de
> postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo
> reproduz�vel.
> >
> >                         _______________________________________________
> >                         R-br mailing list
> >                         [hidden email]
> >
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >                         Leia o guia de postagem
> >                         (http://www.leg.ufpr.br/r-br-guia) e forneça
> >                         código mínimo reproduzível.
> >
> >                         Dados.txt (2K) Download Attachment
> >
> >
> >
> >                         ______________________________________________
> >                         If you reply to this email, your message will
> >                         be added to the discussion below:
> >
> http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659134.html
> >                         To unsubscribe from R-br, click here.
> >                         NAML
> >
> >
> >
> >
> >                 --
> >                 Hélio Gallo Rocha
> >                 IFSULDEMINAS - Câmpus Muzambinho
> >
> >
> >                 _______________________________________________ R-br
> >                 mailing list [hidden email]
> >                 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia
> o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo
> m�nimo reproduz�vel.
> >
> >                 _______________________________________________
> >                 R-br mailing list
> >                 [hidden email]
> >                 https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >                 Leia o guia de postagem
> >                 (http://www.leg.ufpr.br/r-br-guia) e forneça código
> >                 mínimo reproduzível.
> >
> >
> >                 ______________________________________________________
> >                 If you reply to this email, your message will be added
> >                 to the discussion below:
> >
> http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659145.html
> >                 To unsubscribe from R-br, click here.
> >                 NAML
> >
> >
> >
> >
> >         --
> >         Hélio Gallo Rocha
> >         IFSULDEMINAS - Câmpus Muzambinho
> >
> >
> >
> >
> >
> > --
> > Hélio Gallo Rocha
> > IFSULDEMINAS - Câmpus Muzambinho
> >
> > _______________________________________________
> > R-br mailing list
> > [hidden email] <http://user/SendEmail.jtp?type=node&node=4659154&i=2>
> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
>
> _______________________________________________
> R-br mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659154&i=3>
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
> ------------------------------
>  If you reply to this email, your message will be added to the discussion
> below:
> http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659154.html
>  To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
> .
> NAML<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>
>



-- 
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130427/b75df72d/attachment.html>


Mais detalhes sobre a lista de discussão R-br