[R-br] Krigagem
Jônatan
jdtatsch em gmail.com
Sábado Abril 27 22:42:13 BRT 2013
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 m
plaea <- 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 <tcalegario em hotmail.com>
> 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: tcalegario em hotmail.com
> To: r-br em listas.c3sl.ufpr.br
> 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: tcalegario em hotmail.com
> To: r-br em listas.c3sl.ufpr.br
> 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: heliogallorocha em gmail.com
> To: r-br em listas.c3sl.ufpr.br
> 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] <
> ml-node+s2285057n4659134h22 em n4.nabble.com> 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<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: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
> 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.
>
> _______________________________________________
> 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.
>
--
###############################################################
## 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
###############################################################
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130427/ab9ace6e/attachment.html>
Mais detalhes sobre a lista de discussão R-br