基于高光谱遥感技术的积水沉陷区水深反演模型的研究

孟祥来;朱继文;张贺

【摘 要】高光谱遥感是当前内陆水监测的一个重要发展趋势,是有效地捕捉这些水体光谱信息的必要手段.在遥感空间信息技术支持下对内陆积水沉陷区的水体特征进行分析,研究叶绿素a及总悬浮物浓度随着水体深度变化的特征,分析不同叶绿素a浓度、悬浮物浓度以及悬浮物粒径对水体光谱特性的影响特点;同时还对底质对水深信息遥感提取的影响特点进行分析.在此基础上,构建水深遥感反演的定量模型. 【期刊名称】《黑龙江工程学院学报（自然科学版）》 【年(卷),期】2011(025)001 【总页数】3页(P17-19)

【关键词】高光谱遥感;光谱反射率;水深遥感;反演模型 【作 者】孟祥来;朱继文;张贺

【作者单位】黑龙江工程学院测绘工程学院,黑龙江哈尔滨150050;黑龙江工程学院测绘工程学院,黑龙江哈尔滨150050;黑龙江工程学院测绘工程学院,黑龙江哈尔滨150050 【正文语种】中 文 【中图分类】P237

1 叶绿素a浓度与水深变化的关系

利用采水器按照水面以下0.2 m一个间隔，同步分层采集水样。采用以下公式进

行拟合：

式中：C chla为叶绿素a浓度（mg／m3），Z 为水体深度，利用 Matlab工具拟合σ、μ的结果为a1＝12.032，σ＝0.236，μ＝2.477，标准差0.654 9。 2 不同叶绿素a浓度对水体光谱的影响

对应水体中不同含量的叶绿素，水体的光谱曲线在430～700 nm光谱波段会相应出现较明显的差异。图1显示不同叶绿素浓度的水体光谱曲线。从图中可见：水体在波长440 nm处有个吸收峰；蓝光（400～480 nm）反射辐射与叶绿素浓度负相关；在波长520 nm处的辐射值不随叶绿素含量而变化；在波长550 nm处出现反射辐射峰，并随着叶绿素含量增加，反射辐射上升；由于浮游植物分子吸收光后，再反射引起的Raman效应（进行水分子破裂和氧分子生成的光合作用），激发出的能量荧光化的结果，在波长685 nm附近有明显的荧光峰。 图1 不同叶绿素浓度的水体光谱曲线

图2 为沉陷区中心部分采样点，叶绿素a的单位吸收系数随深度有较大的变化，尤其在短波峰值附近上层水体的叶绿素a单位吸收系数要偏大很多，在长波峰值附近的变化较小。上层水体具有较强的光照，但营养盐含量较低，该环境下的浮游植物色素细胞粒径较小，而且所含的辅助色素也较丰富。小的色素颗粒由于对强光的光适应性使得内部色素浓度较低，打包效应较弱，叶绿素a单位吸收系数增大；而且在水体的垂直结构上，因为光适应性，表层藻类细胞内色素的浓度相对深层为低，打包效应的影响也较小，相应地叶绿素a单位吸收系数增大。在叶绿素a浓度最大深度以内的水层，微型藻占有绝对优势；而在该深度以下，微型藻类和大型藻类的比例有所增加，色素打包效应增强。 图2 叶绿素单位吸收系数与水体深度的关系 3 悬浮物浓度与水体深度的关系

数据采集方法同不同水体深度时叶绿素a的采集方法相同。悬浮物浓度与水体深度的关系利用式（2）进行拟合。

式中：Cs为悬浮物浓度（mg／m3），Z表示水体深度。利用Matlab工具拟合σ、μ的结果为a 2＝75.062，a3＝0.653，标准差0.784 6。 4 不同悬浮物浓度对水体光谱的影响

为了解悬浮物浓度和光谱反射率之间的关系，本文计算了两组实验各波段归一化反射率和悬浮物浓度间相关系数：在400～600 nm（低悬浮物浓度为400～584 nm，高悬浮物浓度在400～609 nm）波段范围内，光谱反射率和悬浮物浓度呈负相关，且相关系数的绝对值随波长的增加而增大；在600～900 nm（低悬浮物浓度为585～900 nm，高悬浮物浓度在610～900 nm）波段范围内，光谱反射率和悬浮物浓度呈正相关，相关系数随着波段的增加而增大，在900 nm左右达到最高。其中光谱反射率和悬浮物浓度相关性较高的波段范围为：低悬浮物浓度水体在550～570 nm和750～900 nm处，相关系数的平均值为－0.97和0.96，高悬浮物浓度水体在560～590 nm和740～900 nm处，相关系数的平均值为－0.93和0.94。

由于在进行光谱测量时，天气条件的变化、周围环境的影响及测量角度的变化都会影响反射率数值的大小，为了便于对不同次测量结果进行比较，每条反射率曲线都利用其在可见光范围（450～750 nm）的波段反射率平均值进行归一化处理，归一化反射率计算公式为

式中：LW（λ）和L R（λ）分别为对应于某一波长λ的水面和标准板的反射强度。 因为，FieldSpec Hand Held光谱仪采集的是离散型数据，因此，光谱数据的一阶微分可以用式（4）近似计算。

式中：i＋1，i与i－1为相邻波长，LW（λ＋1）为反射率。

水体光谱反射率和悬浮物浓度呈明显的对数关系，因而计算了悬浮物浓度对数和归一化反射率之间的相关性：低悬浮物浓度水体光谱反射率和悬浮物浓度对数高相关波段位于400～500 nm和620～735 nm处，相关系数的平均值分别为－0.99和0.98；高悬浮物浓度水体相关性高值区位于550～570 nm和720～900 nm处，相关系数的平均值分别为－0.98和0.98。可见悬浮物浓度对数与反射率之间的相关性明显高于原悬浮物浓度的相关性。

一阶微分光谱反映的是反射光谱的斜率，两组实验光谱值都在430 nm、560 nm、730 nm和827 nm附近斜率差别较大。各波段反射率的一阶微分与悬浮物浓度对数的相关系数显示：低悬浮物浓度水体在430～460 nm、490～570 nm 和780～790 nm处呈较高的正相关，相关系数均达到了0.88以上，在720～735 nm、817～834 nm处呈较高的负相关，相关性在－0.85至－0.90。高悬浮物浓度水体光谱的一阶微分与悬浮物浓度对数主要以正相关为主，在590～610 nm和700～735 nm处相关性较高，相关系数均在0.88以上。从整体来看，水体光谱一阶微分与悬浮物浓度对数的相关性不如水体光谱本身与悬浮物浓度对数相关性高，且相关性波动较大。

5 沉陷区底质对水深遥感信息提取的影响

在理论分析上，这里先不考虑大气和表面耀光的作用，由于水体部分含有悬浮质、叶绿素和CDOM等物质，因此，对采用光谱对底质进行提取需要对水体部分进行校正，从水体表面测得的向上辐照度包含水体和底质两部分，即

式中：Ewd（0）是水体产生的向上辐照度，EBu（0）是底质产生的向上辐照度。在水深Z的区域位置

式中：bbd（Z）是水体中向下辐照度的后向散射系数。由于水体的每一深度都会产生向上辐照度，则恰在水面下的向上辐照度

A是底部向下辐照度的后向散射系数，由底部反射到恰在水面下的向上辐照度，计算为

对于没有底部反射的水体来说＝R（∞），则

通过计算，就可以获得水深的信息。在校正过程中，R（∞）一般认为是等于没有底部反射时的水体反射率，通过对临近水域较深的水体进行提取。在现场校正时，通过式（15）来计算。

6 水深遥感反演模型构建 水深与水质参数的关系

将式（17）代入底质修正方程式（13），则可以完成水深信息的底质修正。 7 结束语

本文在计算水深遥感反演的定量模型中，沉陷区积水水体的高光谱遥感数据需要实际测量，目前商用高光谱遥感卫星的空间分辨率有限，无法满足沉陷区这种空间尺度相对较小的内陆水体监测的需要，因此，若想达到实用效果，还需要高空间分辨率的多光谱遥感数据支持。 参考文献

［1］李素菊.利用分析方法建立湖泊水质参数反演算法研究—以安徽巢湖为例

［D］.北京：北京大学，2003：1－20.

［2］武彦斌，彭苏萍，黄明，等.淮南积水沉陷区水深遥感定量分析［J］.中国矿业大学学报，2007，36（4）：537－541.

［3］彭苏萍，王磊，孟召平，等.遥感技术在煤矿区积水塌陷动态监测中的应用—以淮南矿区为例［J］.煤炭学报，2002，27（4）：374－378.

［4］陈清莲，唐军武，王项南，等.东海试验区水体光谱特性现场测量与数据分析［J］.海洋技术，1999，18（3）：21－37.

［5］修鹏，刘玉光，程永存，等.水下高光谱辐照度和辐亮度剖面测量仪［J］.气象水文海洋仪器，2006（2）：7－14.