[R-br] Krigagem

Sábado Abril 27 21:12:05 BRT 2013

Athur,

Infelizmente, desconsidere o e-mail anterior, não o CRM

Quando testei separadamenteos pontos dos fusos 22, 23 e 24, usei para todos
o valor 22, esqueci de mudar os valores dos fusos, nestas linhas:

coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat +zone=*22* +south
+ellps=WGS84"));coords
coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=*22* +south
+ellps=WGS84 +units=m +no_defs"));coords1

ou seja, usando o CRM postado, tem de separar os pontos de fuso 22, os do
fuso 23 e 24 e processar separadamente.

O CRM tá certo, mas tem de processar separadamente os pontos de cada fuso

Quem dá uma luz ???

Abraço

Hélio

Em 27 de abril de 2013 19:39, Hélio Gallo Rocha
<heliogallorocha em gmail.com>escreveu:

> Arthur, fins os testes e não precisou separar os pontos por fuso.
> testei da seguinte forma, fiz a conversão usando o arquivo txt com todos
> os pontos, depois separeis os pontos por fuso, 22,23 e 24 e comparei os
> resultados, não houve diferença.
>
> tai o CRM
>
> library(geoR)
> library(rgdal)
> dados=read.geodata("precipita.txt ",coords.col=1:2,data.col=3,sep =
> "",header=T);dados;n=summary(dados)[[1]];n##o arquivo precipita.txt é o seu
> mesmo, só renomeei
> long= dados$coords[,1] ;long
> lat=  dados$coords[,2] ;lat
> coord=cbind(long,lat);coord
> coords= SpatialPoints(coord,proj4string=CRS("+proj=longlat +zone=22 +south
> +ellps=WGS84"));coords
> coords1=coor1=spTransform(coords, CRS("+proj=utm +zone=22 +south
> +ellps=WGS84 +units=m +no_defs"));coords1
>
>
> Só por curiosidade vou fazer a krigagem usando latlong e utm, pra ver se
> dá o mesmo resultado, tem de dar...
> abraço
> Hélio
>
>
> Em 26 de abril de 2013 14:55, Arthur T. C. [via R-br] <
> ml-node+s2285057n4659145h63 em n4.nabble.com> escreveu:
>
>  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=4659145&i=0>
>> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659145&i=1>
>> 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=4659145&i=2>
>> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659145&i=3>
>> 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=4659145&i=4>
>> To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659145&i=5>
>> 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=4659145&i=6>
>> > 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=4659145&i=7>
>> 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=4659145&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.
>>
>> ------------------------------
>>  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-tp4659130p4659145.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
>

-- 
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/20130427/b0fc1ace/attachment.html>