Buat gambar multiband yang ditransfer dari beberapa file menggunakan gdal atau python?

Diberikan serangkaian waktu dari suatu area (pita tunggal), bagaimana Anda mengelompokkannya dan menggabungkan himpunan bagian menjadi file multiband tunggal (mis. Vrt) menggunakan gdal atau binding gdals untuk python?

Pertama-tama saya akan 'menumpuk' mereka:

gdalbuildvrt -separate -input_file_list my_files.txt my.vrt

di mana my_files berisi daftar file yang ingin Anda gunakan, dalam urutan tertentu. Kemudian pangkas dengan gdal_translate:

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

Saya akan menggunakan gdal_translate dan gdal_merge.py:

Menerjemahkan gambar untuk memotong (mengelompokkannya) menggunakan gdal_translate. Anda dapat menggunakan skrip bash untuk mengotomatisasi. Sesuatu di sepanjang garis di bawah ini.

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

Gunakan gdal_merge.py untuk 'menumpuk' gambar. Di sini kita tidak secara eksplisit mengendalikan urutan tumpukan. Saya percaya ini akan menumpuk secara alfanumerik dalam urutan menaik (0-9, lalu az).

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

Saya menggunakan -o myoutput.tif karena saya tidak tahu apakah myoutput.vrt akan bekerja dengan gdal_merge.py. Saya berasumsi akan, karena ini adalah format yang didukung GDAL, tetapi saya belum pernah mengujinya.