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_spdfnã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
kc.predict Var1 Var2 1439.471584 -749877.155 -712391.6699 1465.329422 -581060.704 -712391.6699 1495.849513 -412244.253 -712391.6699 1532.559175 -243427.802 -712391.6699 1570.512733 -74611.351 -712391.6699 1565.110096 94205.1 -712391.6699 1457.859239 263021.551 -712391.6699 1345.221426 431838.002 -712391.6699 1348.304618 600654.453 -712391.6699 1362.765892 769470.904 -712391.6699 1474.279885 -749877.155 -555353.3091 1500.77349 -581060.704 -555353.3091 1527.140626 -412244.253 -555353.3091 1561.443606 -243427.802 -555353.3091 1521.555604 -74611.351 -555353.3091 1474.766274 94205.1 -555353.3091 1256.024327 263021.551 -555353.3091 1313.511945 431838.002 -555353.3091 1327.413659 600654.453 -555353.3091 1373.914461 769470.904 -555353.3091 1522.508441 -749877.155 -398314.9483 1551.530978 -581060.704 -398314.9483 1535.88269 -412244.253 -398314.9483 1575.61311 -243427.802 -398314.9483 1421.226698 -74611.351 -398314.9483 1321.329377 94205.1 -398314.9483 1335.394278 263021.551 -398314.9483 1159.617269 431838.002 -398314.9483 1302.221988 600654.453 -398314.9483 1409.646779 769470.904 -398314.9483 1572.584981 -749877.155 -241276.5875 1672.958126 -581060.704 -241276.5875 1601.200514 -412244.253 -241276.5875 1503.274484 -243427.802 -241276.5875 1414.150219 -74611.351 -241276.5875 1321.216126 94205.1 -241276.5875 1178.649369 263021.551 -241276.5875 1093.982025 431838.002 -241276.5875 1313.436384 600654.453 -241276.5875 1493.859945 769470.904 -241276.5875 1575.311193 -749877.155 -84238.22665 1597.079043 -581060.704 -84238.22665 1772.743771 -412244.253 -84238.22665 1556.435251 -243427.802 -84238.22665 1375.731907 -74611.351 -84238.22665 1181.369146 94205.1 -84238.22665 932.7270225 263021.551 -84238.22665 908.4123845 431838.002 -84238.22665 1236.646787 600654.453 -84238.22665 1659.9622 769470.904 -84238.22665 1519.472496 -749877.155 72800.13418 1473.881724 -581060.704 72800.13418 1631.115586 -412244.253 72800.13418 1532.545916 -243427.802 72800.13418 1387.758713 -74611.351 72800.13418 1062.71752 94205.1 72800.13418 823.3179681 263021.551 72800.13418 791.6182586 431838.002 72800.13418 1508.027024 600654.453 72800.13418 1978.87457 769470.904 72800.13418 1445.485323 -749877.155 229838.495 1448.321808 -581060.704 229838.495 1478.327233 -412244.253 229838.495 1476.383982 -243427.802 229838.495 1387.620359 -74611.351 229838.495 1024.124504 94205.1 229838.495 856.8627057 263021.551 229838.495 721.0660946 431838.002 229838.495 1101.135134 600654.453 229838.495 2103.231346 769470.904 229838.495 1397.682517 -749877.155 386876.8558 1393.766889 -581060.704 386876.8558 1376.103706 -412244.253 386876.8558 1317.75947 -243427.802 386876.8558 1233.785297 -74611.351 386876.8558 1001.003757 94205.1 386876.8558 862.5257413 263021.551 386876.8558 1206.938437 431838.002 386876.8558 870.5745439 600654.453 386876.8558 1501.733103 769470.904 386876.8558 1362.758142 -749877.155 543915.2166 1349.384967 -581060.704 543915.2166 1317.517255 -412244.253 543915.2166 1258.370748 -243427.802 543915.2166 1159.094384 -74611.351 543915.2166 976.6134194 94205.1 543915.2166 670.2947772 263021.551 543915.2166 711.6220539 431838.002 543915.2166 766.1181506 600654.453 543915.2166 971.9871309 769470.904 543915.2166 1339.219433 -749877.155 700953.5775 1321.56908 -581060.704 700953.5775 1287.411876 -412244.253 700953.5775 1228.404423 -243427.802 700953.5775 1129.792915 -74611.351 700953.5775 965.1359444 94205.1 700953.5775 766.5599213 263021.551 700953.5775 734.5649371 431838.002 700953.5775 694.6837369 600654.453 700953.5775 582.9 769470.904 700953.5775
_______________________________________________ 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.
