[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