作者:CSDN @ _养乐多_
本文将介绍使用 STARFM(Spatial and Temporal Adaptive Reflectance Fusion Model)方法对 MODIS(Moderate Resolution Imaging Spectroradiometer)和 Landsat 遥感影像进行融合,并进行高斯平滑处理的 R 语言代码,该代码和 Google Earth Engne(GEE)平台代码配合使用。
具体而言,代码的主要功能如下:
- 使用StarFM方法融合MODIS和Landsat遥感影像:
- 通过循环迭代处理MODIS时间步骤,并在模拟时填充Landsat输出。
- 选择最近的一对好的MODIS/Landsat日期,一个在模拟日期之前,一个在之后。
- 运行StarFM方法,将MODIS和Landsat图像作为输入,生成通过模拟填充的Landsat图像。
- 对生成的融合图像进行高斯平滑处理:
- 定义高斯核函数,用于生成平滑窗口的权重。
- 将融合图像中的异常值(大于指定阈值的像素)视为噪声,并将其设为NA。
- 使用中值填充将NA值替换为该像素时间序列的中位数。
- 使用卷积运算和高斯核对填充后的图像进行平滑处理。
- 最后,代码保存生成的融合图像和平滑图像,并绘制了融合图像和平滑图像的对比图。
总之,以上代码主要实