[Proj] cs2cs latlong-GB

Francesco Mirabella mirabell at unipg.it
Tue Jan 17 12:05:26 EST 2006


Hi Frank,
thanks for your reply, the multiple ++ were due to cut and paste,
however: thanks, I must say that it works fine but it is funny because it actually works fine if I put long and lat values together.... It goes wrong if I put only long (12.17183333) and then only lat (43.55000000).... You can try, you will get wrong values, both longitude and latitude
It funny I didn't think of putting them together as a couple...

Thanks again,

cheers
Francesco


On Tue, Jan 17, 2006 at 11:27:19AM -0500, Frank Warmerdam wrote:
> On 1/17/06, Francesco Mirabella <mirabell at unipg.it> wrote:
> > Hi to all,
> > I am a newbie on this list, I am triyng to project some lat-long points
> > (international ellipsoid on
> > rome40 datum)( 12.17183333  43.55000000) to Gauss-Boaga Projection
> > (tmerc on
> > +international ellipsoid on rome40). I use:
> >
> >  cs2cs -E +proj=latlong +ellps=intl
> >  ++towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +to +proj=tmerc
> >  +lon_0=15
> >  ++k=0.9996 +x_0=2520000 +ellps=intl
> >   +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68  -v -f %.6f
> >
> >   I have reprojected the values to GB with grass in order to test the
> >   result (the
> >   +result shold be: 2291529.04882672   4825875.89322372 according to
> >   grass)
> 
> Francesco,
> 
> You seem to have doubled '+'es in a couple spots up
> above (++towgs84=..., and ++k).  Is that a reflection of
> what you actually did or some sort of cut and paste error?
> 
> When I do it, things seem to work OK:
> 
> warmerda at gdal2200[208]% cs2cs -E +proj=latlong +ellps=intl
> +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +to +proj=tmerc
> +lon_0=15 +k=0.9996 +x_0=2520000 +ellps=intl
> +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68  -v -f %.6f
> # ---- From Coordinate System ----
> #Lat/long (Geodetic)
> #
> # +proj=latlong +ellps=intl
> # +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
> # ---- To Coordinate System ----
> #Transverse Mercator
> #       Cyl, Sph&Ell
> # +proj=tmerc +lon_0=15 +k=0.9996 +x_0=2520000 +ellps=intl
> # +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
> 12.17183333  43.55000000
> 12.17183333  43.55000000        2291532.166325  4825868.544447 0.000000
> 
> Best regards,
> --
> ---------------------------------------+--------------------------------------
> I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
> light and sound - activate the windows | http://pobox.com/~warmerdam
> and watch the world go round - Rush    | Geospatial Programmer for Rent
> 
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj




More information about the Proj mailing list