[Proj] help with proj.4

zagot maurelio at live.it
Mon Jun 13 09:36:53 EST 2011

Hi to all

if i run this with cs2cs

cs2cs +proj=tmerc +lat_0=0 +lon_0=-3.45233333333333 +k=0.999600 +x_0=1500000
+y_0=0 +ellps=intl +pm=rome +units=m
+to +proj=latlong +datum=WGS84 
1514905.27  5034551.13

i have this result

3d15'42.057"W 45d27'48.522"N 0.00

i need to do this with programmatic but with this code

char proj4Gauss[]="+proj=tmerc +lat_0=0 +lon_0=-3.45233333333333 +k=0.999600
+x_0=1500000 +y_0=0 +ellps=intl +pm=rome +units=m

char projWgs84[]="+proj=latlong +datum=WGS84 ";
projPJ projGauss = pj_init_plus(proj4Gauss);
double x=1514905.27 ;
double y=5034551.13;
double z=0;

double ris = pj_transform(projGauss, projWgs84, 1, 0, &x, &y, &z);

i have this result (in decimals point)

4423479.3363604508      715679.13023441052

waths I wrong? in the source code


sorry for my english


View this message in context: http://osgeo-org.1803224.n2.nabble.com/help-with-proj-4-tp6470395p6470395.html
Sent from the PROJ.4 mailing list archive at Nabble.com.

More information about the Proj mailing list