[R-br] Krigagem
Hélio Gallo Rocha
heliogallorocha em gmail.com
Terça Abril 30 09:38:32 BRT 2013
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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130430/863a1aed/attachment.html>
-------------- Próxima Parte ----------
kc.predict Var1 Var2
1439.471584 -749877.155 -712391.6699
1465.329422 -581060.704 -712391.6699
1495.849513 -412244.253 -712391.6699
1532.559175 -243427.802 -712391.6699
1570.512733 -74611.351 -712391.6699
1565.110096 94205.1 -712391.6699
1457.859239 263021.551 -712391.6699
1345.221426 431838.002 -712391.6699
1348.304618 600654.453 -712391.6699
1362.765892 769470.904 -712391.6699
1474.279885 -749877.155 -555353.3091
1500.77349 -581060.704 -555353.3091
1527.140626 -412244.253 -555353.3091
1561.443606 -243427.802 -555353.3091
1521.555604 -74611.351 -555353.3091
1474.766274 94205.1 -555353.3091
1256.024327 263021.551 -555353.3091
1313.511945 431838.002 -555353.3091
1327.413659 600654.453 -555353.3091
1373.914461 769470.904 -555353.3091
1522.508441 -749877.155 -398314.9483
1551.530978 -581060.704 -398314.9483
1535.88269 -412244.253 -398314.9483
1575.61311 -243427.802 -398314.9483
1421.226698 -74611.351 -398314.9483
1321.329377 94205.1 -398314.9483
1335.394278 263021.551 -398314.9483
1159.617269 431838.002 -398314.9483
1302.221988 600654.453 -398314.9483
1409.646779 769470.904 -398314.9483
1572.584981 -749877.155 -241276.5875
1672.958126 -581060.704 -241276.5875
1601.200514 -412244.253 -241276.5875
1503.274484 -243427.802 -241276.5875
1414.150219 -74611.351 -241276.5875
1321.216126 94205.1 -241276.5875
1178.649369 263021.551 -241276.5875
1093.982025 431838.002 -241276.5875
1313.436384 600654.453 -241276.5875
1493.859945 769470.904 -241276.5875
1575.311193 -749877.155 -84238.22665
1597.079043 -581060.704 -84238.22665
1772.743771 -412244.253 -84238.22665
1556.435251 -243427.802 -84238.22665
1375.731907 -74611.351 -84238.22665
1181.369146 94205.1 -84238.22665
932.7270225 263021.551 -84238.22665
908.4123845 431838.002 -84238.22665
1236.646787 600654.453 -84238.22665
1659.9622 769470.904 -84238.22665
1519.472496 -749877.155 72800.13418
1473.881724 -581060.704 72800.13418
1631.115586 -412244.253 72800.13418
1532.545916 -243427.802 72800.13418
1387.758713 -74611.351 72800.13418
1062.71752 94205.1 72800.13418
823.3179681 263021.551 72800.13418
791.6182586 431838.002 72800.13418
1508.027024 600654.453 72800.13418
1978.87457 769470.904 72800.13418
1445.485323 -749877.155 229838.495
1448.321808 -581060.704 229838.495
1478.327233 -412244.253 229838.495
1476.383982 -243427.802 229838.495
1387.620359 -74611.351 229838.495
1024.124504 94205.1 229838.495
856.8627057 263021.551 229838.495
721.0660946 431838.002 229838.495
1101.135134 600654.453 229838.495
2103.231346 769470.904 229838.495
1397.682517 -749877.155 386876.8558
1393.766889 -581060.704 386876.8558
1376.103706 -412244.253 386876.8558
1317.75947 -243427.802 386876.8558
1233.785297 -74611.351 386876.8558
1001.003757 94205.1 386876.8558
862.5257413 263021.551 386876.8558
1206.938437 431838.002 386876.8558
870.5745439 600654.453 386876.8558
1501.733103 769470.904 386876.8558
1362.758142 -749877.155 543915.2166
1349.384967 -581060.704 543915.2166
1317.517255 -412244.253 543915.2166
1258.370748 -243427.802 543915.2166
1159.094384 -74611.351 543915.2166
976.6134194 94205.1 543915.2166
670.2947772 263021.551 543915.2166
711.6220539 431838.002 543915.2166
766.1181506 600654.453 543915.2166
971.9871309 769470.904 543915.2166
1339.219433 -749877.155 700953.5775
1321.56908 -581060.704 700953.5775
1287.411876 -412244.253 700953.5775
1228.404423 -243427.802 700953.5775
1129.792915 -74611.351 700953.5775
965.1359444 94205.1 700953.5775
766.5599213 263021.551 700953.5775
734.5649371 431838.002 700953.5775
694.6837369 600654.453 700953.5775
582.9 769470.904 700953.5775
Mais detalhes sobre a lista de discussão R-br