环境:Windows10专业版 + PyCharm2022.2.3 + Python3.9.1 + GDAL(GDAL-3.4.3-cp39-cp39-win_amd64)
系列文章:
(一)Python+GDAL实现BSQ,BIP,BIL格式的相互转换
目录
(2)BIL(band interleaved by line format)
(3)BIP(band interleaved by pixel format)
一、BSQ,BIP,BIL格式简介
(插图来源于@CSDN溯水xiangling)
(1)BSQ (band sequential)
像素按波段顺序存储,先保存第一个波段,保存完毕后再保存第二个波段,以此类推。
优点:①便于进行波段间的运算;②便于进行波段间的运算;
(2)BIL(band interleaved by line format)
像素按行存储,先保存第一个波段的第一行,再保存第二个波段的第一行,以此类推。
优点:①像素的空间位置在列的方向上是连续的,既可以形象地表达空间分布特征,又可以反映像素的光谱特征 ;②具有较快的读取速度
(3)BIP(band interleaved by pixel format)
按像元顺序存储,先保存第一个波段的第一个像元,再保存第二个波段的第一个像元,以此类推。
优点:①便于进行像元间的运算;②可以清晰地反映像元的光谱特征;
这三种格式各有优缺点,通常根据具体的应用场景和处理需求来选择合适的格式。例如,BSQ格式适合在读取特定波段时提高效率,因为同一波段的数据连续存储,而BIP和BIL则在需要同时访问多个波段的像素时更为高效。
在实际的数字图像处理中,可能需要将一种格式转换为另一种格式以便于分析或与其他软件兼容。转换这些格式可以通过编程语言结合GDAL库进行转换,或者使用MATLAB编写脚本来实现批量转换;也可以通过数据处理和格式转换的软件或工具。每种格式都有其特定的存储结构和应用场景,在进行格式转换时需要考虑数据的完整性和准确性。由于数据格式的转换可能涉及到大量的数据处理和计算,因此在进行转换时还需要考虑计算机的性能和存储空间等因素。
二、代码实现
# author: jue_chen
# datetime:2022/10/18 20:13
# Attention: 转载,引用请注明出处
from osgeo import gdal
import os
# 设置PROJ_LIB的环境
os.environ['PROJ_LIB'] = "D:\\AppData\\Python\\Setup\\Lib\\site-packages\\osgeo\\data\\proj"
# 原图像8.tif的存储格式为BSQ
srcPath1 = "D:\\Project\\IDEA_Project\\RS01\\src\\rs01\\img\\8.tif"
desPath1 = "D:\\Project\\IDEA_Project\\RS01\\src\\rs01\\img\\8_1.tif"
# BSQtoBIP
gdal.Translate(desPath1, srcPath1, creationOptions=["INTERLEAVE=PIXEL"])
# 原图像8.3_BIL.tif存储格式为BIL
srcPath2 = "D:\\Project\\IDEA_Project\\RS01\\src\\rs01\\img\\8_3_BIL.tif"
desPath2 = "D:\\Project\\IDEA_Project\\RS01\\src\\rs01\\img\\8_3_BSQ.tif"
desPath3 = "D:\\Project\\IDEA_Project\\RS01\\src\\rs01\\img\\8_3_BIP.tif"
# BILtoBSQ
gdal.Translate(desPath2, srcPath2, creationOptions=["INTERLEAVE=BAND"])
# BILtoBIP
gdal.Translate(desPath3, srcPath2, creationOptions=["INTERLEAVE=PIXEL"])
三、实验结果
在ENVI软件里查看,可知转换成功