<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">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.<br>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:<br>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).<br>2) Campionare il DEM in corrispondenza dei punti per avere la z. Salvare il layer in un csv, con xyz negli attributi.<br>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.<br>4) Ricostruire un nuovo DEM mediante il comando di interpolazione; su tale DEM svolgere l'analisi morfologica con QGIS o con SAGA.<br>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).<br>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.<br>Tutti i file sono nello stesso SR: UTM 32N.<br>Grazie per eventuali correzioni o suggerimenti.<br clear="all"></div><div style="font-family:arial,helvetica,sans-serif;font-size:small" class="gmail_default">​Carlo Tersigni​</div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><i>Carlo TERSIGNI</i><div><i><font color="#666666">Via Padre Semeria, 63    00154 - ROMA</font></i></div><div><i><font color="#666666">tel.: 333-4653589<br>e-mail: <a href="mailto:carter.012@gmail.com" target="_blank">carter.012@gmail.com</a></font></i><br></div></div></div></div></div></div>
</div>