科技前沿▏一种顾及辐射畸变的多波束与侧扫声呐通用海底底质分类方法

严俊赵建虎等溪流之海洋人生 2023-05-07

点击上方“溪流之海洋人生”即可订阅哦

声学海底底质分类对海洋环境的科学研究和生态系统管理具有重大意义，已经成为海洋科学一个研究热点。相比于传统的实地采样方法，声学法是目前最有效的大区域海底底质分类的遥感方式。

声学海底底质分类可以通过对多波束声呐(Multibeam echo sounder，MBES)、侧扫声呐(side scan sonar，SSS)、单波束/分裂波束声呐(Single-beam/Split-beam echo sonar，SBES)等常见声呐设备的后向散射回波强度数据进行分析获得。虽然SBES有着成本低、便携性强、易安装等优势，但是采样点密度较低、大面积测量效率低，于是很多时候被MBES和SSS替代。然而MBES与SSS的回波强度数据都受到辐射畸变的严重影响，图像质量较低，因而在进行海底底质分类时必须顾及辐射畸变的影响。一些研究中利用辐射畸变导致的回波强度的角度响应特点进行海底底质分类，并与实际海底情况对比，证明了利用角度响应数据分类的可行性。另一些研究则首先削弱辐射畸变的影响，然后对处理后的MBES和SSS回波强度数据及图像使用BP神经网络、支持向量机、卷积神经网络进行分类并取得了合理的分类结果。

另外，由于SSS图像分辨率较高、纹理特征较丰富，还可利用灰度共生矩阵特征对SSS图像进行海底底质分类。然而这些方法都需要较为复杂的前期预处理或者原始图像质量较好，否则辐射畸变对分类的结果会造成巨大的影响。同时，由于MBES和SSS都是通过换能器发射声波并记录后向散射回波强度数据的方式来形成海底声学图像，且海底底质的变化会直接导致后向散射回波强度的变化，因此MBES和SSS回波强度图像在反映海底底质变化上具有极强的相似性。虽然单独的MBES与SSS底质分类方法研究较多，但两种声呐数据的通用海底底质分类方法研究较少。因此，本文提出了一种顾及辐射畸变的多波束与侧扫声呐通用海底底质分类方法。

一、辐射畸变的影响

MBES与SSS通过发射声波并记录来自海底的后向散射回波强度来获取海底图像以及实现海底底质分类，然而受到换能器波束模式、传播损失、海底散射模式等因素影响后向散射回波强度数据存在严重的辐射畸变，见图1。受到辐射畸变影响严重的区域，辐射畸变对回波强度变化的影响会十分严重，甚至会超过底质变化对回波强度变化的影响，造成底质变化的判断失误．因此底质分类中必须对辐射畸变进行改正。

图1 多波束与侧扫声呐的辐射畸变影响

受到波束主瓣与旁瓣的影响，MBES与SSS的波束模式会导致回波强度随着角度的变化而变化，而SSS的单宽波束模式相对于MBES多波束模式导致回波强度的角度响应更为明显。声波在水中传播时，由于受到球形扩散产生的传播损失和水体的吸收两部分影响，声波能量不断地衰减，在垂航方向上会导致回波强度的斜距响应，进而可以转换为角度响应。声波受到海底散射模式的影响，波束的后向散射回波强度随着波束入射角变化而变化，即为散射模式导致的回波强度角度响应。综上所述，辐射畸变的主要因素都会对垂航方向上的回波强度造成角度响应，因而在海底底质分类的过程中需要处理声呐回波强度的角度响应。

由于SSS声呐高度是不断变化的且无法准确实时获得，因此在沿航方向上也存在由传播损失导致的辐射畸变，通常可以通过后期计算SSS高度的变化来消除。而MBES由于安装在船体上，可以有效地改正传播损失，因而可以忽略沿航方向辐射畸变的影响。

二、顾及辐射畸变的海底底质分类方法

声呐波束模式、传播损失和海底散射模式综合导致了垂航方向上MBES和SSS回波强度的角度响应。基于在同一角度下受到的辐射畸变影响相同这一特征，首先对不同角度下的数据进行归一化；然后，根据归一化后的数据变化趋势与源数据一致的特点，对归一化后的图像进行非监督分类，并通过自适应的迭代方法来确定分类个数；最后利用形态学方法去除分类图像中的噪声。

⒈利用z分数进行角度数据归一化

消除了沿航方向的辐射畸变影响后，垂航方向上的MBES与SSS的回波强度数据在相同波束入射角度下受到的角度响应影响是一样的，因此在这种情况下可认为底

质变化影响是影响强度变化的唯一因素。这就意味着，可以利用同一角度的强度变化直接反映底质的变化。但是，由于辐射畸变导致的回波强度角度响应，MBES与SSS不同角度下的后向散射回波强度明显处于不同强度范围，见图2(a)。通过z分数实现角度数据归一化，保证了不同角度下数据变化的一致性，从而消除了辐射畸变在垂航方向上的影响。反之，如果不进行归一化，回波强度的角度性差异会直接影响分类结果。因而需要首先对不同角度下数据进行归一化。

图2 沿航方向下不同入射角下回波强度序列的归一化

z分数是将多个不同变化区间的参数归算到一个共同变化区间的常用数学统计手段。对于任意同角度下的强度序列B(均值为μ，标准差为σ)，则B中任一强度值b的z分数为：

z＝(b－μ)/σ

不同入射角下的回波强度数据通过z分数归一化后，原来不同强度区间的数据被统一到共同的变化区间中，尤其是同一个底质类型下的回波强度值基本接近，见图2(b)。由于不同底质产生的强度变化趋势仍反映在归一化后的z分数序列中，且保持在同一区间内变化，从而实现了不同角度下底质变化的统一，消除了辐射畸变的影响，进而可以实现全图的海底底质分类。

⒉利用k-means++对归一化后的图像进行分类

在没有先验数据的情况下，可利用非监督分类方法对归一化后的回波强度数据进行分类。相比于其他非监督方法，k-means作为非监督分类的经典方法，具有简单、快速、适用于大数据挖掘等优点，并且顾及到声呐回波强度数据量较大但是数据结构简单，故本文选择k-means算法用于海底底质分类。但k-means同时存在两个基本缺陷，即初始中心变化会导致结果的不唯一性和分类数k的无法确定性。第一个问题可以通过k-means++算法有效地解决；对第二个问题，本文根据海底底质类型的特点给出了自适应的k值迭代确定方法。

⑴k-means++非监督分类

使用k-means++算法对归一化后的z分数图像进行非监督底质分类，具体的步骤如下：

①将z分数图像I_z展成向量V；

②使用k-means++算法对向量V进行k个聚类中心初始化;

③计算V中每个值到每个聚类中心的距离，并将每个值分配到最近中心；

④计算每个聚类的平均参数来建立k个新的聚类中心位置；

⑤重复步骤3～4直到所有的分配不再变化，得到V每个值对应的分布D；

⑥将D还原为二维分类图像I_D。

⑵分类数k 的迭代确定算法

由于k-means分类是将不同的回波强度值分配到不同聚类的中心范围中，每个聚类中心代表该类型底质的中心强度，而每个聚类范围内的回波个数所占比例则代表这种底质在该水域所占比例。在初始k值较大情况下，对占有比例极少(＜5%)的分类，往往是噪声或粗差的影响，应删除对应分类。对于中心强度极为接近(底质简单情况，中心距离＜整体回波强度变化范围的10%；底质复杂情况，中心距离＜整体回波强度变化范围的5%)的两种分类，可以认为是同一种底质，应合并相近分类。综上考虑，在无外部数据下的分类数k 的自适应确定方法步骤如下：

①k默认选择较大的数值(比如k＝7)；

②按照指定的k值对z分数图像进行k-means++非监督分类；

③计算每种聚类类型的中心强度和占有比例；

④删除占有比例小的分类，合并相距中心较近的分类，从而得到新的k值；

⑤重复步骤2～4，直至k值不再变化。

⒊利用形态学方法图像去噪

声呐图像中噪声在底质分类后同样会给底质图像带来影响，可以通过形态学的开闭运算来消除噪声并且连接相同底质。在图像形态学中，图像A开运算指的是利用图像元B对图像A进行先腐蚀后膨胀的过程，主要用于将相近较小片的对象连接在一起；而图像A的闭运算指的是利用图像元B对图像A进行先膨胀后腐蚀的过程，主要用于去除由图像噪声带来的小斑点，通过合理顺序使用这两种运算可以达到需要的目的。由于SSS图像信噪比较高，对SSS底质图像进行先开后闭运算，可以有效地连接图像中的相同底质。而由于MBES图像信噪比较低，对MBES底质图像进行先闭后开的运算，可以有效地去除噪声斑点。

⒋顾及辐射畸变的MBES与SSS通用海底底质分类流程

综上所述，顾及辐射畸变的MBES与SSS通用海底底质分类方法包括: 首先将MBES和SSS数据转换为入射角－回波强度数据，接着利用z分数实现数据归一化，然后通过迭代方法自适应确定分类数k值后利用k-means++实现海底底质分类，最后通过形态学的方法提高底质图像质量，见图3。

图3 海底底质分类流程图

三、实验与分析

为了验证本文方法的有效性，选取2012年在福建近岸相同水域实测的MBES及SSS数据进行实验。MBES数据由EM3002型多波束测深仪(工作频率300kHz，单头最大开角130°，工作时最大采样点数131，测深分辨率为1cm，声强振幅分辨率0.5dB)实测获得。SSS数据则由Edgetech 4100P型拖鱼声呐(工作频率为500kHz单频测量，记录最大斜率为112m，俯角20°，水平波束角为0.5°，垂直波束角为50°)实测获得。整个实验水域由7条MBES测线和4条SSS测线组成，对记录的MBES(*.all)数据、SSS(*.xtf)数据进行解码分别得到原始数据，并且进行常规预处理。其中MBES回波强度为dB值，而SSS回波强度为采集软件归化后的0～2552值。

⒈非监督分类数k的确定

选取实测水域内具有特点的几处区域进行分析，利用文中非监督分类数k迭代确定方法来确定该区域的合理分类数。实验首先验证了其中一个区域海底底质分类数k 的迭代确定过程，以及多次迭代中不同k值对应的海底底质分类结果，见图4。

图4 实测水域分类数k值确定的迭代算法

将图4中每次分类结果对应海底类型的占有比例和中心位置进行统计记录，见表1。

表1 分类数对应的每种类型所占的比例以及中心位置

分类数k	类型1	类型2	类型3	类型4	类型5	类型6	类型7
k＝7	0.1%	0.7%	4.2%	14%	32%	18%	32%
k＝4	4%	30%	33%	33%
k＝3	21%	46%	33%
中心	125	101	61
k＝2	66%	34%
中心	111	62

注：k＝7中类型1、2、3和k＝4中类型1 所占比例过小应剔除，而k＝3中类型1和2中心距离近应该合并

针对该段声呐图像，首先将初始分类数k设为较大值7，k-means++分类后统计每种类型的占有比例得到表1中数据第一行。通过剔除占有比例小于5%的分类将k缩减为4，重新分类后得到表1中第二行。同样剔除小比例分类得到k为3，分类后得到每个类型的比例以及其中心位置(对应灰度声呐图像0～255的值)；由于此时底质情况简单，合并中心距离较近(小于255×10%)的分类得到k＝2。最后在k＝2进行分类后最终可以获得一个稳定合理的分类结果，见图4(d)和表1。

利用实际区域表层沉积物采样得知该区域主要的海底底质为淤泥和砾石，与本文自适应方法最终获得的结果一致，从而证明了自适应分类数确定方法的合理性。

⒉福建水域海底底质分类

利用本文提出的MBES和SSS海底底质分类流程处理测区内MBES和SSS中每条测线，并按照回波强度采样点的实际位置生成地理编码的底质图像。然后，根据每条测线的位置关系将所有测线拼接得到整个区域的海底底质分类图像。MBES测线与SSS测线的区域声呐图像和对应海底底质分类结果，见图5。

图5 福建水域同区域下MBES(a)与SSS(b)海底底质分类结果

通过对比图5(a)和(b)可以看到在实验区域中，MBES与SSS的原始图像在回波强度上差距明显，SSS图像在灰暗变化上要比MBES更加明显；但是两者的底质分类的结果整体基本一致，SSS由于自身位置精度不高与MBES底质分类图像在位置和形状上有细微差异。

对比两种声呐数据的底质分类结果，在整个测区中均匀地选择了10个采样点(P1～P10) 进行分析，见图5。对图5中标定的10个点依次统计其在MBES 和SSS声呐图像中的回波强度值以及对应位置的底质分类结果，见表2，两组强度数据的变化关系统计结果，见图6。

图6 选取点对应的MBES与SSS回波强度对应变化关系

表2 选取图5中点对应位置的MBES与SSS回波强度以及对应分类结果

点名	MBES强度/dB	分类结果	SSS强度×10⁴	分类结果
P1	－23.7	砾石	0.98	砾石
P2	－19.2	砾石	1.12	砾石
P3	－18.4	砾石	1	砾石
P4	－19.4	砾石	1.21	砾石
P5	－33.1	淤泥	0.32	淤泥
P6	－30.4	淤泥	0.38	淤泥
P7	－16.7	砾石	1.02	砾石
P8	－31.6	淤泥	0.22	淤泥
P9	－31.8	淤泥	0.29	淤泥
P10	－17.6	砾石	1.22	砾石

从表2和图6可以看出同名点下MBES和SSS的回波强度的变化关系基本吻合，符合理论规律；细微的差异则由两种声呐系统的测量方式差异导致。通过表2中所有同名点的分类结果完全一致可以得知MBES和SSS的分类结果具有很好的一致性，从而验证了本文方法对于MBES和SSS数据海底底质分类的可行性与通用性。

四、结 论

本文提出了一种顾及辐射畸变的MBES与SSS通用海底底质分类方法，有效地消除了辐射畸变对于海底底质分类中的严重影响，通过数据归一化、非监督分类以及图像去噪等步骤建立了完整的海底底质分类流程。将方法应用于福建水域同区域下实测的MBES与SSS声呐数据，通过合理地判定海底底质类型，实现了从测线到全区域的底质分类图像构建，并通过对比MBES与SSS的海底底质分类结果证明了本文方法在MBES和SSS声呐数据海底底质分类中的可行性和一致性。

【作者简介】本文作者/严俊赵建虎孟俊霞张红梅，分别来自安徽大学资源与环境工程学院、武汉大学测绘学院、武汉大学海洋研究院、安徽建筑大学土木工程学院和武汉大学动力与机械学院。第一作者严俊，1989年出生，男，工学博士，安徽大学资源与环境工程学院，讲师；赵建虎，1970年出生，男，教授，武汉大学海洋研究院，博士生导师。本文为基金项目，国家自然科学基金(41576107、41376109、51804001)、安徽大学博士科研启动经费(J01003270)、安徽省自然科学基金(1808085QE147)、安徽省高校自然科学研究重点项目(KJ2017A038)。文章来自《哈尔滨工业大学学报》（2019年第5期），参考文献略，用于学习与交流，版权归作者及出版社共同拥有，转载也请备注由“溪流之海洋人生”微信公众平台整理。