[Proj] cs2cs latlong-GB

[Frank Warmerdam]

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