Em 27 de abril de 2013 22:42, Jônatan [via R-br] <ml-node+s2285057n4659157h35@n4.nabble.com> escreveu:
Quando trabalho com dados na escala continental eu uso a projeção Azimutal equivalente de Lambert com meridiano centrado na região de interesse, por exemplo para América do Sul em -60. Acho que pode ser uma alternativa.require(sp)require(rgdal)p <- read.table("C:/Users/Jonatan/Downloads/Dados.txt", head = T,blank.lines.skip=T)pll <- SpatialPointsDataFrame(p[,1:2], data=p[3], proj4string=CRS("+proj=longlat +datum=WGS84"))## projecao azimutal de área equivalente de Lambert (meridiano central em -60 para America do Sul)projlaea <- "+proj=laea +lat_0=-15 +lon_0=-60 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs"## dados reprojetados em laea, unidade em mplaea <- spTransform(x=pll, CRSobj=CRS(projlaea))## fazer Krigagem...## reprojetar objeto com resultado para latlon se desejavel (usando ultimo comando)resll <- spTransform(x=resultadoKrigagem, CRSobj=CRS("+proj=longlat +datum=WGS84"))2013/4/26 Arthur Telles Calegario <[hidden email]>
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?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. heheArthur
Date: Fri, 26 Apr 2013 13:09:46 -0300
From: [hidden email]
To: [hidden email]
Subject: Re: [R-br] KrigagemArthurPelo 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 tested=cbind(-17,-50) # estão na zona 22Sdcoord<- SpatialPoints((d), proj4string=CRS("+proj=longlat +datum=WGS84"))coordutm <- spTransform(coord, CRS("+proj=utm +zone=22S"))utmCoord <- SpatialPoints((d), proj4string=CRS("+proj=longlat +zone=22S +datum=WGS84"))Coordutm1 <- spTransform(Coord, CRS("+proj=utm +zone=22S"))utm1a resposta seria W606447.747 e N8120173.351mas em nenhum dos testes deu certo
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] KrigagemArthur.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 temHélioEm 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--
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.
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659134.htmlIf you reply to this email, your message will be added to the discussion below:--
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.--
###############################################################
## Jônatan Dupont Tatsch
## *Pesquisador Associado Recém Doutor## do Programa de Pós Graduação em Meteorologia
## Universidade Federal de Santa Maria (UFSM)
## Santa Maria/RS - Brasil## Telefone: (55) 3301 2084## www.ufsm.br/meteorologia/pos## * Pesquisador colaborador do:###############################################################
## Laboratório de Clima e Biosfera LCB/IAG/USP
## Núcleo de Apoio a Pesquisas em Mudanças Climáticas (INCLINE)
## www.incline.iag.usp.br
_______________________________________________
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.
http://r-br.2285057.n4.nabble.com/R-br-Krigagem-tp4659130p4659157.htmlIf you reply to this email, your message will be added to the discussion below:
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho