APPROFONDIMENTO Creazione di un DTM partendo dalle curve di

APPROFONDIMENTO
Creazione di un DTM partendo dalle curve di livello
Di seguito viene descritto come elaborare un Digital Terrain Model (DTM) partendo da curve di livello vettoriali (per i casi in cui non si avesse a
disposizione il DTM).
Azione
1
Creare un mapset
nuovo
Comando
Se si deve ancora avviare grass si usa
il pulsante crea mapset
Se si è già in un altro mapset
Comando: g.mapset
Si crea un mapset nuovo così si può
impostare la risoluzone al valore che
conviene, senza il rischio di modificare
altre mappe raster già presenti
Opzioni
Scheda Impostazioni name of mapset where to switch scrivo il nome del mapset da creare
Scheda Create si seleziona il quadratino "create
mapset if it doesn't exist”
Permette lavorando dal Mapset A di
vedere anche i contenuti del mapset B
Menu:
2
Collegare agli altri
mapsets
Impostazioni GRASS
ambiente di lavoro accesso
al mapset
Comando: g.mapsets
3
3bis
Visualizzare le
curve di livello
Esempio o note
NON è possibile lavorando da un
mapset A modificare le mappe di un
mapset B.
Creare il mapset “DTM_nuovo”
Chiedere di vedere dal mapset corrente
(“DTM_nuovo”) il mapset “DTM”
NOTA: Da ogni mapset si vedono
sempre di default le mappe del mapset
PERMANET, ma non quelle degli altri
mapset, a meno di specificarlo con
questo comando
Menu:
Scegliere dal menu a tendina la mappa
vettoriale con le curve di livello fatta nell’
esercitazione precedente
Comando: d.vect
In caso si vogliano usare curve di livello in formato vettoriale, occorre prima importarle
Menu: File importa mappa vettoriale formati comuni di importazione
Comando: v.in.ogr
Dispense per il corso di Geomorfologia per il Progetto del Paesaggio, 2013-2014
a cura di Prof. Geol. G.Brancucci, Dott.ssa Arch. P.Salmona, Dott.ssa V.Marin
Geomorfolab http://geomorfolab.arch.unige.it - Scuola Politecnica, Università degli Studi di Genova
4
5
Azione
Visualizzare il
contorno del
bacino del
Bisagno
Impostare la
regione sulla
estensione e
risoluzione
desiderata per il
DTM
Comando
Comando: d.vect
Menu: Impostazione Regione
Scheda Esistente Imposta Region
imposta la regione da questa mappa
vettoriale
Comando: g.region
Più alta è la risoluzione (es. 1,
oppure 5, ecc.) più preciso (e più
pesante) viene il DTM
Comando: v.to.rast
6
Esempio o note
Scegliere dal menu a tendina
bisagno
Menu:
Menu: file conversione tipo di
mappa
vettore a raster
Trasformare le
curve di livello
vettoriali in raster
Opzioni
Attenzione!
Non è fondamentale che le curve
siano continue, ma è necessario
che ogni curva abbia associata
nella tabella la sua quota in
formato numerico!!
Scheda Risoluzione risoluzione 2D della griglia inserire il
numero della risoluzione desiderata
Scheda Richiesto Nome della mappa vettoriale in input:
mappa da trasformare in raster
Nome del raster di output: nome che si
vuole dare alla nuova mappa
Source of raster value attr
Scheda Selection line
Scheda Attributi name of column for 'attr' parameter colonna che contiene le quote
Scheda Esistente imposta la regione da questa mappa
vettoriale bisagno
Scheda Risoluzione risoluzione 2D della griglia 5
Si mette 5 in questo caso per fare poi un
confronto con il DTM "originale"
Scheda Richiesto Input: curve_livello
Output: curve_rast
Source of raster value attr
Scheda Selection line
Scheda Attributi name of column for
'attr' parameter level
NOTA: level è il nome della colonna creata in automatico da GRASS quando si
estraggono le curve di livello dal DTM (come abbiamo fatto noi nella scorsa esercitazione).
Se i file delle curve di livello sono stati importati, verificare prima il nome della colonna in
cui sono riportate le quote e indicare quella.
Dispense per il corso di Geomorfologia per il Progetto del Paesaggio, 2013-2014
a cura di Prof. Geol. G.Brancucci, Dott.ssa Arch. P.Salmona, Dott.ssa V.Marin
Geomorfolab http://geomorfolab.arch.unige.it - Scuola Politecnica, Università degli Studi di Genova
Azione
Comando
Menu: raster calcolatore per
raster mapalgebra
7
Si moltiplica il
valore delle curve
di livello per 1000
Comando: r.mapcalc
Opzioni
Nome per il nuovo raster da creare:
scelgo un nome
Expression Mappa raster delle curve
di livello * 1000
Esempio o note
Nome per il nuovo mappa da creare:
curve_rast_1000
Espressione Curve_rast@DTM_nuovo * 1000
Questo passaggio serve per ottenere valori decimali per le quote del DTM, altrimenti da un pixel all’altro dovrebbe passare almeno
1m
8
Si crea il DTM
Menu: Raster Interpolazione
superfici raster contorni
Comando: r.surf.idw
Menu: raster calcolatore per
raster mapalgebra
9
Si divide il valore
del DTM per 1000
Comando: r.mapcalc
Scheda Richiesto Input: selezionare la mappa con le
curve raster già moltiplicata per 1000
Output: nome del DTM da creare
Nome per il nuovo raster da creare:
scelgo un nome
Expression Mappa raster delle curve
di livello / 1000
Scheda Richiesto Input: curve_rast_1000
Output: DTM_curve_1000
Nome per il nuovo raster da creare:
DTM_curve
Expression Curve_livello_rast@DTM_nuovo / 1000
Questo passaggio serve a riportare le quote del DTM ai valori approssimati al reale
10
Gestire i colori del
DTM
Menu: raster gestisci colori
tabella dei colori
Comando: r.colors
Scheda Richiesto Input: selezionare la mappa da modificare
Colors Tipo tavola colore scelgo
una tabella
Cliccare equalizzazione
dell'istogramma
Scheda Richiesto Input: DTM_curve
Colors Tipo tavola colori elevation
Cliccare equalizzazione
dell'istogramma
Dispense per il corso di Geomorfologia per il Progetto del Paesaggio, 2013-2014
a cura di Prof. Geol. G.Brancucci, Dott.ssa Arch. P.Salmona, Dott.ssa V.Marin
Geomorfolab http://geomorfolab.arch.unige.it - Scuola Politecnica, Università degli Studi di Genova