Python з GDAL та PostGIS як заміна ArcGIS

Кілька разів тому я працював із щоденними супутниковими знімками температури поверхні ділянки з розміром Аляски.

Мені довелося обробити близько тисячі зображень. Я завантажував їх, мозаїкував, проектував, обрізав до межфайлів, витягував кошти з обрізаних растерів і розміщував засоби на графіку температури часу. Графік часу та температури був кінцевим результатом. Я все це робив в ArcGIS, і хоча я створив сценарій ArcPy, який допоміг автоматизувати частину роботи, це все ще було кропітко.

Тепер я хочу створити програму, яка робить те саме, але на 100% автоматично. Я просто хочу ввести в програму координати кутів, щоб визначити цікаву область і отримати кінцевий вихід, який є графіком температури та часу.

Як ви вважаєте, це можливо за допомогою Python з бібліотекою GDAL та PostGIS як растровим сховищем даних?

Так - це коротка відповідь. Я думаю, вам будуть потрібні NumPy та SciPy як частина вашого рішення Python. Погляньте на модуль scipy.ndimage, коли мова заходить про обчислення засобів (це зональні засоби?), Оскільки це буде набагато швидше, ніж робити це за допомогою лише NumPy. Також, використовуючи багатопроцесорний модуль Python , ви отримаєте значне збільшення швидкості. Однак зменшіть кількість підпроцесів, оскільки обробка растру може бути трудомісткою, і ви можете уповільнити роботу машини, а не прискорити її. Я пропоную на процес менше 1 процесу, ніж ядра на вашій машині.