[R-br] Krigagem
Paulo Justiniano
paulojus em leg.ufpr.br
Domingo Abril 28 21:12:12 BRT 2013
Uma possivel alternativa é manter as coordenadas lat-log
e usar funcoes que calculam distancias no geoide
como a fields:::rdist.earth
On Sat, 27 Apr 2013, Hélio Gallo Rocha wrote:
> 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]> 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]> 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]
> > 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-tp4659130p4659154.html
> To unsubscribe from R-br, click here.
> NAML
>
>
>
>
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
>
Mais detalhes sobre a lista de discussão R-br