QGIS2ThreeJs

QGis2ThreeJS e' un comodo plugin di QGis per creare dei modelli digitali dal terreno che possono essere visualizzati in modo interattivo su Web mediante la libreria ThreeJS

Per la prova ho preso i dati DEM ripresi da Aster (http://earthexplorer.usgs.gov/)

Non avendo delle immagini georiferite da spalmare sul DEM ho ripreso una Geotiff del rilievo a falsi colori da questo indirizzo http://geodati.fmach.it/gfoss_geodata/SRTM-Italy/

Una volta caricati i dati come raster su QGis si apre il plugin QGis2ThreeJS  scegliendo il tema del DEM ed il tema dell'immagine

deve essere selezionata anche l'esagerazione verticale (a tentativi ho trovato il valore di 0.02)

Una volta scelto il nome del file del progetto html si clicca Run e si apre in automatico il browser e vengono caricati i dati generati dal plugin

Time Machine su Linux

Mi e' capitato di dover trasferire un backup Time Machine da una macchina Mac ad un disco di una macchina Linux piu' che altro per usare ShotWell al posto di foto ... e sono nati i problemi

Su Ubuntu collegando il disco esterno formattato da Mac, questo viene correttamente ed automaticamente montato (almeno su Ubuntu) e compare nella directory /media. Il problema e' che l'organizzazione dei file non e' immediatamente chiara

le directory principali sono
/Backups.backupdb
/.HFS+ Private Directory Data

dentro /Backups.backupdb si trovano directory con il nome coincidente con la data nella quale sono stati fatti i vari backup, esiste anche una directory Latest che indica dove e' contenuto l'ultimo backup

Se si entra in queste directory si trova la struttura dell'albero come normale su un sistema unix
Per le fotografie il path sara' simile a questo
cd /media/luca/Backup1T/Backups.backupdb/MacBook\ Pro\ di\ Luca/Latest/Macintosh\ HD/Users/lucainnocenti/

qui pero' non si trovano pero' i file. Facendo un lista dei file si osserva che (ls -la)

-r--r--r-- 2067021 root 2075999     0 gen  3 17:06 Music
-r--r--r-- 2067033 root 2076001     0 gen  3 17:06 Pictures
-r--r--r-- 1513259 root 2075617     0 gen  3 17:06 .pip
-r--r--r-- 1332692 root 2076002     0 gen  3 17:06 Public

non esiste una directory Pictures.
Per trovare i file il trucco e' prendere nota del numero subito prima di root (in questo caso 2067033 che per essere precisi e' un inode) ed andare nella directory

/media/luca/Backup1T/.HFS+ Private Directory Data/dir_2067033

qui si trovano dei file di immagini e delle sottodirectory. Per copiare i file si puo' fare direttamente una copia (cp) mentre le directory sono di nuovo degli inode-r--r--r-- 1069123 root 2031303       0 gen  3 17:06 2000

-r--r--r-- 1069205 root 2031304       0 gen  3 17:06 2001
-r--r--r-- 1069248 root 2031305       0 gen  3 17:06 2002
-r--r--r-- 1069326 root 2031306       0 gen  3 17:06 2003
-r--r--r-- 1069891 root 2031307       0 gen  3 17:06 2004

per ottenere i file delle foto dell'anno 2001 dovremo cambiare directory annotando l'inode 1069205
e copiare con il comando

cp -R /media/luca/Backup1T/.HFS+ Private Directory Data/dir_1069205/* /home/luca/copia/

poco lineare ma funziona

Configurare traffico client su Tomato VPN

Questo post fa seguito al precedente dove sono indicate le istruzioni per configurare OpenVPN su Tomato

Con la configurazione standard i client si connettono alla VPN ma non si vedono reciprocamente.

Pingando i vari client direttamente sul router questi rispondono ma pingando un client verso l'altro questo non risponde

Per instradare il traffico tra i vari client si deve andare nelle impostazioni avanzate e settare Manage Client Specific Option-> Allow Client<->Client

Fotogrammetria da terra per geologia (2)

Dopo questo post ho cercato di rendere la procedura di estrazione dei piani automatizzata ed ho trovato dentro CloudCompare un plugin denominato RANSAC Shape Detection che fa esattamente cio' (il plugin e' presente solo dalla versione 2.6.1 e si attiva solo se e' selezionato un tema di nuvola di punti)

Il programma, basato su questo articolo, estrae tutte le feature possibili dall'immagine. Nell'albero a sinistra le feature di maggior peso (per esempio i piani maggiormente rappresentativi come quello di cava) si trovano in alto per scendere al maggiore dettaglio)

Piano di faglia

Piano di faglia e fronte di cava

Come si vede dalle immagini  sono stati correttamente selezionati il piano di cava e lo specchio di faglia e il piano di stratificazione

Stratificazione

Stromboli 3D

Applicazione di quanto visto in questo post a dati ripresi dall'alto

Ho ripreso delle foto che mi sono state passate circa 7 anni fa da un amico riguardanti il vulcano Stromboli scattate da un elicottero con una panoramica sulla Sciara del Fuoco

Il modello di Google Earth con il modello digitale del terreno e una immagine ottica ad alta risoluzione

Le fotografie sono state elaborate con Agisoft Photoscan. E' molto evidente il gradino morfologico sulla sinistra della Sciara ed alcuni strutture delle strutture piu' piccole

Da Kinect a Stampante 3D per paleontologia

Grazie ad un amico sono riuscito ad utilizzare una stampante a 3D.
Lo scopo era quello di, partendo da un fossile, effettuare una scansione 3D con Kinect e poi ottenere nuovamente un modello fisico mediante la stampante 3D

Kinect non e' molto adatto a questo scopo perche' ha un precisione sull'asse verticale di 1 mm e non per modelli piccoli si perde molta risoluzione

Inoltre, dato che la plastica per la stampa 3D e' costosa, mi e' stato concessa una prova su un modello molto ridotto ed alla risoluzione peggiore
La stampa e' durata una 40 (video in tempo estremamente accelerato)

Fotogrammetria da terra per geologia

Visti i due precedenti post (1 e 2) ho voluto provare a mettere alla prova il sistema su un problema geologico

Vicino a casa esiste questa vecchia cava di Macigno (o Pietra Serena)

il problema che mi sono posto e' stato quello di misurare l'orientazione dei due piani che emergono dal fronte di cava soltanto con l'utilizzo delle fotografie

Nella foto sottostante sono individuati i due angoli diedri da misurare (probabilmente si tratta di faglie)

L'area si trova all'interno di una proprieta' privata e comunque il fronte di cava non e' raggiungibile per una misura diretta

Mediante l'utilizzo di un comune telefono cellulare (Moto G) ho scattato una decina di foto dell'affioramento da una distanza calcolato mediante Carta Tecnica Regionale di circa 100 m

Sempre dalla stessa fonte si puo' calcolare (un po' grossolanamente ma la precisione non e' fondamentale) che la cava e' orientata lungo la direzione 60-240° di azimuth

Mediante VisualFSM sono riuscito a calcolare in modo non supervisionato la nuvola dei punti 

e tramite una opportuna rotazione si puo' stimare che l'angolo tra la superficie di faglia e il piano di cava e' di circa 30°. A questo punto si puo' stimare che l'orientazione del piano incognito e' di circa 30-210°

La stessa elaborazione puo' essere fatta , con risultato graficamente piu' gradevoli, anche con AGI Soft Photoscan