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: heliogallorocha@gmail.com To: r-br@listas.c3sl.ufpr.br 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] <ml-node+s2285057n4659130h98@n4.nabble.com> 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] 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 -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.

Olá, Curiosidade: Arthur, onde vc viu que (-17,-50) (WGS84?) seria W606447.747 e N8120173.351 (Qual ellps?)? Segundo http://www.spatialreference.org/ref/epsg/29182/ a coordenada (-51, -19.065) (WGS84) seria (em utm, ellps=GRS67): (500000, 7891973.267179) em R temos: require(rgdal) coo1 <- SpatialPoints(matrix(c(-51, -19.065), 1), proj4string=CRS("+proj=longlat +zone=22 +south +ellps=WGS84")) spTransform(coo1, CRS("+proj=utm +zone=22 +south +ellps=GRS67 +units=m +no_defs")) att. Elias. On 04/26/2013 01:09 PM, Hélio Gallo Rocha wrote:

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] < ml-node+s2285057n4659134h22@n4.nabble.com> 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<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>

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.

Olá Elias, inverter não adiantou, mas seu ultimo comando funcionou perfeitamente, não analisei direito qual a diferença, mas fez... Uma pergunta E para os pontos que não estão na zona 22,pois tem pontos nas zonas 22, 23 e 24??? pode converter todos usando a 22? vou testar Em 26 de abril de 2013 14:09, Elias Teixeira Krainski [via R-br] < ml-node+s2285057n4659140h8@n4.nabble.com> escreveu:

ah... basta usar (-50,-17) em vez de (-17,-50)

On 04/26/2013 01:09 PM, Hélio Gallo Rocha wrote:

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=4659140&i=0>> 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> <http://user/SendEmail.jtp?type=node&node=4659134&i=0>

To: [hidden email] <http://user/SendEmail.jtp?type=node&node=4659134&i=1> <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> <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> <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> <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> <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> <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> <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%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> <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>

_______________________________________________ R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4659140&i=1>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=4659140&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.

------------------------------ 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-tp4659130p4659140.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

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: tcalegario@hotmail.com To: r-br@listas.c3sl.ufpr.br 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. heheArthur Date: Fri, 26 Apr 2013 13:09:46 -0300 From: heliogallorocha@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Krigagem ArthurPelo 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")) Coordutm1 <- 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] <ml-node+s2285057n4659134h22@n4.nabble.com> 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] To: [hidden email] 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]> 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] 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 -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho _______________________________________________ R-br mailing list [hidden email] 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] 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 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 -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.

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

Caro Hélio, Não sou especialista na área, mas acho que não é bom separar as coordenadas de cada zona, projetar considerando cada zona de referência e depois juntá-las novamente para fazer a análise. Pois, dessa forma, as distâncias entre locais de uma zona a locais de outra estarão erradas. Att. Elias. On Sat, 2013-04-27 at 21:12 -0300, Hélio Gallo Rocha wrote:

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@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@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] To: [hidden email] 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] To: [hidden email] 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] To: [hidden email] 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]> 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] To: [hidden email] 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]> 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] 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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

______________________________________________ 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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.

Uma possivel alternativa é manter as coordenadas lat-log e usar funcoes que calculam distancias no geoide como a fields:::rdist.earth On Sat, 27 Apr 2013, Hélio Gallo Rocha wrote:

Concordo Elias. Esse é um problema que não conheço a resolução, de ter pontos em diversos fusos, por isso perguntei se alguém tinha uma luz. Teria de ver qual o algoritimo que os programas de topografia fazem esse cálculo...não sei.

O Arthur vai ter de fazer a krigagem usando latlong mesmo, mas isso não é problema, já que ainda são dados em um plano cartesiano.

Hélio

Em 27 de abril de 2013 21:35, Elias Teixeira Krainski [via R-br] <ml-node+s2285057n4659154h42@n4.nabble.com> escreveu: Caro Hélio,

Não sou especialista na área, mas acho que não é bom separar as coordenadas de cada zona, projetar considerando cada zona de referência e depois juntá-las novamente para fazer a análise. Pois, dessa forma, as distâncias entre locais de uma zona a locais de outra estarão erradas.

Att. Elias.

On Sat, 2013-04-27 at 21:12 -0300, Hélio Gallo Rocha wrote:

...

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 <[hidden email]> 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         +[hidden email]> 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]                 To: [hidden email]                 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]                 To: [hidden email]                 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]                 To: [hidden email]                 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]> 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]                         To: [hidden email]                         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]> 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]                                 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                                                                                                                         --                         Hélio Gallo Rocha                         IFSULDEMINAS - Câmpus Muzambinho                                                                         _______________________________________________ R-br mailing list [hidden email] 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]                         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                                                                                                 ______________________________________________                         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                                                                                 --                 Hélio Gallo Rocha                 IFSULDEMINAS - Câmpus Muzambinho                                                 _______________________________________________ R-br                 mailing list [hidden email]                 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]                 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.                 NAML                                         --         Hélio Gallo Rocha         IFSULDEMINAS - Câmpus Muzambinho        

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

_______________________________________________ R-br mailing list [hidden email] 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] 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-tp4659130p4659154.html To unsubscribe from R-br, click here. NAML

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

Boa sugestão Jônatan, Ja usei essa opção recomendada por um colega. Porém, eu uso lat_0 = ponto médio da latitude lon_0 = ponto médio da longitude este ponto médio calculado com base apenas no menor e menor valor de cada coordenada dos dados. Assim, a região de estudo será centrada em (0,0). Outra alternativa é usar funções que calculam distâncias euclieanas a partir de coordenadas longlat. O Sudipto Banerjee tem um artigo avaliando o efeito dessas coias (usar distancias em graus ou em unidades de distância) em modelos geoestatísticos. Na página dele http://www.biostat.umn.edu/~sudiptob/ a partir de um link em 'softwares', ver em 'Other R Programs' algumas funções para isso. Att. Elias. On Sat, 2013-04-27 at 22:42 -0300, Jônatan wrote:

+lat_0=-15 +lon_0=-60

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 <heliogallorocha@gmail.com>

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

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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 ###############################################################

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] < ml-node+s2285057n4659167h45@n4.nabble.com> 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

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

Prezados, Estou começando a usar or programa R e preciso fazer curva de sobrevivência (Kaplan-Meier) sobre o acompanhamento de pacientes ao longo do tempo. Gentileza, se alguém tiver conhecimento ou algum tutorial sobre essa análise, o mais detalhado possível, agradeceria. Atenciosamente, -- Labcor <http://www.labcor.com.br/> Ana Luiza Cassin Monitora de Pesquisa Monitor Research Tel.com.: +55(31)3313-1251 Tel.cel.: +55(31)9308-0072 Labcor www.labcor.com.br <http://www.labcor.com.br/>

#Oi segue abaixo um exemplo, espero ter ajudado! #Alexandro #modelo tempo=c(40,30,20,10,5,45,35,25,15,5) censura=c(1,1,1,1,1,1,1,0,1,1) sexo=c("m","f","m","f","m","f","m","m","m","f") doente=factor(c(1,1,0,0,0,0,1,1,1,1)) idade=c(20,30,40,50,60,25,35,45,55,65) dados=data.frame(tempo,censura,sexo,doente,idade);dados require(survival) .Survfit <- survfit(Surv(tempo, censura) ~ doente + sexo, conf.type="log", conf.int=0.95, type="kaplan-meier", error="greenwood", data=dados) summary(.Survfit) plot(.Survfit, col=1:4, lty=1:4, mark.time=TRUE) legend("bottomleft", legend=c("0.f","1.f","0.m","1.m"), title="doente + sexo", col=1:4, lty=1:4, bty="n")#INTERACAO

________________________________ De: Ana Luiza Cassin <pesquisaclinica@labcor.com.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 6 de Maio de 2013 15:35 Assunto: [R-br] Survival

Prezados,

Estou começando a usar or programa R e preciso fazer curva de sobrevivência (Kaplan-Meier) sobre o acompanhamento de pacientes ao longo do tempo. Gentileza, se alguém tiver conhecimento ou algum tutorial sobre essa análise, o mais detalhado possível, agradeceria.

Atenciosamente,

--

Ana Luiza Cassin Monitora de Pesquisa Monitor Research Tel.com.: +55(31)3313-1251 Tel.cel.: +55(31)9308-0072   Labcor www.labcor.com.br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.

Prezados, Alexadro e Daniel, agradeço. Com relação a importar planilha com os dados, queria saber se é possível ser do excell, pois tentei com os comandos mas nao permitiu. Grata, Ab Ana Em 06/05/2013 19:37, Alexandro (Yahoo) escreveu:

#Oi segue abaixo um exemplo, espero ter ajudado! #Alexandro

#modelo tempo=c(40,30,20,10,5,45,35,25,15,5) censura=c(1,1,1,1,1,1,1,0,1,1) sexo=c("m","f","m","f","m","f","m","m","m","f") doente=factor(c(1,1,0,0,0,0,1,1,1,1)) idade=c(20,30,40,50,60,25,35,45,55,65) dados=data.frame(tempo,censura,sexo,doente,idade);dados

require(survival) .Survfit <- survfit(Surv(tempo, censura) ~ doente + sexo, conf.type="log", conf.int=0.95, type="kaplan-meier", error="greenwood", data=dados) summary(.Survfit) plot(.Survfit, col=1:4, lty=1:4, mark.time=TRUE) legend("bottomleft", legend=c("0.f","1.f","0.m","1.m"), title="doente + sexo", col=1:4, lty=1:4, bty="n")#INTERACAO

------------------------------------------------------------------------ *De:* Ana Luiza Cassin <pesquisaclinica@labcor.com.br> *Para:* r-br@listas.c3sl.ufpr.br *Enviadas:* Segunda-feira, 6 de Maio de 2013 15:35 *Assunto:* [R-br] Survival

Prezados,

Estou começando a usar or programa R e preciso fazer curva de sobrevivência (Kaplan-Meier) sobre o acompanhamento de pacientes ao longo do tempo. Gentileza, se alguém tiver conhecimento ou algum tutorial sobre essa análise, o mais detalhado possível, agradeceria.

Atenciosamente, -- Labcor <http://www.labcor.com.br/> Ana Luiza Cassin Monitora de Pesquisa Monitor Research Tel.com.: +55(31)3313-1251 Tel.cel.: +55(31)9308-0072 Labcor www.labcor.com.br <http://www.labcor.com.br/>

-- Labcor <http://www.labcor.com.br/> Ana Luiza Cassin Monitora de Pesquisa Monitor Research Tel.com.: +55(31)3313-1251 Tel.cel.: +55(31)9308-0072 Labcor www.labcor.com.br <http://www.labcor.com.br/>

o rgdal possui dependencias externas que voce tem que instalar. Vaja pa página e documentacao do pacote se estiver usando linux (com apt-get) apt-get --build-dep rgdal On Mon, 6 May 2013, Arthur Telles Calegario wrote:

ow Galera,Comecei a ver esse CMR contudo não consigo instalar o  pacote rgdal. Li nele que é necessário instalar um tal de Gdal. É isso mesmo? Vcs podem me adiantar algo? Agradeço novamente pelo auxílio na resolução do problema. Arthur.

_______________________________________________________________________________________________________________________________________ Date: Tue, 30 Apr 2013 17:24:02 -0300 From: heliogallorocha@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Krigagem

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]> Caro JônatanFiz 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]> 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]> 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]> 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]> 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] To: [hidden email] 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] To: [hidden email] 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] To: [hidden email] Subject: Re: [R-br] Krigagem

ArthurPelo 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]> 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] To: [hidden email] 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]> 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] 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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

[IMAGE] Dados.txt (2K) Download Attachment

_______________________________________________________________________________________________________________________________________ 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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.

OK, atente para as dependencias externas On Mon, 6 May 2013, Arthur Telles Calegario wrote:

Hm.Não uso linux. Mas vou procurar algo na net pra tenar instalar.

Date: Mon, 6 May 2013 15:17:17 -0300 From: paulojus@leg.ufpr.br To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Krigagem

o rgdal possui dependencias externas que voce tem que instalar. Vaja pa página e documentacao do pacote

se estiver usando linux (com apt-get)

apt-get --build-dep rgdal

On Mon, 6 May 2013, Arthur Telles Calegario wrote:

...

ow Galera,Comecei a ver esse CMR contudo não consigo instalar o  pacote rgdal. Li nele que é necessário instalar um tal de Gdal. É isso mesmo? Vcs podem me adiantar algo? Agradeço novamente pelo auxílio na resolução do problema. Arthur.

_____________________________________________________________________________________________________________________________________ __ Date: Tue, 30 Apr 2013 17:24:02 -0300 From: heliogallorocha@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Krigagem

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 f oi 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]> Caro JônatanFiz 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]> 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]> 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]> 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]> 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] To: [hidden email] 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] To: [hidden email] 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] To: [hidden email] Subject: Re: [R-br] Krigagem

ArthurPelo 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]> 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] To: [hidden email] 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]> 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] 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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

[IMAGE] Dados.txt (2K) Download Attachment

_____________________________________________________________________________________________________________________________________ __ 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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

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

-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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???d igo m???nimo reproduz???vel.

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.