注：本程序适合任何所有四个波段的国产卫星数据主要分为如下几步：
1、 GF4号数据预处理
2、 MOD09地表反射率数据制作
3、 百分比匹配查找暗像元
4、 气溶胶反演
5、 非暗像元进行气溶胶插值
6、 精确大气校正
第一步：
1、LC8转换存储格式以及波段选择减少数据量，GF4转换存储格式以及波段选择，因为GF4第一波段是全色波段，无用。
2、基于参考影像的RPC正射校正，DEM是全球200m分辨率DEM。从图中看几何校正效果明显
3、基于控制点的RPC正射校正，DEM是全球200m分辨率DEM，20/50,小于半个象元，满足需求
4、矢量裁剪
5、定标为表观反射率，用插件添加必要信息
第二步：
1、数据准备：MOD09，按照年份存储好（本次以两年数据为例讲解，
每年两景数据，演示镶嵌以及取最小的功能）
2、使用插件进行处理（需要授权文件），在一个就是可以不进行重投影，建议
不进行掩膜提取（因为ENVI5.3的shp裁剪接口不稳定）
3、通过对比可以取最小
第三步：
1、数据准备（高分四号表观反射率，地表反射率数据库（相同地区）
正弦投影）
2、百分比匹配（阈值的选择）
3、云，水体，城市地区就会变空值
视频地址教程地址：
注：6个步骤全
https://www.bilibili.com/video/av14679327/
反演结果图
参考文献
孙林, 于会泳, 傅俏燕,等. 地表反射率产品支持的GF-1PMS气溶胶光学厚度反演及大气校正[J]. 遥感学报, 2016, 20(2):216-228.
xxin blog , 版权所有丨如未注明 , 均为原创丨本网站采用BY-NC-SA协议进行授权 , 转载请注明高分四号气溶胶反演（视频教程）！
