viernes, 24 de febrero de 2017

Processing taudem:peukerdouglas en GNU/Linux Debian

La instalación y ejecución de algoritmos TauDEM en Windows fue considerada en un post anterior. En Linux, la instalación puede ser automatizada utilizando el script obtenido de esta página. A pesar de que el nombre del script sugiera que sólo puede instalarse en Ubuntu, si miran dentro del script verán que fue desarrollado inicialmente para Debian y es de validez universal para cualquier sistema derivado de Debian.



Después de instalado, averiguamos cuales son los parámetros a introducir en el código para el algoritmo taudem:peukerdouglas:


Los parámetros númericos asumidos a continuación son los que sugiere la Toolbox para este comando. Por tanto, se desarrolló el código de prueba siguiente:

1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
import processing

layer = iface.activeLayer()

path = processing.runalg('taudem:peukerdouglas', 
                         layer,      #ELEVATION_GRID <ParameterRaster>
                         0.400000,   #CENTER_WEIGHT <ParameterNumber>
                         0.100000,   #SIDE_WEIGHT <ParameterNumber>
                         0.050000,   #DIAGONAL_WEIGHT <ParameterNumber>
                         None)       #STREAM_SOURCE_GRID <OutputRaster>

STREAMSOURCEGRID = QgsRasterLayer(path['STREAM_SOURCE_GRID'],
                                        'STREAMSOURCEGRID')

QgsMapLayerRegistry.instance().addMapLayer(STREAMSOURCEGRID)

Cuando se ejecuta en la Python Console de QGIS con este ráster DEM:


se obtiene lo siguiente:


Si se hace zoom, con 50 % de transparencia, a una zona montañosa del ráster se puede observar su patrón de drenaje.


Mas detalles sobre la ejecución del TauDEM algoritmo pueden ser obtenidos de aquí.

No hay comentarios: