[R] converting coordinates from utm to longitude / latitude
Werner Wernersen
pensterfuzzer at yahoo.de
Wed Aug 20 00:08:56 CEST 2008
It would be nicer to convert directly the entire shapefile object to long/lat coordinates but if that is not possible, I will convert the other points to UTM.
Hence, I am playing around with rgdal.
library(rgdal)
SP <- SpatialPoints(cbind(32.29252, -0.3228500),
proj4string=CRS("+proj=longlat"))
spTransform(SP, CRS("+proj=utm +zone=36"))
> spTransform(SP, CRS("+proj=utm +zone=36"))
SpatialPoints:
coords.x1 coords.x2
[1,] 421274.4 -35687.37
Coordinate Reference System (CRS) arguments: +proj=utm +zone=36 +ellps=WGS84
This result corresponds with what I get when using convUL() but my map of that area in UTM coordinates does not extend to the negative.
An external program converts the point to x=420994 y=9964407 which also seems correct with respect to the map. Fore sure, I am using the function wrongly somehow. Can anyone give me a hint?
That's very much appreciated!
Thanks,
Werner
----- Ursprüngliche Mail ----
Von: Werner Wernersen <pensterfuzzer at yahoo.de>
An: r-help at stat.math.ethz.ch
Gesendet: Dienstag, den 19. August 2008, 20:28:29 Uhr
Betreff: converting coordinates from utm to longitude / latitude
Hi,
is there a function in R to convert data read with read.shape and which is originally in UTM coordinates into longitude / latitude coordinates?
I found the convUL() function from the PBSmapping package but I have no idea how I could apply that to the read.shape object.
Many thanks,
Werner
______________________________
fügt über einen herausragenden Schutz gegen Massenmails.
http://mail.yahoo.com
__________________________________________________
Do
ragenden Schutz gegen Massenmails.
http://mail.yahoo.com
More information about the R-help
mailing list