[gdal-dev] projection with EPSG transformation

[Jean-Claude Repetto]

Le 11/11/2010 12:26, Michele Beneventi a écrit :
> Hi list,
>
> I have to reproject with "gdalwarp" some tiff file from "+proj=tmerc
> +lat_0=0.0 +lon_0=9.0 +k=0.9996 +x_0=500000.0 +y_0=0.0 +ellps=WGS84
> +datum=WGS84 +units=m" (EPSG:32632) to "+proj=tmerc +lat_0=0.0
> +lon_0=9.0 +k=0.9996 +x_0=1500000.0 +y_0=0.0 +ellps=intl +units=m"
> (EPSG:3003).
> To enhance the precision (<  4m) I would like to include in the command
> line also the information about the local tranformation  (EPSG:1662).
>
> How can I include this information by the  command parameters?

Hello,

You must add a +togs84 parameter. Something like :

+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl 
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +units=m +no_defs

(These values are for EPSG:1660, not EPSG:1662).

Jean-Claude