[R-br] Krigagem

Sábado Abril 27 19:39:33 BRT 2013

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]<http://user/SendEmail.jtp?type=node&node=4659145&i=0>
> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659145&i=1>
> 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]<http://user/SendEmail.jtp?type=node&node=4659145&i=2>
> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659145&i=3>
> 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]<http://user/SendEmail.jtp?type=node&node=4659145&i=4>
> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659145&i=5>
> 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]<http://user/SendEmail.jtp?type=node&node=4659145&i=6>
> > 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]<http://user/SendEmail.jtp?type=node&node=4659134&i=0>
> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659134&i=1>
> 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]<http://user/SendEmail.jtp?type=node&node=4659134&i=2>
> > 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] <http://user/SendEmail.jtp?type=node&node=4659130&i=0>
> 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<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble:email.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble:email.naml-instant_emails%21nabble:email.naml-send_instant_email%21nabble:email.naml>
>
>
>
>
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
> _______________________________________________ R-br mailing list [hidden
> email] <http://user/SendEmail.jtp?type=node&node=4659134&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.
>
> _______________________________________________
> R-br mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659134&i=4>
> 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<http://r-br.2285057.n4.nabble.com/attachment/4659134/0/Dados.txt>
>
>
> ------------------------------
>  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<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble:email.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble:email.naml-instant_emails%21nabble:email.naml-send_instant_email%21nabble:email.naml>
>
>
>
>
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
> _______________________________________________ R-br mailing list [hidden
> email] <http://user/SendEmail.jtp?type=node&node=4659145&i=7>
> 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=4659145&i=8>
> 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<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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130427/a8129393/attachment.html>