[Proj] converting WGS84 to NAD27

Anna Hentzel wibge at gmail.com
Sun Oct 4 15:22:03 PDT 2009


Hi,

I'm trying to convert WGS84 coordinates to UTM and NAD27. I'm using
the route-me library, which is a thin wrapper on proj4. It just
creates the PJ objects using pj_init_plus on the strings, and calls
pg_transform. The transformation to UTM works correctly, but the NAD27
transformation doesn't change the coordinates. Does anyone know what I
might be doing wrong?

CLLocationCoordinate2D coord = {60, -119};
int zone = ((coord.longitude + 186) / 6);

RMProjection* latLonProj = [[RMProjection alloc]
initWithString:@"+proj=latlong +ellps=WGS84 +datum=WGS84"];
RMProjection* utmProj = [[RMProjection alloc] initWithString:
[NSString stringWithFormat:@"+proj=utm +zone=%d +ellps=WGS84", zone]];

RMTransform* trans = [[RMTransform alloc] initFrom:latLonProj To:utmProj];

CLLocationCoordinate2D transCoord = [trans projectForward:coord AtZoom:9];
NSLog(@"UTM Coord %f, %f %d" , transCoord.latitude, transCoord.longitude, zone);

RMProjection* nad27Proj = [[RMProjection alloc]
initWithString:@"+proj=latlong +ellps=WGS84 +datum=NAD27 +no_defs"];
trans = [[RMTransform alloc] initFrom:latLonProj To:nad27Proj];
transCoord = [trans projectForward:coord AtZoom:9];
NSLog(@"NAD27 Coord %f, %f" , transCoord.latitude, transCoord.longitude);

**output**
UTM Coord 6653097.435298, 388455.958023 11
NAD27 Coord 60.000000, -119.000000


Thanks,
Anna



More information about the Proj mailing list