海洋学报

环境变量	莫兰指数	Z值	P值
水深	0.511	8.853	0.001
蓝光波段	0.631	10.967	0.001
绿光波段	0.501	8.954	0.001
红光波段	0.301	5.531	0.001

环境变量

莫兰指数

Z值

P值

水深

0.511

8.853

0.001

蓝光波段

0.631

10.967

0.001

绿光波段

0.501

8.954

0.001

红光波段

0.301

5.531

0.001

环境变量	莫兰指数	Z值	P值
水深	0.511	8.853	0.001
蓝光波段	0.631	10.967	0.001
绿光波段	0.501	8.954	0.001
红光波段	0.301	5.531	0.001

环境变量

莫兰指数

Z值

P值

水深

0.511

8.853

0.001

蓝光波段

0.631

10.967

0.001

绿光波段

0.501

8.954

0.001

红光波段

0.301

5.531

0.001

方法	均方根误差/m	平均绝对误差/m	决定系数
对数比值模型	2.067	1.608	0.797
普通克里金模型	1.894	1.487	0.845
随机森林模型	1.635	1.058	0.888
空间自相关随机森林模型	1.338	0.998	0.923

方法

均方根误差/m

平均绝对误差/m

决定系数

对数比值模型

2.067

1.608

0.797

普通克里金模型

1.894

1.487

0.845

随机森林模型

1.635

1.058

0.888

空间自相关随机森林模型

1.338

0.998

0.923

方法	均方根误差/m	平均绝对误差/m	决定系数
对数比值模型	2.067	1.608	0.797
普通克里金模型	1.894	1.487	0.845
随机森林模型	1.635	1.058	0.888
空间自相关随机森林模型	1.338	0.998	0.923

方法

均方根误差/m

平均绝对误差/m

决定系数

对数比值模型

2.067

1.608

0.797

普通克里金模型

1.894

1.487

0.845

随机森林模型

1.635

1.058

0.888

空间自相关随机森林模型

1.338

0.998

0.923

顾及空间自相关特征的机器学习水深反演方法研究

PDF下载

王鑫 ¹^,⁴ , 贝祎轩 ¹ , 陈卓 ³ , 张凯 ¹^,²^,^*

海洋学报 | 论文 2022,44(11): 159-169

收起

海洋学报 | 论文 2022, 44(11): 159-169

顾及空间自相关特征的机器学习水深反演方法研究

全屏

王鑫¹^,⁴, 贝祎轩¹, 陈卓³, 张凯¹^,²^,^*

作者信息

1.山东科技大学测绘与空间信息学院，山东青岛 266590

2.自然资源部第二海洋研究所，浙江杭州 310012

3.中兵勘察设计研究院有限公司，北京 100053

4.广州三海海洋工程勘察设计有限公司，广东广州 510220

王鑫（1996-），男，江苏省扬州市人，主要从事海洋水深反演研究。E-mail：wx849445406@qq.com

通讯作者:

张凯（1983-），男，副教授，主要从事海洋测绘相关研究。E-mail：zk0773@163.com

Retrieving shallow bathymetry by integrating spatial autocorrelation features with machine learning

Xin Wang¹^,⁴, Yixuan Bei¹, Zhuo Chen³, Kai Zhang¹^,²^,^*

Affiliations

1. College of Geodesy and Geomatics, Shandong University of Science and Technology, Qingdao 266590, China

2. Second Institute of Oceanography, Ministry of Natural Resources, Hangzhou 310012, China

3. China Ordnance Industry Survey and Geotechnical Institute Co., Ltd., Beijing 100053, China

4. Guangzhou Sanhai Marine Engineering Surveying and Designing Co., Ltd., Guangzhou 510220, China

出版时间: 2022-11-01 doi: 10.12284/hyxb2022033

文章导航

摘要

收起

基于多光谱影像的水深反演方法是获取近岸水深信息的高效手段，然而反演精度低一直是制约其广泛应用的瓶颈。本文聚焦于实测水深与多光谱数据自身的空间自相关特性，提出在机器学习框架下将学习样本的空间自相关特征与统计互相关特征相结合，以提高水深反演精度。西沙北岛海域的实验结果表明：在实测数据量较小的情况下，相比传统机器学习，顾及自相关特征的新方法可获得18%的精度提升；而当实测数据量充足时，精度提升可达到27%。结果表明，将数据源的空间自相关特征融入机器学习算法中，可显著提升多光谱水深反演结果的精确性，进而为浅海海洋研究提供有效数据支撑。

关键词

水深反演 / 随机森林 / 机器学习 / 空间自相关性

Abstract

收起

Retrieving shallow water depth based on multispectral satellite imagery is highly cost-effective. However, the extensive application of satellite-derived bathymetry has been restricted by its low prediction accuracy. To improve about the accuracy of the retrieved bathymetry, spatial autocorrelation features within the in situ depth measurements and the multi-spectral image are focused in this research. To this end, we develop a machine learning method combining with spatial autocorrelation features and statistical intercorrelation features of learned samples. The experimental results of Xisha Beidao show that compared with the traditional machine learning, the accuracy of the new method is improved by 18% when the number of in situ depths is small. On the contrary, when the number of in situ depths is large, an improvement of 27% in root mean square error is achieved. This demonstrates that incorporating the spatial autocorrelation features of data sources into the machine learning can significantly improve the prediction accuracy, and then provide effective data support for shallow ocean research.

Key words

bathymetry retrieval / random forest / machine learning / spatial autocorrelation

引用本文

王鑫, 贝祎轩, 陈卓, 张凯. 顾及空间自相关特征的机器学习水深反演方法研究. 海洋学报, 2022 , 44 (11) : 159 -169 . DOI: 10.12284/hyxb2022033

Xin Wang, Yixuan Bei, Zhuo Chen, Kai Zhang. Retrieving shallow bathymetry by integrating spatial autocorrelation features with machine learning[J]. Haiyang Xuebao, 2022 , 44 (11) : 159 -169 . DOI: 10.12284/hyxb2022033

正文

收起

1 引言

收起

随着我国海洋强国战略的逐步实施，近岸水深数据的需求显著增加。因此，如何高效、准确地获取近岸水深信息成为当前研究热点。目前，水深测量主要依靠船载声呐测量^[1-2]与机载激光雷达测量^[3-4]。其中，船载多波束测深系统横向覆盖宽度达深度的3～4倍，可高效获取海底的准确水深信息，但在浅水区效率较低且存在搁浅风险。机载激光雷达测量系统可高效测量浅水区域水深信息，但其广泛应用受到测量成本高和飞行区域合法性的限制。相比之下，利用遥感卫星的多光谱影像反演近岸水深信息^[5-6]，因其高效、高时空分辨率、低成本的特点，表现出显著的应用价值。

自20世纪80年代Lyzenga^[7]提出双层流理论水深反演模型以来，遥感卫星技术的进步有力推动了多光谱遥感影像反演水深理论的发展^[8]，所形成的反演方法大致可分为3类：理论解释模型、半理论半经验模型、统计模型。其中，理论解释模型基于水光场辐射传输方程^{[7, 9-11]}，通过建立光学传感器接收到的辐亮度与水深和底质反射的解析表达式来反演水深。这类方法理论严密，但所需水体光学参数较多且获取困难，模型构造复杂。而半理论半经验模型通过结合理论模型和经验参数来实现水深反演，其中以Stumpf对数比值^[12]为代表的方法，将辐射亮度表示为深水区与海底反射辐亮度之和，建立其与水深之间的解析表达式来预测水深，得益于模型构造简单与物理机制清晰的优势被广泛应用^[13-16]。然而对数比值法特征维度较少，导致反演精度受到模型表达能力的限制^[17]。近年来，随着计算机技术的迅速发展，基于机器学习算法的统计模型逐渐成为水深反演最热门的前沿研究领域^[18-21]。该方法不需要考虑水深遥感的物理机制，而是从水深与图像辐亮度值之间的统计关系出发，自发学习数据之间的本质联系来建立模型。凭借其在解决多变量、非线性复杂问题等方面的优势，被引入统计水深反演模型中，取得了优于对数比值法的反演效果。然而受卫星影像信噪比低和实测数据量少的制约，多光谱水深反演精度提升有限，成为限制该方法广泛应用的瓶颈。因此，通过深入挖掘数据潜在信息提高水深反演精度，是多光谱水深反演方法研究的一个发展方向。

上述方法通过挖掘影像光谱值与水深之间的函数关系与统计互相关特征进行水深反演。除此之外，水深信息本身的空间统计自相关特征亦可用于估计未采样点位上的水深信息，如基于稀疏水深数据进行空间插值（如反距离加权法^[22]、克里金法^[23-24]）估计未采样位置的水深信息。因此，在上述传统水深反演方法基础上，通过挖掘水深信息的空间自相关特征，有望有效提高水深反演精度。对此，Su等^[25]提出回归克里金法（Regression Kriging, RK），将空间相关性特征纳入了Stumpf对数比值模型中，有效提升了水深反演精度。然而，RK和Stumpf模型皆属于线性模型，在描述复杂非线性映射关系时表达能力有限。同时，RK基于平稳假设的前提在实际观测环境中也常常难以得到有效满足^[26]。因此，该方法的效果在环境各向异性显著的区域明显下降。

针对上述问题，本文以提高多光谱水深反演精度为目标，利用机器学习在解决多变量与非线性问题方面的优势，通过融合数据源的空间自相关特征与统计互相关性来研究高精度水深反演方法，以期在现有研究基础上，进一步提高预测水深的准确性。

2 方法

收起

2.1 对数比值模型（Stumpf模型）

对数比值法是代表性的半理论半经验水深反演模型^[12]，该方法以蓝、绿波段反射率的对数比值作为反演因子，其模型表达式为

(1)

$ d = \beta_{0} + \beta_{1} \cdot \frac{{\ln \left[ {{{n}} \cdot R_{w}\left( {\lambda _{{\rm{b}}}} \right)} \right]}}{{\ln \left[ {{{n}} \cdot R_{w}\left( {\lambda _{{\rm{g}}}} \right)} \right]}} = \beta_{0} + \beta _{1 }\cdot x \text{，} $

式中，d为反演出的水深；β₀和β₁是通过回归方程得到的系数；n为缩放比例因子；R_w(λ_b)和R_w(λ_g)分别是蓝、绿波段的反射率。

该方法具有构造简单，对环境因素干扰（如大气、水体和海底反射率的变化）不敏感等优势。但作为线性拟合模型，该方法在实际应用中仍然受到复杂水质和海底底质变化的显著影响，难以准确描述复杂海区反射率与水深呈现出的非线性关系^[27]。同时，未能综合利用多个波段的光谱测量信息亦是该方法的一个缺点。

2.2 随机森林

机器学习方法可以对不同变量之间的复杂映射机制进行有效学习和表达。其中，随机森林（Random Forest，RF）模型凭借其优秀的非线性回归和泛化能力、多变量特征融合能力以及运算速度快等优势^[28-30]，得到广泛应用。用于多光谱反演水深时，RF模型通过并行处理策略将一组回归决策树组合生成学习器，并利用未抽取样本来泛化模型误差，综合分析每个模型的预测结果。通过上述学习过程，RF模型可以获得优于对数比值法的水深反演精度，其算法流程如图1所示。

光学水深遥感反演的物理基础是光对水体的穿透能力，受水体的漫衰减影响，光在水体中的穿透深度与水体特性、不同谱段的漫衰减系数密切相关性，可通过机器学习发掘其中显著而复杂的函数关系。RF模型的数学表达式为

(2)

$ d={\displaystyle \sum _{k=1}^{N}[h({\boldsymbol{X}}, {{\boldsymbol{\textit{Θ}}}}_{k})]},\;\;k=1,2,\cdot \cdot \cdot,N $

式中，h（）代表随机森林算法的函数关系式；N为决策树的总数；输入向量X为波段光谱值；Θ_k为样本通过模型训练得到的第k棵树的参数向量；d为综合所有决策树预测结果的平均值所得到的反演水深。

2.3 空间自相关随机森林

2.3.1 空间自相关随机森林模型

传统RF模型通过分析多光谱数据与实测水深数据之间的统计互相关特征预测水深，但忽视了数据自身的空间自相关特征，故未能充分利用有效信息进行学习和表达^[31]。为此，本文在训练模型过程中，以实测水深及对应的波段光谱值作为基础因/自变量前提下，将待测点周围具有高度空间自相关性的已知点所对应的属性值作为自相关变量纳入RF的学习框架中，以此实现统计互相关特征与空间自相关特征的融合。其算法流程图如图2所示。

2.3.2 空间自相关与最佳空间间隔

空间自相关特征在遥感图像分析中的应用十分广泛，常用于变化检测与分类优化^[32-33]。空间自相关性指的是研究区域内要素的属性值之间潜在的相互依赖性，可用来衡量地理现象空间聚集程度。研究表明，预测点周围具有较高自相关性的已知点所包含的属性信息可用来提高预测的精度^[34]，理论上，数据之间的统计相关性越强，预测精度提升幅度越大。对此，全局莫兰指数（Moran’s I）^[35]常用于度量研究数据的聚集性于自相关程度。莫兰指数越大，表明研究对象的聚集性越强，自相关程度越高。其计算公式为

(3)

$ I = \frac{{N\displaystyle\sum\limits_{i = 1}^N {\displaystyle\sum\limits_{j = 1}^N {{\boldsymbol{w}}_{ij}\left( {x_{i} - \overline x} \right)({x_{j} - \overline x} )} } }}{{\left( {\displaystyle\sum\limits_{i = 1}^N {\displaystyle\sum\limits_{j = 1}^N {\boldsymbol{w}_{ij}} } } \right)\displaystyle\sum\limits_{i = 1}^N {{{\left( {{x}_{i} - \overline x} \right)}^2}} }} \text{，} $

式中，N为要素的总数目；x_i和x_j分别为第i和第j个要素的光谱值与水深值；$\overline x $为上述属性值的均值；w_ij为要素i和j的空间权重值，这里为基于距离的邻接矩阵。

由于Moran’s I基于空间随机性分布的零假设，因此需要利用Z得分和P值来判断假设是否成立，并检验自相关指数的显著性^[36]。其中，P值是由已知分布的曲线得出的面积近似值，而Z值按以下形式计算：

(4)

$ Z = \frac{{I - E\left[ I \right]}}{{\sqrt {V\left[ I \right]} }} \text{，} $

式中，E[I]和V[I]分别为Moran’s I指数的期望值和方差。

受空间距离的影响，变量的自相关程度仅在一定范围内较为显著，若超过该距离，则认为发生了空间变异。为确定局部最佳自相关的空间间隔，采用半方差来体现空间变异程度^[37]，半方差越大，自相关越弱，其计算公式为

(5)

$ r(h) = \frac{1}{{2n}}\sum\limits_{i = 1}^n {\mathop {[Z(x_{i}) - Z(x_{i} + h)]}\nolimits^{\text{2}} } \text{，} $

式中，h为样本点之间的距离；n为距离h范围内成对样本点的数量；Z为样本点所对应的不同波段光谱值与水深值。

2.3.3 空间自相关变量

为了获取训练样本的自相关特征信息，需要建立局域子窗进行提取，具体步骤如下：

（1）计算训练样本的水深与不同波段值的半方差函数和莫兰指数，用于确定搜索窗口初始尺寸。

（2）根据计算窗口中待测点到m个已知点的二维欧式距离进行排序，得到衡量自相关程度的距离向量D_n，并将遍历所有待测点得到的n个距离向量，合成为空间自相关距离矩阵（EDF）。

（3）由于训练数据分布的不均匀性，D_n中元素数量k应取所有距离向量的最少元素数目，且需保证不等于0，以满足向量组成矩阵的条件；此外，为在不影响预测精度条件下，提升模型训练效率，还需对搜索窗口大小进行调整并重复步骤（2），将得到的EDF作为表征数据空间自相关特征的补充变量，EDF的具体计算公式计算为

(6)

$ d_{ij} = \sqrt {{{(x_{i} - x_{j})}^2} + {{(y_{i} - y_{j})}^2}} \;\; (i=1, 2, \cdots, n, j=1, 2, \cdots,m), $

式中，d_ij为i和j两点二维欧式距离；（x_i，y_i）和（x_j，y_j）分别为已知点和待测点坐标。

(7)

$ {\boldsymbol{D}}_{n} = \left[ {dn_{1},dn_{2},\cdots,dn_{k}} \right] （k\leqslant m\text{，}k\ne 0） \text{，} $

(8)

$ {\boldsymbol{EDF}} = \left[ {{\boldsymbol{D}}_{1};{\boldsymbol{D}}_{2};\cdots;{\boldsymbol{D}}_{n}} \right] \text{，} $

式中，D_n为距离向量；EDF为n个距离向量的集合。

以每个距离向量中k个元素为前提，按其自相关程度的强弱顺序，依次获取每个已知点所对应的属性值并将其与EDF共同组成空间自相关变量X₂，与基础变量X₁构成反演因子。

(9)

$ X_{1} = {\text{\{ }}r{_{\text{1}}},r_2,r{_{\text{3}}\} } \text{，} $

(10)

$ X_2 = {\text{\{ }}r_{1k},r_{2k},r_{3k},{\textit{z}}_k,{\boldsymbol{EDF}}\} \text{，} $

式中，r₁、r₂、r₃为预测点对应的蓝、绿、红波段光谱值比值，即B_blue/B_green，B_green/B_red，B_blue/B_red；r_1k、r_2k、r_3k、z_k分别为k个已知点的光谱值比值与水深值。

2.3.4 空间自相关随机森林模型训练/验证

基于上述步骤，空间自相关随机森林（Spatial Autocorrelation Random Forest，SARF）模型训练与验证过程存在一定差异，下面列出关键步骤说明。

（1）对水深点个数为a的训练数据而言，所用的水深点既充当待测点，也作为已知点参与模型训练，即按给定的搜索半径计算每个当前水深点到周围邻近点的距离，用以构成大小为a×k的自相关矩阵EDF。将当前点水深值作为监督学习的因变量，与空间自相关变量、基础自变量输入RF模型中训练，得到训练参数。

（2）对于验证数据（数量为b）来说，每个水深点都为待测点，故需将训练数据中的所有水深点作为已知点来构建自相关矩阵（大小为b×k），并将待测点对应的空间自相关变量与基础自变量输入训练好的模型中，即可求得待测点水深。

3 实验区域与数据源

收起

3.1 实验区域概况

本文研究区域为北岛周边浅海水域。北岛位于南海西沙群岛海域七连屿的中部，面积约为3.6 km²，岛礁呈长条形（图3a）。作为典型的珊瑚岛，其底质由珊瑚、砂和贝屑组成，水下地形特征复杂。由于远离大陆且受人类活动影响小，岛屿周边水质清澈，适合多光谱水深反演。

3.2 数据源

本研究使用的原始影像为高分辨率的WorldView-2 卫星影像（图3b），拍摄于2017年3月11日11时。多光谱分辨率为1.84 m，包含蓝、绿、红以及近红外4个波段。使用的数据来自于2013年4月测量得到的北岛1:2 000陆域及水下地形测量图，对其进行矢量化与筛选后得到1 700个水深点（图3b）。

3.3 数据预处理

数据的获取过程中存在时间、空间以及仪器测量带来的误差，因此在利用遥感影像与水深数据进行定量水深研究前需要对这两类数据源进行预处理。对于Worldview-2影像，首先通过辐射定标和大气校正得到真实地物反射率，并对影像进行几何校正以及云雾和陆地部分的掩膜，使用Hedley法^[38]消除受到太阳光、海面波浪和光线入射角度等因素造成的太阳耀斑。对于水深数据而言，将其坐标转换到与影像同一坐标系UTM WGS-84下，并进行地理配准，以提取水深点对应的像元值。由于遥感影像采用的是卫星过境时瞬时海面深度，而实测水深数据是以理论深度基准面为基准的稳态水深，因此还需要根据潮位预报表提供的潮高进行潮汐改正。

本文中水深分布区如图4所示，主要集中于0～6 m浅水区域。为模拟小样本浅海水深反演这一常见应用场景，同时确保数据具有代表性，以1 m为间隔分层抽样，共选取150个水深数据作为训练样本，而剩余1 550个水深点则作为验证样本。为定量比较4种反演模型效果，本文选取决定系数（R²），均方根误差（RMSE）和平均绝对误差（MAE）进行精度评价。

4 实验结果分析

收起

4.1 变量的空间自相关性

通过计算训练样本的Moran’s I指数、Z值和P值对水深与波段光谱值的空间自相关性加以评估（表1）。以蓝光波段为例，Moran’s I指数为0.631，说明该波段光谱值呈现出较强的聚类型式，而其Z值（10.967）和P值（0.001）则表明可以拒绝空间随机性分布的零假设，并通过了置信度为99%的显著性检验。表1中结果显示，训练数据的水深、蓝、绿波段皆具有显著的自相关性，而红光波段自相关性较弱。

4.2 搜索窗口尺寸确定

研究区域变量的全局Moran’s I及半方差随空间间隔的变化如图5所示，通过计算波段值的半方差，并除以最大值归一化后与全局Moran’s I对比发现，随着空间间隔变大，全局Moran’s I减小，而半方差值呈上升趋势，这表明变量的空间变异程度增加，而空间自相关性逐渐减小。值得注意的是，二者相交所产生的交点反映了空间自相关与空间变异的平衡。对于水深及蓝、绿、红波段而言，对应的最佳空间间隔分别为175 m、156 m、152 m和76 m，这表明上述变量在该局域内具有显著空间自相关性。基于上述观察，本文将初始搜索窗口尺寸定为175 m，并根据距离向量的最小元素数原则，将搜索窗口尺寸扩大为220 m，用以获取空间自相关变量并构建SARF水深反演模型。

4.3 精度评价

为验证实验方法效果，本文首先将反演结果与传统方法进行了对比（表2）。可以看到，受限于模型表达能力，Stumpf模型反演精度最低，RMSE为2.067 m，R²仅为0.797，相关性也最低。普通克里金（Ordinary Kriging, OK）模型受限于空间插值方法对样本点密度等要求，在小样本实测水深数据集的条件下，其预测精度亦不甚理想（RMSE为1.894 m）。RF模型得益于综合多个波段光谱信息，RMSE为1.635 m，决定系数达到0.888，相比前两种模型预测结果更准确。在这4种模型中，SARF模型表现效果最好，RMSE仅为1.338 m，相对于RF模型预测精度提升18%，同时相关性参数也提高到0.923。

4.3.1 误差分析

为了比较4种模型的反演精度和优劣性，根据水深反演图（图6）对不同方法的预测结果与真实水深进行对比分析。图6中黑线是斜率为1的指标线，黑线上方表示预测水深值大于实测水深值，下方则表示预测值偏小；红线是通过最小二乘拟合的趋势线，用以衡量与指标线的拟合程度；为体现预测值聚集区间，采用色棒来标识，聚集程度越高，则数值越大且颜色越浅。

Stumpf模型的水深反演效果如图6a所示，在水深浅于10 m的区域，反演结果与实测水深符合较好。但是受限于线性表达方式，出现反演水深值小于0的不合理情况，且对于深于10 m的数据，水深反演结果产生了明显的偏差。很显然，该偏差现象产生的根源在于线性模型无法准确拟合数据的非线性特征；类似地，作为非线性局部拟合方法，OK模型较为准确地预测了0～8 m区域水深（图6b），而受制于8 m以深实测数据量的减少，无法真实还原深水区情况。虽然该方法顾及了水深的自相关特征信息，但未考虑水深与光谱值之间的相关性，因此只能通过增加实测数据量来弥补模型的缺陷。对比前两种模型，图6c和图6d中机器学习预测精度提升效果明显，较为完整地还原了研究区域内真实水深分布情况。但从细节上比较发现，RF模型受训练过程中低信噪比与实测数据量不足影响，指标线上方出现大量离散点且与趋势线相距较远，并影响了预测结果可靠性；与之相比，SARF模型具有明显的抗离群值效果以及更高的反演精度。

残差分布图分别给出了不同方法在表达水深−波段值耦合机制过程中的差异，从而通过残差的分布形式来对4种模型的细节差异进行比较（图7）。图7a显示，Stumpf模型的残差在10 m以深区域呈现出明显的递增趋势，预测值与真实值存在较大差异；由图7b可知，OK模型整体残差值较大，特别是在8 m以深区间，大量的离散点明确显示出模型自身的缺陷；而RF模型在浅水区域的效果较好（图7c），残差值在0两侧均匀分布，但在水深超过10 m的区域，其预测结果受到离群值干扰；与前者相比，无论是从残差整体分布形式，亦或是深水区数据受自身统计值波动影响所产生的离群点数量，SARF模型都表现出更明显的优势（图7d）。

4.3.2 不同训练样本数量下反演结果的精度对比

训练样本数量的变化可显著影响机器学习方法的预测效果。通过计算训练样本数为500时SARF模型和RF模型的残差标准差（图8a）可以看到，RF模型残差标准差（1.295）远高于SARF模型（1.012），因此后者残差离散程度较小，分布更为集中；其次，通过观察0.05显著性水平下残差分布区间发现，RF模型残差主要分布于[1.25, 1.35]，而SARF模型残差分布区间集中在[0.97, 1.05]，后者的残差分布幅度更小，说明预测值与实测值拟合效果更好。除此之外，通过观察训练样本占全部样本比例下的均方根误差（图8b）发现，SARF模型的误差大小及其下降速率始终小于RF模型，且随着训练数据量的增加，该优势愈加明显；当训练样本占比达到60%时，相比于RF模型，SARF模型的误差降低约27%，该结果进一步验证了将数据源的空间自相关特征纳入机器学习框架中，可显著提升预测结果的准确性。

本文基于SARF模型使用350个训练数据得到北岛海域反演水深，并将其与实测陆上高程点生成反演水深与陆域地形图（图9a）。可以看到，珊瑚岛周围海底地形分布表现出复杂的特征。近岸浅水区域海底主要由沙质沉积物构成，水深变化较为平衡。而在离岸较远的区域，受海流的影响，岛屿南北两侧呈现出截然不同的海底地貌细节特征（图9b）。显然，海底地形信息表现出较强的空间相关性，从而验证了利用空间相关性辅助水深反演的有效性。

5 结论与讨论

收起

随着多光谱卫星反演水深技术的发展，如何在小样本现场实测数据的前提下，有效提高预测精度一直是研究热点。为此，本研究基于机器学习的框架，提出了结合数据源自相关特征的空间自相关随机森林水深反演模型，并在北岛海域开展多光谱水深反演实验。通过对比改进方法与对数比值、普通克里金、随机森林方法的实验结果，得到以下结论：

（1）Stumpf模型通过建立波段反射率与水深之间的数学关系进行预测，具有模型构造简单、对实测数据量要求较低等优势，因此得到了广泛应用。然而在水深反演过程中，受到来自深度、水质以及底质等因素变化影响，反射率与水深之间的线性关系通常不成立，而Stumpf这类线性模型受限于模型的表达能力，渐渐无法满足日益增长的高精度遥感测深需求。对此，以RF模型为代表的机器学习方法凭借其优异的非线性映射能力，可有效提升水深反演精度。

（2）实验表明，以OK模型为代表的空间插值方法仅仅利用实测水深的自相关性特征信息，即取得与Stumpf模型相似的水深预测精度。这一结果客观上印证了本文通过引入空间自相关特征信息以提高水深反演精度这一思路的可行性。对此，本文在机器学习框架下引入空间自相关特征进行建模，提出了空间自相关随机森林模型，通过深入挖掘不同数据源之间、相邻像元内部观测值的有效信息，将输入变量的统计互相关性与空间自相关性两类特征进行融合，得到了更准确的预测模型。该方法比单独使用两类特征的反演模型更具优势。SARF模型在输入反演因子时引入两类相关性变量的信息，提高了预测过程中的信噪比，具有比传统机器学习更出色的抗噪能力。同时，SARF模型在建模过程中充分利用训练样本的观测值，通过建立训练样本与测试样本之间的联系，从而提高有效信息的利用率。除此之外，将数据源的空间自相关特征融入机器学习算法中，有效减弱了空间非平稳性因素、极浅区域高辐射亮度值以及深水区统计值波动的影响，可显著提升不同实测数据量条件下的多光谱水深反演结果的精度。

（3）鉴于SARF模型在遥感测深研究中取得了显著的精度优势，基于相同思路，将自相关性特征引入其他机器学习模型中亦可望取得类似的提升效果。需要注意的是，本研究区域内，水深主要集中于0～6 m浅水区，而考虑到不同研究区数据源的差异，将其应用于更多地区还需要进一步探讨。基于现有方法，未来的研究方向可考虑挖掘额外的环境信息，利用机器学习融合多变量的优势，增加更多环境相关特征（如水体指数、叶绿素浓度）并融合不同时序的多光谱影像数据源进行研究。

基金

收起

山东省自然科学基金(ZR2020MD084)
国家自然科学基金重点基金(41930535)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

刘经南, 赵建虎. 多波束测深系统的现状和发展趋势[J]. 海洋测绘, 2002, 22(5): 3−6.

Liu Jingnan, Zhao Jianhu. Status and development tendency of multi-beam bathymetric system[J]. Hydrographic Surveying and Charting, 2002, 22(5): 3−6.

赵建虎, 欧阳永忠, 王爱学. 海底地形测量技术现状及发展趋势[J]. 测绘学报, 2017, 46(10): 1786−1794.

Zhao Jianhu, Ouyang Yongzhong, Wang Aixue. Status and development tendency for seafloor terrain measurement technology[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1786−1794.

Richard Z, Daniel I, David R, et al. Habitat classification of temperate marine macroalgal communities using bathymetric LiDAR[J]. Remote Sensing, 2014, 6(3): 2154−2175.

杨必胜, 梁福逊, 黄荣刚. 三维激光扫描点云数据处理研究进展、挑战与趋势[J]. 测绘学报, 2017, 46(10): 1509−1516.

Yang Bisheng, Liang Fuxun, Huang Ronggang. Progress, challenges and perspectives of 3D LiDAR point cloud processing[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1509−1516.

蒋兴伟, 林明森, 张有广, 等. 海洋遥感卫星及应用发展历程与趋势展望[J]. 卫星应用, 2018(5): 10−18.

Jiang Xingwei, Lin Mingsen, Zhang Youguang, et al. Development and prospect of ocean remote sensing satellite and its application[J]. Satellite Application, 2018(5): 10−18.

Hodul M, Bird S, Knudby A, et al. Satellite derived photogrammetric bathymetry[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2018, 142(8): 268−277.

Lyzenga D R. Passive remote sensing techniques for mapping water depth and bottom features[J]. Applied Optics, 1978, 17(3): 379−383.

马毅, 张杰, 张靖宇, 等. 浅海水深光学遥感研究进展[J]. 海洋科学进展, 2018, 36(3): 5−25.

Ma Yi, Zhang Jie, Zhang jingyu, et al. Progress in shallow water depth mapping from optical remote sensing[J]. Advances in Marine Science, 2018, 36(3): 5−25.

Lyzenga D R, Malinas N P, Tanis F J. Multispectral bathymetry using a simple physically based algorithm[J]. IEEE Transactions on Geoscience & Remote Sensing, 2006, 44(8): 2251−2259.

陈启东, 邓孺孺, 秦雁, 等. 广东飞来峡库区水深遥感[J]. 中山大学学报(自然科学版), 2012, 51(1): 122−127.

Chen Qidong, Deng Ruru, Qin Yan, et al. Water depth extraction from remote sensing image in Feilaixia Reservoir[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2012, 51(1): 122−127.

Figueiredo I N, Pinto L, Gonalves G. A modified Lyzenga’s model for multispectral bathymetry using Tikhonov regularization[J]. IEEE Geoscience & Remote Sensing Letters, 2016, 13(1): 53−57.

Stumpf R P, Holderied K, Sinclair M. Determination of water depth with high-resolution satellite imagery over variable bottom types[J]. Limnology and Oceanography, 2003, 48(1/2): 547−556.

Liang Jian, Zhang Jie, Ma Yi. A spatial resolution effect analysis of remote sensing bathymetry[J]. Acta Oceanologica Sinica, 2017, 36(7): 102−109.

陈本清, 杨燕明, 罗凯. 基于高分一号卫星多光谱数据的岛礁周边浅海水深遥感反演[J]. 热带海洋学报, 2017, 36(2): 70−78.

Chen Benqing, Yang Yanming, Luo Kai. Retrieval of island shallow water depth from the GaoFen-1 multi-spectral imagery[J]. Journal of Tropical Oceanography, 2017, 36(2): 70−78.

陈琛, 马毅, 张靖宇. GF-1 WFV图像经验模分解的光谱保真性与水深遥感探测[J]. 海洋学报, 2018, 40(4): 51−60.

Chen Chen, Ma Yi, Zhang Jingyu. Spectral fidelity and water depth remote sensing detection of EMD of GF-1 WFV images[J]. Haiyang Xuebao, 2018, 40(4): 51−60.

Xia H, Li X, Zhang H, et al. A Bathymetry mapping approach combining log-ratio and semianalytical models using four-band multispectral imagery without ground data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(4): 2695−2709.

王燕红, 陈义兰, 周兴华, 等. 基于多项式回归模型的岛礁遥感浅海水深反演[J]. 海洋学报, 2018, 40(3): 121−128.

Wang Yanhong, Chen Yilan, Zhou Xinghua, et al. Research on reef bathymetry using remote sensing based on polynomial regression model[J]. Haiyang Xuebao, 2018, 40(3): 121−128.

Wang Yanjiao, Zhang Peiqun, Dong Wenjie, et al. Study on remote sensing of water depths based on BP artificial neural network[J]. Marine Science Bulletin, 2007, 9(1): 26−35.

王锦锦, 马毅, 张靖宇. 基于模糊隶属度的多核SVR遥感水深融合探测[J]. 海洋环境科学, 2018, 37(1): 130−136.

Wang Jinjin, Ma Yi, Zhang Jingyu. Multiple kernel support vector regression based on fuzzy membership for remote sensing water depth fusion detection[J]. Marine Environmental Science, 2018, 37(1): 130−136.

朱金山, 纪轩禹, 宋珍珍. 基于支持向量机和BP神经网络的水深反演研究[J]. 测绘与空间地理信息, 2019, 42(6): 11−14.

Zhu Jinshan, Ji Xuanyu, Song Zhenzhen. Water depth inversion based on support vector machine and BP neural network[J]. Geomatics & Spatial Information Technology, 2019, 42(6): 11−14.

温开祥, 李勇, 王华, 等. 基于遥感和机器学习的内陆水体水深反演技术[J]. 热带地理, 2020, 40(2): 314−322.

Wen Kaixiang, Li Yong, Wang Hua, et al. Estimating inland water depth based on remote sensing and machine learning technique[J]. Tropical Geography, 2020, 40(2): 314−322.

刘光孟, 汪云甲, 张海荣, 等. 空间分析中几种插值方法的比较研究[J]. 地理信息世界, 2011, 9(3): 41−45.

Liu Guangmeng, Wang Yunjia, Zhang Hairong, et al. Comparative study of several lnterpolation methods on spatial analysis[J]. Geomatics World, 2011, 9(3): 41−45.

印兴耀, 刘永社. 储层建模中地质统计学整合地震数据的方法及研究进展[J]. 石油地球物理勘探, 2002, 37(4): 423−430.

Yin Xingyao, Liu Yongshe. Methods and development of integrating seismic data in reservoir model-building[J]. Oil Geophysical Prospecting, 2002, 37(4): 423−430.

何红艳, 郭志华, 肖文发. 降水空间插值技术的研究进展[J]. 生态学杂志, 2005, 10(15): 1187−1191.

He Hongyan, Guo Zhihua, Xiao Wenfa. Review on spatial interpolation techniques of rainfall[J]. Chinese Journal of Ecology, 2005, 10(15): 1187−1191.

Su H, Liu H, Wu Q. Prediction of water depth from multispectral satellite imagery-the regression Kriging alternative[J]. IEEE Geoscience & Remote Sensing Letters, 2015, 12(12): 2511−2515.

Brunsdon C, Fotheringham A S, Charlton M E. Geographically weighted regression: A method for exploring spatial nonstationarity spatial[J]. Geographical Analysis, 2010, 28(4): 281−298.

Zhang K, Wang X, Wu Z, et al. Improving statistical uncertainty estimate of satellite-derived bathymetry by accounting for depth-dependent uncertainty[J]. IEEE Transactions on Geoscience and Remote Sensing, 2021, PP(99): 1−9.

Breiman L. Random forests[J]. Machine Learning, 2001, 45(1): 5−32.

Ham J, Chen Y, Crawford M M, et al. Investigation of the random forest framework for classification of hyperspectral data[J]. IEEE Transactions on Geoscience & Remote Sensing, 2005, 43(3): 492−501.

Stumpf A, Kerle N. Object-oriented mapping of landslides using random forests[J]. Remote Sensing of Environment, 2011, 115(10): 2564−2577.

邱耀炜, 沈蔚, 纪茜. 随机森林模型在遥感水深反演中的应用[J]. 海洋技术学报, 2019, 38(5): 98−103.

Qiu Yaowei, Shen Wei, Ji Qian. Satellite-derived bathymetry using random forest model[J]. Journal of Ocean Technology, 2019, 38(5): 98−103.

朱钟正, 苏伟. 基于局部空间统计分析的SPOT 5影像分类[J]. 遥感学报, 2011, 15(5): 957−972.

Zhu Zhongzheng, Su Wei. The analysis of the classification of SPOT 5 image based on local spatial statistics[J]. National Remote Sensing Bulletin, 2011, 15(5): 957−972.

张涛, 方宏, 韦玉春, 等. 顾及空间自相关性的高分遥感影像中建设用地的变化检测[J]. 自然资源学报, 2020, 35(4): 212−225.

Zhang Tao, Fang Hong, Wei Yuchun, et al. Detection of the construction land change in fine spatial resolution remote sensing imagery coupling spatial autocorrelation[J]. Journal of Natural Resources, 2020, 35(4): 212−225.

Behrens T, Schmidt K, Viscarra R, et al. Spatial modelling with Euclidean distance fields and machine learning[J]. European Journal of Soil Science, 2018, 69(5): 757−770.

Anselin L. Local indicator of spatial association-LISA[J]. Geographical Analysis, 1995, 27(2): 93−115.

张松林, 张昆. 全局空间自相关Moran指数和G系数对比研究[J]. 中山大学学报(自然科学版), 2007, 46(4): 93−97.

Zhang Songlin, Zhang Kun. Comparison between general Moran’s index and Getis-Ord general G of spatial autocorrelation[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2007, 46(4): 93−97.

Carr J R, De Miranda F P. The semivariogram in comparison to the co-occurrence matrix for classification of image texture[J]. IEEE Transactions on Geoscience & Remote Sensing, 1998, 36(6): 1945−1952.

Hedley J D, Harborne A R, Mumby P J. Technical note: Simple and robust removal of sun glint for mapping shallow-water benthos[J]. International Journal of Remote Sensing, 2005, 26(10): 2107−2112.

2022年第44卷第11期

PDF下载

132

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2022033

接收时间：2021-10-11
首发时间：2026-02-01
出版时间：2022-11-01

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2021-10-11
修回日期：2021-12-03

基金

山东省自然科学基金(ZR2020MD084)

国家自然科学基金重点基金(41930535)

作者信息

1.山东科技大学测绘与空间信息学院，山东青岛 266590

2.自然资源部第二海洋研究所，浙江杭州 310012

3.中兵勘察设计研究院有限公司，北京 100053

4.广州三海海洋工程勘察设计有限公司，广东广州 510220

通讯作者:

张凯（1983-），男，副教授，主要从事海洋测绘相关研究。E-mail：zk0773@163.com

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2022033

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT