[Talk-it] [OT] convertire shape da CTR a OSM/Google

emmexx emmexx a tiscalinet.it
Mar 11 Ott 2011 17:21:46 BST


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



Maggiori informazioni sulla lista Talk-it