Antecedentes
Estoy comenzando un proyecto que se hará íntegramente a través de Python, salvo quizás al final para la visualización. Hasta ahora lo he hecho con ArcGIS y arcpy, pero nos gustaría pasarlo a una plataforma de código abierto.
He estado buscando y no he podido averiguar cómo encaja exactamente QGIS con GDAL/OGR (y otros proyectos de código abierto como GRASS). Parece que QGIS incluye GDAL/OGR pero tiene muchos añadidos propios. Estoy tratando de averiguar cuál es la ruta más fácil para realizar un montón de operaciones simples de alto nivel.
Por el momento lo que estoy haciendo con arcpy es la conversión entre rasters y puntos, la extracción de valores raster en puntos existentes, la extracción de tablas y geometría a NumPy, la hidrología, la búsqueda de intersecciones y la distancia más cercana entre características. Sólo hay una operación que imagino que habría que rehacer bastante (insertar puntos cada x metros a lo largo de un conjunto de vectores).
Mi pregunta
¿Es la interfaz Python de QGIS suficiente para manejar esto, o hay alguna ventaja en el uso de GDAL/OGR y otros paquetes? Me gustaría que fuera lo más sencillo y legible posible, ya que mucha otra gente lo utilizará y modificará cuando hayamos terminado.
