[R-br] Krigagem
Elias Teixeira Krainski
eliaskrainski em yahoo.com.br
Sábado Abril 27 21:35:03 BRT 2013
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
> <heliogallorocha em gmail.com> 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
> +s2285057n4659145h63 em n4.nabble.com> 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
> R-br em listas.c3sl.ufpr.br
> 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.
Mais detalhes sobre a lista de discussão R-br