[R-br] Krigagem

Domingo Abril 28 21:42:04 BRT 2013

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 <heliogallorocha em gmail.com>

> 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
>
> _______________________________________________
> 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/20130428/76bb9221/attachment-0001.html>