基于PCNN与混合加权的多波束声呐图像镶嵌方法.pdf

第4 2卷 第4期2 0 2 3年8月 山东科技大学学报(自然科学版)J o u r n a l o f S h a n d o n g U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y(N a t u r a l S c i e n c e)V o l.4 2 N o.4A u g.2 0 2 3D O I:1 0.1 6 4 5 2/j.c n k i.s d k j z k.2 0 2 3.0 4.0 0 1 文章编号:1 6 7 2-3 7 6 7(2 0 2 3)0 4-0 0 0 1-1 0基于P C N N与混合加权的多波束声呐图像镶嵌方法曲 萌1,高兴国1,王 峰2,闫循鹏2,阳凡林2,3(1.山东电力工程咨询院有限公司,山东 济南 2 5 0 0 1 3;2.山东科技大学 测绘与空间信息学院,山东 青岛 2 6 6 5 9 0;3.自然资源部 海洋测绘重点实验室,山东 青岛 2 6 6 5 9 0)摘 要:多波束声呐图像是进行海底底质分类与目标识别的主要数据源之一。原始反向散射强度经过改正处理后仍存在残余误差,需要进行图像镶嵌以削弱残差的影响,构建大区域声呐图像。本研究针对多波束声呐图像镶嵌中的拼接线消除与信息综合问题,提出一种在N S S T域内结合P C NN与混合加权的图像镶嵌方法。首先,对单条带声呐图像进行N S S T分解,再对重叠区低频项按照混合加权进行镶嵌,高频项按照P C NN进行镶嵌。实验结果表明,在拼接线消除方面,镶嵌后条带拼接线两侧声强均值差由-6.1 8降为-1.2 9 d B;在信息综合方面,镶嵌结果的信息熵相对于混合加权法提升2.5%,空间频率提升3 1.6%。本研究方法能够在削弱拼接线的同时,有效综合相邻条带纹理信息,实现多波束声呐图像的镶嵌。关键词:多波束测深系统;图像镶嵌;混合加权;非下采样剪切波变换;脉冲耦合神经网络中图分类号:P 7 1 7;P 2 2 9.1 文献标志码:A收稿日期:2 0 2 2-0 9-0 3基金项目:国家自然科学基金重点项目(4 1 9 3 0 5 3 5);山东科技大学科研创新团队支持计划项目(2 0 1 9 T D J H 1 0 3)作者简介:曲 萌(1 9 9 4),男,山东潍坊人,工程师,硕士,主要从事海洋工程勘测相关研究工作.王 峰(1 9 9 7),男,安徽广德人,博士研究生,主要从事多波束声散射图像处理与应用研究,本文通信作者.E-m a i l:s d u s t w f 1 6 3.c o mM u l t i b e a m s o n a r i m a g e m o s a i c m e t h o d b a s e d o n P C N N a n d b l e n d i n g w e i g h t e d f u s i o nQU M e n g1,GAO X i n g g u o1,WANG F e n g2,YAN X u n p e n g2,YANG F a n l i n2,3(1.S h a n d o n g E l e c t r i c P o w e r E n g i n e e r i n g C o n s u l t i n g I n s t i t u t e C o.L t d,J i n a n 2 5 0 0 1 3,C h i n a;2.C o l l e g e o f G e o d e s y a n d G e o m a t i c s,S h a n d o n g U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y,Q i n g d a o 2 6 6 5 9 0,C h i n a;3.K e y L a b o r a t o r y o f O c e a n G e o m a t i c s,M i n i s t r y o f N a t u r a l R e s o u r c e s,Q i n g d a o 2 6 6 5 9 0,C h i n a)A b s t r a c t:M u l t i b e a m s o n a r i m a g e s a r e o n e o f t h e m a i n d a t a s o u r c e s f o r s e a f l o o r s e d i m e n t c l a s s i f i c a t i o n a n d t a r g e t d e t e c t i o n.H o w e v e r,t h e r e a r e s t i l l r e s i d u a l e r r o r s i n b a c k s c a t t e r i n t e n s i t y a f t e r c o r r e c t i n g,a n d i m a g e m o s a i c s a r e n e e d e d t o r e d u c e t h e i m p a c t o f r e s i d u a l e r r o r s t o b u i l d a l a r g e-a r e a s o n a r i m a g e.T h i s p a p e r p r o p o s e d a n i m a g e m o s a i c m e t h o d c o m b i n i n g p u l s e c o u p l e d n e u r a l n e t w o r k(P C NN)a n d b l e n d i n g w e i g h t e d m e t h o d i n t h e n o n-s u b s a m p l e d s h e a r l e t t r a n s f o r m(N S S T)d o m a i n t o e l i m i n a t e s e a m l i n e s a n d f u s e i n f o r m a t i o n.F i r s t l y,t h e s i n g l e s o n a r i m a g e w a s d e c o m p o s e d b y N S S T a n d t h e l o w-f r e q u e n c y t e r m s w e r e f u s e d a c c o r d i n g t o t h e b l e n d i n g w e i g h t e d m e t h o d w h i l e t h e h i g h-f r e q u e n c y t e r m s w e r e f u s e d a c c o r d i n g t o P C NN.T h e r e s u l t s s h o w t h a t i n t e r m s o f s e a m l i n e e l i m i n a t i o n,t h e m e a n d i f f e r e n c e o f b a c k s c a t t e r o n b o t h s i d e s o f t h e s e a m l i n e s a f t e r t h e m o s a i c r e d u c e d f r o m-6.1 8 t o-1.2 9 d B.I n t e r m s o f i n f o r m a t i o n f u s i o n,t h e i n f o r m a t i o n e n t r o p y i n t h e m o s a i c r e s u l t w a s i n c r e a s e d b y 2.5%c o m p a r e d w i t h t h e b l e n d i n g w e i g h t e d m e t h o d,a n d t h e s p a t i a l f r e q u e n c y w a s i n c r e a s e d b y 3 1.6%.T h e p r o p o s e d m e t h o d c a n e f f e c t i v e l y f u s e t h e t e x t u r e i n f o r m a t i o n o f a d j a c e n t s t r i p s w h i l e w e a k e n i n g t h e s e a m l i n e s,t h u s r e a l i z i n g t h e m o s a i c o f m u l t i b e a m s o n a r i m a g e s.K e y w o r d s:m u l t i b e a m e c h o s o u n d e r s y s t e m;i m a g e m o s a i c;b l e n d i n g w e i g h t e d f u s i o n;N S S T;P C NN山东科技大学学报(自然科学版)2 0 2 3年第4期多波束测深系统(m u l t i b e a m e c h o s o u n d e r s y s t e m,MB E S)能够获取海底的反向散射强度(b a c k s c a t t e r,B S),利用反向散射强度可以生成声呐图像并用于海底底质分类与海底目标识别1-2。受测量中海洋环境和仪器自身误差的影响,原始数据存在明显的系统性误差3。经过传播损失改正、角度响应改正等处理后,声强数据中的系统性误差明显减小4。但由于改正不完善,在条带间仍存在由残差引起的声强差异,表现在声呐图像上为重叠区域明显的拼接线,对海底底质分类与目标识别造成干扰,影响图像质量和应用。此外,由于测量时声波入射角不同,重叠区域在不同条带声呐图像中呈现不同的纹理特征。因此,在应用声呐图像之前,需对声呐图像进行镶嵌以消除重叠区域差异,并尽量保留更多细节信息,对提高声呐图像质量具有重要意义。声呐图像处理需进行图像配准和图像镶嵌。图像配准的主要目的是消除图像的几何畸变。多波束测深系统通常采用“平均声强法”对回波序列进行采样,得到的声强点与测深点一一对应,经过波束归位后,采样点具有较高的几何精度,可认为不存在几何畸变。图像镶嵌的主要目的是消除图像灰度畸变,实现多图像的无缝拼接。在光学图像拼接中,常采用的方法有平均法、混合加权法(b l e n d i n g w e i g h t e d f u s i o n,BWF)5、泊松法6等。其中,混合加权法由于计算简便、效率高等优势,在图像拼接中应用广泛。目前,对于多波束声呐图像相邻条带的拼接研究相对较少。孙文川等7提出一种基于加权最小二乘估计的声强数据归一化方法,以消除重叠区域与非重叠区域的强度差异。赵建虎等8提出一种通过滑动窗口逐p i n g确定拼接线的拼接方法,但该方法需要进行相邻条带p i n g的匹配。朱正任等9提出一种基于多条带最小二乘拟合的多波束声呐图像残差处理方法,有效削弱了残差对声呐图像的影响。上述方法主要针对相邻条带拼接线的消除问题,在一定程度上削弱了条带间声强的残差,但未顾及重叠区域纹理信息存在的差异,易造成信息损失。在光学图像处理中,常基于图像融合实现信息的综合,并以多尺度变换作为图像融合框架,如离散小波变换、非下采样轮廓波变换、非下采样剪切波变换(n o n-s u b s a m p l e d s h e a r l e t t r a n s f o r m,N S S T)等。在变换域内,对不同尺度的图像分量按照不同的规则进行融合,可以得到比较理想的效果。巩稼民等1 0在N S S T域内结合脉冲皮质模型(s p i k i n g c o r t i c a l m o d e l,S CM)与引导滤波实现了红外图像与可见光图像的融合。曹义亲等1 1提出一种基于N S S T的脉冲耦合神经网络与稀疏表示相结合的多元遥感卫星图像融合方法。在声呐图像中,相关研究主要是多波束声呐图像与侧扫声呐图像在小区域的信息综合与互补1 2-1 3,并未考虑声呐图像拼接对图像信息综合的影响。综上所述,目前多波束声呐图像镶嵌方法未能实现大区域无缝拼接和信息综合的有效统一。为实现多波束声呐图像的无缝拼接,同时顾及重叠区域纹理差异,实现图像信息综合,本研究基于脉冲耦合神经网络(p u l s e c o u p l e d n e u r a l n e t w o r k,P C NN)和混合加权法,在N S S T变换域内进行多波束声呐图像的镶嵌,消除图像中存在的灰度畸变,实现声呐图像的拼接和信息综合。图1 多波束声呐图像镶嵌流程图F i g.1 F l o w c h a r t o f t h e m u l t i b e a m s o n a r i m a g e m o s a i c1 N S S T域内多波束声呐图像镶嵌方法 针对多波束声呐图像镶嵌问题,为了实现相邻条带的无缝拼接,本研究提出一种在N S S T域内结合P C-NN与混合加权的多波束声呐图像镶嵌方法(如图1)。首先对多波束反向散射强度数据进行传播损失、角度响应等系统误差的改正,再按照地理坐标对两条带图像进行配准。多波束声呐图像镶嵌过程在N S S T的基础上进行,主要步骤包括源图像N S S T变换、低频图像融合、高频图像融合与N S S T逆变换。其中,低频图像融合2曲 萌等:基于P C NN与混合加权的多波束声呐图像镶嵌方法采用混合加权方法,以削弱图像间的拼接线。高频图像融合采用P C NN,充分保留图像的细节纹理信息。最后通过N S S T逆变换完成多波束声呐图像的镶嵌。图2 非下采样剪切波变换(N S S T)过程示意图F i g.2 D i a g r a m o f n o n-s u b s a m p l e d s h e a r l e t t r a n s f o r m(N S S T)p r o c e s s1.1 N S S T图像变换图2展示了图像N S S T变换的过程,主要包括多尺度分解和方向局部化。图中的方框表示频域图像。首先使用上采样得到的非下采样金字塔(n o n-s u b s a m p l e d p y r a m i d,N S P)滤波器组对频域图像进行多尺度分解,得到一个低频和一个高频子带系数1 4,对低频子带可继续进行分解得到更小尺度的低频和高频系数。然后对所有高频子带进行方向局部化,方向局部化采用剪切波滤波器(s h e a r l e t f i l t e r,S F)。在频域内将剪切波滤波器窗口作用于分解得到的高频子带,获得不同方向子带1 5。如图2所示,经过N S S T变换后得到一个低频子带系数和若干不同方向的高频子带系数。N S S T可以避免图像的下采样操作,同时避免了伪吉布斯现象,具有良好的平移不变性。1.2 低频系数的混合加权融合低频系数包含纹理、结构等信息,能够反映声呐图像之间的强度差异。为了削弱图像中存在的灰度畸变,同时充分保留低频系数的有效信息,实现条带间的无缝拼接,采用混合加权法对低频系数进行拼接。混合加权法是根据重叠区域的宽度,以每个像素点到重叠区边界的距离进行定权,使得图像在重叠区边界平滑过渡,削弱声强差异,表达式为5:I(x,y)=Il(x,y),(x,y)Il;W1(x,y)Il(x,y)+W2(x,y)Ir(x,y),(x,y)(IlIr)Ir(x,y),(x,y)Ir。;(1)式中:Il和Ir分别表示左图像与右图像;x和y分别表示像素点在图像坐标系中的坐标;W1和W2是两条带重叠区域像素点(x,y)的权重函数,且W1+W2=1。权重函数为:W1(x,y)=Xr-xXr-Xl,W2(x,y)=x-XlXr-Xl。(2)式中,Xl和Xr分别表示重叠区域的左边界和右边界坐标。图3 脉冲耦合神经元结构F i g.3 P u l s e c o u p l e d n e u r o n a l s t r u c t u r e1.3 高频系数的P C N N融合高频系数反映图像的突变特性,含有丰富的细节信息。本研究将P C NN方法引入多波束声呐图像镶嵌领域,通过P C NN点火次数可有效获取图像的纹理等细节信息,实现多波束声呐图像的细节信息综合。P C NN是由若干神经元组成的单层二维神经网络,每一个神经元都与图像中的一个像素点对应。单个神经元由接 收 域、调 制 域 和 脉 冲 发 生 器 三 部 分 组成1 6,如图3所示。其中,F为反馈输入,L为连接输入,U为内部活动项,Y为输出部分,W为权阵,为连接强度,E为活动阈值。传统的P C NN模型存在参数众多、迭代计算较慢等缺陷,因此衍生了较多简化的P C NN模型。其中,S P C NN(s i m p l i f i e d P C NN)是由C h e n等1 7在S CM基础上提出的,并提供了参数的自动设置方法,S P C NN3山东科技大学学报(自然科学版)2 0 2 3年第4期模型的迭代过程为:Fi j(n)=Si j,(3)Li j(n)=VLk lNi jWi j,k lYk l(n-1),(4)Ui j(n)=e-FUi j(n-1)+Fi j(n)(1+Li j(n)。(5)Yi j(n)=1,Ui j(n)Ei j(n-1);0,其他。(6)Ei j(n)=e-EEi j(n-1)+VEYi j(n)。(7)式中:Fi j表示反馈输入项,Li j表示连接输入,Ui j表示内部活动项,Yi j表示脉冲输出,Ei j表示动态阈值,n为迭代次数,下标i、j、k、l均表示像素坐标位置,VL和VE分别为连接输入和动态阈值的幅值,F和E分别为反馈输入和动态阈值的时间衰减常数。式(3)将外部激励Si j作为反馈输入。式(5)为脉冲调制部分,通过将反馈输入和连接输入进行耦合,得到内部活动项Ui j,由于神经元的漏电积分特性,调制项Ui j随着每次迭代按照F进行指数衰减。将Ui j和Ei j进行比较,当Ui j大于阈值时,该神经元点火,即Yi j=1;否则神经元不点火,Yi j=0,得到的二维图Y称为点火图。如果神经元点火,动态阈值按VE跃升,如果不点火,则按照E进行指数衰减1 6。表1 S P C N N参数自动设置1 7T a b l e 1 A u t o m a t i c s e t t i n g o f S P C NN p a r a m e t e r s参数值VL1FF=l o g(1/(I)=(Sm a x/S)-1)/6VLVEVE=e-F+1+6VLEE=l nVE(S(1-e-3F(1-e-F)-1+6VL)-1S P C NN模 型 中 的5个 参 数 分 别 是VL、F、VE、E,均可进行自适应计算,如表1所示。其中,(I)表示图像的标准差,Sm a x和S 分别表示最大像素值和归一化O t s u阈值,其余参数同S P C NN模型一致。连接强度是S P C NN模型中一个重要的参数,反映了各神经元模型之间的连接强度,表征了某 一 神 经 元 受 其 邻 近 神 经 元 的 影 响 大小1 7。传统P C NN算法中,所有神经元均采用相同的连接强度,未考虑不同像素点的贡献差异。实验采用改进的拉普拉斯能量ML和拉普拉斯能量SML定义连接强度1 8:ML(i,j)=2I(i,j)-I(i-1,j)-I(i+1,j)+2I(i,j)-I(i,j-1)-I(i,j+1),(8)SML=1p=-11q=-1ML(i+p,j+q)2。(9)式中:i和j表示像素坐标,I为高频子带系数,拉普拉斯算子步长设置为1,ML值累加窗口大小为33。权阵W按照欧式距离倒数进行计算:W=0.7 0 7 1 0.7 0 71 0 1 0.7 0 7 1 0.7 0 7 。(1 0)将高频系数Hl,d1和Hl,d2作为S P C NN模型的反馈输入项F。分别对各层各方向的对应高频系数按照S P C NN模型进行计算,得到累计点火次数:Tl,d(n)=Tl,d(n-1)+Yl,d(n)。(1 1)式中:Yl,d(n)表示l层d方向的第n次点火次数;Tl,d(n)表示l层d方向的前n次迭代累计点火次数,即该像素点在迭代计算过程中被激发的总次数。以点火次数作为判决算子,按照点火次数的大小进行融合:Hl,d=Hl,d1,Tl,d1Tl,d2;Hl,d2,Tl,d1Tl,d2。(1 2)式中:Hl,d1和Hl,d2分别表示条带1与条带2的高频系数,Hl,d为融合后的高频系数。4曲 萌等:基于P C NN与混合加权的多波束声呐图像镶嵌方法1.4 N S S T逆变换经过低频信息融合和高频信息融合后,实现了两图像的信息综合,消除了低频信息中的拼接线并保留高频特征信息。最后将两幅图像融合后的低频系数L和一系列高频系数Hl,d进行N S S T逆变换,得到两幅图像镶嵌后的结果。2 实验与分析2.1 实验数据本研究采用多波束测深系统实测数据进行声呐图像镶嵌实验。实验数据采集于2 0 1 1年3月8日,来源于2 0 1 2年浅水调查会议资料1 9。调查区域位于新西兰惠灵顿附近海域,平均水深约2 0 m(图4),图4(b)中X、Y表示高斯投影后的平面坐标。使用K o n g s b e r g EM 2 0 4 0多波束测深系统以平均声强方式获取的实验数据进行验证。多波束测深系统的工作频率为3 0 0 k H z,波束宽度为1 1,波束开角为1 2 0,采样模式为等距采样。图4 实验区域概况F i g.4 E x p e r i m e n t a l a r e a2.2 评价指标为准确、客观地评价信息综合效果,选择信息熵、标准差、空间频率用于定量分析,从不同角度反映信息综合效果。上述指标值越大,表明图像的质量越高。信息熵作为一种图像质量评价指标,常被用来衡量图像中包含信息的丰富程度。信息熵越大,图像包含的信息量越多,计算公式为:IE=-Lk=0Pkl g(Pk)。(1 3)式中:IE表示信息熵,L表示待测图像像素灰度最大值,Pk表示像素值为k的像素个数。标准差通过对图像进行统计,得到像素的离散程度。标准差越大,图像的灰度值离散程度越高,视觉效果越好,计算公式为:SD=1nn1(xi-x)2。(1 4)式中:SD表示标准差,n表示像素个数,x表示像素值,x表示像素值均值。空间频率是基于梯度的图像质量指数,能够有效衡量图像的梯度分布,用来估计图像的纹理和细节2 0。SF越大,表明图像具有丰富的边缘和纹理,计算公式为2 0:SF=CF2+RF2。(1 5)其中:5山东科技大学学报(自然科学版)2 0 2 3年第4期RF=1MNMi=1Nj=2(F(i,j)-F(i,j-1)2,CF=1MNMi=2Nj=1(F(i,j)-F(i-1,j)2。(1 6)式中:RF与CF分别表示图像的空间列频率与行频率,M和N表示图像的行列数。2.3 实验结果与分析为验证本研究方法的有效性,对调查区域的7个条带数据进行镶嵌实验。首先对数据进行预处理,主要包括:对原始多波束反向散射强度进行解析,并结合导航、姿态等辅助数据对各条带进行波束归位,将各条带数据纳入统一的坐标框架下;对各条带数据进行传播损失、入射角效应和中央异常改正等各项声强误差处理,以削弱系统误差的影响;对图像进行旋转,使得航迹线方向大致与图像纵轴平行,利用图像蒙版操作提取相邻条带重叠区域,完成预处理。再使用本研究方法对各条带进行镶嵌,主要过程为:首先对各条带重叠区域图像进行N S S T分解,分解层数设置为2,各层方向数设置为1 6,得到高低频子带图像;对低频子带图像,沿着图像横轴方向,按照混合加权法进行拼接;对高频子带图像,按照P C NN进行融合,迭代次数为2 0 0,其余参数按照表1进行计算。对该调查区域的7个条带数据进行镶嵌,镶嵌前后的结果如图5所示。图5 实验区多波束声呐图像镶嵌结果F i g.5 M o s a i c r e s u l t s o f m u l t i b e a m s o n a r i m a g e i n e x p e r i m e n t a l a r e a从图5(a)中可以看出,尽管经过误差改正处理,条带间仍存在系统性残差,表现为明显的拼接线,该现象影响声呐图像的进一步使用。经本研究方法镶嵌后,如图5(b)所示,能有效削弱图像的灰度畸变,无明显拼接痕迹,且底质和纹理特征的边界清晰,能够较好地保留图像的纹理信息,实现多条带的无缝拼接。为了说明本研究方法的有效性,选取实验区域内具有代表性的相邻两条带进行分析,两条带重叠区域大、强度差异明显。采用平均法、混合加权法和本研究方法分别对条带重叠区域图像进行镶嵌处理,实验结果如图6所示。从图6(a)可以看出,原始数据存在的残差造成了重叠区域边界明显的拼接线。图6(b)中,经平均法镶嵌后,两条带总体声强较一致。但由于只对重叠区域数据进行平均化,削弱了部分随机噪声,导致重叠区域更加平滑,与非重叠区域间的纹理细节存在一定差异,因此镶嵌结果仍存在较明显的拼接线。图6(c)中,混合加权法总体上效果较好,条带间声强过渡均匀,无明显的拼接线。但在加权计算过程中,损失了地物地貌的细节信息,使图像存在一定程度的模糊。图6(d)中,本研究方法在低频系数上利用混合加权法消除了拼接线,并对高频系数进行信息提取,保留了细节信息,有效实现了多波束声呐图像的镶嵌。下面分别从拼接效果和信息综合效果两方面分析本研究方法的有效性。1)拼接效果分析6曲 萌等:基于P C NN与混合加权的多波束声呐图像镶嵌方法为了定量评价本研究方法在消除条带间拼接线方面的有效性,首先在声呐图像镶嵌结果图(图6)中沿垂直航迹线方向取一横向截面(黄色直线),该截面跨越了相邻条带的拼接线且底质类型基本一致;然后分别计算拼接线两侧一定范围内声强的统计特征,如表2所示。图6 多波束声呐图像镶嵌效果对比F i g.6 M o s a i c e f f e c t c o m p a r i s o n o f m u l t i b e a m s o n a r i m a g e s表2 拼接线两侧声强统计特征T a b l e 2 S t a t i s t i c a l c h a r a c t e r i s t i c s o f b a c k s c a t t e r i n t e n s i t y o n b o t h s i d e s o f t h e s e a m l i n ed B声强接缝左侧均值接缝右侧均值接缝左侧标准差接缝右侧标准差原图像-3 1.8 7-2 5.6 92.0 51.2 9平均法-3 2.1 4-2 9.3 31.3 61.2 6混合加权法-3 2.1 4-3 0.9 91.3 61.2 1本研究方法-3 2.1 4-3 0.8 51.3 61.3 6图7显示了截面处不同镶嵌方法的拼接效果。红色实线与黄色实线分别表示拼接线两侧的声强均值。在图7(a)中,原始图像条带之间在重叠区域存在明显的声强差异,约为-6.1 8 d B。在图7(b)平均法处理结果中,拼接线两侧均值相差2.8 1 d B,条带边界处呈明显的阶梯状。此外,如表2所示,平均法处理结果中拼接线两侧的声强标准差相差0.1 d B,表明重叠区经过平均化处理后,削弱了随机噪声影响,与非重叠区域存在一定差异。因此,在镶嵌结果图6(b)中有较明显的拼接痕迹。在图7(c)混合加权方法中,拼接线两侧均值相差1.1 5 d B,比平均法降低了5 9%,表明该方法能够有效削弱拼接线两侧声强差异,实现声强的均匀过渡。该方法拼接线两侧声强标准差相差0.1 5 d B,说明加权处理削弱了随机噪声的影响。在图7(d)本研究方法中,拼接线两侧声强均值相差1.2 9 d B,总体上无明显阶跃,表明对N S S T分解后的低频系数进行混合7山东科技大学学报(自然科学版)2 0 2 3年第4期加权处理能实现声强的无缝拼接。此外,如表2所示,本研究方法拼接线两侧声强标准差基本一致,表明对N S S T分解后的高频系数进行P C NN融合能保持声呐图像的原始细节信息,避免加权处理弱化图像的纹理信息。图7 拼接线两侧声强比较F i g.7 C o m p a r i s o n o f b a c k s c a t t e r i n t e n s i t y o n b o t h s i d e s o f t h e s e a m l i n e2)信息综合效果分析通过对声呐图像进行镶嵌,实现了相邻条带重叠区域的信息融合。为定量评价3种方法在信息综合方面的有效性,对图像的信息熵IE、标准差SD和空间频率SF进行计算,结果如表3所示。表3 声呐图像信息综合效果评价T a b l e 3 E v a l u a t i o n o f f u s i o n r e s u l t s评价指标平均法混合加权法本研究方法SD6 2.2 96 2.2 36 2.8 1IE3.1 73.1 23.2 5SF1 8.8 41 7.2 92 2.7 6从表3可以看出,3种镶嵌方法的标准差相近,本研究方法实验结果标准差稍高,与平均法和混合加权法相比,分别提升0.8%和0.9%。本研究方法实验结果的信息熵最高(3.2 5),平均法次之,混合加权法最低。主要原因是混合加权法在计算过程中未考虑具体区域的纹理信息,仅按照像素点位置加权,可能会损失部分纹理信息。与其他两种方法相比,本研究方法的信息熵分别提升2.5%和4.2%,表明该方法在保持信息丰富程度上具有一定优势。本研究方法的空间频率为2 2.7 6,分别比平均法的1 8.8 4、混合加权法的1 7.2 9提升了2 0.8%和3 1.6%,表明镶嵌结果图中具有丰富的边缘信息,图像质量更高。从图6中的源图像和拼接后声呐图像中提取子区域(黑色矩形),如图8所示。从图8(a)和8(b)中可以看出,条带之间的纹理信息差异明显,不同条带的重叠区域的细节信息不同,且呈较强的互补性。主要原因是采用多波束声呐测量同一区域时,不同条带的波束入射角和方位角不同,导致返回的信号存在差异。经过8曲 萌等:基于P C NN与混合加权的多波束声呐图像镶嵌方法镶嵌处理后,平均法(图8(c)和混合加权法(图8(d)的镶嵌结果均存在一定的模糊。本研究方法的镶嵌结果(图8(e)纹理特征明显、边界清晰,能够保持更丰富的细节信息。图8 声呐图像镶嵌信息综合结果F i g.8 I n f o r m a t i o n f u s i o n r e s u l t s o f s o n a r i m a g e s m o s a i c3 结论本研究针对多波束声呐图像镶嵌中存在的拼接线消除与信息综合问题,提出一种在N S S T域内结合P C NN和混合加权法的图像镶嵌方法,并利用实测多波束海底反向散射强度数据进行实验分析。实验结果表明,本研究方法能够有效削弱条带间的拼接线并消除灰度畸变,同时能够保留较为全面的纹理细节信息,可为声呐图像的镶嵌任务提供参考。但当重叠区域较窄、海底底质较为单一且无明显纹理特征时,该方法存在一定局限性,需要进一步研究。参考文献:1 HE R KU L K,P E T E R S ON A,P A E K I V I S.A p p l y i n g m u l t i b e a m s o n a r a n d m a t h e m a t i c a l m o d e l i n g f o r m a p p i n g s e a b e d s u b-s t r a t e a n d b i o t a o f o f f s h o r e s h a l l o w sJ.E s t u a r i n e,C o a s t a l a n d S h e l f S c i e n c e,2 0 1 7,1 9 2:5 7-7 1.2 阳凡林,朱正任,李家彪,等.利用深层卷积神经网络实现地形辅助的多波束海底底质分类J.测绘学报,2 0 2 1,5 0(1):7 1-8 4.YAN G F a n l i n,Z HU Z h e n g r e n,L I J i a b i a o,e t a l.S e a f l o o r c l a s s i f i c a t i o n b a s e d o n c o m b i n e d m u l t i b e a m b a t h y m e t r y a n d b a c k-s c a t t e r u s i n g d e e p c o n v o l u t i o n n e u r a l n e t w o r kJ.A c t a G e o d a e t i c a e t C a r t o g r a p h i c a S i n i c a,2 0 2 1,5 0(1):7 1-8 4.3 韩冰,冯成凯,印萍,等.基于小批量K均值预分类的多波束反向散射强度角度影响改正J.山东科技大学学报(自然科学版),2 0 2 0,3 9(3):3 3-4 2.HAN B i n g,F E NG C h e n g k a i,Y I N P i n g,e t a l.C o r r e c t i o n o f m u l t i-b e a m b a c k s c a t t e r i n t e n s i t y a n g l e b a s e d o n M i n i B a t c h K-m e a n s p r e-c l a s s i f i c a t i o nJ.J o u r n a l o f S h a n d o n g U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y(N a t u r a l S c i e n c e),2 0 2 0,3 9(3):3 3-4 2.4 S CH I ME L A C G,B E AU D O I N J,P A R NUM I M,e t a l.M u l t i b e a m s o n a r b a c k s c a t t e r d a t a p r o c e s s i n gJ.M a r i n e G e o p h y s i-c a l R e s e a r c h,2 0 1 8,3 9:1 2 1-1 3 7.5 黄樟钦,母召,岑陈,等.一种多路视频的3 6 0度实时拼接技术J.西安电子科技大学学报,2 0 2 0,4 7(2):1-8.HUAN G Z h a n g q i n,MU Z h a o,C E N C h e n,e t a l.3 6 0-d e g r e e r e a l-t i m e s t i t c h i n g t e c h n o l o g y f o r m u l t i-c h a n n e l v i d e oJ.J o u r-n a l o f X i d i a n U n