Tutorial de interpolación triangular geoespacial con Python, Scipy, Geopandas y Rasterio

Bajo el concepto de “Python geoespacial aplicado” hemos desarrollado algunos procedimientos / tutoriales en Python de algunas tareas comunes de análisis espacial realizadas en software GIS de escritorio. El objetivo no es reinventar la rueda, sino explorar las herramientas y bibliotecas de Python actuales que pueden crear, analizar y representar datos espaciales vectoriales y ráster.

La interpolación triangular es uno de los varios tipos de interpolación disponibles tanto en Python como en software GIS, sin embargo, la ventaja de trabajar con Python es que la interpolación es una función en la que puede obtener el valor interpolado en un punto específico mientras que en el software GIS es necesario para crear un ráster para luego muestrear valores a partir del ráster (.. hasta donde sabemos).

Hemos creado un tutorial con un procedimiento completo en Python para importar puntos con elevación como atributo, crear una función de interpolación triangular y generar dos salidas espaciales: un ráster geoespacial interpolado en formato TIFF y un shapefile con atributo de elevación para otro conjunto de puntos. El tutorial utiliza varias bibliotecas de Python como Matplotlib, Rasterio, Geopandas, Scipy.

Datos de entrada
Puede descargar los datos de entrada en este link: