ArcGIS'in yerine GDAL ve PostGIS'li Python

Birkaç kez önce, Alaska büyüklüğünde bir alanın günlük kara yüzey sıcaklığı uydu görüntüleri ile çalıştım.

Yaklaşık bin görüntü işlemek zorunda kaldım. Onları indirdim, mozaikledim, yansıttım, şekil dosyası sınırlarına kırptım, kırpılan rasterlerden çıkarılan araçları aldım ve araçları bir zaman-sıcaklık grafiği üzerine çizdim. Zaman-sıcaklık grafiği nihai çıktıydı. Tüm bunları ArcGIS'de yaptım ve işin bir kısmını otomatikleştirmeye yardımcı olan bir ArcPy betiği oluştursam da, yine de titizdi.

Şimdi aynı şeyi yapan bir program oluşturmak istiyorum, ancak% 100 otomatik olarak. İlgilenilen alanı tanımlamak ve zaman-sıcaklık grafiği olan son çıktıyı elde etmek için programa köşelerin koordinatlarını girmek istiyorum.

Python'u GDAL kütüphanesi ve PostGIS'i raster veri deposu olarak kullanmak mümkün mü?

Evet - kısa cevap. Sana ihtiyacım olacağını düşünüyorum numpy ve scipy Python çözümün bir parçası olarak. Hesaplama araçları söz konusu olduğunda scipy.ndimage modülüne bir göz atın (bu bölgesel araçlar mı?) Bu sadece NumPy ile yapmaktan çok daha hızlı olacaktır. Ayrıca, Python çoklu işlem modülünü kullanarak önemli bir hız kazancı elde edersiniz. Ancak, tarama işlemi işlemci açısından yoğun olabileceğinden ve makineyi hızlandırmak yerine yavaşlatabileceğiniz için alt işlemlerin sayısını azaltın. Makinenizdeki çekirdeklerden 1 daha az işlem öneririm.