[R-sig-Geo] UTM problem

Werner W. pensterfuzzer at yahoo.de
Wed Oct 13 15:37:00 CEST 2010


Hi,

this questions is not really about R but it would be great if somebody could 
point me to a forum or so where I can get help with this.

I only know a few basics about GIS and have a problem with some UTM coordinates 
from a survey.
Coordinates are given as, e.g.
3258843     32780
which are supposedly zone 36N on system WGS84 (according to the documentation).
If I convert it to long/lat using gdal:
> require(rgdal)
> SP <- SpatialPoints(cbind(X=3258843, Y=32780),proj4string=CRS("+proj=utm 
>+zone=36 +ellps=WGS84"))
> spTransform(SP, CRS("+proj=longlat +ellps=WGS84")) # x1=longitude x2=lat
SpatialPoints:
           X         Y
[1,] 57.0489 0.2706471
Coordinate Reference System (CRS) arguments: +proj=longlat +ellps=WGS84 

The point should be close to Kampala which has the coordinates 0.3°, 32.633333° 
and UTM 36N 459199 33160 on WGS84. My data's easting seems quite off. 

Has anybody an idea what I am doing wrongly? 

Thanks so much,
  Werner





More information about the R-sig-Geo mailing list