[QGIS-it-user] Integrazione Lidar DEM

[Carlo Tersigni]

Salve. Sto facendo una analisi geomorfologica per una ferrovia in Africa.
E' stato fatto un rilievo Lidar (risoluzione 15 metri) lungo una fascia di
circa 1,5 km a cavallo del tracciato, con la restituzione della nuvola di
punti in formato xyz o csv. Ho scaricato i DEM per l'area vasta, con
risoluzione 30 metri, dal sito USGS, per poter comprendere nell'analisi
anche i bacini di una certa importanza.
Per poter lavorare in maniera omogenea e poter aumentare l'accuratezza,
vorrei integrare il DEM con il rilievo Lidar. Ho pensato di procedere in
questo modo:
1) Costruire una maglia di punti distanti 30 metri, larga circa 30 km e
lunga circa 120 km (in totale più di 10 milioni di punti).
2) Campionare il DEM in corrispondenza dei punti per avere la z. Salvare il
layer in un csv, con xyz negli attributi.
3) Fare un merge tra il layer dei punti Lidar e quello dei punti
ricampionati dal DEM, ottenendo un unico layer con un addensamento di punti
nella fascia a cavallo del tracciato.
4) Ricostruire un nuovo DEM mediante il comando di interpolazione; su tale
DEM svolgere l'analisi morfologica con QGIS o con SAGA.
C'è però qualcosa che non va in questo piano di lavoro. Fino al
campionamento dei punti va tutto bene, ma in seguito il processo è
lentissimo e difficoltoso. Creare due nuove colonne di attributi in cui
riversare le coordinate x e y nel file di campionamento, in modo da avere
un i tre valori x, y e z nella tabella attributi, copiando le coordinate $x
e $y è praticamente impossibile; comunque risulta difficile copiare anche
l'attributo y in una apposita colonna (dato che la posizione ed il nome
sono errati).
Cosa c'è che non va in questo piano di lavoro? Qualcuno sa dirmi dove è
sbagliato o sa suggerirmi un percorso più efficiente, per poter integrare
il DEM con i dati del Lidar? Ho una buona esperienza con il calcolatore dei
campi per i vettori, un po' meno con il calcolatore per i DEM.
Tutti i file sono nello stesso SR: UTM 32N.
Grazie per eventuali correzioni o suggerimenti.
​Carlo Tersigni​

-- 
*Carlo TERSIGNI*
*Via Padre Semeria, 63    00154 - ROMA*

*tel.: 333-4653589e-mail: carter.012 a gmail.com <carter.012 a gmail.com>*
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/qgis-it-user/attachments/20180402/01e47a9a/attachment.html>