python大气校正_Sentinel-2卫星影像的大气校正方法

4 结果分析

为了初步分析三种大气校正方法的结果, 选取10 m分辨率(B2, B3, B4, B8)及20 m分辨率(B5, B6, B7, B8a)8个波段, 并选择冬小麦、裸土、水泥路和油菜地4种地物的像元实测值、大气上层表观反射率及三种模型的校正结果进行分析。其中冬小麦和油菜地为混合像元, 首先通过第2.2节的方法计算植被覆盖度, 得到混合像元中植被与土壤的比例, 并结合实测地物光谱计算出混合光谱, 结果如图4所示。

由图4可见, Sentinel-2卫星影像经过三种大气模型校正后, 其光谱曲线与实测地物的光谱曲线的变化趋势一致, 有较高的拟合度, 蓝波段Band 2经大气校正后相比于大气上层表观反射率TOA明显下降, 红边、近红外波段明显上升。在冬小麦和油菜地的光谱曲线中, Band 2的三种大气校正结果高于实测地物的反射率, 误差较大。Sen2cor方法与SMAC的误差约为0.04, 6S模型的误差约为0.06;

图4

Fig. 4

图4 各类地物大气校正前后的反射光谱曲线。(a)冬小麦; (b)油菜地; (c)裸土; (d)水泥路面Fig. 4 Reflectance spectral curves of different ground objects before and after atmospheric correction. (a) Winter wheat; (b) rape field; (c) bare soil; (d) cement pavement

Band 3~Band 5的大气校正结果与实测值拟合效果较好, 误差均小于0.02; 在大于705 nm的Band 6~Band 8a波段中, 实测地物的反射率高于三种大气校正结果, 误差范围为0.01~0.04。在土壤光谱曲线中, 大气校正结果相比于大气上层表观Band 2反射率明显下降, Band 4~Band 8a波段的反射率明显增加, Band 2的三种大气校正结果高于实测地物的反射率, 其余波段的实测地物反射率均高于三种大气校正