[Gvsig_italian] Sovrapposizione layer .shp - Scarto di 500 km
e-paul
paolo a caporro.net
Dom 19 Dic 2010 23:33:58 PST
A questo punto... non ci sto più capendo nulla!
Scusate ma, come ho premesso, sono autodidatta e alle prime armi.
Dunque, ricapitolando;
- come faccio a sapere quale è il CRS associato al mio quadro di unione
della CTR? Il file qu_10u33.prj contiene quanto segue:
GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
quindi siamo in WGS 1984, ma in che zona?
Il fatto che nel nome del file ci sia un "33" mi fa pensare che siamo
nella zona 33, ma questa informazione è contenuta nel file?
- Il file ISTAT dei Comuni effettivamente fa riferimento alla zona 32N.
Il file "com2010_s.prj" contiene quanto segue:
PROJCS["ED_1950_UTM_Zone_32N",GEOGCS["GCS_European_1950",DATUM["D_European_1950",SPHEROID["International_1924",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",9.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]
- riproiezione al volo vuol dire che i layer vengono correttamente
rappresentati indipendentemente dal CRS adottato per la vista che sto
utilizzando? Devo essere io comunque a specificare la trasformazione
necessaria o il software individua i CRS di ciascun layer autonomamente?
Il 17/12/2010 10.22, Antonio Falciano ha scritto:
> Il 17/12/2010 9.20, Pietro d'Orio ha scritto:
>> Io opererei in questo modo:
>>
>> 1. Apri una vista con codice EPSG 23030 e ci carichi i dati originali Istat
>> 2. Apri il processo di riproiezioni a partire da 23030 al codice EPSG
>> della tua CTR
>> 3. Chiudi la vista, ne fai una nuova con il codice EPSG della tua CTR e
>> carichi i nuovi confini riproiettati.
Alla fine ho aeguito un po' una via di mezzo tra le varie proposte che
ho letto qui nella lista:
- ho creato una nuova vista con CRS 32633
- ho caricato il quadro di unione CTR (non dovrebbero essere necessarie
riproiezioni, se ho capito bene, quindi non ne applico alcuna)
- ho caricato lo shapefile dei comuni ISTAT riproiettando al volo (cioè
quando carico il file con trasformazione da WGS_84 a EPSG 32633
I layer sono ancora distanti 500 km :-(
Se invece riproietto (tirando a indovinare) il layer dei comuni ISTAT da
WGS84 a EPSG 32634 si sovrappongono correttamente.
Qualcuno di voi potrebbe essere disposto a fare la stessa operazione di
sovrapposizione per verificare i risultati? in caso posso inviare i
files necessari in email privata.
Grazie ancora a tutti e scusate per il protrarsi di questa discussione!
Paolo Caporro
More information about the Gvsig_italian
mailing list