大气校正6S与FLAASH计算与对比

（1）FLAASH大气校正

大气校正的目的是消除大气和光照等因素对地物反射的影响，获得地物反射率和辐射率、地表温度等真实物理模型参数，用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响，消除大气分子和气溶胶散射的影响。

FLAASH大气校正的特点：支持传感器种类多；FLAASH采用了MODTRAN4+辐射传输模型，该算法精度高；通过影像像素光谱上的特征来估计大气的属性，不依赖遥感成像时同步测量的大气参数数据；可以有效地去除水蒸气/气溶胶散射效应，同时基于像素级的校正，矫正 目标像元和邻近像元交叉辐射的“邻近效应”；对由于人为抑止而导致波谱噪声进行光谱平滑处理。作为结果，除了真实地表反射率外，还可以得到整幅图像内的能见度、卷云与薄云的分类影像、水气含量数据。

（2） ENVI5.1下的FLAASH大气校正步骤

在Toolbox中打开FLAASH工具：/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction。启动FLAASH Atmospheric Correction Module Input Parameters面板。

1） Input Radiance Image：选择辐射定标结果数据，要求为BIL存储格式；
2） Output Reflectance File：设置输出路径和文件名；
3） Output Directory for FLAASH Files：设置其他文件输出目录；
4） Scene Center Location：自动获取；
5) Sensor Type：Landsat-8 OLI；Sensor Altitude:自动读取；Pixel Size:自动读取；
6) Flight Time：在头文件里查找。单击Ok。

（1）基于6S模型大气校正

Landsat 8 地表反射率生产基于6S（version 1.1）模型，采用基于查找表逐像元计算的方式，6S模型需要的输入参数来自于和Landsat 8卫星过境时间同步的Terra MODIS大气产品数据。其中水蒸汽采用MODIS MOD05产品，臭氧数据来自MODIS MODO7 Total Ozone 产品，气溶胶光学厚度来自MODIS MOD04 Aerosol Optical Thickness产品。针对Landsat-8 OLI传感器，利用其大气校正的影响因素太阳天顶角、水汽、臭氧、气溶胶、以及高程等五个输入参数建立了一个五维度的大气校正系数查找表，该查找表适用于所有的Landsat-8 OLI数据的地表反射率反演。利用处理后的MODIS大气参数影像、DEM数据，以及Landsat8数据的成像几何参数，基于查找表逐像元插值计算该点的大气校正系数，最后再根据大气校正系数得到其地表反射率图像。其中所采用的插值算法为多维线性插值法。

（2）输入的各个参数的介绍

1）需要输入的为几何参数igeom，其中有8情况，当选择igeom=0时，需要手动输入太阳高度角、方位角，卫星高度角、方位角以及影像获取月份与日期。剩下的7种为内置的一些卫星参数，如AVHRR、TM、GOES等，大部分都是需要输入日期，行列号等。

2）大气模式idatm的选择，6S里内置了0-6共其中大气模式，依次为无大气吸收，热带，中纬度夏季，中纬度冬季，亚寒带夏季，亚寒带冬季以及美国标准62大气模式，或者选择自定义大气廓线，选择7的话，需要输入海拔，压强，温度，水汽浓度，臭氧浓度。选择8需要输入水汽与臭氧浓度，大气廓线使用美国62。

3）选择气溶胶模式iaer，内置了7中标准气溶胶模式，为0-无气溶胶模式，1-大陆型，2-海洋型，3-城市型，5-沙漠背景的shettle模型，6-生物燃烧模型，7-平流层模型。或者自定义自己的气溶胶模式，当iaer=4时，需要输入各组分的百分比：粉状尘，水溶性物质，海洋性物质，煤烟。8为最多四种模型联合运算的一个正态分布模型。9为改进型的一种伽马分布模型。10为Junge谱幂分布模型。11为选择太阳光度计的测量值进行计算，12为读取已经设定好的参数进行计算。一般情况下如果不想做气溶胶的模式研究，就直接用刚开始介绍的7个标准模式。

4）为气溶胶浓度，关于这个由两种输入方式，气溶胶光学厚度taer55或者可见度v。可见度v的单位为千米，如选择光学厚度taer55，可见度v的值就为0，taer55为550纳米处的气溶胶光学厚度。在气溶胶模式iaer中如果选择无气溶胶，可见度v的值就要为-1。其中气溶胶光学厚度与可见度有一个公式可以转换，为了快速计算，也可以使用4.8除光学厚度得到可见度，但这只是一个约值。

5）目标海拔高度xps的计算，如果xps>=0，那么表示目标在水平面上。xps<0表明你知道地面的高程，所以输入具体的值，不过前面加个负号。

6）传感器Xpp高度的输入，xpp=-1000 表明传感器为星载传感器，xpp=0 表明传感器在地面上-100

7）光谱条件iwave的输入，其中6S内置了很多传感器的光谱条件，如MODIS，TM，AVHRR，SPOT，NOAA等。也可以使用自定义的光谱条件。其中光谱响应函数的计算是以2.5微米为间隔。

由于6S所支持的光谱区间为0.25-4，代入上式iinf=(0.25-.25)/0.0025+1.5 =1.5=1，isup=(4-.25)/0.0025+1.5=1501.5=1501。所以6S所能处理的所有光谱曲线就位1-1501条。以TM数据蓝光波段为例，其光谱范围为0.430-0.560，代入上式可以得到iinf=(0.430-.25)/0.0025+1.5=73，isup=(0.56-.25)/0.0025+1.5=125，在TM的光谱响应曲线的就可以看做前72个数值为0，后1376个也为0值。在6S之中，计算各种参数的时候，尤其是在计算光谱响应的起始与结束值时，定义的为整型，但计算结果可能会为一个小数，这就要考虑截断了，FORTRAN中并没有使用四舍五入的情况，而是直接取整。对于6S没有内置的一些光谱响应也就可以自己按照这个规则来写一个Fortran文件。

8）地表特性的选择，这是一个比较复杂的选项，可以选择地面作为朗伯体或者各向异性来构建BRDF。inhomo选择0代表地表是考虑地面为朗伯体，选择1代表地表各向异性。如果inhomo=0进入第二步的选择，再次选择0表明没有方向效应的影响，然后输入roe（表面反射比）。选择1表明有方向效应的影响，然后选择内置的一些BRDF模型来进行计算。如果inhomo=1，那么接着输入igrou1（roc）（目标的反射）、igrou2（roe）（周围环境反射）、radius r（表面半径）。

9）rapp大气校正模型， 如果rapp<-1那么就不进行大气改正，这就为遥感之中的正演，即为通过地表反射率来求得表观反射率。如果rapp>0 程序将重新赋值表观反射率，使用输入的这个辐射亮度值进行表观反射率的转换。

为了评价和验证两种模型对Landsat 8遥感影像大气辐射校正的效果, 对FLAASH 和6S 模型校正后的图像与辐射定标求表观反射率的图像分别进行典型地物光谱曲线和NDVI分析。

典型地物光谱分析

分别对校正以后的影像中选取植被、水体、建筑物、空地四种典型地物进行光谱曲线分析。典型地物具体如表：

具体操作步骤
点击File—>open—>选LC81190422014347_Radiance_Reflectance.dat
点击Region of Interest（ROI）tool
如下图所示选择四种典型地物，然后保存。如图。

在ENVI主菜单中，选择File>Open，打开所保存的样本文件。
回到Spectral Library Builder面板中，选择Options>from ROI/EVF from input file，选择高光谱文件作为波谱来源。如图：

点击Plot即可看到光谱图：

分别将6S和Flaash校正后的四钟典型地物的波普数据保存输出并且分别绘制曲线图进行比较：

经大气校正后植被反射率在近红外波段上升明显，对于植被的大气校正，6S大气校正的效果要好于Flaash大气校正；对于水体，Flaash大气校正后的结果要好于6S大气校正的结果，这是因为Flaash的大气校正是基于整幅影像的，通过图像上的黑暗像元来获取整幅影像的气溶胶厚度，而6S是逐一像元校正的。因此对于不同用途的影像我们应该使用合适的大气校正方法。

更多遥感知识，关注我的微信公众号：