[Proj] convert from LatLong to transverse mercator LUREF

[Norman Vine]

Oscar van Vlijmen writes:
> 
> >> someone could help me to solve my problem. I would like to convert geodetic
> >> coordinates to a transverse mercator projection for Luxembourg.
> 
> >> I have the following coordinates:
> 
> >> Latitude: 49d34'17".60287N
> >> Longitude: 5d55'50".69323E
> 
> >> And want to convert them into:
> 
> >> X (North) 70910.00
> >> Y (East)  62935.00
> 
> > Ulrich,
> 
> >  cs2cs +proj=latlong +ellps=intl +to +proj=tmerc +x_0=80000 +y_0=100000  +lon_0=6d10 +lat_0=49d50 +ellps=intl +k=1
> > 5d55'50.69 49d34'17.60
> > 62934.93        70909.91 0.00
> 
> > Seems like pretty good agreement.
> 
> 
> I can't believe PROJ or cs2cs are that inaccurate!


Huh,  What inaccuracy ???


# Luxembourg 1930 / Gauss
<2169> +proj=tmerc +lat_0=49.83333333333334 +lon_0=6.166666666666667 +k=1.000000  \
             +x_0=80000 +y_0=100000 +ellps=intl +towgs84=-193,13.7,-39.3,-0.41,-2.933,2.688,0.43 +units=m +no_defs  <>


$ proj +proj=tmerc +lat_0=49.83333333333334 +lon_0=6.166666666666667 +k=1.000000 +x_0=80000 +y_0=100000 +ellps=intl +to
wgs84=-193,13.7,-39.3,-0.41,-2.933,2.688,0.43 +units=m +no_defs
5d55'50.69323 49d34'17.60287
62935.00        70910.00


$ cs2cs +proj=latlong +ellps=intl +to +proj=tmerc +x_0=80000 +y_0=100000  +lon_0=6d10 +lat_0=49d50 +ellps=intl +k=1
5d55'50.69323 49d34'17.60287
62935.00        70910.00 0.00


HTH

Norman