[R-br] Krigagem

Hélio Gallo Rocha heliogallorocha em gmail.com
Terça Abril 30 17:24:02 BRT 2013


Feito, grato....

Pra mim tá resolvido

Vou postar mais uma dúvida neste capitulo mesmo, já que deva ter vários
interessados

Teria como fazer um clip para obter apenas os dados de uma área limite, mas
que não seja uma área regular.
Anexei a bordadura limite do exemplo do Arthur em latlong.

vai o CRM de conversão da borda de latlong pata lambert:
pb <- read.table("k:/2013/r/arthur/borda.txt", head = T,blank.lines.skip=T)
; pb
pllb <- SpatialPointsDataFrame(pb[,1:2],
data=pb[3],proj4string=CRS("+proj=longlat +datum=WGS84"));pllb
cllb <- data.frame(lon = round(mean(range(pb$long)),3), lat =
round(mean(range(pb$lat)),3));cllb   #média
## projecao azimutal de área equivalente de Lambert (meridiano central em
 -60 para America do Sul)
projlaeab <- "+proj=laea +lat_0=-15.65 +lon_0=-45.135 +x_0=0 +y_0=0
+a=6370997 +b=6370997 +units=m +no_defs"; projlaeab
## dados reprojetados em laea, unidade em m
plaeab <- spTransform(x=pllb, CRSobj=CRS(projlaeab)) ;      plaeab
tabela=edit(data.frame(plaeab))
write.table(tabela,file="k:\\2013/r/arthur/borda_lamb.txt",sep =
'\t',col.names = NA)

##para ler o arquivo em lambert na krigagem:
borda=read.table("borda_lamb.txt",header=T);borda
borda=cbind(borda$long,borda$lat);borda


Também fiz a krigagem usando latlong , ficou bastante parecido, com
pequenas diferenças, mas o melhor método de ajuste para latlong foi da
Verossimilhança e para lambert foi quadrados minimos ponderado

Hélio


Em 30 de abril de 2013 12:09, Jônatan [via R-br] <
ml-node+s2285057n4659201h99 em n4.nabble.com> escreveu:

> 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 <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659201&i=0>
> >
>
>> 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] <[hidden email]<http://user/SendEmail.jtp?type=node&node=4659201&i=1>
>> > 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.
>>> 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=4659201&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=4659201&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-tp4659130p4659201.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/8654c64c/attachment-0001.html>
-------------- Próxima Parte ----------
long	lat	
-38.13025	-9.23351
-42.06612	-9.37441
-43.10941	-11.03425
-44.32118	-10.99939
-45.00786	-12.0727
-47.50827	-12.32703
-46.235	-14.06068
-47.19565	-15.30697
-52.14593	-15.53839
-50.58539	-17.47698
-47.55472	-19.44578
-45.01159	-22.06289
-42.40095	-21.23248
-41.51141	-20.31369
-41.0361	-19.28989
-41.55565	-18.50874
-41.30522	-17.51506
-39.2054	-17.45865
-38.88361	-15.85488
-38.90264	-13.64869
-38.18265	-12.56051
-38.57341	-11.38175
-38.2983	-11.04811
-39.17153	-10.25997
-38.13025	-9.2335


Mais detalhes sobre a lista de discussão R-br