Kilka razy temu pracowałem z codziennymi zdjęciami satelitarnymi temperatury powierzchni ziemi obszaru o wielkości Alaski.
Musiałem przetworzyć około tysiąca zdjęć. Ściągnąłem je, mozaikowałem, rzutowałem, przycinałem do granic kształtów, wyciągałem środki z przyciętych rastrów i narysowałem je na wykresie czas-temperatura. Wykres czas-temperatura był końcowym wyjściem. Zrobiłem to wszystko w ArcGIS i chociaż stworzyłem skrypt ArcPy, który pomógł zautomatyzować część pracy, wciąż był żmudny.
Teraz chcę zbudować program, który robi to samo, ale w 100% automatycznie. Chcę jedynie wprowadzić do programu współrzędne narożników, aby zdefiniować obszar zainteresowania i uzyskać końcową generowaną moc wyjściową, jaką jest wykres czas-temperatura.
Czy uważasz, że jest to możliwe przy użyciu Pythona z biblioteką GDAL i PostGIS jako magazynem danych rastrowych?
źródło