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