1. 引言

草型湖泊具有多种生态系统功能，对维持生态平衡、保护物种多样性和维持农业发展具有重要作用 [1] [2] 。在草型湖泊中，植物群落结构复杂 [3] ，水域和沿岸地带生长水生植物，包括沉水植物、挺水植物、浮叶植物和藻类等 [4] 。

由于现场和卫星获取到的遥感反射率(Remote sensing reflectance, Rrs)数据存在很多不确定性。因此，在草型湖泊中进行遥感反射率产品的质量控制，获取准确的 $R_{rs}$尤为重要 [5] 。Wei等(2016)提出了一个Rrs质量控制系统，可以准确识别出“好”与“不好”Rrs光谱。该系统可应用于现场测量和MODIS等卫星获取的Rrs产品的质量控制 [6] 。Qing等(2022)基于全球海域的光谱漫衰减系数，提出了光学分类与质量控制模型(OC_QC模型)，并将模型应用于(Medium Resolution Imaging Spectrometer Instrument, MERIS)和(Moderateresolution Imaging Spectroradiometer, MODIS)遥感数据 [7] 。质量控制系统具有潜在的应用价值，然而目前还没有对内陆湖泊水体遥感反射率进行相关的质量控制研究。

国际上目前已有很多基于实测和卫星遥感反射率光谱的光学分类研究 [8] 。例如Wei等(2022)利用k-means聚类将海洋、沿海和内陆水域划分为23个类别，并检验了可见光红外成像辐射计(VIIRS)产品的可用性 [9] 。此外，还有一些研究使用了机器学习 [10] 、密度分割分类 [11] 、迭代聚类ISODATA [12] [13] 和分类树 [14] 等方法对水体进行分类。

因此，本文的主要目的是：(1) 构建草型湖泊实测遥感反射率水体的质量控制与光学分类模型；(2) 分析水体光学类型及特征。

2. 方法

2.1. 水体光学类型库构建

本研究参考了wei等 [6] 的研究方法：(1) 首先对遥感反射率光谱数据集进行归一化；(2) 利用间隙法确定归一化遥感反射率光谱nRrs数据的聚类数k；(3) 其次使用kmeans聚类法进行聚类； (4) 确定乌梁素海水体类型，最后提取每种水体类型nRrs光谱在离散波段(13个波段)上的平均值、最大值、最小值作为参考光谱的质心及上下边界值，构成质量控制与光学分类的核心。

2.2. R_rs光谱质量控制

通过比较目标光谱 $R_{rs}^{*}$与对应水类型的上、下边界来计算QA得分。其中 $C{\lambda }_i$为波长特定得分，N为总波长数。

$C_{tot}=\frac{C\left({\lambda }_1\right)+C\left({\lambda }_2\right)+\cdots +C\left({\lambda }_N\right)}{N}$(1)

2.3. 水体光学分类

通过波长将目标光谱 $n{R}_{rs}^{*}$与参考光谱 $n{R}_{rs}^{ref}$进行匹配，为目标光谱分配水类型。其中 $a$为参考光谱与目标光谱形成的夹角。

$cosa=\frac{{\displaystyle {\sum }_{i=1}^N\left[n{R}_{rs}^{*}\left({\lambda }_i\right)\cdot n{R}_{rs}^{ref}\left({\lambda }_i\right)\right]}}{\sqrt{{{\displaystyle {\sum }_{i=1}^N\left[n{R}_{rs}^{*}\left({\lambda }_i\right)\right]}}^2{\displaystyle {\sum }_{i=1}^N{\left[n{R}_{rs}^{ref}\left({\lambda }_i\right)\right]}^2}}}$(2)

3. 数据

乌梁素海(108˚43′~108˚57′, 40˚36′~41˚03′)位于中国内蒙古自治区巴彦淖尔市，是世界上为数不多的干旱地区草型湖泊之一，具有生物多样性和生态功能，对维护黄河中上游重要的涵养、储水和调水 [15] 起着重要的作用(图1)。本研究在2015~2021年乌梁素海水域开展了14次野外调查，共获取了595组遥感反射率光谱(Rrs)数据(图2)以及每个站点所对应的水体透明度和水草高度数据，本文从高光谱库中提取了13个波长(412、443、490、510、530、560、620、665、681、705、740、783、865 nm)的遥感反射率数据，这些波长大都可以在卫星传感器中获取。

3.1. 遥感反射率数据

本研究中使用地物光谱仪ASD Field Spec 4，采集的遥感反射率波段范围在400 nm~900 nm之间，采样间隔为1.4 nm，光谱测量采用水面以上方法 [16] 。遥感反射率(R_rs)计算公式如下：

$R_{rs}=\left[L_{sw}-r\ast L_{SK}\left(\lambda \right)\right]{\rho }_p\left(\lambda \right)/\pi \ast L_p\left(\lambda \right)$(3)

式中， $L_{sw}$为测量的水体辐亮度， $L_{SK}\left(\lambda \right)$为天空辐亮度， $L_p\left(\lambda \right)$为标准版辐亮度，r为海气界面反射率， ${\rho }_p\left(\lambda \right)$为标准版的反射率。

3.2. 透明度及水生植被高度数据

水体透明度数据通过直径为30 cm的黑白塞氏圆盘放入水中直到磁盘不再可见来确定。水草高度由铅球标尺读取，由水平面至水草高度间的距离确定水草高度数值。水体遥感反射率在不同的透明度和不同淹没水生植被高度的影响下有不同的光谱响应 [17] 。

4. 结果

4.1. 构建水体光学类型库

对 $n{R}_{rs}$采用间隙法确定最优聚类数，聚类簇的数量k最初设置为20，间隙曲线随着簇k数量的增加而递增，但从某个k开始上升趋势明显平缓(图3(a))，在簇数接近9的位置表明了有适当的集群数量 [18] ，因此确定最优聚类数为9。利用k-means非监督算法对 $n{R}_{rs}$光谱进行聚类，得到9种水体光学类型和9种平均归一化中心光谱(图3(b))。

4.2. 水体光学类型及特征

根据聚类后的 $n{R}_{rs}$结果进行监督分类，调整个别误分类结果，发现尽管这些光学水体光谱很复杂，但水体透明度以及水生植被高度对各类水体光谱有显著贡献(图4)，最终确定了9种光学水体类型数据集(图5)，1~2类为纯水体类型、3~8类为沉水植被类型、第9类为挺水类型，随着透明度和水生植被高度不同而呈现出不同的光谱特征。最后提取了每种聚类在13个波段上的平均值、最大值、最小值作为参考光谱的中心光谱、上下边界，为质量控制和水体光学分类模型提供参照依据，9种水体类型的平均归一化遥感反射率光谱特征(图6)，各类水体光谱间存在显著差异，表1描述了水体类型主要特点。

5. 讨论

基于本文的质量控制与光学分类建模使用了实测数据13个波段(412、443、490、510、530、560、620、665、681、705、740、783、865 nm)的遥感反射率值，包括大多数传感器的波段，模型是否适用于遥感应用还需评估波段完整性，考虑不同传感器的波段完整性是否影响质量控制与分类还需进一步分析，例如Sentinel-2 MSI和Landsat-8 OLI影像的波段与模型的拟合性。本文使用的光学分类方法与传统的水色指数方法(Forel-Ule Index, FUI)和Spyakos等人 [19] 的内陆水体分类方法有很大差异，FUI主要根据水体颜色及水质进行分类，由于乌梁素海为草型湖泊，因此不适用于本研究区，而内陆水体光学分类方法根据水的吸收系数和生物光学特性进行了描述。本文的研究方法更适用于复杂的草型湖泊遥感反射率光谱中。

6. 结论

本研究通过乌梁素海草型湖泊全面获取的现场测量遥感反射率光谱R_rs数据集构建了质量控制与光学分类方法，期望获取准确的遥感反射率数据，与此同时更好地了解草型湖泊水体光学类型变化，并为未来的研究提供一个框架，得到如下重要结论：

(1) 对于质量控制，如人为造成的原位测量误差、影像质量的误差等有问题的数据，质量控制系统都可以很好地识别出来。

(2) 草型湖泊的不同水体生物光学参数对光谱类型都有很大影响。如水体透明度、水草高度的不同都会影响水体反射率的呈现，这为光学分类提供了重要依据。

(3) 基于实测数据研究存在成本高、时效性差等问题，我们的遥感反射率质量控制和光学分类方法还可以应用于遥感当中，需要注意的是大气校正有关的不确定性对反射信号的影响。

基金项目

内蒙古师范大学科研创新项目“复杂内陆湖泊水体遥感反射率光学分类与质量控制(CXJJS22129)”。

参考文献