Creazione del DTM con GRASS 6.0.2
Da Wikipedia, l'enciclopedia libera.
A partire dalle curve di livello è possibile creare un modello digitale del terreno mediante il comando:
v.surf.rst
dopo aver settato la regione di lavoro ad una risoluzione di 5m.
Si riportano i valori dei parametri settati per il territorio (estensione di circa 30 kmq) con zona pianeggiante pari a circa il 40%.
ten= 40 (valori bassi sono indicati per regioni pianeggianti con pochi punti sparsi) smooth= 0,3 dmin= 2 dmax= 100 segmax= 80 npm= 300
Una volta ottenuto il DTM desiderato si possono effettuare ulteriori elaborazioni quali analisi del soleggiamento con r.los, affinamento del DTM stesso aggiungendo l'elevazione degli edifici. Per fare questo si sono utilizzati i dati delle schede relative agli edifici catalogate in un file vettoriale; fra i dati memorizzati vi è anche la componente "altezza media edificio". Innanzitutto è necessario convertire il vettoriale degli edifici in raster assegnado alle celle il valore trovato nel campo "altezza_media":
v.to.rast in=edifici_vect out=edifici_rast col=altezza_media
In questo modo è stato creato un raster in cui ogni cella riporta il valore dell'altezza media degli edifici; altra operazione da fare è convertire tutti i valori NULLs del raster in valori zero.
r.null edifici_rast null=0
mediante r.mapcalc "aggiungiamo" la mappa raster degli edifici al DTM:
r.mapcalc "dtm_edifici=if(edifici_rast > 0, dtm+edifici_rast, dtm)"
Con quest'ultima operazione abbiamo ottenuto il dtm definitivo "estruso" in corrispondenza degli edifici di una quantità pari all'altezza media degli stessi.
