TINITALY DEM: COME SCARICARE LA REGIONE SICILIA E USARLO IN QGIS

Una buona notizia (grazie a Gimmi Gis) la sezione di Pisa dell’Istituto Nazionale di Geofisica e Vulcanologia, (alias INGV, grazie mille!!!), ha messo a disposizione qui il Modello Digitale di Elevazione di tutta l’Italia, alla risoluzione di 10 metri pixel, con il nome di Tinitaly DEM.

Tinitaly DEM è disponibile in formato Geotiff nel sistema di proiezione UTM WGS 84 zona 32 (EPSG: 32632), con licenza aperta CC BY 4.0 (leggi qui) che permette di condividere e di modificare il dato, purché venga indicata una menzione di paternità adeguata, e si indichi se sono state fatte delle modifiche.

Dal sito è possibile cliccare su una cella del reticolo e scaricare il DEM relativo, ma la procedura sarebbe parecchio lunga se volessimo scaricare, per esempio, i tasselli di una intera regione (per esempio la Sicilia ha 27 tasselli). Per velocizzare il processo di download ho realizzato uno script che permette di scaricare tutti i tasselli che ricadono in una determinata area.

1. Tarquini S., Isola I., Favalli M., Mazzarini F., Bisson M., Pareschi M. T., Boschi E. (2007). TINITALY/01: a new Triangular Irregular Network of Italy, Annals of Geophysics 50, 407 – 425.
2. Tarquini S., Vinci S., Favalli M., Doumaz F., Fornaciai A., Nannipieri L., (2012). Release of a 10-m-resolution DEM for the Italian territory: Comparison with global-coverage DEMs and anaglyph-mode exploration via the web, Computers & Geosciences 38, 168-170. doi: doi:10.1016/j.cageo.2011.04.018

A seguire uno script bash per scaricare una porzione del Tinitaly DEM, in particolare per scaricare l’intera Sicilia.

Script bash per scaricare i tasselli della Sicilia

#!/bin/bash
# scarica tutti i tasselli della Sicilia TINITALY DEM 10m/px CC BY 4.0
### requisiti ###
# wget https://www.gnu.org/software/wget/
# 7z https://www.7-zip.org/
# ogr2ogr https://gdal.org/programs/ogr2ogr.html
# GDAL https://gdal.org/
# aborruso https://twitter.com/aborruso
### requisiti ###
folder="$(cd "$(dirname "${BASH_SOURCE[0]}")" && pwd)"
mkdir -p ./rawdata
mkdir -p ./processing
# cancella file
rm -r "$folder"/rawdata/*
rm -r "$folder"/processing/*
rm -r "$folder"/out*.txt
# crea lista nomi celle per la Sicilia cod_reg= 19
ogr2ogr -f CSV -dialect sqlite -sql "SELECT ID FROM Reg01012020_g_WGS84 r, './Quadro_unione/_qu_retiled_32.shp'._qu_retiled_32 u where cod_reg= 19 and st_intersects (u.geometry, r.geometry) = 1" output.csv ./Reg01012020_g/Reg01012020_g_WGS84.shp
# togli prima riga
tail -n +2 output.csv >out1.txt
# elimina _10 finale
sed -e 's/_10//' out1.txt >lista.txt
# fai un ciclo per ogni tassello della Sicilia
for i in $(cat lista.txt)
do
# vai nella cartella rawdata
cd "$folder"/rawdata/
# scarica il file zip
wget http://tinitaly.pi.ingv.it/data/"$i"_s10/"$i"_s10.zip
# vai nella cartella processing
cd "$folder"/processing/
# unzippa il file
7z x "$folder"/rawdata/"$i"_s10.zip
# sposto file
mv "$folder"/processing/"$i"_s10/"$i"_s10.tif "$folder"/processing/
# cancella le cartelle inutili
rmdir "$folder"/processing/"$i"_s10
done
# vai nella cartella processing e crea lista dei file tif
cd "$folder"/processing/
ls *.tif >list.txt
# mergia tutti i file e tasselali
gdal_merge.py -o merged.tif -co "TILED=YES" -co "BLOCKXSIZE=512" -co "BLOCKYSIZE=512" -co "COMPRESS=DEFLATE" -co "ZLEVEL=9" -co "BIGTIFF=YES" -init -9999 -a_nodata -9999 -n -9999 -ot Float32 –optfile list.txt
# crea piramidi
gdaladdo -r nearest –config COMPRESS_OVERVIEW DEFLATE –config GDAL_TIFF_OVR_BLOCKSIZE 512 "$folder"/processing/merged.tif 2 4 8 16

Per la sola Sicilia :

NOTE FINALI: Lo script bash scarica tutti i tasselli, ma la magia è di GDAL che permette di mergiare, tassellare e piramidare rendendo facile la navigazione (in QGIS) di un DEM pesante circa 2GB.

Per poter usufruire al meglio dello script occorre scaricare lo shapefile ISTAT regioni (Reg01012020_g_WGS84.shp) e il quadro di unione della Tinitaly DEM (vedi nei riferimenti, gentilmente concesso dal dott. Tarquini dopo richiesta di Andrea Borruso tramite email, presto sarà messo nel sito).

Per chi usasse la OSGeo4W 64b shell ho realizzato lo stesso script, lo trovate qui (leggere i requisiti nello script).

Riferimenti

RINGRAZIAMENTI

  • I ragazzi di INGV per tutto il lavoro svolto;
  • Gimmi GIS per la divulgazione sui social e il suo ottimo articolo;
  • Andrea Borruso, perché gli voglio bene e mi aiuta sempre.

video demo


Se il blog post vi è piaciuto cliccate su ‘Mi piace’, grazie!!!
if you liked the blog post click on ‘Like’, thank you !!!

SE IL POST VI È STATO UTILE CONTRIBUITE A MANTENERLO AGGIORNATO PAYPAL


2 pensieri su “TINITALY DEM: COME SCARICARE LA REGIONE SICILIA E USARLO IN QGIS

Rispondi

Inserisci i tuoi dati qui sotto o clicca su un'icona per effettuare l'accesso:

Logo di WordPress.com

Stai commentando usando il tuo account WordPress.com. Chiudi sessione /  Modifica )

Google photo

Stai commentando usando il tuo account Google. Chiudi sessione /  Modifica )

Foto Twitter

Stai commentando usando il tuo account Twitter. Chiudi sessione /  Modifica )

Foto di Facebook

Stai commentando usando il tuo account Facebook. Chiudi sessione /  Modifica )

Connessione a %s...

Questo sito utilizza Akismet per ridurre lo spam. Scopri come vengono elaborati i dati derivati dai commenti.