[R-br] Krigagem
Jônatan
jdtatsch em gmail.com
Terça Abril 30 12:09:14 BRT 2013
Os coordenadas de referência estão definidas como x_0=0 +y_0=0, pq? Tem que
definir como cll (ver no ultimo email que enviei).
2013/4/30 Hélio Gallo Rocha <heliogallorocha em gmail.com>
> Caro Jônatan
> Fiz a transformação de latlong para lambert como orientou, fiz a krigagem(
> no caso o melhor modelo foi por quadrados mininos ponderados).
> Para facilitar a verificação,usei l=10 no loci.
> Gerei um txt ,(em anexo) das coordenadas em lambert e o kc$predict para
> poder voltar a latlong, como disse se quiser, mas para deixar o CRM
> completo, dai a dificuladade, olha o CRM da volta:
> out_krig=read.table("k:/2013/r/arthur/graficos/out_krig_lambert.txt", head
> = T,blank.lines.skip=T) ; summary(out_krig) ;class(out_krig) ;out_krig#OK
> out_krig_spdf <- SpatialPointsDataFrame(out_krig[,2:3],
> data=out_krig[1],proj4string=CRS("+proj=laea + x_0=0 +y_0=0 +a=6370997
> +b=6370997+units=m +no_defs"))# +datum=WGS84"))
> out_krig_spdf
>
> não deu certo...
>
> Hélio
>
>
>
> Em 28 de abril de 2013 21:42, Jônatan [via R-br] <
> ml-node+s2285057n4659167h45 em n4.nabble.com> escreveu:
>
>> Hélio segue CMR da comparação do método que sugeri com outros (incluindo
>> dicas do Elias e do Paulo).
>>
>> ## comparacão da distancia entre 2 pontos (no dominio espacial dos dados)
>> por diferentes metodos
>> require(sp);require(rgdal);require(fields) # para comparacao
>>
>> ## Funcoes da dica do Elias
>> source("http://www.biostat.umn.edu/~sudiptob/Software/distonearth.R")
>>
>> ## dados
>> 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"))
>>
>> ## centro aproximado
>> cll <- data.frame(lon = round(mean(range(p$Longitude)),3), lat =
>> round(mean(range(p$Latitude)),3))
>>
>> ## adicionando ponto a uma distancia conhecida
>> loc <- rbind(cll, cll + 1)
>> locll <- SpatialPoints(loc, proj4string=CRS("+proj=longlat +datum=WGS84"))
>>
>> ## pontos e centro
>> plot(pll,col="gray")
>> points(cll,col=2,pch=20)
>>
>> ## indicando os 2 pontos
>> text(coordinates(locll),labels=rownames(loc))
>>
>> ## projecao azimutal de área equivalente de Lambert
>> projlaea0 <- "+proj=laea +lat_0=YY +lon_0=XX +x_0=0 +y_0=0 +a=6370997
>> +b=6370997 +units=m +no_defs"
>> projlaea <- gsub("XX", cll$lon, projlaea0)
>> projlaea <- gsub("YY", cll$lat, projlaea)
>>
>> ## dados reprojetados em laea, unidade em m
>> loclaea <- spTransform(x=locll, CRSobj=CRS(projlaea))
>>
>> ## distancia em km a partir das coordenadas em laea
>> res1 <- dist(coordinates(loclaea))/1000
>>
>> ## aprox grosseira, 1 grau ~ 111 km
>> res2 <- sqrt(111^2+111^2)
>>
>> ## distancia estimada pelo fields
>> res3 <- rdist.earth(loc, miles=F)
>>
>> ## Uma das Funcoes do Sudipto Banerjee
>> res4 <- rdistearth(loc)
>>
>> ## diferencas
>> res1 - rbind(res1, res2, unique(res3[res3>0]), unique(res4[res4>0]))
>>
>> Quanto as duvidas sobre as propriedades da proj. LAEA:
>> http://en.wikipedia.org/wiki/Lambert_azimuthal_equal-area_projection
>>
>>
>>
>>
>> 2013/4/28 Hélio Gallo Rocha <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659167&i=0>
>> >
>>
>>> 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] <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659167&i=1>
>>> > 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.
>>>> 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
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659167&i=2>
>>> 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=4659167&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.
>>
>> ------------------------------
>> 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-tp4659130p4659167.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
>
> _______________________________________________
> 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/20130430/238371b1/attachment.html>
Mais detalhes sobre a lista de discussão R-br