Python修复Landsat ETM+条带：GDAL与OpenCV方案对比

Alice

发布时间：2026-10-13| 2人浏览过

原创

Landsat 7 ETM+影像因扫描线校正器故障产生条带，影响数据可用性。本文介绍两种Python修复方案：利用GDAL的FillNodata函数进行波段插值，以及使用OpenCV的inpaint算法进行像素重建。对比显示，GDAL处理速度快但质量一般，OpenCV修复效果更优但耗时较长。文章提供完整代码示例及参数说明，帮助开发者根据需求选择合适工具。

GDAL修复Landsat ETM+影像条带

利用GDAL库修复Landsat ETM+卫星影像中的条带问题

由于Landsat 7 ETM+卫星的扫描线校正器发生故障，后续获取的影像数据中出现了明显的平行条带噪声。如图所示，这些条带均匀分布在整个影像中，严重影响了数据的可视化和分析质量。

以下是使用GDAL库进行条带修复的具体代码实现：


def gdal_repair(tif_name, out_name, bands):
  """
    tif_name(string): 源影像名
    out_name(string): 输出影像名
    bands(integer): 影像波段数
  """
  # 打开影像文件
  tif = gdal.Open(tif_name)
  
  # 根据文件类型获取对应的驱动程序
  driver = gdal.GetDriverByName('GTiff')
  
  # 根据指定文件的驱动程序，使用现有数据集创建新的可写数据集
  # 所有支持创建新文件的驱动程序都支持该`CreateCopy()`方法，   # 但仅`Create()`部分支持该方法
  # CreateCopy的第一个参数为目标文件名，第二个参数为源数据集
  # 第三个参数的值是`0`或`1`，值是`0`。即使无法将原始数据准确地转换为目标数据，程序仍将执行
  new_img = driver.CreateCopy(out_name, tif, 0)
 
  for i in tqdm(range(1, bands)):
    # 分别对每个波段处理
    band = new_img.GetRasterBand(i)
    
    # 使用FillNodata对条带部分进行插值
    gdal.FillNodata(targetBand = band, maskBand = band, maxSearchDist = 15, smoothingIterations=0)
    
    # 将修复好的波段写入新数据集中
    new_img.GetRasterBand(i).WriteArray(band.ReadAsArray())

经过GDAL修复处理后的影像效果如图所示：

OpenCV修复Landsat ETM+影像条带

利用OpenCV库修复Landsat ETM+卫星影像中的条带问题

以下是使用OpenCV库进行条带修复的具体代码实现：


def cv2_repair(tif_name):
  # 读取tif影像
  tif_data = gdal_array.LoadFile(tif_name).astype('float32')

  # 获取掩膜
  mask = tif_data.sum(axis=0)
  mask = (mask == 0).astype(np.uint8)
  
  bands = tif_data.shape[0]

  res = []
  for i in tqdm(range(bands)):
    # cv.Inpaint(src, inpaintMask, dst, inpaintRadius, flags)
    # src：源图像，可以是8位、16位无符号整型和32位浮点型1通道或者8位无符号3通道
    # inpaintMask：掩膜，8位无符号整型
    # dst：和源图像具有一样大小的输出
    # inpaintRadius：算法考虑的每个已修复点的圆形邻域的半径     # flags：修复算法类型，可选cv2.INPAINT_NS和cv2.INPAINT_TELEA
    
    repaired = cv2.inpaint(tif_data[i], mask, 3, flags=cv2.INPAINT_TELEA)
    res.append(repaired)

  return np.array(res)

经过OpenCV修复处理后的结果影像如图所示：