raster_aguaLa semana pasada recibí la entrega del proyecto de Fin de Master GIS con Python de dos alumnos de la Agencia de Medio Ambiente y Agua de la Junta de Andalucía.  Han montado un proyecto muy interesante que cubre una de sus necesidades del día a día. La gestión de riesgo del incendios forestales es parte del núcleo de su actividad principal y la lectura e interpretación de datos meteorológicos es, como se puede entender, una rama fundamental para su desempeño diario. Pues bien, nos han planteado una herramienta que lee de forma automática los ficheros de datos meteorológicos, los une con los puntos de estaciones de medida y hace una análisis de interpolación con uso incluso de técnicas de Machine Learning, como al Cluster and Outliers Analysis.

logo_Agencia-Medio-Ambi

La herramienta se implementa sobre un script de Python en ArcPy. Y para usarla se puede:

 – Instalar un Addin

 – Lanzarla directamente desde la ToolBox desarrollada a tal efecto

La Toolbox y los archivos que utiliza se pueden descargar en esta url y más abajo os dejamos el código Python completo.

Para utilizar la herramienta debemos lanzar la Toolbox o el Addin y nos va a pedir los siguientes parámetros.

Máscara: sería el raster que sirve para definir la máscara, tamaño de pixel del raster de salida y snap a ajustas. Adjuntamos un ráster de referencia llamado malla100.

Workspace: para definir el entrono de trabajo

Estaciones: Fichero .txt con el identificador y coordenadas de la red de estaciones

Valores DC: Fichero .txt con los valores DC de un día determinado. El DC es un índice de sequía muy empleado para definir el riesgo de incendios forestales. Su cálculo lo tenemos implementado en un desarrollo Visual Basic, donde introduciendo datos de la red de estaciones se calcula el índice para un día determinado, dando como resultado el archivo .txt  a introducir en esta casilla.

Leyenda: opción para aplicar al mapa de salida una leyenda definida. Adjuntamos una leyenda de referencia llamada Legend.lyr

Borrar shapes intermedio: opción para mantener o borrar los .shp intermedios que genera la herramienta.Tool_Python_2

Desde el punto de vista de la programación es interesante, por ejemplo:

-La lectura del fichero de texto

– La generación de Outlayers con el módulo de Análisis espacial

Os animamos a utilizar esta interesante herramienta y a compartir con nosotros vuestras impresiones.

Además, el script completo queda como sigue.