Tạo hình ảnh multiband được đặt lại từ nhiều tệp bằng gdal hoặc python?

Đưa ra một chuỗi thời gian của một khu vực (một băng tần), làm thế nào bạn tập hợp chúng và kết hợp các tập hợp con thành một tệp đa cấp (ví dụ: vrt) bằng cách sử dụng các ràng buộc của gdal hoặc gdals cho python?

Trước tiên tôi sẽ 'xếp' chúng:

gdalbuildvrt -separate -input_file_list my_files.txt my.vrt

trong đó my_files chứa danh sách các tệp bạn muốn sử dụng, theo thứ tự cụ thể. Sau đó cắt với gdal_translate:

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

Tôi sẽ sử dụng gdal_translate và gdal_merge.py:

Dịch các hình ảnh để cắt (tập hợp chúng) bằng gdal_translate. Bạn có thể sử dụng tập lệnh bash để tự động hóa. Một cái gì đó dọc theo dòng dưới đây.

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

Sử dụng gdal_merge.py để 'xếp chồng' hình ảnh. Ở đây chúng tôi không rõ ràng kiểm soát thứ tự ngăn xếp. Tôi tin rằng những thứ này sẽ xếp chồng chữ và số theo thứ tự tăng dần (0-9, sau đó là az).

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

Tôi đã sử dụng -o myoutput.tif vì tôi không biết liệu myoutput.vrt sẽ hoạt động với gdal_merge.py. Tôi cho rằng nó sẽ như vậy, vì nó là định dạng được hỗ trợ GDAL, nhưng tôi chưa bao giờ thử nghiệm nó.