[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