1. ESTARFM模型
模型参考：Zhu, Xiaolin, et al. "An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions." Remote Sensing of Environment 114.11 (2010): 2610-2623.
步骤参考：https://xiaolinzhu.weebly.com/open-source-code.html
我自己的步骤也不知道是不是正确的，因此写博客供大家参考，不正确的希望大家指正
1.1 IDL版本
注意：Landsat Collection 2 的转化系数为0.0000275-0.2，而MODIS09为0.0001。我不知道要不要考虑比例系数的问题，又觉得把代码的整形和浮点型会出毛病，所以我把L8C2的像元值*1000转变为整形。
(1)影像预处理，包括MODIS(MRT重投影、拼接)、Sentinel-2需要进一步投影转换、镶嵌等，最后为ENVI波段合成、裁剪等。
(2)Resize
参数设置：降尺度采样 30*30 方法： Bilinear
为了节省时间，我把MODIS影像先按L8大小大致裁剪
（3） Registration image to image
L8选择了B5(NIR)波段 对应的MODIS选择B2
因为我的L8和MODIS影像并不是在同一天合成，因此只能大致校正
（4）Clip
裁剪就很简单了，保证L8与MODIS边界即可，我比较喜欢ArcMAP的掩膜提取
（5）RUN
打开ENVI+IDL依次放入数据即可
1.2 PYTHON版本
朱博士团队已经写得很详细了，以后再补充详细步骤
唯一补充的就是 不要打开QQ拼音！不要打开QQ拼音！！不要打开QQ拼音！！！