こちらではgdal_merge.pyをつかって分割されたGeotiffをマージしました。
今回ここでやりたいことは、分割されているそれぞれのGeotiffに複数のバンドが格納されており、またそれぞれでデータタイプが異なるため、１枚のGeotiffにマージしつつ、バンドをごとのGeotiffを作成することです。
#gdalbuildvrtでvrtファイルを作成する
# -b でバンドを指定
gdalbuildvrt OUTPUT_b01.vrt -b 1 INPUT1.tif INPUT2.tif INPUT3.tif
gdalbuildvrt OUTPUT_b02.vrt -b 2 INPUT1.tif INPUT2.tif INPUT3.tif
gdalbuildvrt OUTPUT_b03.vrt -b 2 INPUT1.tif INPUT2.tif INPUT3.tif
#gdal_translateでvrtファイルをtifに変換する
#band1をByteタイプでエクスポート
gdal_translate -co COMPRESS=LZW -co BIGTIFF=YES -ot Byte -a_nodata 0 OUTPUT_b01.vrt OUTPUT_b01.tif
#band2をInt16タイプでエクスポート
gdal_translate -co COMPRESS=LZW -co BIGTIFF=YES -ot Int16 -a_nodata 0 OUTPUT_b02.vrt OUTPUT_b02.tif
#band3をFloat32タイプでエクスポート
gdal_translate -co COMPRESS=LZW -co BIGTIFF=YES -ot Float32 -a_nodata 0 OUTPUT_b03.vrt OUTPUT_b03.tif