[Gvsig_italian] Proiezione del DEM della Regione Autonoma della Sardegna

Mar 25 Maggio 2010 11:10:31 PDT

Ciao
queste sono le informazioni supplementari che mi avete chiesto, spero vi possano essere utili per aiutarmi a risolvere il problema

1) in che CRS e' 
definito il DEM? ED50 o ROMA40 o altro?
   - Come già detto il CRS del DEM è EPSG 4022, "unknow datum based upon the International 1924 ellipsoid
    proiettato: no
    Semiasse maggiore: 6378388.0
    Inverse flattening: 297.0
    Meridiano d'origine: Greenwich
    Longitudine del meridiano: 0.0

2) perche' non usi EPSG:3003 
per la vista?
    - Ho necessità di lavorare su CRS Gauss-Boaga Fuso Ovest

3) 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.
   - Perdonate l'ignoranza, ma non so come ricavare il gvsig.log, mi potreste indicare come fare?

3)di che entita' e' lo sfasamento tra DEM e 
CTR?
    - Circa 74m circa verso N355

4)> 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?
    - Apparentemente si

Grazie, Beppe

________________________________
Da: Antonio Falciano <afalciano a yahoo.it>
A: gvsig_italian a lists.gvsig.org
Inviato: Mar 25 maggio 2010, 17:56:45
Oggetto: Re: [Gvsig_italian] Proiezione del DEM della Regione Autonoma della Sardegna

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
_______________________________________________
Gvsig_italian mailing list
Gvsig_italian a lists.gvsig.org
https://gvsig.org/lists/mailman/listinfo/gvsig_italian

-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/gvsig-italian/attachments/20100525/280b4c0d/attachment.html>