gdalまたはpythonを使用して複数のファイルからサブセット化されたマルチバンド画像を作成しますか？

あるエリア（単一バンド）の時系列を前提として、それらをサブセット化し、gdalまたはpythonのgdalsバインディングを使用してサブセットを単一のマルチバンドファイル（例：vrt）に結合する方法は？

最初にそれらを「スタック」します。

gdalbuildvrt -separate -input_file_list my_files.txt my.vrt

my_filesには、使用するファイルのリストが特定の順序で含まれています。次に、gdal_translateでトリミングします。

gdal_translate -projwin ulx uly lrx lry my.vrt my.tif

私はgdal_translateとgdal_merge.pyを使用します：

gdal_translateを使用して、画像を変換してトリミング（サブセット化）します。bashスクリプトを使用して自動化できます。以下の行に沿って何か。

for f in *.tif;do gdal_translate -projwin ulx uly lrx lry "$f" "$f".cropped.tif ; done

gdal_merge.pyを使用して画像を「スタック」します。ここでは、スタックの順序を明示的に制御していません。これらは昇順（0-9、次にaz）で英数字順にスタックすると思います。

gdal_merge.py -separate -o myoutput.tif *.cropped.tif

myoutput.vrtがgdal_merge.pyで機能するかどうかわからないため、-o myoutput.tifを使用しました。GDALでサポートされている形式であるため、そうなると思いますが、テストしたことはありません。