基于GEE的Landsat 7ETM+条带填补

项目简介

该项目使用Google Earth Engine (GEE)平台，对Landsat 5和Landsat 7卫星影像进行预处理与影像填补操作。主要功能包括影像的选取、波段处理、缺失影像的填补以及最终影像的导出。代码中的核心功能是通过空间回归方法对Landsat 7和Landsat 5影像进行时序配准与插值填补。

主要步骤

1. 定义感兴趣区域 (ROI)：
    - 使用四个坐标点绘制矩形多边形，确定需要处理的地理区域。

2. 选择影像数据集：
    - 使用Landsat 5和Landsat 7数据集 (`ee.ImageCollection`)，分别选择时间范围为2010年5月1日至2010年5月31日之间的影像。

3. 影像预处理：
    - 定义影像预处理函数，将选择的波段进行裁剪、波段重命名，并转换像元值范围，以便与其他数据进行处理和分析。

4. 影像填补：
    - 使用空间回归方法，计算两个影像集之间的关系，通过滑动窗口对像元进行回归分析并填补缺失区域。主要步骤包括：
      - 选择匹配的填补影像；
      - 计算回归模型，生成比例和偏移量；
      - 处理填补失败的情况，使用默认值填补。

5. 显示与导出结果:
    - 显示原始Landsat 7、Landsat 5影像和填补后的Landsat 5影像，并将结果导出为影像文件，便于后续分析。

核心代码说明

- 预处理函数 `preprocessImage`:
  该函数处理Landsat影像集中的波段并应用裁剪区域。通过乘以系数与增加偏移量将影像值缩放至0-1范围。

- 影像填补函数 `GapFill`:
  这是实现影像插值填补的主要函数，输入Landsat 7影像集和Landsat 5影像集，并对它们进行空间回归分析，最终生成无缺失的填补影像。

- 导出函数 `Export.image.toDrive`:
  使用GEE的导出功能将填补后的影像导出为文件，以便本地存储和处理。

如何运行

1. 打开Google Earth Engine的代码编辑器，并将上述代码粘贴到编辑器中。
2. 修改`startDate`和`endDate`参数以更改需要处理的影像时间范围。
3. 运行代码，预处理并显示影像，同时将填补