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