[Gfoss] Problema conversione coordinate

[Antonio Falciano]

chri ha scritto:
> Ho provato a salvare lo shapefile con le coordinate Gauss Boaga e
> convertirlo con:
> 
> ogr2ogr -f "ESRI Shapefile" -s_srs "roma40.srs" -t_srs "EPSG:4326"
> wgs84_shp IP000GM.shp
> 
> dove roma40.srs è
> 
> +proj=tmerc +lat_0=0 +lon_0=15 +k=0.999600 +x_0=2520000 +y_0=0
> +ellps=intl +units=m
> +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
> 
> e tutto va alla grande. Non capisco perché il seguente frammento di codice:
> 
> import osgeo.osr
> import osgeo.ogr
> gb = osgeo.osr.SpatialReference()
> gb.ImportFromEPSG(3004)
> wgs84 = osgeo.osr.SpatialReference()
> wgs84.ImportFromEPSG(4326)
> tr = osgeo.osr.CoordinateTransformation(gb, wgs84)
> tr.TransformPoint(2365622, 5118141, 0)
> 
> non è equivalente. Qualche suggerimento?

Non è equivalente perchè EPSG:3004 è pari a

> +proj=tmerc +lat_0=0 +lon_0=15 +k=0.999600 +x_0=2520000 +y_0=0
> +ellps=intl +units=m
> +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68

meno +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68

cioè

> +proj=tmerc +lat_0=0 +lon_0=15 +k=0.999600 +x_0=2520000 +y_0=0
> +ellps=intl +units=m

Ciao
Antonio

PS: Scusa la risposta tardiva.