Il 10/11/2011 05:51 PM, Federico Cozzi scrisse:
> Io ho provato il tool della regione. Ho messo in input le tue coordinate:
> 9,170837078374765 45,525406361655726
> e ho ricevuto in output
> 9,17050223531423 45,5260796970209
> (differenze alla terza / quarta cifra decimale)
>
> Non so se la conversione sia giusta, ma qualcosa รจ successo...
>
9.17047473,45.52613180
Questo e' il valore che mi hanno passato nel file kml.
> Tu prova a convertire con ogr2ogr in due passi:
> 1. da EPSG:4265 a EPSG:3003 (nessun cambio di ellissoide: da Roma1940
> geografico a Roma1940 / Monte Mario 1)
> 2. da EPSG:3003 a WGS84 con i parametri di Luca
Ho provato con:
ogr2ogr -a_srs "EPSG:4265" -t_srs "EPSG:3003" ...
seguito da:
ogr2ogr -a_srs "gaussboaga1.prj" -s_srs "latlong_wgs84.prj" -t_srs
"EPSG:4326"
segnala un problema. Si aspetta latitudine e longitudine (???).
Riepilogo:
- la prima trasformazione mi genera uno shape che ha coordinate cartesiane.
Layer name: mioshape
Geometry: Line String
Feature Count: 185
Extent: (1503803.618361, 5029013.003523) - (1520914.343248, 5041637.401461)
Layer SRS WKT:
PROJCS["Monte Mario / Italy zone 1",
GEOGCS["Monte Mario",
DATUM["Monte_Mario",
SPHEROID["International_1924",6378388,297]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",9],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",1500000],
PARAMETER["false_northing",0],
UNIT["Meter",1]]
- la seconda trasformazione (a questo punto mi aspettavo che
applicandola ad un file in coordinate cartesiane fosse la cosa corretta)
invece si aspetta coordinate geografiche.
La mia confusione aumenta. O forse e' errata la sintassi dei comandi. O
Luca mi ha mandato dei file non corretti (ma in realta' ne ho trovati di
simili in rete).
Anche questa genera punti in coord. cartesiane.
Layer name: CICLOMILANO RILIEVO
Geometry: Line String
Feature Count: 185
Extent: (1503803.618361, 5029013.003523) - (1520914.343248, 5041637.401461)
Layer SRS WKT:
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_1984",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]]
Se a questo punto do il comando:
ogr2ogr -f "KML"
ottengo una segnalazione di coordinate fuoriscala.
ciao
maxx
_______________________________________________
Talk-it mailing list
Talk-it@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-it
