<div dir="ltr"><div><div><div><div><div><div><div>gentilissimo,<br></div>non voglio entrare troppo nel merito del tuo lavoro, però un consiglio che ti voglio dare è quello di lavorare in ambiente raster anziché vettoriale (come mi sembra di aver capito che tu stia lavorando).<br></div><br></div>converti la tua nuvola di punti in una griglia di risoluzione 15X15 metri e poi fai un resample del dem, dividendo il pixel in 4 (hai questa fortuna che le due risoluzioni siano una multipla dell'altra).<br></div>e vedrai che qualunque analisi di tipo raster sarà molto più performante (soprattutto con grass).<br><br></div>facci sapere.<br><br></div>un saluto,<br></div>francesco<br><div><div><div><div><br><br><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">Il giorno 2 aprile 2018 17:02, Carlo Tersigni <span dir="ltr"><<a href="mailto:carter.012@gmail.com" target="_blank">carter.012@gmail.com</a>></span> ha scritto:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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.<span class="m_7365132470317388722HOEnZb"><font color="#888888"><br clear="all"></font></span></div><span class="m_7365132470317388722HOEnZb"><font color="#888888"><div style="font-family:arial,helvetica,sans-serif;font-size:small" class="gmail_default">​Carlo Tersigni​</div><br>-- <br><div class="m_7365132470317388722m_-2907647920433773039gmail_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>
</font></span></div>
<br>______________________________<wbr>_________________<br>
QGIS-it-user mailing list<br>
<a href="mailto:QGIS-it-user@lists.osgeo.org" target="_blank">QGIS-it-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/qgis-it-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/qgis-it-user</a><br>
<br></blockquote></div><br></div></div>