Alguém poderia me explicar como faço a converção de utm para lat lon, estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão Tito Conte
Fala Tito ai vai CRM corrigida: library(rgdal) #1.1 buscando dados em utm dados.utm=read.table("k:/2013/r/transforma/dados.utm.txt", head = T,blank.lines.skip=T) #1.2 Transformando em "SpatialPointsDataFrame" dados.utm.spdf <- SpatialPointsDataFrame(dados.utm[,3:4], data=dados.utm[2],proj4string= CRS("+proj=utm +zone=23+ x_0=0 +y_0=10000000 +south +datum=WGS84")) #1.3 transformando utm em longlat dados.longlat <- spTransform(dados.utm.spdf , CRS("+proj=longlat +datum=WGS84"))#+ellps=WGS84"))#+lat_0=0 +lon_0=0")) dados.longlat abraço Hélio 2013/6/3 TitoConte [via R-br] <ml-node+s2285057n4659529h67@n4.nabble.com>
Alguém poderia me explicar como faço a converção de utm para lat lon,
estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão
Tito Conte
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4659529&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-Coordenadas-UTM-para-Lat-Lon-tp465952... 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
Vc precisa informar a projeção no formato Proj4 na função project() On 03/06/13 19:08, Tito Conte wrote:
Alguém poderia me explicar como faço a converção de utm para lat lon,
estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão
Tito Conte
Veja se isso funciona: UTM.to.Lat_Long<-function(coord){ b=6356752.3142 a=6378137 e=0.0818191909 eisq=0.0067394968 k0=0.9996 ei=(1-(1-e*e)^(1/2))/(1+(1-e*e)^(1/2)) C1=3*ei/2-27*ei^3/32 C2=21*ei^2/16-55*ei^4/32 C3=151*ei^3/96 C4=1097*ei^4/512 polo=coord[,1] zone=coord[,2] lon=coord[,3] lat=coord[,4] Corrected.Northing=ifelse(polo=="N",lat,10000000-lat) East.Prime=500000-lon Arc.Length=lat/k0 mu=Arc.Length/(a*(1-e^2/4-3*e^4/64-5*e^6/256)) phl=mu+C1*sin(2*mu)+C2*sin(4*mu)+C3*sin(6*mu)+C4*sin(8*mu) C=eisq*cos(phl)^2 T=tan(phl)^2 N=a/(1-(e*sin(phl))^2)^(1/2) R=a*(1-e*e)/(1-(e*sin(phl))^2)^(3/2) D=East.Prime/(N*k0) F1=N*tan(phl)/R F2=D*D/2 F3=(5+3*T+10*C-4*C*C-9*eisq)*D^4/24 F4=(61+90*T+298*C+45*T*T-252*eisq-3*C*C)*D^6/720 LF1=D LF2=(1+2*T+C)*D^3/6 LF3=(5-2*C+28*T-3*C^2+8*eisq+24*T^2)*D^5/120 delta.lon=(LF1-LF2+LF3)/cos(phl) zone.cm=6*zone-183 raw.lat=180*(phl-F1*(F2+F3+F4))/pi LAT=ifelse(polo=="N",raw.lat,-1*raw.lat) LON=zone.cm-delta.lon*180/pi lon.deg=as.integer(LON) lon.min=abs(as.integer(60*(LON-lon.deg))) lon.sec=3600*(abs(LON)-abs(lon.deg)-lon.min/60) lat.deg=as.integer(LAT) lat.min=abs(as.integer(60*(LAT-lat.deg))) lat.sec=3600*(abs(LAT)-abs(lat.deg)-lat.min/60) df=data.frame(coord,LON,LAT,lon.deg,lon.min,lon.sec,lat.deg,lat.min,lat.sec) return(df) } coord=data.frame(polo=sample(c("N","S"),100,r=T),zone=sample(0:23,100,r=T),lon=sample(0:90000,100,r=T),lat=sample (0:90000,100,r=T)) UTM.to.Lat_Long(coord) Em 4 de junho de 2013 09:21, Elias Teixeira Krainski < eliaskrainski@yahoo.com.br> escreveu:
Vc precisa informar a projeção no formato Proj4 na função project()
On 03/06/13 19:08, Tito Conte wrote:
Alguém poderia me explicar como faço a converção de utm para lat lon,
estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão
Tito Conte
______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.
-- *Luciano F. Sgarbi* * Mestre em Ecologia e Evolução Dep. Ecologia, Instituto de Ciências Biológicas I* *Universidade Federal de Goiás, campus II, Goiânia-GO, Brasil*
Obrigado a todos! Tito Conte 2013/6/4 Luciano F. Sgarbi <luciano.f.sgarbi@gmail.com>
Veja se isso funciona:
UTM.to.Lat_Long<-function(coord){
b=6356752.3142 a=6378137 e=0.0818191909 eisq=0.0067394968 k0=0.9996 ei=(1-(1-e*e)^(1/2))/(1+(1-e*e)^(1/2)) C1=3*ei/2-27*ei^3/32 C2=21*ei^2/16-55*ei^4/32 C3=151*ei^3/96 C4=1097*ei^4/512 polo=coord[,1] zone=coord[,2] lon=coord[,3] lat=coord[,4] Corrected.Northing=ifelse(polo=="N",lat,10000000-lat) East.Prime=500000-lon Arc.Length=lat/k0 mu=Arc.Length/(a*(1-e^2/4-3*e^4/64-5*e^6/256)) phl=mu+C1*sin(2*mu)+C2*sin(4*mu)+C3*sin(6*mu)+C4*sin(8*mu) C=eisq*cos(phl)^2 T=tan(phl)^2 N=a/(1-(e*sin(phl))^2)^(1/2) R=a*(1-e*e)/(1-(e*sin(phl))^2)^(3/2) D=East.Prime/(N*k0) F1=N*tan(phl)/R F2=D*D/2 F3=(5+3*T+10*C-4*C*C-9*eisq)*D^4/24 F4=(61+90*T+298*C+45*T*T-252*eisq-3*C*C)*D^6/720 LF1=D LF2=(1+2*T+C)*D^3/6 LF3=(5-2*C+28*T-3*C^2+8*eisq+24*T^2)*D^5/120 delta.lon=(LF1-LF2+LF3)/cos(phl) zone.cm=6*zone-183 raw.lat=180*(phl-F1*(F2+F3+F4))/pi LAT=ifelse(polo=="N",raw.lat,-1*raw.lat) LON=zone.cm-delta.lon*180/pi lon.deg=as.integer(LON) lon.min=abs(as.integer(60*(LON-lon.deg))) lon.sec=3600*(abs(LON)-abs(lon.deg)-lon.min/60) lat.deg=as.integer(LAT) lat.min=abs(as.integer(60*(LAT-lat.deg))) lat.sec=3600*(abs(LAT)-abs(lat.deg)-lat.min/60)
df=data.frame(coord,LON,LAT,lon.deg,lon.min,lon.sec,lat.deg,lat.min,lat.sec) return(df) }
coord=data.frame(polo=sample(c("N","S"),100,r=T),zone=sample(0:23,100,r=T),lon=sample(0:90000,100,r=T),lat=sample (0:90000,100,r=T))
UTM.to.Lat_Long(coord)
Em 4 de junho de 2013 09:21, Elias Teixeira Krainski < eliaskrainski@yahoo.com.br> escreveu:
Vc precisa informar a projeção no formato Proj4 na função project()
On 03/06/13 19:08, Tito Conte wrote:
Alguém poderia me explicar como faço a converção de utm para lat lon,
estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão
Tito Conte
______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.
--
*Luciano F. Sgarbi* * Mestre em Ecologia e Evolução Dep. Ecologia, Instituto de Ciências Biológicas I* *Universidade Federal de Goiás, campus II, Goiânia-GO, Brasil*
_______________________________________________ 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.
participantes (4)
-
Elias Teixeira Krainski -
Hélio Gallo Rocha -
Luciano F. Sgarbi -
Tito Conte