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 _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo