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