[Proj] Differences in reprojection between using the cs2cs app and using the api
Adri CS
acsantome at gmail.com
Wed Feb 11 12:07:59 EST 2015
Hi,
I built Proj4 4.8.0 for Windows with MinGW. I'm trying to reproject some
WGS84 lat/long coordinates (degrees) from Brazil to UTM Sirgas 2000. The
UTM zone would be 23S.
**) Using the cs2cs app:
*cs2cs.exe -f %.6f +init=epsg:3273 +proj=latlong +to +init=epsg:31938
input.txt > output.txt*
**The input file would contain coords like this: -45.491596 -23.570147
769.750000
And the output would be: 449834.145105 7393276.784821 769.750000
/***********************************************/
**) Using the api (I omit the error checking):
*projPJ wgs = pj_init_plus("+init=epsg:32723 +proj=latlong");*
*projPJ sirgas = = pj_init_plus("+init=epsg:31983");*
*double lat = -45.491596;*
*double lon = -23.570147;double alt = 769.750000;*
*lat *= DEG_TO_RAD;*
*lon *= DEG_TO_RAD;*
*pj_transform(wgs, sirgas, 1, 1, &lon, &lat, &alt);*
This outputs: 449834.127632 7393276.821296 769.750000
In both cases, the left part of the returned UTM easting and northing
coordinates are the same, but the numbers after the decimal point are
different.
I'm only transforming the input lat/lon coordintates when using the api,
since if I transform them before calling cs2cs, it will output wron coords
(the easting will have 7 digits before the decimal points).
What's happening here? The only difference I can see it's the
transformation of the lat/lon coordinates from degrees to radians: I'm
multypling directly while in the cs2cs source, there's a call to the
*dmstor* function (which I can't call directly, it seems).
Thanks for your time & cheers!
Adri.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20150211/a80c4ea6/attachment.htm
More information about the Proj
mailing list