Re: [R-br] coordenadas geograficas
Natalia e colegas, bom dia! Ainda não estou muito familiarizado com o R, mas acho que este exemplo que rodei irá ajudar... Não utilizei o project(), mas sim o spTransform()... ############################################################# #install.packages("rgdal", dep=T) #install.packages("sp", dep=T) require(sp); require(rgdal) ### long coord.x = c(-54.607845, -54.599831, -54.430241, -54.258459, -54.344408, -54.333258, -53.607359, -53.589245, -53.546762, -53.105123, -52.972362, -52.107579, -52.089572, -51.940453, -51.718151, -51.263728, -50.741531, -50.734130, -50.667779, -50.662230, -49.986102, -49.970018, -49.901474, -49.911543, -49.678260, -49.654858, -49.610676, -49.599658, -49.511570, -49.200709, -49.222508, -49.317621, -49.311839, -49.275882, -48.581791, -48.578971, -48.582604, -48.568999, -48.528527, -48.411012, -48.251622, -48.027824, -48.029886, -48.023501, -48.590796, -48.594522, -48.643601, -49.175230, -49.555033, -49.567207, -49.689784, -49.941462, -49.942167, -50.571750, -50.591636, -50.553035, -50.718453, -51.208345, -51.216170, -51.292329, -51.239319, -51.258597, -51.283372, -51.411480, -53.551638, -53.552753, -53.643791, -53.835953, -53.842011, -53.858677, -53.955544, -54.099576, -54.108952, -54.429627, -54.438295, -54.593660, -54.593683, -54.619090, -54.617907, -54.607845) ### lat coord.y = c(-25.446263, -25.434436, -25.159459, -24.357597, -24.140630, -24.116248, -22.951327, -22.936230, -22.890090, -22.621923, -22.570479, -22.516685, -22.521237, -22.567335, -22.669007, -22.669047, -22.961387, -22.962502, -22.897375, -22.895686, -22.897501, -22.910601, -22.999147, -23.051110, -23.165054, -23.210043, -23.851684, -23.857080, -23.925858, -24.344464, -24.362727, -24.566630, -24.663615, -24.693608, -24.671609, -24.685526, -25.051554, -25.051240, -25.100784, -24.980119, -24.978739, -25.221284, -25.226189, -25.230964, -25.976982, -25.977087, -25.955943, -26.001844, -26.237326, -26.233284, -26.188456, -26.008879, -26.008876, -26.003790, -26.022134, -26.055436, -26.244779, -26.297107, -26.313987, -26.404261, -26.628383, -26.625825, -26.653990, -26.717283, -26.292744, -26.292051, -26.253072, -25.971930, -25.687925, -25.680270, -25.646984, -25.495677, -25.495742, -25.695430, -25.691293, -25.592352, -25.544262, -25.457555, -25.452072, -25.446263) pr.df <- data.frame(coord.x,coord.y); head(pr.df) ### criar um objeto "SpatialPoints" pr.spgeo <- pr.df; head (pr.spgeo) coordinates(pr.spgeo) <- c('coord.x','coord.y') ### definir projeção e datum de entrada proj4string(pr.spgeo) <- CRS('+proj=longlat +datum=WGS84') class(pr.spgeo) head(pr.spgeo) plot(pr.spgeo) ### lista parâmetros disponíveis projInfo("proj") projInfo("datum") projInfo("ellps") ### transforma para utm ### atenção para o uso do parâmetro '+south' ### ter claro que o Paraná utiliza mais de uma Zona UTM e maior parte do estado em 22S pr.sputm <- spTransform(pr.spgeo, CRS('+proj=utm +zone=22 +south')); head(pr.sputm) plot(pr.sputm) ############################################################# Esse exemplo assumiu o WGS84 como datum padrão de entrada e saída. Para datums (data) diferentes tem que dar uma estudada mais a fundo, porque tem poucas opções pré-configuradas no pacote que usei (como pode ser visto usando <projInfo("datum")>). -- ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00
----------------------------------------------------------------------
2012/10/29 Natalia Martins <nsmbarreto@gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=nsmbarreto@gmail.com>
Walmes, bom dia. Estou tentando transformar minhas coordenadas geograficas (das estações metereologicas do PR) em utm, procurando no help do R encontrei:
data(state) s <- project(state.center, "+proj=merc")
mas eu não estou conseguindo entender qual proj eu devo usar, no link
http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html
até encontri mas quando tento fazer aparece o seguinte erro
s <- project(state.center, "+proj=utm") Erro em project(state.center, "+proj=utm") : elliptical usage required
se vc puder me ajudar...
Muito obrigada!!
Em 28 de outubro de 2012 11:11, Walmes Zeviani <walmes@ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br>> escreveu:
Sim, para calculo de distâncias você tem que ter coordenadas decimais e
não exagesimais (º ' '').
==========================================================================
Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
==========================================================================
--
Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743
Mensagem Verificada Pelo AntiSpam UFPR
Muito obrigada pessoal, deu certo.. Em 31 de outubro de 2012 09:19, Eder Comunello <ecomunel@gmail.com>escreveu:
Natalia e colegas, bom dia!
Ainda não estou muito familiarizado com o R, mas acho que este exemplo que rodei irá ajudar...
Não utilizei o project(), mas sim o spTransform()...
#############################################################
#install.packages("rgdal", dep=T) #install.packages("sp", dep=T) require(sp); require(rgdal)
### long coord.x = c(-54.607845, -54.599831, -54.430241, -54.258459, -54.344408, -54.333258, -53.607359, -53.589245, -53.546762, -53.105123, -52.972362, -52.107579, -52.089572, -51.940453, -51.718151, -51.263728, -50.741531, -50.734130, -50.667779, -50.662230, -49.986102, -49.970018, -49.901474, -49.911543, -49.678260, -49.654858, -49.610676, -49.599658, -49.511570, -49.200709, -49.222508, -49.317621, -49.311839, -49.275882, -48.581791, -48.578971, -48.582604, -48.568999, -48.528527, -48.411012, -48.251622, -48.027824, -48.029886, -48.023501, -48.590796, -48.594522, -48.643601, -49.175230, -49.555033, -49.567207, -49.689784, -49.941462, -49.942167, -50.571750, -50.591636, -50.553035, -50.718453, -51.208345, -51.216170, -51.292329, -51.239319, -51.258597, -51.283372, -51.411480, -53.551638, -53.552753, -53.643791, -53.835953, -53.842011, -53.858677, -53.955544, -54.099576, -54.108952, -54.429627, -54.438295, -54.593660, -54.593683, -54.619090, -54.617907, -54.607845)
### lat coord.y = c(-25.446263, -25.434436, -25.159459, -24.357597, -24.140630, -24.116248, -22.951327, -22.936230, -22.890090, -22.621923, -22.570479, -22.516685, -22.521237, -22.567335, -22.669007, -22.669047, -22.961387, -22.962502, -22.897375, -22.895686, -22.897501, -22.910601, -22.999147, -23.051110, -23.165054, -23.210043, -23.851684, -23.857080, -23.925858, -24.344464, -24.362727, -24.566630, -24.663615, -24.693608, -24.671609, -24.685526, -25.051554, -25.051240, -25.100784, -24.980119, -24.978739, -25.221284, -25.226189, -25.230964, -25.976982, -25.977087, -25.955943, -26.001844, -26.237326, -26.233284, -26.188456, -26.008879, -26.008876, -26.003790, -26.022134, -26.055436, -26.244779, -26.297107, -26.313987, -26.404261, -26.628383, -26.625825, -26.653990, -26.717283, -26.292744, -26.292051, -26.253072, -25.971930, -25.687925, -25.680270, -25.646984, -25.495677, -25.495742, -25.695430, -25.691293, -25.592352, -25.544262, -25.457555, -25.452072, -25.446263)
pr.df <- data.frame(coord.x,coord.y); head(pr.df)
### criar um objeto "SpatialPoints" pr.spgeo <- pr.df; head (pr.spgeo) coordinates(pr.spgeo) <- c('coord.x','coord.y')
### definir projeção e datum de entrada proj4string(pr.spgeo) <- CRS('+proj=longlat +datum=WGS84') class(pr.spgeo) head(pr.spgeo) plot(pr.spgeo)
### lista parâmetros disponíveis projInfo("proj") projInfo("datum") projInfo("ellps")
### transforma para utm ### atenção para o uso do parâmetro '+south' ### ter claro que o Paraná utiliza mais de uma Zona UTM e maior parte do estado em 22S pr.sputm <- spTransform(pr.spgeo, CRS('+proj=utm +zone=22 +south')); head(pr.sputm) plot(pr.sputm)
#############################################################
Esse exemplo assumiu o WGS84 como datum padrão de entrada e saída. Para datums (data) diferentes tem que dar uma estudada mais a fundo, porque tem poucas opções pré-configuradas no pacote que usei (como pode ser visto usando <projInfo("datum")>).
--
================================================ Éder Comunello
Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]
Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00
----------------------------------------------------------------------
2012/10/29 Natalia Martins <nsmbarreto@gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=nsmbarreto@gmail.com>
Walmes, bom dia. Estou tentando transformar minhas coordenadas geograficas (das estações metereologicas do PR) em utm, procurando no help do R encontrei:
data(state) s <- project(state.center, "+proj=merc")
mas eu não estou conseguindo entender qual proj eu devo usar, no link
http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html
até encontri mas quando tento fazer aparece o seguinte erro
s <- project(state.center, "+proj=utm") Erro em project(state.center, "+proj=utm") : elliptical usage required
se vc puder me ajudar...
Muito obrigada!!
Em 28 de outubro de 2012 11:11, Walmes Zeviani <walmes@ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br>> escreveu:
Sim, para calculo de distâncias você tem que ter coordenadas decimais e
não exagesimais (º ' '').
==========================================================================
Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
==========================================================================
--
Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743
Mensagem Verificada Pelo AntiSpam UFPR
_______________________________________________ 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.
-- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743
participantes (2)
-
Eder Comunello -
Natalia Martins