How to read and decently plot a DEM with Python

I worked for long time to build  GISes (e.g. JGrass, uUdig and so on), an adventure that was supported and eventually taken in charge by Hydrologis. After accepting the reality the QGIS won the battle, I am not still satisfied with it and always looking to alternatives. One that I practice a little is to use as vehicle for mapping Jupyter and Jupyterlab. Therefore,  copying  around, I was able to produce decent maps. You can find one of the result of my search in the figure below.

Nothing really original from me: I copied from the bests I found and adapted to my desiderata. Clicking on the Figure above you access the Jupyter notebook where all the procedure is explained. Aa I wrote,  it is very much a work in progress. Therefore any suggestion is welcomed.

Updated material on Horton Machine

The Horton Machine is our set of tools for terrain analysis actually maintained by Hydrologis. They since 2000 and something took care of keeping my tools alive on various platforms: GRASS, JGrass, uDig and now on gvSIG. We also have a porting of the Horton machine that works under GEOframe but is the same codebase.

Silvia Franceschi and Andrea Antonello gave a tutorial at this year gvSIG conference and here it is the material.

You can find the same material also at this osf project. Soon they will follow some theoretical lectures.

The Horton Machine: JGrasstools resources for DEM treatment and delineation (in Italian)

Ecco dunque la documentazione e vari esempi dell'utilizzo della Horton Machine, ovvero come si analizzano i dati digitali del terreno con i Jgrasstools/STAGE

• Il tutorial (ancora in bozza - 38Mb)
• Il vecchio manuale di riferimento della Horton Machine (in Inglese) per capire qualcosa di quello che si sta facendo ed avere un po' di bibliografia.

E le varie presentazioni:

• p-JGrassTools (jgt) (ovvero: scaricare e installare STAGE, iniziare ad usare QGIS, creare una location in QGIS, usare il plugin per legende con scale colori continue
•  p-SpatialToolbox (con: un'introduzione ai jgrasstools, una spiegazione di come utilizzare i comandi dei jgt attraverso lo spatial toolbox di STAGE e di come utilizzare Mosaic12 per unire le tavole del DTM)
• p-HortonMachine (ovvero come utilizzare i comandi HortonMachine dei jgt per fare un analisi idro-geomorfologica)

•  Qui trovate alcuni DEM per esercitarvi e per confrontare, nel caso, i vostri risultati:  il rio Valpiana (100 Mb). Contiene anche una relazione tipo fatta sullo stesso bacino.

Quanto sopra è il bagaglio che serve per il mio corso di idrologia. A questo proposito

• Un'altra relazione di esempio può essere trovata qui.
• Qui si trovano degli script di R che che possono essere utili nell'elaborare alcuni risultati dei jgrasstools: 
• come graficare la distribuzione dei valori di una mappa
• come estrarre la tabella attributi di uno shapefile

Per la vostra curiosità: qui potete avere qual è il costo delle nozioni apprese (costi del software a parte) in una iniziativa commerciale.

Open Topography

The site

http://www.opentopography.org/

hosts free LIDAR data, and tools contained in OpenTopography for their treatment. Data are usually provided in form of points clouds that can be subsequently treated with their tools to obtain DEMs, and in already treated DEMs. They also provide a ciberinfrastructure for the storage and delivery of those data that can be of interest for those who wants to manage professionally a data storage service.
Woth to note is that you can use the jgrasstools for manage them.