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. _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione [email protected] http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it.
