python 大气校正_基于6s模型的遥感影像大气校正方法

基于6s模型的遥感影像大气校正工程化实现

目的：针对光学遥感影像(Landsat-8、Sentinel-2、GF-1、GF-2)的大气校正自动化实现方法，大多数是使用ENVI的FLASSH模块进行大气校正，但是该模块工程化起来比较麻烦，虽然现在老美提供了FLASSH模块的调用接口。但是当需要对大规模数据进行处理时，尤其是在linux环境下，ENVI对环境配置特别麻烦，此外该模块的处理速度相对较慢。因此在实际工作中选用了开源的6s模型，来进行大气校正批处理。

一、大气校正的作用及原理

1.1 几个重要的概念

数字量化值(Digital Number: DN)：用来描述还没有校准到具有意义单位的像素值。是对卫星接收到辐射进行量化的一个数值，大小与量化深度有关，比如8bit量化，Max(DN)=255，无单位，无实际意义。

辐射亮度(Radiance)：是指单位面积、单位波长、单位立体角的辐射通量，辐射亮度值有单位标识，如W/cm2.μm.sr(瓦特/平方厘米.微米.球面度)，将未定标的DN图像转换成辐射亮度图，只需要一个线性转换公式即可，需要用到偏移和增益两个参数，一般在元数据文件中可以获取.

表观反射率(Top of Atmosphere Reflectance)：又叫做卫星上反射率，大气表层反射率，是未经大气校正的，卫星上接收到的反射率。

地表反射率：是地物真实反射率，由表观反射率去处大气影响后得到的反射率。可以通过辐射传输模型、统计模型等去除大气影响。

地表反射率：是地物真实反射率，由表观反射率去处大气影响后得到的反射率。可以通过辐射传输模型、统计