珠海一号高光谱卫星的于桥水库水质参数反演初步研究

殷子瑶，李俊生*，范海生，高 敏，谢 娅

1.中国科学院空天信息创新研究院，数字地球重点实验室，北京 100094 2.中国科学院大学，北京 100049 3.珠海欧比特宇航科技股份有限公司，广东 珠海 519080

引 言

卫星遥感技术在内陆水体水质监测中具有重要作用，可以大范围、较快速、低成本的实现主要水质参数的反演[1-2]。不过，内陆水体水质监测也对卫星遥感技术提出了较高要求：小型水体需要较高的空间分辨率；快速变化的水质要求较高的时间分辨率；复杂多变的水体光学特性要求较高的光谱分辨率。因此，同时满足高空间、高时间、高光谱分辨率的卫星数据在内陆水体监测中具有巨大的优势，但是单颗卫星难以达到这一指标，只能通过星座的方式实现。珠海欧比特宇航科技股份有限公司于2018年4月26日成功发射了4颗珠海一号高光谱卫星，2019年9月19日又成功发射了4颗珠海一号高光谱卫星。珠海一号高光谱卫星具有一致的参数设置，在400～1 000 nm的波谱范围内具有32个波段，光谱分辨率达到2.5 nm；成像幅宽150 km，空间分辨率10 m；8颗卫星组网形成的高光谱卫星星座可以实现2.5 d重返周期[3]。因此，珠海一号高光谱卫星星座同时具有了高光谱、高空间、高时间分辨率，在内陆水体水质监测中具有重要潜力。不过，由于珠海一号高光谱卫星投入运行不久，利用其数据进行的研究还较少，并且目前研究方向主要集中在地物分类，关于内陆水体的研究很少，并且缺少星地同步实验检验其应用效果。

以2018年11月22日获取的覆盖于桥水库珠海一号高光谱影像为遥感数据源，基于星地同步实验获取的水质参数数据，开展珠海一号的水质参数反演建模和评价工作，探索珠海一号高光谱卫星在水质参数反演方面的潜力。

1 研究区与研究数据

1.1 研究区概况

于桥水库(39°99′71″—40°07′44″N，117°43′09″—117°68′26″E)是天津市的主要供水源之一，坐落在天津市蓟州区城东，流域总面积超过2 000 km2，是国家级重点大型水库[4-5]。近年来随着经济的发展，于桥水库水环境状况恶化，富营养化问题较为严重[6]，而且地表径流携带了较多悬浮泥沙和有色可溶性有机物，导致于桥水库水体光学特性比较复杂，水质参数遥感监测的难度较大。通过分析于桥水库中总悬浮物浓度(total suspended matter,TSM)、水体透明度(secchi disk depth,Zsd)的分布特征，可以预测及了解于桥水库水质的发展趋势，为相关部门进行水资源规划提供支撑。

1.2 遥感数据

使用的遥感数据为珠海一号高光谱数据，该数据的辐射量化级别为10 bit，辐射动态范围为1 024级，数据级别为L1A级，即DN值，具体波段参数信息如表1所示。获取了与实测数据同步的一景珠海一号高光谱影像，成像日期为2018年11月22日。

表1 珠海一号高光谱卫星波谱基本信息Table 1 Basic spectral information of Zhuhai-1 hyperspectral satellite

1.3 卫星同步水面实验

2018年11月22日在于桥水库开展水面实验，布设了13个采样点，如图1所示。现场利用赛克盘法测量了水体透明度，并采集水样，送到实验室内利用烘干称重法测量悬浮物浓度。

图1 2018年11月22日珠海一号高光谱彩色合成图像及水面实验采样点位Fig.1 Hyperspectral color synthetic image of Zhuhai-1 and sampling points of water surface experiment on November 22,2018

2 实验部分

2.1 珠海一号高光谱图像定量化预处理

基于卫星遥感图像的水质参数定量反演，要对图像进行定量化预处理，首先计算了原始图像各波段的信噪比；然后进行辐射定标，将图像数码值转为大气层顶辐亮度；再进行大气校正，将大气层顶辐亮度转为地面的反射率。

同时采用FLAASH(fast line-of-sight atmospheric analysis of spectral hypercubes)模型的大气校正方法，对珠海一号高光谱影像进行大气校正，获取了影像的地表反射率，其中用到的能见度可以根据当天MODIS气溶胶反演产品换算得到。水质参数反演一般使用离水反射率或遥感反射率，它与地表反射率存在差别，需要在地表反射率的基础上进一步校正天空光反射率的影响。因此在地表反射率的基础上采用了一种减去近红外(NIR)波段最小值来校正天空光反射的影响，得到离水反射率[7]；离水反射率除以π，便可最终获得水色遥感中常用的遥感反射率[Rrs(λ)]。

2.2 水质参数反演建模方法

纯水以及水中的悬浮颗粒和可溶解物质等对太阳辐射的吸收和散射及水体的状态决定了水体的光谱特征。太阳辐射在水中的传输和分布主要受悬浮颗粒物、浮游植物、溶解性有机物等物质的影响，并且传感器接收的离水辐射的差异主要取决于水体中不同组分对光的吸收和散射。因此，用遥感技术定量监测水质参数的本质是水体中不同组分含量、浓度的差异在一定波谱范围内引起的反射率差异[8]的度量。

常用的水体悬浮物浓度和透明度反演算法主要可以分为两大类，即解析/半解析法和经验/半经验法。本工作主要是初步检验珠海一号高光谱在悬浮物和透明度反演中的效果，因此使用了比较简单的经验/半经验方法进行反演建模分析。在统计学中，常见的回归方法主要包括单自变量的线性、非线性回归，多自变量的多元回归、主成分回归，以及在各种自变量数据转换之后的相关分析等。本研究尝试对各波段及波段组合进行分析，选择最优方法建立基于珠海一号数据的于桥水库悬浮物浓度、透明度反演模型。

根据13个采样点的经纬度位置，提取珠海一号高光谱离水反射率图像中对应的数值，组成数据对，用于水质参数反演建模和精度评价。为了减少采样点与图像象元的匹配误差，需要保证采样点附近的水体图像是比较均匀的。因此，要求每个采样点对应的图像像素为中心的3×3象元的变化率小于40%[9-10]。

利用最终确定的实测和图像数据对，对悬浮物和透明度反演的常用波段和波段组合进行统计分析，构建适用于珠海一号高光谱图像的于桥水库悬浮物和透明度反演模型。

2.3 水质参数反演模型精度评价方法

由于同步采样点比较少，只有13个，不宜再划分为训练样本和检验样本，因此采用每次留一点法(leave-one-out cross validation,LOOCV)进行精度评价[11-12]。该方法的主要思想是每次保留1个采样点的值，用其余12个点的采样值进行建模，然后把模型应用于保留的这个点；这样迭代13次，就得到13个点对应的反演结果，从而进一步计算反演精度。与其他验证方法相比，LOOCV方法不会受到建模和验证样本的分离的影响，消除了潜在异常值对算法精度的影响，更适合于小样本数据的验证[10-13]。在此基础上，采用统计参数对水质参数的反演精度进行检验，包括平均相对误差(mean relative error,MRE)和均方根误差(root mean square Error,RMSE)二者公式如式(1)和式(2)

(1)

(2)

式中，Y是水质参数的反演值，X是水质参数的实测值，n为样本数。

3 结果与讨论

3.1 水质参数反演模型构建结果

利用大气校正后的遥感反射率与同步采样点的悬浮物浓度和透明度值进行单波段/波段组合的Pearson相关性分析和显著性检验，并利用LOOCV法建模，发现Rrs(684)/Rrs(540)与悬浮物浓度的相关系数最高，Rrs(656)/Rrs(556)与透明度的相关系数最高。据此，建立了于桥水库悬浮物浓度和透明度的反演模型：

(3)

(4)

式(3)中，TSM为于桥水库的悬浮物浓度，Rrs(684)/Rrs(540)为珠海一号684 nm波段与540 nm波段遥感反射率的比值。式(4)中，Zsd为于桥水库的水体透明度，Rrs(656)/Rrs(556)为珠海一号656 nm波段与556 nm波段遥感反射率的比值，建模的结果如图2所示。

3.2 水质参数反演模型精度评价结果

根据留一验证法，对构建的悬浮物浓度和透明度模型进行了精度验证(图3)。验证发现悬浮物反演模型的平均相对误差为8.6%，均方根误差为1.0 mg·L-1，透明度反演模型的平均相对误差为11.7%，均方根误差为18.2 cm。两种水质参数实测值与反演值的散点都较均匀的分布在1∶1线的两侧，说明构建的模型可以用来反演于桥水库的水体悬浮物浓度和透明度。

图3 于桥水库悬浮物浓度和透明度模型的验证Fig.3 Model validation of suspended matters concentration and water clarity in Yuqiao Reservoir

3.3 于桥水库水质参数反演结果

利用于桥水库悬浮物浓度和水体透明度的反演模型，应用于珠海一号高光谱遥感反射率图像，得到了于桥水库悬浮物浓度和水体透明度的空间分布图(图4)。从图中可以看到，于桥水库的悬浮物浓度空间分布呈现出明显的北高南低的特点，变化范围6～20 mg·L-1，而水体透明度呈现出北低南高的特点，变化范围为40～170 cm。

图4 2018年11月22日珠海一号高光谱影像反演的于桥水库水体悬浮物浓度和透明度分布图Fig.4 Distribution of suspended matters concentration and water clarity in Yuqiao Reservoir retrieved from Zhuhai-1 hyperspectral satellite;data on Nov. 11,2018

于桥水库悬浮物浓度北高南低的分布特点与于桥水库的水深分布有关。于桥水库南部是古河道，水体相对比较深，北部水体相对较浅[14]。水体越深，底泥越不容易再悬浮[15]，水体表层的悬浮物浓度一般会越低，透明度也会越高。

4 结 论

高光谱卫星具有较多波段，可以从中选择最适用于研究区水质参数反演建模的波段。使用珠海一号高光谱卫星数据，以于桥水库为研究区，定量反演了于桥水库的悬浮物浓度、水体透明度，对反演数据进行了验证，并绘制了于桥水库水体悬浮物浓度和透明度的空间分布图。通过研究得到如下结论：

(1)珠海一号高光谱遥感影像Rrs(684)/Rrs(540)波段组合是于桥水库悬浮物浓度反演建模的最佳波段组合，模型的决定系数R2为0.5，平均相对误差为8.6%，均方根误差为1.0 mg·L-1。

(2)以珠海一号高光谱遥感影像Rrs(656)/Rrs(556)波段组合是于桥水库透明度反演建模的最佳波段组合，模型的决定系数R2为0.7，平均相对误差为11.7%，均方根误差为18.2 cm。

(3)于桥水库的水体悬浮物浓度空间分布呈现北高南低的特点，水体透明度的空间分布则与之相反，呈现出北低南高的分布特征，这种分布特征的形成应该是于桥水库北部水较浅，南部水较深导致的。

(4)基于星地同步实验的于桥水库水质参数反演建模和检验结果初步说明了珠海一号高光谱卫星在定量反演内陆水体水质参数方面具有潜力，未来需要开展更多的星地同步实验进一步完善预处理和反演模型。