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 <[email protected]> > 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] < > [email protected]> 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 > > _______________________________________________ > R-br mailing list > [email protected] > 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 [email protected] 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.
