Python dengan GDAL dan PostGIS sebagai pengganti ArcGIS

Beberapa kali yang lalu, saya bekerja dengan citra satelit suhu permukaan tanah harian suatu daerah dengan ukuran Alaska.

Saya harus memproses sekitar seribu gambar. Saya mengunduhnya, mosaik, diproyeksikan, dipotong ke batas-batas shapefile, diekstraksi berarti dari raster terpotong, dan diplot sarana pada grafik suhu-waktu. Grafik waktu-suhu adalah hasil akhir. Saya melakukan semua ini di ArcGIS, dan walaupun saya membuat skrip ArcPy yang membantu mengotomatisasi sebagian pekerjaan, itu masih melelahkan.

Sekarang saya ingin membangun program yang melakukan hal yang sama, tetapi secara otomatis 100%. Saya hanya ingin memasukkan ke program koordinat dari sudut-sudut untuk menentukan bidang yang diminati dan mendapatkan hasil akhir yang dihasilkan yang merupakan grafik suhu-waktu.

Apakah Anda pikir ini mungkin menggunakan Python dengan perpustakaan GDAL dan PostGIS sebagai penyimpan data raster?

Ya - adalah jawaban singkatnya. Saya pikir Anda akan membutuhkan NumPy dan SciPy sebagai bagian dari solusi Python Anda. Lihatlah modul scipy.ndimage ketika datang ke cara menghitung (apakah ini zona berarti?) Karena ini akan jauh lebih cepat daripada melakukannya hanya dengan NumPy. Selain itu, dengan menggunakan modul multiprosesor Python , Anda akan mendapatkan peningkatan kecepatan yang signifikan. Namun, pertahankan jumlah subproses menjadi lebih rendah karena pemrosesan raster bisa menjadi prosesor yang intensif dan Anda bisa memperlambat mesin daripada mempercepatnya. Saya sarankan 1 proses lebih sedikit dari core pada mesin Anda.