Python con GDAL y PostGIS como reemplazo de ArcGIS

Hace algunas veces, trabajé con imágenes satelitales diarias de la temperatura de la superficie terrestre de un área del tamaño de Alaska.

Tuve que procesar alrededor de mil imágenes. Los descargué, mosaicé, proyecté, recorté los límites del archivo de forma, extraje los medios de los rásteres recortados y tracé los medios en un gráfico de tiempo y temperatura. El gráfico de tiempo-temperatura fue el resultado final. Hice todo esto en ArcGIS, y aunque creé un script de ArcPy que ayudó a automatizar una parte del trabajo, todavía fue minucioso.

Ahora quiero construir un programa que haga lo mismo, pero 100% automáticamente. Simplemente quiero ingresar al programa las coordenadas de las esquinas para definir el área de interés y obtener el resultado final generado, que es el gráfico de tiempo y temperatura.

¿Crees que esto es posible usando Python con la biblioteca GDAL y PostGIS como un almacén de datos ráster?

Sí, es la respuesta corta. Creo que necesitará NumPy y SciPy como parte de su solución Python. Eche un vistazo al módulo scipy.ndimage cuando se trata de calcular medios (¿son estos medios zonales?), Ya que será mucho más rápido que hacerlo solo con NumPy. Además, al usar el módulo de multiprocesamiento Python , obtendrá una ganancia de velocidad significativa. Sin embargo, mantenga baja la cantidad de subprocesos, ya que el procesamiento ráster puede requerir un uso intensivo del procesador y podría terminar ralentizando la máquina en lugar de acelerarla. Sugiero 1 proceso menos que los núcleos en su máquina.