[Proj] ED50_UTM to ED50_TM (Not sure!)

[Hermann Peifer]

On 23/09/2010 10:47, Jean-Claude REPETTO wrote:
>
> cs2cs +proj=tmerc +lat_0=0 +lon_0=3 +k=0.9996 +x_0=500000 +y_0=0
> +ellps=intl +towgs84=-89.5,-93.8,-123.1,0.000,0.000,-0.156,1.2 +units=m
> +no_defs +to +proj=tmerc +lat_0=0 +lon_0=0 +k=0.9996 +x_0=500000 +y_0=0
> +ellps=intl +towgs84=-89.5,-93.8,-123.1,0.000,0.000,-0.156,1.2 +units=m
> +no_defs
> 391401.4 6381482.7
> 570890.80	6380649.35 0.00
>
> The result is not too far from yours.

As Hamid also asked for "...using proj", here an alternative option:

echo "391401.4 6381482.7" |
invproj +proj=tmerc +lat_0=0 +lon_0=3 +k=0.9996 \
   +x_0=500000 +y_0=0  +ellps=intl |
proj +proj=tmerc +lat_0=0 +lon_0=0 +k=0.9996 \
   +x_0=500000 +y_0=0 +ellps=intl
570890.80       6380649.35

Hermann