0

我正在尝试按照本教程将一组 .tif 图像转换为一个多波段 .tif 。我没有收到任何错误,但是每次我运行脚本时,当我有比这更多的输入 .tifs 时,它都会生成一个三波段图像。我已经尝试了 -r 标志,它没有改变任何东西。

我有:

wdir="/Users/<mydir>"
creoleLA="$wdir/data/s2l2a_tiffs/T15RVP"
dhaka="$wdir/data/s2l2a_tiffs/T45QZG"
easthouston="$wdir/data/s2l2a_tiffs/T15RTN"
kolkata="$wdir/data/s2l2a_tiffs/T45QXF"
murcia="$wdir/data/s2l2a_tiffs/T30SXH"
wuhan="$wdir/data/s2l2a_tiffs/T50RKU"
#=============================================
# Main script
#=============================================
read -p 'ROI: ' roi
indir=${!roi}
outdir="$wdir/data/s2l2a_tiffs"

ls ${indir}*"_B"*.tif > "$outdir/btif_${roi}.txt"
gdalbuildvrt -separate -overwrite -input_file_list "$outdir/btif_${roi}.txt" "$outdir/S2L2A_${roi}.vrt"
gdal_translate -strict -ot uint16 "$outdir/S2L2A_${roi}.vrt" "$outdir/S2L2A_${roi}_mb.tif"

rm "$outdir/S2L2A_${roi}.vrt"
rm "$outdir/btif_${roi}.txt"

"$outdir/btif_${roi}.txt"像这样的 GeoTIFF 文件路径的纺织品在哪里:

/Users/<mydir>/data/s2l2a_tiffs/T15RVP_20200831T164849_B01_60m.tif
/Users/<mydir>/data/s2l2a_tiffs/T15RVP_20200831T164849_B02_60m.tif
...

我正在处理 Sentinel-2 图像并使用 OSX 11.6。

4

0 回答 0