[gdal-dev] Trasformation from proj 4

zagot maurelio at live.it
Sun Jun 5 12:57:25 EDT 2011


Hi to all,

i'm new and i have this code in proj 4

        String proiezioneCassini=""+proj=cass +ellps=bessel
+lat_0=43.3182575  +lon_0=11.33221063888889";
               
        projUV spazio;
        projPJ projCassini;

        spazio.u=32
        spazio.v=48
       
        projCassini=pj_init_plus(proiezioneCassini);

        spazio=pj_inv(spazio, projCassini);
        pj_free(projCassini);
       
        String proiezioneGauss=" +proj=tmerc +lat_0=0
+lon_0=-3.45233333333333 +k=0.999600 +x_0=1500000 +y_0=0 +ellps=intl
+pm=rome +units=m +no_defs +towgs84=-225,-65,9";
        String proiezioneWgs84="+proj=latlong +datum=WGS84";
       
        projPJ projOrigine=pj_init_plus(proiezioneWgs84);
       
        projPJ projDestinazione=pj_init_plus(proiezioneGauss);
       
        double x=spazio.u;
        double y=spazio.v;
        double z=0;
       
        int ris=pj_transform(projOrigine, projDestinazione, 1, 0, &x, &y,
&z);
        pj_free(projOrigine);
        pj_free(projDestinazione);

anyone can help me to translate this code in GDal?

at lest can you help me how can I use projUV concept in GDal

please help me

thanks

sorry for my english


--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/Trasformation-from-proj-4-tp6442626p6442626.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.


More information about the gdal-dev mailing list