[Proj] help with proj.4

zagot maurelio at live.it
Mon Jun 13 07:36:53 PDT 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
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
+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
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68";

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

thanks

sorry for my english

m.



--
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