こちらでは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

参考：gdal_translate, gdalbuildvrt