Python con GDAL e PostGIS in sostituzione di ArcGIS

Poche volte, ho lavorato con immagini satellitari giornaliere della temperatura della superficie terrestre di un'area delle dimensioni dell'Alaska.

Ho dovuto elaborare circa mille immagini. Li ho scaricati, mosaicati, proiettati, ritagliati ai confini dello shapefile, estratto i mezzi dai raster ritagliati e tracciato i mezzi su un grafico tempo-temperatura. Il grafico tempo-temperatura era l'output finale. Ho fatto tutto questo in ArcGIS, e anche se ho creato uno script ArcPy che ha aiutato ad automatizzare una parte del lavoro, è stato comunque scrupoloso.

Ora voglio costruire un programma che faccia la stessa cosa, ma automaticamente al 100%. Voglio semplicemente inserire nel programma le coordinate degli angoli per definire l'area di interesse e ottenere l'output finale generato che è il grafico tempo-temperatura.

Pensi che sia possibile utilizzare Python con la libreria GDAL e PostGIS come archivio dati raster?

Sì - è la risposta breve. Penso che avrai bisogno di NumPy e SciPy come parte della tua soluzione Python. Dai un'occhiata al modulo scipy.ndimage quando si tratta di calcolare i mezzi (sono questi mezzi zonali?) In quanto sarà molto più veloce rispetto a farlo con solo NumPy. Inoltre, utilizzando il modulo multiprocessing di Python , otterrai un significativo guadagno di velocità. Tuttavia, mantenere basso il numero di sottoprocessi poiché l'elaborazione raster può richiedere molto tempo al processore e si potrebbe finire per rallentare la macchina anziché accelerarla. Suggerisco 1 processo in meno rispetto ai core sulla tua macchina.