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@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