ArcGISの代替としてのGDALおよびPostGISを使用したPython

ほんの少し前に、私はアラスカのサイズの地域の毎日の地表面温度衛星画像を扱いました。

約1000枚の画像を処理する必要がありました。それらをダウンロードし、モザイク化し、投影し、シェープファイルの境界にクリップし、クリップされたラスターから平均を抽出し、平均を時間-温度グラフにプロットしました。時間と温度のグラフが最終的な出力です。これらはすべてArcGISで行いましたが、作業の一部を自動化するのに役立つArcPyスクリプトを作成しましたが、それでもまだ骨の折れる作業でした。

同じことを行いますが、100％自動的に実行するプログラムを作成したいと思います。関心のある領域を定義し、時間と温度のグラフである最終的な出力を生成するために、プログラムにコーナーの座標を入力したいだけです。

これは、GDALライブラリとPostGISをPythonでラスターデータストアとして使用することで可能だと思いますか

はい-短い答えです。Pythonソリューションの一部として、NumPyとSciPyが必要になると思います。これがNumPyだけで計算するよりもはるかに速いので、平均を計算するとき（これらのゾーン平均ですか？）については、scipy.ndimageモジュールを見てください。また、Python マルチプロセッシングモジュールを使用すると、速度が大幅に向上します。ただし、ラスター処理はプロセッサを集中的に使用する可能性があり、マシンを高速化するのではなく低速化する可能性があるため、サブプロセスの数を減らしてください。あなたのマシンのコアよりも1プロセス少ないことをお勧めします。