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