[Gvsig_italian] Proiezione del DEM della Regione Autonoma della Sardegna

[Antonio Falciano]

Il 25/05/2010 15.34, VenerioG ha scritto:
>
> Ciao a tutti!
> Devo eseguire un lavoro in un comune della Sardegna utilizzando gvSIG OADE
> 2010 Beta2. Dopo aver caricato nella vista i CTR ho cercato di caricare
> un'estrazione dal DEM della Regione Sardegna contenente l'area di interesse
> e, visto che la proiezione del raster non coincide con quello della vista e
> del CTR (Roma 1940 Gauss Boaga Ovest), compare una finestra di opzioni che
> mi chiede se voglio proiettare il raster. Scelgo di proiettarlo nella
> proiezione della vista e clicco su "Accetta" nella finestra "Aggiungi
> layer". A questo punto compare una nuova finestra ("livello da
> riproiettare") nel quale compare il nome del file origine (w001001.adf), la
> proiezione di origine (4022, "unknow datum based upon the International 1924
> ellipsoid"), la proiezione di destinazione (102094, Roma 1940 Gauss Boaga
> ovest, ovvero la proiezione della vista) e seleziono su "livello" la voce

1) in che CRS e' definito il DEM? ED50 o ROMA40 o altro?
2) perche' non usi EPSG:3003 per la vista?

> "Creazione di file", clicco su "accetta" e nella finestra "seleziona file"
> gli attribuisco un nome (DEM_Ulas) assegnandogli come tipo di file GeoTIFF,
> clicco su "salva" e a questo punto il programma va in crash lanciando il
> seguente lunghissimo messaggio di errore:
>
> http://osgeo-org.1803224.n2.nabble.com/file/n5098475/Messaggio_di_errore.jpg

non va bene... ci puoi inviare il file gvsig.log?
il crash forse e' dovuto al fatto che il CRS della vista non e' definito
con un codice EPSG, ma con uno ESRI... se mi mandi il gvsig.log si
dovrebbe isolare il problema.

> Lo stesso problema si verifica anche se scelgo di aprire in memoria il file
> anzichè creare il file GeoTiff.
> Invece il programma non crasha se decido di non riproiettare il layer; però
> dal confronto tra il DEM e i CTR è rilevabile un certo sfasamento tra i due.

di che entita' e' lo sfasamento tra DEM e CTR?

> Informazioni del dataset DEM utilizzato
> File: w001001.adf
> Dimensione: 3,31 MB (3.467.568 bytes)
> Larghezza X Altezza: 1427.0 X 562.0
> Formato: adf
> Georeferenziando: si
> Numero di bande: 1
> Tipo di Dato: Float
>
> Coordinate geografiche
> Alto sinistra: 1535603.8788478847, 4411157.680373412
> Inferiore destra: 1549873.8788478847, 4405537.680373412
> Alto destra: 1549873.8788478847, 4411157.680373412
> Inferiore sinistro: 1535603.8788478847, 4405537.680373412

Ma queste coordinate sono gia' in Gauss-Boaga Fuso Ovest, o sbaglio?

> Dimensione pixel X: 10.0
> Dimensione pixel Y: -10.0
> Rotazione in X: 0.0
> Rotazione in Y: 0.0
>
> Origine
> Band 1: Type=Float, ColorInterp=Undefined
>
> Proiezione
>    PROJCS["unnamed",     GEOGCS["Unknown datum based upon the International
> 1924 ellipsoid",
>        DATUM["Not_specified_based_on_International_1924_ellipsoid",
>          SPHEROID["International 1924",6378388,297,
>            AUTHORITY["EPSG","7022"]],         AUTHORITY["EPSG","6022"]],
>        PRIMEM["Greenwich",0,         AUTHORITY["EPSG","8901"]],
>        UNIT["degree",0.01745329251994328,         AUTHORITY["EPSG","9122"]],
>        AUTHORITY["EPSG","4022"]],     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["METERS",1]],

Il CRS/proiezione del layer deve essere definito con cura, altrimenti
gvSIG (o meglio la libreria GDAL che gestisce le riproiezioni) non sara'
in grado di individuare una trasformazione compatibile con i CRS di
origine e di destinazione.
Nella documentazione si afferma che sono utilizzate solo le
trasformazioni EPSG (no manuale, no composte, no griglie) e quindi e'
opportuno che tutti i CRS siano riferiti a tale db.
A livello di precisione, naturalmente solo con le conversioni si puo'
ottenere qualcosa di assolutamente accurato.

ciao
Antonio

-- 
Antonio Falciano
http://www.linkedin.com/in/antoniofalciano