Python com GDAL e PostGIS como uma substituição do ArcGIS

Há poucas vezes, trabalhei com imagens diárias de satélite da temperatura da superfície de uma área do tamanho do Alasca.

Eu tive que processar cerca de mil imagens. Fiz o download deles, em mosaico, projetados, cortados nos limites do arquivo de forma, extraímos os meios dos rasters cortados e plotei os meios em um gráfico de temperatura e tempo. O gráfico tempo-temperatura foi a saída final. Fiz tudo isso no ArcGIS e, apesar de criar um script do ArcPy que ajudou a automatizar uma parte do trabalho, ainda era trabalhoso.

Agora eu quero criar um programa que faça a mesma coisa, mas 100% automaticamente. Eu apenas quero inserir no programa as coordenadas dos cantos para definir a área de interesse e obter a saída final gerada, que é o gráfico de temperatura e tempo.

Você acha que isso é possível usando o Python com a biblioteca GDAL e o PostGIS como um armazenamento de dados raster?

Sim - é a resposta curta. Eu acho que você precisará do NumPy e SciPy como parte da sua solução Python. Dê uma olhada no módulo scipy.ndimage quando se trata de meios de cálculo (são meios zonais?), Pois isso será muito mais rápido do que com apenas o NumPy. Além disso, usando o módulo de multiprocessamento Python , você obterá um ganho de velocidade significativo. No entanto, mantenha o número de subprocessos baixo, pois o processamento de varredura pode exigir muito do processador e você pode acabar diminuindo a velocidade da máquina, em vez de acelerá-la. Sugiro 1 processo a menos do que núcleos em sua máquina.