Acho que o correto seria informar a região na projeção UTM, não na longlat.
Assim, você teria:
Coord <- SpatialPoints(cbind(dados$coord_x,dados$coord_y),
proj4string=CRS("+proj=longlat +datum=WGS84"))
utm <- spTransform(Coord, CRS("+proj=utm +zone=22S"))
Me parece que essa coordenada fica na zona 23S...
Elias T. Krainski
>________________________________
> De: Samuel Carvalho <[email protected]>
>Para: r-br <[email protected]>
>Enviadas: Terça-feira, 15 de Maio de 2012 10:32
>Assunto: [R-br] tranformar LatLong para UTM
>
>
>Bom dia caros membros do r-br,
>Alguém poderia auxiliar na transformação de coordenadas coletadas em LatLong
>para UTM?
>Segue um cmr seguido de um print de parte dos dados (adaptado do post do
>Marcelo Costa em 20/09/2011)
>
>
>
>dados[1,]
>coord_x coord_y
>-45.33402 -23.31597
>Coord <- SpatialPoints(cbind(dados$coord_x,dados$coord_y),
>proj4string=CRS("+proj=longlat +zone=22S +datum=WGS84"))
>utm <- spTransform(Coord, CRS("+proj=utm"))
>
>
>Apesar do comando ter rodado, os resultados diferem quando comparado a um site
>de transformacao de coordenadas que encontrei na web
>(_http://home.hiwaay.net/~taylorc/toolbox/geography/geoutm.html)
>Outro ponto é que estas coordenadas constituem como uma variável do dataframe
>"dados" e com a transformação sugerida acima eu nao consigo incorporar os
>resultados como uma nova coluna no dataframe
>
>
>Grato pela atenção
>Samuel
>
>
>
>
>====================================
>Samuel P. C. Carvalho
>Mestre em Ciências Florestais [UFLA]
>Doutorando em Recursos Florestais [ESALQ/USP]
>=============================================
>_______________________________________________
>R-br mailing list
>[email protected]
>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
>mínimo reproduzível.
>
>_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
