[R-br] Krigagem
Hélio Gallo Rocha
heliogallorocha em gmail.com
Domingo Abril 28 14:55:29 BRT 2013
Caro Jônatan,
Nunca usei lambert
duas dúvidas, talvez a resposta da primeira já resolva:
no CRM: +lat_0=-15 +lon_0=-60, centro da projeção
quais os limites máximos para que não haja distorções? de forma rápida
achei que a escala máxima seria de 1:500.000.
Num papel de 50cmX50cm teríamos uma escala de 1:3.000.000, para cobri a
área do Arthur
No Caso do problema do Arthur, que temos:
long lat minima -38.13 -9.24 máxima -52.14 -22.06 amplitude 14.01
12.82
multiplicando por 111Km cada grau, teríamos um quadrado de 1500kmx1500km,
aproximadamente
usando lat_0=-15 para as coordenadas
( long -52.14) e (lat -15.54 )
a latitude -15.54 maior que -15
resultou em lambert
(841253, -75134.7) 1575.2
seria isso mesmo?
ou colocaria lat_0=-25, já que a maior latitude do exemplo do Arthur é
-22.06?
Hélio
Em 27 de abril de 2013 22:42, Jônatan [via R-br] <
ml-node+s2285057n4659157h35 em 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 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 <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659157&i=0>
> >
>
> 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=4659157&i=1>
>> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659157&i=2>
>> 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=4659157&i=3>
>> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659157&i=4>
>> 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=4659157&i=5>
>> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659157&i=6>
>> 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=4659157&i=7>
>> > 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=4659157&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.
>>
>> _______________________________________________
>> R-br mailing list
>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659157&i=9>
>> 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] <http://user/SendEmail.jtp?type=node&node=4659157&i=10>
> 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-tp4659130p4659157.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/20130428/2eb085b0/attachment.html>
Mais detalhes sobre a lista de discussão R-br