Grazie, esaurientissimo!
chalten
chalten
Nelle istruzioni per la creazione dei file tif è riportato che bisognerebbe inserire Proiezione e Datum relativi alla mappa su cui saranno attivati i DEM:
"For raster maps, you will have to use the correct projection and datum. Replace the projection string (
+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
) with that which is used for your map. If there is a mismatch between the map and the DEM data QLandkarte GT will tell you and display the map's string, too."
Dovresti cambiare quindi la riga di comando inserendo questi parametri:
+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0
A mio avviso però è inutile, il programma dovrebbe convertire al volo Proiezione e Datum dei DEM e sovrapporli esattamente sulla mappa base, anche se nella conversione probabilmente viene introdotto qualche metro di errore nelle coordinate. Ho fatto un paio di prove ed infatti rilevo qualche metro di quota di errore.
Le differenze che rilevi di poche decine di metri sono fisiologiche considerando che si tratta di DEM con griglia di 90 mt, mentre le IGM sono molto più dettagliate.
Comunque, se stai usando le IGM online e converti i file .hgt con questa riga di comando non avrai il messaggio di errore perchè i file .hgt originali sono già con proj=latlong
gdalwarp -r cubic c:\percorsodem\*.hgt c:\percorsodem\dem.tif
francesco