使用遥感手段获取水体信息时，由于大气分子和气溶胶等大气成分吸收与散射的影响，卫星传感器接收到的水体离水辐射信号一般不多于大气层顶总辐射亮度的10%^[1]，有效剔除大气影响是水色遥感定量反演的关键问题。目前水色遥感的大气校正方法主要有暗像元法^[2]、经验线性法^[3]以及直方图均衡化^[4]等基于图像特征的方法，Gordon提出的基于NIR的一类水体标准大气校正方法、改进的NIR-SWIR混合校正方法^[5]以及辐射传输模型法等。对于二类水体，在近红外波段清洁水体反射率为0的假设不成立，且某些传感器没有设置短波红外波段，暗像元法和NIR-SWIR等方法并不适用。辐射传输模型方法依据电磁波在大气中的辐射传输原理来计算地物反射率，依据同步大气参数，如能见度、气溶胶参数等，可在不同大气条件下以较高的精度获得真实地表反射率。其中6S模型^[6]考虑了目标物的海拔高度、地表非均匀状况和气体对辐射的吸收影响，对分子和气溶胶散射作用的计算使用近似和逐次散射算法，是一种运行速度较快，操作相对简单，应用较为广泛的大气校正方法。

水体光谱信息微弱，辐射探测的误差较大，6S辐射传输模型在陆地大气校正中精度较高，但对于水体的校正精度仍在改进和研究当中^{[7,8,9,10,11]}。Vermote等^[6]认为6S大气校正模型在水色遥感中只适用于做敏感性研究。另外，内陆湖泊水体的定量遥感反演往往需要高空间分辨率遥感数据，而当前针对高分数据的大气校正算法大多应用于陆地区域^[12,13,14]，在水体区域讨论较少。例如，刘佳等^[15]的研究表明，除水体外，6S和FLAASH大气校正的结果相差不大，其中偏差最大的在蓝光波段，二者具有较好的可比性。因此，需要进一步研究6S大气校正模型在基于高分影像的水色遥感中的参数敏感性和精度。

6S辐射传输模型输入参数的敏感性研究表明，气溶胶类型以及组分是影响遥感反演水体反射信号精度的重要因素之一^[16]，而气溶胶模式的假定是反演气溶胶光学厚度和大气校正的主要误差来源之一^[17]。由于内陆湖泊水体上空气溶胶特性复杂多变，使用默认参数会导致较大误差，因此根据水域特点进行参数敏感性分析，优化与选择大气参数，对于改进6S大气校正模型的区域应用非常重要。

以典型内陆湖泊水体太湖为研究区，使用2016年4月29日的高分一号宽幅相机影像（GF-1/WFV）和同步采样的地面实测数据，进行6S参数敏感性分析并遴选适用于太湖水体的大气参数。针对WFV影像特征，逐像元计算观测几何，考虑分区气溶胶类型，使用分区暗像元和Spline插值确定气溶胶光学厚度（Aerosol Optical Depth，AOD），提出了一种基于6S辐射传输模型的逐像元大气校正方法，并与FLAASH大气校正进行了结果比较。研究结果对改进6S模型在湖泊水体的区域应用具有重要参考意义。

太湖是我国第三大淡水湖泊，位于119°55.3′~120°53.6′E，30°56′~31°33′N，地处江苏和浙江省交界处和我国经济发展最具活力的长三角地区。太湖面积2 427.8 km²，平均水深2.12 m，属于典型的亚热带东南季风气候区，夏季高温多雨，冬季温和少雨。近30 a来由于该地区经济持续快速的增长导致环境问题显著，诸如蓝藻水华频发、水环境质量下降等，水质状态为中富营养化到重富营养化^[18]。2016年4月29日在太湖地区进行了野外采样，采样时间从上午9点至下午4点，17个有效样点的分布如图1所示。

水体的野外光谱测量使用美国ASD(Analytical Spectral Devices)野外光谱辐射仪(ASD FieldSpec.)。为减少水体镜面反射和船体自身阴影的影响，测量采用二类水体水面以上光谱测量方法^[19]，测量时距离水面0.5 m。测量的数据有标准灰板(反射率为30%)、水体和天空光的光谱辐亮度信息，每个对象采集10条光谱数据，选取每个样点重复测量的遥感反射率的中值作为该点的反射率。在测量得到的光谱曲线中，选择350~900 nm的数据计算遥感反射率^[19]。利用GF-1/WFV的光谱响应函数和地面实测光谱反射率值，按照式（1），模拟计算WFV各个波段的等效遥感反射率：

其中:Ri代表传感器第i通道的光谱反射率，λsi是通道i的起始波长，λei是通道i的终止波长，Rrs(λ)是波长λ处的实测光谱反射率，间隔为1 nm，φi(λ)是波段i在波长λ处的光谱响应函数，来自中国资源卫星应用中心网站。计算的遥感反射率和模拟的等效遥感反射率如图2所示。

获取了与野外实验同步的2016年4月29日太湖地区WFV3传感器影像，卫星过境时间为北京时间11:26，云覆盖率为3%，影像质量良好。对该影像分别进行几何纠正、影像裁剪和辐射定标等预处理，作为大气校正的输入数据。以研究区地形图为基准，校正模型为二次多项式，内插方法为最邻近内插法，纠正后的误差控制在1个像元范围以内。对几何校正后的影像进行裁剪，得到覆盖太湖水体的影像。

根据绝对定标系数和辐射定标公式，对GF-1/WFV图像进行定标，将卫星载荷观测值DN转化成辐亮度L。WFV数据的辐射定标公式为：L(λ)= Gain×DN+ Offset，其中L(λ)为转换后辐亮度，Gain 为波段增益，DN 为卫星载荷观测值，Offset为偏移量。WFV相机辐射定标参数，包括4 个波段的增益、偏移量参数及定标公式，均从中国资源卫星应用中心网站获取。

对定标后的辐亮度数据进行表观反射率的计算，计算公式如下：

其中:ρTOA为表观反射率；L为辐亮度；d为日地相对距离，其值在1左右，随日期略微变动，取日地平均距离为1；θs为太阳天顶角；E_λ为波段平均太阳辐射值，可根据GF-1/WFV传感器的光谱响应函数和对应区间的太阳光谱函数来计算，本文采用刘佳等^[15]的计算结果。

为了评价大气校正的效果，使用均方根误差（RMSE，Root Mean Square Error）和平均相对误差（ARE，Average Relative Error）作为评判标准来检验不同大气校正方法得到的结果值与ASD实测的遥感反射率是否一致。计算方法如下：

其中:y是ASD实测的遥感反射率，y’是遥感图像经过大气校正后得到的遥感反射率，n是样点数目。

大气组分中除水汽（H₂O）和臭氧（O₃）外一般变化不大。本文对比了两种标准大气模式：中纬度夏季（H₂O和O₃含量分别为2.93 g/cm²和0.319 cm-atm）、中纬度冬季（H₂O和O₃为0.853 g/cm²、0.395 cm-atm），以及3种自定义的大气模式：包括2016年4月29日的全球地基气溶胶监测网AERONET太湖站(图1中绿色点)的观测结果A1（H₂O和O₃为1.75 g/cm²、0.317 cm-atm），已有文献中使用的大气模式A2(H₂O和O₃为1.54 g/cm²、0.34 cm-atm)^[10]和A3(H₂O和O₃为2.93 g/cm²、0.293 cm-atm)^[11]。

使用中心像元的观测几何、太湖地区常用的大陆型气溶胶^[11,17]类型，2016年4月29日AERONET太湖站观测的气溶胶浓度。由于AERONET上无法提供550 nm处的AOD，因此使用AERONET提供的440、500、675和870 nm 4个通道处的AOD值，进行二项式拟合^[20]，插值得到550 nm处的AOD，计算结果为1.002 3。

对比5种不同的水汽和臭氧含量的大气模式对大气校正的结果影响，计算结果表明，不同的H₂O和O₃含量对蓝、绿、红波段的校正结果影响不大，校正后的遥感反射率的极差在0.07%以内；对近红外波段的校正影响略大，校正后的遥感反射率最大偏差为0.11%。由于大气模式对大气校正的结果影响有限，本研究使用中纬度夏季的大气模式进行后续分析。

6S中预定义了7种缺省的气溶胶模式，均由4种基本类型的气溶胶粒子（沙尘类、海洋类、水溶类、煤烟类）按不同比例混合组合。本文对比分析了3种标准气溶胶模式，包括大陆型（沙尘性70%，水溶性29%，海洋性0，煤烟性1%）、城市型（沙尘性17%，水溶性61%，海洋性0，煤烟性22%）、海洋型（沙尘性0，水溶性5%，海洋性95%，煤烟性0），以及已有文献中春夏季太湖使用的3种不同的气溶胶模式^[10,21]，包括M1（0、50%、40%、10%）、M2（10%、56%、10%、24%）、M3（48%、50%、0、2%）。

基于中心像元卫星观测几何，标准的中纬度夏季大气模式，AERONET插值获得的AOD（1.002 3），使用6种不同粒子组成的气溶胶模式对WFV数据进行大气校正。任意选择太湖水体中3个样点，对比不同气溶胶模式校正后的遥感反射率如图3所示。

由图3看出，不同气溶胶模式校正后的遥感反射率相差较大，说明不同的粒子比例组合对大气校正的结果影响较大，这与已有的研究结论一致^[16,17]。从光谱形态来看，海洋型、M3和大陆型校正后的光谱曲线与二类水体光谱相似：蓝波段由于水体吸收值较低，且蓝绿波段反射率低于近红波段。城市型、M1和M2型对蓝波段校正结果较差。与实测光谱（图2）对比来看，海洋型气溶胶校正后的各波段反射率值最接近实测光谱，B1波段反射率在0.02~0.03附近,且B4波段反射率低于0.02，显著优于其他几种模式的校正结果。

使用已有文献中常用的大陆型气溶胶^[11,17]，中纬度夏季大气模式，中心像元的观测几何，分析不同的AOD对大气校正结果的影响。AOD范围从0.5~1.2以0.1递增，不同AOD值校正后任意3个样点的遥感反射率如图4所示。

由图4看出，不同AOD的变化对波长较短的波段影响较大，校正后反射率在蓝波段和绿波段的反射率最大差异达 1.06% 和0.64 %，而在近红外波段的反射率差异仅为0.09%。使用不同AOD值进行大气校正得到遥感反射率差异较大，对于面积较大的二类湖泊水体，必须考虑AOD空间分布的不均匀性。

对于整景WFV影像和裁剪后的太湖影像，分别获取4个角点（左上、右上、左下、右下）的观测几何（表1）。由表看出，对于原始影像，太阳天顶角和方位角极差为2.96°和10.29°，卫星天顶角极差达到17.63°。对于裁剪后的太湖区域，太阳天顶角和方位角极差为0.84°和3.39°，卫星天顶角极差为2.64°。

假定大气模式为中纬度夏季、大陆型气溶胶，AOD为1.002 3，使用表1中4个角点以及影像中心像元的观测几何，分别进行6S大气校正。遥感反射率计算结果表明，使用中心像元观测几何替代各像元观测几何进行校正，会导致4个点在蓝、绿、红和近红波段分别出现0.14%、0.13%、0.10%、0.24%的最大误差。对于宽视场的WFV数据，观测几何的差异会产生校正误差，逐像元计算观测几何可以改进大气校正的效果。对于太湖4个角点，使用中心像元观测几何替代样点像元的观测几何，6S校正出的遥感反射率的最大偏差为0.12%。

以上分析表明，6S辐射传输模型的不同参数设置对大气校正结果影响不同，其中大气模式中水汽和臭氧含量的变化，对校正后的水体遥感反射率影响不大。气溶胶模式对大气校正影响最为明显，海洋型气溶胶校正结果更接近实测光谱。不同AOD值对短波反射率影响较大。观测几何对幅宽较大的WFV影像大气校正结果存在一定影响。

假设陆地表面为均匀朗伯体、大气垂直均匀变化的条件下，卫星传感器所接受的大气层顶表观反射率可表示为:

其中:ρTOA是大气层顶表观反射率；ρ0是大气的路径辐射项等效反射率；ρs为地表二向反射率；μv=cos(θv)，μs=cos(θs)，θv和θs分别为观测天顶角和太阳天顶角；ϕ为太阳与传感器之间的相对方位角；S为大气下界的半球反射率；T(μs)和T(μv)分别为向下和向上整层大气透过率。

6S模型的输入参数包括卫星几何参数、大气模式、气溶胶模式、光谱条件和地表反射率。输出参数为校正系数xa，xb和xc，再使用公式(6)就可以计算地表反射率。

由于WFV是宽幅相机，其幅宽为800 km，图像边缘像元与中心点像元的观测几何角度存在较大差异，因此需要计算每个像元的观测角度。整景影像中心像元的太阳天顶角和方位角、卫星天顶角和方位角从卫星自带参数中查找（其中xml文件中提供的天顶角其实为高度角）。

由图像4个角点的经纬度进行插值处理得到每个像元的经纬度，再结合卫星过境时间，按下式逐像元计算太阳天顶角θ和太阳方位角A^[22]：

其中:φ为像元的地理纬度，δ表示太阳赤纬角，t表示太阳时角。

根据相机视场角插值获得观测天顶角。WFV影像的xml文件中提供了整景影像中心像元的传感器的扫描角度，将其插值到整幅影像的不同位置。根据下式中传感器扫描角与观测天顶角的相互关系^[23]，逐像元计算观测天顶角的值。

其中:扫描角为α，观测天顶角为θ，R0为地球半径，h为卫星高度。由于卫星方位角对定量应用影响不大，直接使用中心像元的卫星方位角值。

气溶胶来源复杂且其类型处于不断地变化中，不同的气溶胶类型对辐射的影响不同^[24]。由于太湖面积较大，其不同湖区气溶胶类型差异较大^[11,25]，目前尚未有适用的气溶胶类型的确定方法。本文参照金鑫等^[26]的思路对太湖进行分区的气溶胶模型建立，即将太湖按照辐亮度的差异分为北太湖和其他湖区两个部分。其中北太湖受到沿岸陆地影响，考虑为大陆型气溶胶，其他湖区可考虑为海洋型气溶胶。

参考已有文献中的方法^[27]，使用分区暗像元和Spline 插值的方法计算太湖各个区域的AOD。传统暗像元大气校正算法认为研究区域上空的AOD呈均匀分布状态，但对于气溶胶类型复杂的二类湖泊水体, 暗像元算法的均匀性假设将不再适用。本研究将太湖划分为9 个子区域, 每个子区域内根据近红外波段最小值确定暗像元，再计算该像元的AOD，最后使用Spline 插值方法获取整个太湖的AOD。图5给出了太湖地区的9个子区域划分，查找出的各区暗像元的位置和计算出的AOD值，以及Spline插值后太湖全区的AOD。

基于分区暗像元和Spline插值方法确定的太湖AOD最大值为0.935 9,最小值为0.570 8,平均值为0.690 9，标准差为0.027 8，其数值范围与太湖多年平均气溶胶光学厚度（AOD）保持一致^[28,29]。图5体现了太湖AOD的空间异质性：北部竺山湾和贡湖湾AOD值最大，东太湖AOD值最小，整个区域AOD分布差异较大，使用统一的AOD值进行大气校正会造成一定误差。

目前常用的AOD获取方法有MODIS气溶胶产品、AERONET观测和暗像元反演等。其中MODIS气溶胶产品分辨率较低、容易受到云的影响和出现空值。2016年4月29日AERONET太湖站（图1）观测数据计算出的550 nm处的AOD为1.002 3，代表站点所在的北部湖区的AOD值。使用未分区的暗像元方法反演出的太湖AOD为0.614 7，其结果的真实性低于分区暗像元的反演结果。

为了提高计算效率，首先利用6S模型计算不同观测几何（太阳天顶角、太阳方位角和卫星天顶角）、气溶胶类型和AOD情况下的校正系数xa,xb和xc，建立针对这5个因素的查找表。设定6S模型的基本参数：几何条件为自定义，时间为4月29日；大气模式为中纬度夏季；气溶胶模式为用自定义；目标物高为-0；卫星观测高度为-645；使用GF-1/WFV的传感器波段范围450~520 nm、520~600 nm、630~690 nm、760~900 nm。假设地表均一且为朗伯体，无方向效应，选择6S自带的湖水类型作为地表模型。

变化参数的取值范围为：太阳天顶角从0°~80°以10°的步长变化，共有9个值；太阳方位角从0°~360°，以20°的步长变化，共19个值；卫星天顶角从0°~80°以10°的步长变化，共有9个值；卫星方位角为头文件中自带的282.283度。气溶胶类型分别取值为1（大陆型）、2（海洋型）、3（城市型）；AOD从0.5~1.2以0.1的步长变化，共8个值。利用IDL将上述参数组合代入6S模型循环计算，得到GF-1/WFV 4个波段在不同太阳天顶角、太阳方位角、卫星方位角、气溶胶类型和AOD下的大气校正系数xa,xb和xc的查找表。输出内容依次为太阳天顶角、太阳方位角、卫星方位角、气溶胶类型、AOD以及对应的校正参数。°′″

根据裁剪后影像的观测角度范围，缩小查找表范围，即查询低于影像输入参数最小值和高于最大值的第一条记录。逐像元计算的观测角度、气溶胶类型和AOD，对查找表进行插值^[30]，计算当前像元对应波段的大气校正系数xa,xb和xc。根据6S反射率计算公式，计算该像元的地表反射率。逐像元6S大气校正的过程如图6所示。

使用中纬度夏季大气模式，逐像元计算观测几何、气溶胶类型和AOD，对太湖水体进行6S大气校正。校正前后的太湖WFV影像（真彩色合成图）以及水体光谱曲线（影像中红框处的点）如图7所示。由校正前后的影像图看出，经过大气校正之后反射率的遥感影像图像更清晰。由光谱曲线看出，与TOA相比，校正后反射率值普遍降低，尤其降低了蓝绿波段的大气影响。校正后的遥感反射率曲线体现了湖泊水体的特征：绿波段（Band2）的反射率较高，而在蓝波段（Band1）和近红外波段（Band4）由于水体吸收的影响反射率较低。含藻水体在红波段（Band3）的0.67 μm附近有一个窄的吸收峰^[16]，在0.7 μm附近有一个荧光峰。4个波段经过大气校正之后，反射率显著降低，说明大气校正减弱了空气中水汽、臭氧和气溶胶等粒子对水体信号的干扰，恢复了水体本身的信号。

使用ENVI 5.0软件的FLAASH模块进行大气校正，首先将定标后的大气顶部反射率转化成BIL格式，使用比例系数10将输入光谱数据的单位转换成μW/(cm²×sr×nm)。FLAASH大气校正的参数为：软件自动判别影像的中心经度和中心纬度，成像日期和时间由传感器自带的文件中获取；设定传感器高度为645 km，地面海拔为0 km，像元大小为16 m，大气模型为中纬度夏季。气溶胶类型也使

用海洋型，AOD使用AERONET插值后获取的1.002 3，并参照文献[31]中的公式转换成能见度，结果为8.96 km。

比较17个样点的实测光谱、FLAASH校正结果和逐像元6S大气校正结果，结果如图8所示。由图看出，6S逐像元计算后精度与FLAASH的校正效果大体相当，其中6S逐像元校正后反射率与实测光谱更为接近。由于近红外波段水面遥感反射率较低，且该波段传感器信噪比较差，其大气校正精度显著低于其他几个波段。

基于实测光谱，计算6S逐像元和FLAASH大气校正后所有样点在各个波段的均方根误差和平均相对误差，结果如表2所示。由表看出，6S逐像元校正后所有波段的校正精度均优于FLAASH方法，与FLAASH方法相比，4个波段的ARE分别降低了1.84%、7.78%、4.79%以及17%，所有波段的ARE平均值降低了7.85%；RMSE的结果也均优于FLAASH，降低了0.01以内。结果表明，逐像元6S校正方法显著优于基于统一参数输入的FLAASH校正结果。

两种校正结果的相对误差如图9所示。结果表明, 使用统一参数输入的FLAASH校正得到的遥感反射率与实测的遥感反射率之间的平均相对误差大部分大于经过分区和逐像元计算后的6S校正得到的遥感反射率。可见, 基于不同湖区的气溶胶类型分区和逐像元的大气校正能够提高大气校正的精度。

太湖上空气溶胶组分复杂，较常使用的是大陆型气溶胶^[11]、城市型和大陆型的气溶胶^[17]，或者城市型和工业型气溶胶^[32]，但还没有确定的气溶胶类型。Higurashi等^[33]使用0.8的Angstrom指数α分界线来区分气溶胶类型，2016年4月29日AERONET观测的α系数为1.25，对应的是吸收性的“黑炭气溶胶”或者非吸收性的“硫酸盐气溶胶”。这是因为太湖受到长三角地区工业城市排放的大气污染物（如非吸收性的硫酸盐粒子）的影响^[34]，且快到夏季从海洋上空输送的非吸收性硫酸盐粒子至太湖上空。虽然可以依据气溶胶指数α来判定气溶胶的类型，但不同类型气溶胶粒子的AOD 值在大小上并不是严格区分的，即使是同一类型下的气溶胶，其AOD和气溶胶指数α也存在较大的偏差，这给判断气溶胶的组分带来困难。不同类型粒子比只是反映气溶胶类型的一个参数，还需要结合气溶胶的各种光学特性参数，辅助以其他分析手段进行综合处理^[24]。

对于6S模型，当不考虑气溶胶分区时，使用海洋型气溶胶的校正结果最好（图3）。对于FLAASH方法，假设AOD为1.002 3（参照文献[31]中的公式转换成能见度8.96 km），分别使用城市型或海洋型气溶胶类型进行大气校正，与实测光谱比较表明，海洋型气溶胶的结果优于城市型。一般认为，空间动力学直径>1 μm的海盐气溶胶粒子随海洋气团进入陆地的过程中，在重力和湍流的作用下，其浓度会以指数形式衰减，从而迅速从大气中去除，其空间尺度大约为10 km^[35]。Willison等^[36]研究指出，海洋气溶胶对离海岸100 km以外的大气气溶胶粗、细粒子的质量贡献分别为10%~20%和2%~3%。而本文的结果表明，基于AERONET观测值的大气模式和气溶胶浓度条件，海洋型气溶胶更适用于春季太湖GF-1/WFV影像的大气校正。这可能是由于4月底5月初太湖地区降雨量逐渐增多，相对湿度的增加，大气中水溶性的细粒子增加^[37]；另外受到季风气候影响，气团从海洋向大陆运动，从海洋上空带来的大量的细粒子气溶胶^[38]。也有文献表明海洋型气溶胶适用于沿海的太湖地区^[39]，本文采用海洋型气溶胶用于6S大气校正，结论与之一致。

对于太湖地区大气校正，已有的研究中大多针对的是中低空间分辨率的遥感影像，如MODIS^[40]、TM^[11]和HJ1^[26]等，且这些研究较少考虑不同气溶胶粒子比对大气校正结果的影响。随着高空间分辨率遥感影像应用的普及，使用高分影像进行内陆水体的水色遥感应用将日益广泛。目前针对高分影像水色遥感产生了大量的大气校正算法^{[12,13,14,15]},但这些算法在水体区域的应用效果讨论较少。本文对太湖地区GF1影像的大气校正结果表明，海洋型气溶胶可以用于春季太湖数据，6S逐像元模型可以改进校正精度。本文研究对于太湖水体的大气校正可以提供一定的应用参考。

使用太湖地区的GF-1 WFV影像和同步实测光谱数据，基于6S大气校正参数的敏感性分析结果，逐像元计算了观测几何，使用分区气溶胶类型，使用分区暗像元和Spline插值计算AOD，提出了一种6S逐像元大气校正方法。实验结果表明，气溶胶模式对大气校正的影响最大，逐像元计算观测几何和气溶胶参数大气校正精度有所改进。与FLAASH方法相比，逐像元的6S模型校正后4个波段的平均相对误差分别降低了1.84%、7.78%、4.79%和17%，4个波段的平均值降低了7.85%。本文提出的6S逐像元模型对于水体信息的大气校正精度有所改进。

下一步需要使用多期星地同步数据对该校正方法进行进一步验证，构建针对太湖地区的高精度大气校正算法。同时，还需进一步扩展本算法的应用范围，针对不同的湖泊、河流和水库等水体信息选取适宜的大气参数，构建校正参数查找表，改进6S辐射传输模型在水色遥感大气校正中的适用性。