基于无人机高光谱遥感的河湖水环境探测-莱森光学

一、引言

河湖水环境监测是人类一直以来高度重视的环境问题。相对于传统监测手段，遥感技术具有快速、大面积同步观测、周期性等特点，对于获取长期、大范围河湖水环境的时空变化具有显著优势。无人机遥感平台搭载高光谱传感器可以获得高空间、高时间、高光谱分辨率的遥感数据，利用该数据可实现狭长河道的长时间精准观测，对于水环境的持续性遥感监测以及紧急重点排查具有重要意义。本研究基于实测和无人机高光谱遥感反射率数据计算水体颜色参量并反演水质参数，利用Hueangle对水体进行分类，通过水体颜色参量和水质参数反演结果，分析上海市崇明岛河湖小水体的水体颜色变化，进而识别河湖水环境中的疑似污染水体。

二、研究区概况及数据源

2.1研究区概况

本研究选取上海市崇明岛、青浦区大莲湖作为无人机高光谱水环境监测的主要研究区域。崇明岛有2条市级河道（南横引河、北横引河）贯通南北形成环岛运河。作为中国第一大生态岛，岛内河道水体的水环境状况直接影响崇明岛的可持续发展和岛内人民的生活质量。大莲湖位于青浦区西部，北部为淀山湖，南部联通黄浦江，是上海市黄浦江上游重要的水源保护区。

a）研究区域分布（b）大莲湖采样点及无人机飞行区域

C）崇明岛采样点及无人机飞行区域

图1研究区域及采样点

2.2实测数据

海市青浦区大莲湖有8个采样点，崇明区河道有86个采样点，其中无人机高光谱数据覆盖采样点（包含同步测量和一小时间隔内测量）35个。使用地物光谱仪根据水面以上测量方法，测量并计算得到采样点水体的遥感反射率数据。使用便携式浊度计现场测量浊度。根据《地表水和污水监测技术规范》采集和保存水样，带回实验室采用分光光度法测量叶绿素a浓度（Chl-a）、有色溶解有机物吸收系数（CDOM）、总氮（TN）和总磷（TP），对于内陆水体选择440nm处的吸收系数ag（440）作为CDOM的代表。采用称重法测量水体总悬浮物浓度（TSM），实测采样点各类水质参数浓度分布如表1所示。

表1采样点水质参数浓度分布

本研究的无人机平台为多旋翼无人机高光谱系统，于2019年5月、9月、11月以及2020年5月在崇明岛河道进行了无人机飞行，共获取了11条崇明区重点关注河道的无人机高光谱数据。

三、高光谱数据处理

3.1遥感反射率测定

本研究基于反射率基法发展了一种无人机高光谱遥感反射率测定方法。该方法在无人机飞行河湖区域铺设面积1m2的标准反射率为20%~30%的漫反射参考板，无人机经过该区域时Rref(λ)和下行辐亮度Lref(λ)，获得太阳总辐照度。并根据水面以上光谱测量方法同步测量参考板上方的天空光辐亮度Lsky(λ)

3.2数据几何矫正及降噪处理

本文采用快速傅里叶变换（FFT）去除高光谱影像的条带噪声，通过计算高光谱影像每列（条带噪声纵向分布）的条带度量值（Si）评估条带噪声去除的效果。采用对噪声敏感的最小噪声分离法（MNF)实现光谱维噪声的去除，利用局部均值和局部标准差法计算高光谱影像的信噪比和噪声标准差，检验高光谱影像的降噪效果。

图2无人机高光谱影像降噪流程及噪声评估

图2（a）显示，利用FFT去除条带后，整体Si有了明显降低，低于条带充分去除标准Si=0.005，证明该方法可以有效去除高光谱数据的条带噪声。图2（b）显示，MNF方法在光谱维上有效剔除了大量噪声，得到了相对平滑的光谱曲线。图2（c），（d）显示，MNF去除噪声后信噪比、噪声标准差有大幅度改善，信噪比均值提升128%，噪声标准差也显著降低，平均值3.42×10-4（图2（d））；FFT去除条带噪声后，信噪比、噪声标准差未显著改善。

四、研究方法

4.1水体颜色参量计算

根据国际照明委员会（CIE）制定的CIE-XYZ颜色标准系统以及定量的描述颜色方法，基于本研究采集的无人机高光谱数据，计算河湖水体的Hueangle。除此之外，还采用给定光谱的主导波长来表示。AVW因为包含了紫外和近红外波段，对于颜色更蓝或更红的极端水域，相对Hueangle具有更高的变化范围。故为了充分利用所获得的高光谱信息来准确量化内陆水体水色，本研究从水色主导波长的角度利用无人机高光谱全波段数据（400～800nm）计算AVW，计算公式为:

4.2水体色相分类

实地调查发现，采样区域水环境恶化主要表现为大量浮萍漂浮的绿色水体和污染源附近的黄棕色污染水体。参考设定Hueangle阈值对内陆河湖水体的分类方法。从实测数据中选出19条不同类型水体（5条浮萍漂浮的绿色水体、7条一般水体和7条污染源附近黄棕色水体）的高光谱反射率曲线计算Hueangle。浮萍漂浮的水体呈绿色，具有植物光谱特征。污染源附近黄棕色水体主要组分为非色素颗粒物或有色溶解有机物，有色溶解有机物中含有的腐殖酸和富里酸随着浓度的升高会使水体呈现黄褐色。根据所选不同水体的Hueangle值将水体分为绿色异常水体（Hueangle≤218°）、一般水体（218°≤Hueangle≤225°）、黄棕色异常水体共3类（Hueangle≥225°）（图3）。

图3不同类型水体光谱曲线及对应Hueangle

4.3水质参数

高光谱遥感反演通过对同步实测的遥感反射率与水质参数（Chl-a，TSM，CDOM，浊度，TN，TP）进行统计分析，将不同水质参数选择相关性最高的波段或波段组合作为自变量，通过指数、幂函数、线性、多项式、偏最小二乘回归等方法构建反演模型（表2），通过Ｒ2和ＲMSE对模型的反演效果进行评估。

表2水质参数反演模型

参考综合营养指数（TLI）计算方法，基于Chl-a、TN、TP浓度估算河湖水体的TLI指数。通过水体的TLI指数，将水体营养状态分为贫营养（TLI＜30）、中营养（30≤TLI≤50）、轻度富营养（50＜TLI≤60）、中度富营养（60＜TLI≤70）、重度富营养（70＜TLI）5种状态。

五、结果讨论

5.1无人机-实测遥感反射率验证

遥感反射率整体测定误差表明：无人机高光谱计算的遥感反射率，各波段ε在7.7%～27.9%之间（平均值13.34%），ＲMSE为0.0034sr-1～0.0062sr-1之间（平均值0.0046sr-1），Ｒ在0.63～0.93之间（平均值0.83），整体误差较低。由于该定标方法没有准确去除光照强度变化对遥感反射率计算的影响，因此一小时间隔内定标点的测定精度略低于同步定标精度（图4）。考虑到大多数无人机由于载重局限，仅搭载一个向下探测的传感器，在光照条件较为稳定均匀的情况下，利用该测定方法获得的遥感反射率精度较高，可以推广应用。

图4无人机-实测Rrs对比验证

5.2水体颜色参量与水质参数反演效果评估

如图，无人机高光谱反演的水体颜色参量基本与原位观测高光谱反演值一致（Hueangle：Ｒ2=0.97，ＲMSE=0.86°；AVW：Ｒ2=0.93，ＲMSE=2.01nm）。相对于多光谱数据反演的水体颜色参量需要通过实测高光谱校正，无人机高光谱数据反演的水体颜色参量更加准确且避免了校正带来的误差。根据相关分析构建的各类水质参数反演模型（图5（b）），从与实测数据的对比结果来看具有较高的反演精度。

图5无人机高光谱反演水体颜色参量、水质参数精度验证

5.3河段水体色相分类

利用25条绿色异常水体、33条一般水体和17条黄棕色异常水体的原位观测高光谱数据对Hueangle水体分类方法进行验证，验证精度分别为68%，91%和82%，证明该方法具有较高的分类精度。根据Hueangle阈值对2019-2020年崇明岛重点河湖水体进行分类。图6（a）显示，2019年11月份的南横引河河段整体都为黄棕色异常水体，结合实际采样调查分析，南横引河作为崇明岛重要的航运河道，船只航行直接造成水体浊度增高，且该河段临近崇明岛外围港口，长期停靠在河道两侧的船只产生的生活污水也会对水体颜色造成影响。同一时间段的七效港河（图6（b））大部分为一般水体，城镇居民生活区一侧的河道有明显的黄棕色异常。2020年5月份崇明岛的河湖水体颜色相对2019年11月份有明显变化，七效港河的黄棕色异常水体面积有明显的减小（图6（d））。南横引河整体的黄棕色异常现象消失转为一般水体，但临近夏季丰水期，开始出现绿色异常水体（图6（c））。通过分析2019年11月和2020年5月崇明岛重点河段Hueangle分类结果，发现崇明岛河湖在枯水期有较多的黄棕色异常水色现象。

图6上海市崇明岛河湖水体Hueangle分类

产生该现象的原因可能与崇明“闸控型”水系有关，在非引排期且雨水量不太充沛的时候，河流之间近乎封闭，导致水动力严重不足，再加上企业、居民生活废水的排放，使得河湖水体颜色出现明显异常。另外南横引河作为主要航运河道，在11月和5月的水体颜色变化明显，怀疑该河道受人为因素影响较为严重，需加强监管和治理。

5.4多参数水环境分析

图7（a）-（i）显示了崇明岛北港东岸转河，水体颜色参量和水质参数基于无人机高光谱数据定量反演的结果，图7（j）显示北港东岸转河左右两侧为黄棕色异常水体。各类水质参数浓度都较高，其中（d）-（e）显示TN，TP浓度达到Ⅳ类水体的标准（1≤TN≤1.5mg/L，0.2≤TP≤0.3mg/L）。TLI指数（图7（i））表明该区域的黄棕色异常水体为轻度富营养化水体。河道中间的绿色异常水体TLI指数也显示为轻度富营养化，部分水质参数（Chl-a，TSM，CDOM，TP）显示高浓度特点（图7（a）-（c），7（e）），但该区域的浊度、TN（图7（h），7（e）)相对黄棕色异常水体较低，分析原因在于该异常区域水流相对缓慢，泥沙等非色素颗粒物沉淀，水体透光性较好，浊度较低，浮游植物光合作用增强，再加上附近多处沟渠和隐蔽排污口不定期排放高营养盐（TP主导）污染物导致出现藻华现象。而河道两侧的异常水域现场调查时正在排放污水，排污过程导致水体浊度升高，水中非色素颗粒物增多，水体呈黄棕色。

图7上海市崇明岛重点观察河道疑似污染水体识别

结果显示Hueangle和AVW均能够有效地对异常颜色水体进行甄别，此外，Hueangle分类可以有效地划分不同的水体颜色类型，用于判别异常水体的空间位置和水色异常类型，从而丰富了水环境遥感探测信息。结合水质参数定量反演，也有助于对疑似污染水体进一步分析，对于突发性河湖水体污染的快速监测提供了一定辅助支持，同时为现场观测提供了有效信息。

六、结果讨论

本研究发展了一种无人机高光谱遥感反射率标定方法，通过该方法计算的高光谱遥感反射率各波段平均ε为13.34%，ＲMSE平均为0.0046sr-1，Ｒ平均为0.83。可以在光照稳定的情况下计算较为准确的水体遥感反射率。后续仍需研究如何去除光照变化计算遥感反射率。通过FFT和MNF方法对高光谱数据的条带噪声和光谱维噪声进行了去除，高光谱整体信噪比提升了128%，有效提高了数据质量。利用高光谱数据反演水体Hueangle，AVW和水质参数，有效利用了高光谱数据丰富的光谱信息。根据Hueangle能够区分不同颜色水体的特点，通过设定Hueangle阈值的方法对水体进行分类，应用该分类方法对上海市崇明岛河湖水体颜色的时空变化进行分析。同时利用水体颜色参量以及水质参数反演结果对上海市崇明岛的重点观察河段进行疑似污染水体的识别。以上结果表明高光谱数据在水体颜色和水质参数反演应用方面的潜力，结合无人机高时效性、高空间覆盖度的特点，对于辅助河湖水环境监测提供了强有力的支持。