（一）Python+GDAL实现BSQ，BIP，BIL格式的相互转换-编程知识

（一）Python+GDAL实现BSQ，BIP，BIL格式的相互转换

环境：Windows10专业版 + PyCharm2022.2.3 + Python3.9.1 + GDAL(GDAL-3.4.3-cp39-cp39-win_amd64)

系列文章：

（一）Python+GDAL实现BSQ，BIP，BIL格式的相互转换

（二）BSQ,BIL,BIP存储格式的相互转换算法

（三）单波段图像的伪彩色合成：密度分割（含介绍OpenCV中的Mat类）

（四）图像的%2线性拉伸

（五）图像的标准假彩色合成

（六）图像的直方图均衡化

（七）图像的均值滤波

（八）图像的中值滤波

（九）图像的高斯低通滤波

（十）图像的梯度倒数加权平滑

（十一）图像的罗伯特梯度锐化

（十二）图像的Sobel梯度锐化

（十三）图像的拉普拉斯梯度锐化

一、BSQ，BIP，BIL格式简介

（1）BSQ (band sequential)

（2）BIL(band interleaved by line format)

（3）BIP(band interleaved by pixel format)

二、代码实现

三、实验结果

①BSQtoBIP

②BILtoBSQ

③BILtoBIP

一、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"])