杜雨春子,青 松,曹萌萌,袁瑞强,顺布日,郝艳玲

(1：内蒙古师范大学地理科学学院,呼和浩特 010022)

(2：内蒙古大学生态与环境学院,呼和浩特 010021)

沉水植物在浅水湖泊湿地发挥着不可忽视的生态功能，如吸收氮磷等富营养化物质，抑制浮游植物生长，为无脊椎动物、鱼、水禽提供栖息地，为浮游动物提供庇护所，促进底层沉积物的淤积等. 因此实时监测湖泊沉水植物，对于水资源调查研究、开发保护、水环境生态修复具有重要意义[1-4].

传统监测沉水植物分布与生长状况的现场调查法和潜水测量法耗时耗力、监测覆盖范围小，而且影响水体环境[4]. 遥感技术能够大面积同步获取地物信息，所以在水生植物时空分布监测方面具有传统方法不可比拟的优势[5]. 在多光谱应用方面，许多学者利用GF，Landsat系列(包括 MSS、TM、ETM+)，SPOT和HJ-CCD等遥感数据对水生植物进行监测、初步分类，提取和制图，并获得了很好的研究结果[6-9]；其中，Dogan等[6]利用Quickbird卫星数据绘制了浅湖沉水植物覆盖范围，并指出冠层水深是影响大规模遥感监测沉水植物覆盖度和丰度的主要因素. Rotta等[10]探讨了利用SPOT-6估算内陆水域沉水植物高度(SAV)的可行性，并指出在深度大于5 m，高度为1.5 m的情况下，沉水植物反射率的灵敏度低. 近几年来，随着高光谱遥感技术的应用与发展，使遥感由定性分析向定量或半定量的转化成为可能，高光谱遥感具有空间分辨率和光谱分辨率高的特点[11]. 目前，国内外学者使用实测高光谱数据对内陆湖泊的沉水植物进行种类识别、分类、覆盖度、生物量评估及沉水植物结构特征与光谱信息关系等方面的定量研究[12-13]. Pande-Chhetri等[14]利用机载CASI高光谱传感器，对浅水湖泊(深度<3 m)沉水植物进行了遥感分类. Zou等[15]在富营养化湖中，利用实测高光谱数据研究了覆盖度与沉水植物光谱特征之间的关系，并指出沉水植物光谱反射率随覆盖度的减少而下降，水深是影响沉水植物覆盖度的主要因素. 由此可见水深这一重要因素一直以来直接影响着沉水植物的遥感监测[16-19]. 由于水体环境的复杂性，目前对沉水植物光谱特征及沉水植物群落冠层在水下深度(冠层水深)对沉水植物光谱影响的定量研究还需进一步深入探讨.

本文以内蒙古巴彦淖尔市境内乌梁素海沉水植物群落为研究对象，基于地物光谱仪测定的沉水植物光谱数据，开展如下研究：1) 分析沉水植物与其他地类光谱特征差异; 2) 探讨冠层水深对沉水植物群落光谱特征的影响; 3) 建立冠层水深遥感反演模型；4)利用Landsat-8 OLI数据反演得到乌梁素海冠层水深分布图.

1 数据和方法

1.1 研究区概况

乌梁素海(40°36′～41°03′N，108°43′～108°57′E)(图1)位于内蒙古巴彦淖尔市境内，是中国的八大淡水湖之一，也是河套灌区唯一的排涝湖. 湖泊呈南北长、东西窄的狭长形态，湖岸线长130 km. 湖泊蓄水量2.5亿～3亿 m3，水深为0.5～3.0 m，平均水深1.5 m，水域面积273.32 km2. 所在地区的多年平均气温为7.3℃，多年平均降雨量为224 mm，蒸发量为1502 mm，全年无霜期为152 d，湖水于每年11月初结冰，次年4月份解冻，冰封期为5个月. 乌梁素海还是横跨欧亚大陆鸟类栖息和迁移的重要场所，是黄河中上游重要的保水、蓄水和调水基地，也是全球范围内荒漠半荒漠化地区极为少见的具有很高生态效益的大型多功能草型湖泊，是地球上同一纬度最大的自然湿地，在中国湿地、荒漠及动物物种3大生态系统保护中均居于重要地位. 由于长期受农田退水、工业废水等的影响，湖水被严重污染，湖区内挺水植物、沉水植物等大型水生植物疯长，导致乌梁素海成为重度富营养化草型湖泊[20-23].

乌梁素海水生植被生长茂盛，在全湖范围内广泛分布，植被种类大致可以分为挺水植被、沉水植被和漂浮藻类3类. 其中，挺水植被主要有芦苇、宽叶香蒲，芦苇为优势种，4月开始发芽出青，9月份开始衰退枯黄，挺水植物占湖面的1/2，主要分布在湖的中部、西岸和北部，芦苇外围生长着少量香蒲，面积不大；沉水植被有龙须眼子菜、菹草、狐尾藻、茨藻、轮藻、竹叶眼子菜等，龙须眼子菜和穗花狐尾藻为沉水植被优势种，5月中旬开始萌芽，10月开始逐渐衰退，主要分布在明水区，但是南部明水区分布较少. 漂浮藻类主要为黄苔，黄苔由水绵、双星藻、转板藻3个属的藻类组成，其丝状绿藻最初生长于水下，缠绕在沉水植被上、或依附在底泥上生长，当聚集到一定程度后会形成团块，随着光合作用产生的气泡飘浮并出露水面，受到强光照射部分死亡而呈现出黄色，被称为“黄苔”，5月初“黄苔”开始出现，7-8月份“黄苔”大面积暴发，主要分布在中部和东部[24-26].

图1 乌梁素海地理位置及实测站位

1.2 野外实测

本研究分别于2018年7月16日、8月17日、9月2日、9月18日、11月5日和2019年6月17日在乌梁素海开展了6次为期2日的野外调查，且时间与Landsat-8卫星过境时间同步，共获得有效数据378组，包括挺水植物(30组)、沉水植物(228组)、漂浮藻类(8组)和开阔水体(112组)的遥感反射率光谱、冠层水深/水深和现场照片(图1和图2).

光谱测量采用美国 ASD 公司生产的 ASD Field Spec 4 地物光谱仪，其测量波段范围为350～2500 nm，光谱分辨率为3 nm，光纤视场为25°，350～1000 nm内采样间隔为1.4 nm，1001～2500 nm 内采样间隔为1.1 nm. 光谱测量方法采用唐军武等提出的水面以上水体光谱测量方法[27]. 遥感反射率测量在无风、天气晴朗时进行，且每个样点需要测量标准灰板、水体和天空光的辐射亮度，剔除异常数据之后，取均值作为该点的光谱数据. 为了避开水体对太阳直射反射和船体阴影等影响，需要采用一定的观测角度，仪器观测平面与太阳入射平面夹角为90～145°之间，且仪器与水体之间的距离为1 m. 用卷尺现场测量站点沉水植物群落从水面到植物冠层的深度，即冠层水深(WDC). 冠层水深呈正态分布，其平均值和标准差分别为42.73和25.00 cm. 用塞氏盘测量水体透明度，利用手持GPS定位仪记录每个站点的位置.

图2 挺水植物、沉水植物、漂浮藻类、开阔水体照片

1.3 遥感影像数据及其预处理

本文选用2018年9月2日和2019年6月17日的无云覆盖的2景Landsat-8 OLI影像数据(http://earthexplorer.usgs.gov/). Landsat-8卫星发射于2013年2月11日，重访周期为16 d，搭载2个传感器，OLI(Operational Land Image)和TIRS (Thermal Infrared Sensor). OLI有9个波段，其中全色波段的分辨率为15 m，其余8个波段的分辨率为30 m[28]，9个波段的中心波长分别为443 nm(Coastal)、483 nm(Blue)，561 nm(Green)、655 nm(Red)、865 nm(NIR)、1609 nm(SWIR1)、2210 nm(SWIR2)、592 nm(Pan)和1373 nm(Cirrus). 由于乌梁素海遍布挺水植物、沉水植物、漂浮藻类，水环境复杂，精确的气溶胶校正存在困难. 因此OLI影像只做了大气瑞利校正，在此基础上进行WDC的遥感反演. 参考Quinten Vanhellemont & Kevin Ruddick[28-29]进行瑞利校正，步骤如下：

首先，利用公式(1)进行遥感影像的辐射定标：

LTOA=ML·DN+AL

(1)

式中，LTOA为辐射亮度(W/(m2·μm·sr))，DN值为遥感影像像元亮度值，ML(multiplicative factor，gain)为增益系数(W/(m2·μm·sr))，AL(additive factor，offset)为偏移量(W/(m2·μm·sr))，分别可从OLI数据头文件中获取.

其次，获取大气层顶反射率：

(2)

式中，ρTOA为大气层顶反射率(sr-1);d为以天文单位表示的日地距离(AU)；F0为波段平均的大气层外太阳辐照度(W/(m2·μm))；θ0为太阳高度角(°).

大气瑞利校正公式为：

(3)

式中，ρr是瑞利反射率(sr-1)；τr是瑞利光学厚度(无单位)；Pr是瑞利散射相位函数；θv是观测天顶角(°)；Δφ是相对太阳和传感器间高度角(°).

另外，本文试图利用瑞利校正后的OLI数据进行冠层水深的反演. 因此，利用OLI数据的波段响应函数，将实测光谱反射率积分到OLI波段上，波段积分运算见公式(4)：

(4)

式中，Rrs(λ)(sr-1)表示实测遥感反射率；Fi表示OLI数据的第i波段的光谱响应函数.

1.4 方法

1.4.1 光谱特征分析 基于实测遥感数据，分析沉水植物与漂浮藻类和挺水植物之间的光谱差异，辨识沉水植物与其他水生植物光谱区分特征；研究沉水植物光谱随冠层深度的变化特征.

1.4.2 相关分析 基于沉水植物光谱特征，进行冠层水深与沉水植物光谱之间的相关分析，确定相关性最大的单波段和波段比. 相关分析同时考虑了皮尔逊(线性)和斯皮尔曼(单调性)相关分析法.

1.4.3 回归模型建立 在1.4.1和1.4.2工作的基础上，建立单波段/波段比与冠层水深之间的回归模型，模型分别为：

(5)

(6)

1.4.4 精度评估 采用均方根误差(root mean square error，RMSE)，平均相对误差(mean relative percentage error，MRPE)以及Pearson相关系数评价反演模型精度，计算公式为：

(7)

(8)

式中，xipre和ximea分别表示通过模型计算得到的冠层水深和实际测量得到的冠层水深；n为采样点数.

2 结果

2.1 沉水植物与挺水植物和漂浮藻类光谱差异特征

图3为0深度沉水植物(WDC=0)、挺水植物和漂浮藻类实测光谱，由于1349～1449、1750～1979和2350～2500 nm段噪音大，因此剔除. 由图3A可见，沉水植物(WDC=0)和挺水植物、漂浮藻类的光谱曲线表现出典型的植被光谱特征，在可见光部分的蓝波段(480 nm为中心)和红光波段(680 nm为中心) 附近有较强的吸收，形成两个吸收谷，在540～560 nm附近形成了明显的反射绿峰；在670～720 nm之间反射率增高最快，出现“红边”；波长>700 nm时，其光谱反射率增强，形成了水生植物在近红外波段的高反射区.

三类植物光谱特征差异性表现在：挺水植物、沉水植物(WDC=0)和漂浮藻类在波段602、1272、1662和2223 nm光谱差异较大(图3B). 图4和表1中对比了挺水植物、沉水植物(WDC=0)和漂浮藻类在这4个特征波段处的取值范围. 可以看出，挺水植物反射率光谱，在近红外和短波红外波段范围内的1272、1662和2223 nm附近分别有一个明显的反射峰，这是挺水植物区别于沉水植物和漂浮藻类的重要波段. 沉水植物在近红外高反射率的最大值都在0.099以下，而挺水植物都在0.152以上；漂浮藻类的近红外高反射率在0.070～0.109之间，这与沉水植物的光谱波段有交叉，变化趋势比较一致，而且沉水植物和漂浮藻类在1350 nm之后反射率逐渐减少，均趋于0. 虽然沉水植物和漂浮藻类的波谱形状相似，但是在550～690 nm漂浮藻类的反射率值却高于沉水植物(WDC=0)的反射率值，尤其在602 nm处漂浮藻类的反射率值最低为0.046，而沉水植物的反射率值最高不到 0.021. 所以在550～690 nm漂浮藻类和沉水植物(WDC=0)的光谱曲线具有明显的分层现象，光谱特征差异明显.

图3 表层沉水植物、漂浮藻类、挺水植物实测光谱(a)和平均光谱(b)对比

图4不同特征波段沉水植物(WDC=0)、漂浮藻类和挺水植物光谱反射率箱型图：(a)为602 nm处；(b)为1272 nm处； (c)为1662 nm处；(d)为2223 nm处

表1 表层沉水植物、漂浮藻类和挺水植物实测光谱参数

图5 不同冠层水深沉水植物群落的实测光谱曲线

2.2 不同冠层水深的沉水植物光谱特征

图6 810 nm处不同冠层水深沉水植物和开阔水体光谱反射率箱型图

图5为不同冠层水深的沉水植物实测光谱曲线. 可以看出，光谱反射率随冠层水深增加而下降，在近红外区域700～900 nm尤为突出，这是因为水体通过吸收和散射造成的强烈光衰减对沉水植物光谱产生明显干扰[16]. 随着冠层水深的增加，沉水植物光谱在700～900 nm波段范围内出现以714 nm和810 nm为中心的2个反射峰，这主要是由于这2个波段之间正好是一个明显的水吸收带[31]. 当冠层离水面的距离为50 cm时水体与沉水植物群落光谱反射率非常接近，但在810 nm处存在细微的差异. 图6对比了810 nm处不同冠层水深沉水植物和开阔水体的光谱范围. 可以看出，随着冠层水深的增加，沉水植物反射率的集中区域逐渐降低变窄，并且逐步趋于0，而且51～160 cm深度范围的沉水植物光谱反射率与开阔水体的光谱反射率区域非常接近(图5e和图5f)，这表明冠层水深会影响沉水植物的光谱反射率.

2.3 冠层水深反演模型

2.3.1 冠层水深与沉水植物光谱反射率的相关分析 通过对不同冠层水深的沉水植物光谱特征分析，可以发现沉水植物群落的光谱反射率随冠层水深的增加而降低，在700～900 nm 波段范围内差异最为明显. 由于乌梁素海各月沉水植物物候特征随季节变化明显，所以分别对2018年7月、8月、9月、11月和2019年6月沉水植物群落的光谱反射率与冠层水深进行相关性分析(图7). 在700～735 nm波段附近，沉水植物群落光谱反射率与冠层水深呈显著负相关，是光谱反射率对冠层水深变化最敏感的波段. 这与Han等[32]的研究结果一致. 其中7月沉水植物群落的光谱反射率与冠层水深皮尔逊最高的相关系数位于715 nm处，r=-0.77，斯皮尔曼相关系数为 -0.91. 8月皮尔逊最高的相关系数位于707 nm处，r= -0.75，在760 nm处斯皮尔曼相关系数最高，ρ=-0.88. 9月皮尔逊最高的相关系数位于719 nm处，r= -0.64，在819 nm处斯皮尔曼相关系数最高，ρ=-0.75. 11月皮尔逊最高的相关系数位于967 nm处,r=-0.91，在954 nm处斯皮尔曼相关系数最高，ρ=-0.98. 2019年6月皮尔逊最高的相关系数位于715 nm处，r=-0.64，在733 nm处斯皮尔曼相关系数最高，ρ=-0.96. 所有月份的皮尔逊最高的相关系数位于720 nm处，r=-0.64，在733 nm处斯皮尔曼相关系数最高，ρ=-0.84. 从8月以后的相关性可以看出皮尔逊相关系数与斯皮尔曼相关系数所对应的特征波段差异很大，这是因为在707 nm之后沉水植物光谱反射率与冠层深度可能存在非线性关系. 而且6、7、8、9月为沉水植物群落生长旺季，沉水植物群落全部呈绿色，植株距水面近，近红外光谱反射率最高. 11月为沉水植物群落衰退期，沉水植物群落“失绿”呈黄褐色，植株距水面远，沉水植物群落近红外光谱反射率较低. 而在700～900 nm波段附近，是光谱反射率对冠层水深变化最敏感的波段，这可能也是影响沉水植物群落与冠层水深相关性的因素[33]. 通过迭代算法对400～1350 nm(5 nm间隔)波段区间的实测反射率逐一两两比值，与沉水植物群落冠层水深进行相关性计算，发现在表2波段比处r和ρ最高.

表2 各月波段比与冠层水深相关系数

图7 不同月份沉水植物光谱反射率与冠层水深之间的单波段相关系数(其中圆点和三角形分别代表皮尔逊和斯皮尔曼相关系数最优的波段位置)

表3 单波段/波段比冠层水深反演模型

图8 各月冠层水深实测值与反演值对比

图8为各月单波段/波段比冠层水深反演模型的实测冠层水深与反演冠层水深对比图，从图可以看出，波段比的反演精度要高于单波段，这是因为波段比值可在一定程度上消除不同时间和空间上水表面光滑度和大气对电磁波的影响，并部分地减少其他物质的干扰[34]，而单波段却不能. 在单波段的情况下，幂函数的反演精度要比线性函数高，幂函数的决定系数R2>0.5，均方根误差<16.18 cm，而线性函数的决定系数R2>0.4，均方根误差<18.87 cm，这是因为在700～735 nm波段附近，沉水植物群落光谱反射率与冠层水深呈显著负相关，而且在707 nm之前两者之间的线性关系较好，在707 nm之后的近红外波段反射率与冠层水深存在很好的非线性关系. 在波段比的情况下，幂函数与线性函数的决定系数R2非常接近，而且R2>0.70，均方根误差<13.70 cm. 7月波段比反演模型的线性函数和幂函数的决定系数分别为0.87和0.84，均方根误差分别为10.39和12.70 cm，平均相对误差分别为22.06%和22.41%；8月波段比的线性函数和幂函数的决定系数分别为0.79和0.76，均方根误差分别为9.23和10.07 cm，平均相对误差分别为36.99%和36.43%；9月波段比的线性函数和幂函数决定系数分别为0.72和0.70，均方根误差分别为13.36和13.33 cm，平均相对误差分别为28.75%和27.06%；11月波段比的线性函数和幂函数的决定系数分别为0.91和0.80，均方根误差分别为4.05和5.11 cm，平均相对误差分别为13.55%和14.75%；6月波段比的线性函数和幂函数的决定系数分别为0.95和0.81，均方根误差分别为6.59和12.40 cm，平均相对误差分别为15.59%和14.66%；以上所有月份波段比的的线性函数和幂函数的决定系数分别为 0.70 和0.71，均方根误差分别为13.58和13.70 cm，平均相对误差分别为31.42%和27.57%. 通过各月的RMSE、MRPE和散点图可以看出大部分反演值位于1∶1线附近，与实测值有较好的一致性，反演精度较高，单波段模型最好的冠层水深反演精度集中在10～30 cm，而波段比模型最好的冠层水深反演精度集中在10～60 cm.

沉水植物具有净化和改善水质的作用[17]，通过现场调查发现，在乌梁素海沉水植物分布区域水体透明度高，悬浮物浓度和叶绿素浓度低，对沉水植物光谱的影响小，而水面到植物冠层的深度对沉水植物光谱影响较大，所以本文重点探讨冠层水深对沉水植物光谱的影响.

2.3.3 模型推广 本文试图利用瑞利校正后的OLI数据进行冠层水深的反演. 首先，基于波段响应函数，将实测光谱反射率积分到Landsat-8 OLI波段上，见公式(4)；然后，进行了OLI单波段反射率/波段比值与实测冠层水深之间的相关性分析，见图9. OLI第5波段(中心波段为865 nm)的沉水植物光谱反射率与冠层水深的相关性最大，其皮尔逊相关系数为r=0.54，斯皮尔曼相关系数为ρ=0.73. 在波段比Rrs(Costal)/Rrs(NIR)处沉水植物光谱反射率与冠层水深的皮尔逊相关系数最大，r=0.67. 在波段比Rrs(Blue)/Rrs(NIR)处斯皮尔曼相关系数最大，ρ=0.76. 这是因为植物的叶绿素对蓝光的强吸收和对近红外波段(850～880 nm)高反射导致的[35-37]. 因此，将Rrs(Costal)/Rrs(NIR)和Rrs(Blue)/Rrs(NIR)作为自变量，利用LOOCV方法建立OLI波段比的冠层水深反演模型，其中以Rrs(Costal)/Rrs(NIR)为自变量的波段比线性函数模型为y=14.295x+17.534，R2=0.44，以Rrs(Blue)/Rrs(NIR)作为自变量的波段比幂函数模型为y=23.845x0.716，R2=0.49.

图9 Landsat-8 OLI单波段、波段比沉水植物光谱反射率与冠层水深的相关系数:(a)为单波段,(b)为波段比(其中圆点和三角形分别代表皮尔逊和斯皮尔曼相关系数最优的波段位置)

图10为OLI冠层水深反演值与实测值对比图. 由图10可知，OLI冠层水深反演结果与实测值之间有良好的一致性，其中波段比的幂函数反演结果最好. 波段比的线性反演模型的决定系数R2为0.44，均方根误差为18.60 cm，平均相对误差为63.67%，幂函数反演模型的决定系数R2为0.49，均方根误差为18.17 cm，平均相对误差为40.05%. 通过135个RMSE、MRPE和散点图10可以看出大部分反演的冠层水深位于1∶1线附近，表明本文使用的Landsat-8 OLI幂函数反演模型获得了较好的结果.

图10 Landsat-8 OLI冠层水深反演值与实测值对比

2.3.4 Landsat-8 OLI波段比模型反演乌梁素海冠层水深 将OLI波段比幂函数模型应用于瑞利校正后的OLI影像(2018年9月2日和2019年6月17日)，图11对比了OLI反演的冠层水深和实测值. 从散点图可以看出，反演冠层水深和实测冠层水深具有较好的一致性，R2为0.67. 然而，OLI数据只做了瑞利校正，未做气溶胶校正，因此反演冠层水深与实测值之间存在一定的偏差.

因此，以瑞利校正后Landsat-8 OLI的Band2/Band5反射率为自变量建立了冠层水深反演模型，模型为y=57.68x0.98，R2=0.70. 图12 中对比了大气瑞利校正后Landsat-8 OLI影像数据波段比模型反演得到的冠层水深和实测值. 由图12 可以看出冠层水深散点大部分都位于1∶1线附近，R2为0.70，RMSE为22.34 cm，MRPE为38.66%. 因此利用此模型得到了2018年9月2日和2019年6月17日乌梁素海沉水植物冠层水深分布图，见图13. 由图13可以看出，6月与9月沉水植物冠层水深差异明显，从整体来看，9月中部和东部的沉水植物冠层水深较小，而6月较大，这是因为9月沉水植物仍为生长旺季，植株距水面近，而6月沉水植物刚刚进入生长旺季，植株距水面较远. 从局部来看，9月北部沉水植物冠层水深较大，而6月却较小；中部9月沉水植物平均冠层水深较小，而6月平均冠层水深较大；东部9月沉水植物冠层水深较小，而6月冠层水深较大. 这种冠层水深空间分布特征与现场船舶走航观测结果大体一致.

图11 Landsat-8 OLI影像数据波段比幂函数模型 反演冠层水深的反演值与实测值对比

图12 大气瑞利校正后Landsat-8 OLI影像数据波段 比模型反演冠层水深的反演值与实测值对比

图13 乌梁素海沉水植物冠层水深分布(手动提取了挺水植物)

3 讨论与结论

本文利用实测水生植物遥感光谱数据，分析了沉水植物与挺水植物和漂浮藻类光谱特征差异及沉水植物光谱随冠层水深的变化特征，得出以下结论：

1)沉水植物和漂浮藻类在可见光和近红外的光谱反射率明显低于挺水植物；挺水植物在短波红外1662 nm和2223 nm附近分别有一个反射峰，这是挺水植物区别于沉水植物和漂浮藻类的重要波段. 0深度沉水植物(WDC=0)与漂浮藻类的光谱反射率非常接近，但是在550～690 nm有明显差异，漂浮藻类的反射绿峰要高于沉水植物，并具有明显分层现象. 因此，可以利用绿波段和短波红外波段的光谱特征来区分挺水植物、沉水植物和漂浮藻类.

2)沉水植物群落的光谱反射率随冠层水深的增加而降低，在700～900 nm 波段范围内变化最为明显，且在700～735 nm波段附近，沉水植物群落光谱反射率与冠层水深呈显著负相关. 8月以后的r与ρ所对应的特征波段差异很大，在707 nm之前沉水植物群落光谱反射率与冠层水深的线性关系较好，在707 nm之后两者存在很好的非线性关系.

3)利用LOOCV方法，以单波段/波段比为自变量建立了沉水植物群落冠层水深反演模型，其中波段比反演模型要优于单波段反演模型，波段比反演模型的决定系数R2>0.70，均方根误差<13.70 cm，平均相对误差<28%. 反演精度较好，适用于10～60 cm沉水植物冠层水深的反演. 利用冠层水深与湖水深度，可以绘制沉水植物高度图，与覆盖度相结合可以估算沉水植物生物量，该模型可以有效地推导出冠层水深，这可以为遥感研究沉水植物提供有效依的科学依据.

4)利用波段积分运算将实测光谱反射率扩展应用到Landsat-8 OLI波段上，建立了OLI单波段/波段比冠层水深反演模型，其中波段比幂函数反演模型效果最好，其决定系数R2为0.49，均方根误差为18.17 cm，平均相对误差为40.05%. 利用大气瑞利校正后Landsat-8 OLI影像数据波段比模型得到了2018年9月2日和2019年6月17日乌梁素海沉水植物冠层水深分布图. OLI遥感反演的冠层水深具有一定的不确定性，原因主要为：1)实测点位与遥感数据空间尺度不一致. 现场实测为点位数据，而OLI为面数据(空间分辨率30 m×30 m)，在900 m2内乌梁素海沉水植物水深不均一，这是造成反演结果存在误差的重要原因. 2)由于乌梁素海水环境复杂，精确的大气校正实现困难，本文所用的OLI数据只做了瑞利校正，未做气溶胶校正，这也是导致反演误差的重要因素之一. 因此，在以后的研究中可以结合无人机高光谱对乌梁素海沉水植物冠层水深的遥感反演模型进行验证.