一种基于图像匹配的吊放式声纳方位修正方法
【专利摘要】本发明公开了一种基于图像匹配的吊放式声纳方位修正方法，其特征在于，包括：步骤1、对声纳图像进行图像栅格处理：步骤2、对图像栅格处理的输出进行二值分割；步骤3、选择参考模板，所述参考模板为声纳图像时间序列上的任意一帧，或者对时间序列上连续多帧进行统计后获得；步骤4、根据相似性测度原则，计算二值分割后图像和参考模板的相关矩阵；步骤5、根据相关矩阵判定当前图像的偏移信息；步骤6、根据所述偏移信息对二值分割后图像进行修正。本发明的方法对吊放式声纳方位进行有效修正。
【专利说明】
一种基于图像匹配的吊放式声纳方位修正方法
技术领域
[0001] 本发明涉及水下主动探测技术领域，具体涉及一种基于图像匹配的吊放式声纳方 位修正方法。
【背景技术】
[0002] 吊放式声纳广泛应用于直升机、舰船和码头，可用于监视、侦测海洋环境、海洋生 物和水下人类活动。由于平台和成本限制，吊放式声纳往往没有安装指向罗经和姿态伺服 器，受涌浪影响，声纳在水下必然存在旋转和摆动。缺少姿态补偿时，声纳图像将存在左右 抖动的现象。通过图像匹配技术，获取当前图像航向偏移量，可有效修正吊放式声纳图像左 右抖动的现象。
[0003] 时效性和匹配精度是图像匹配算法的两个重要指标。目前的图像匹配方法基本上 分3类：基于区域的匹配方法，基于特征的匹配方法，基于解释的匹配方法。基于区域的方 法中，交叉相关方法时效性较高，但匹配精度较差；互信息交互方差的方法匹配精度较高， 时效性较差。基于特征的方法，时效性较高，但匹配精度依赖于特征提取的质量。基于解释 的图像匹配方法建立在对图像正确解译的基础上，从而对于图像的质量和先验知识的要求 较高。考虑吊放式声纳除了航向变化，同时存在纵横摇变化，其导致波束照射至水底或水 面，引起声纳图像混响背景剧烈变化。结合不同方法的优势，本文考虑使用基于二值分割的 图像特征匹配方法，并通过图像分割，保证较高的匹配精度和时效性。
【发明内容】
[0004] 本发明要解决的技术问题是现有吊放式声纳因水下姿态不稳定引起的声纳图像 左右抖动的问题。为了解决上述问题本发明提供一种基于图像匹配的吊放式声纳方位修正 方法。
[0005] 本发明提供的基于图像匹配的吊放式声纳方位修正方法，包括： 步骤1、对声纳图像进行图像栅格处理：
々（·）为图像分害U函数，为声纳图像以（r，c)为左上点，尺寸为#乂成勺第i 行，第J列区域子图，为经图像分割处理后的图像的第i行、第^/列像素值； 步骤2、对图像栅格处理的输出进行二值分割：
为二值分割后图像的第i行，第^/列像素值，出为二值化判决域值，皿'<3/7为均值函 数，为方差函数，J为修正值，J的取值大小在〇与声纳图像信噪比之间，信噪比越大，乂 的取值越大； 步骤3、选择参考模板，所述参考模板为声纳图像时间序列上的任意一帧，或者对时间 序列上连续多帧进行统计后获得； 步骤4、根据相似性测度原则，计算二值分割后图像和参考模板的相关矩阵； 步骤5、根据相关矩阵判定当前图像的偏移信息； 步骤6、根据所述偏移信息对二值分割后图像进行修正。
[0006] 进一步，所述参考模板为对声纳图像时间序列上连续多帧进行统计后获得的情况 下，参考模板第i行，第j列的像素值为+ +』；I为时间序列上任 意帧的第i行，第j列的像素值。
[0007] 进一步，二值分割后图像和参考模板的相关矩阵的计算公式为
;f为模板图像，为模板的第ar i行、 第U列像素值；c为模板匹配输出相关矩阵，4"为相关矩阵的第4于、第λ列相关值； (ck, cr, cc)含义为c的最大值为ck，ck在c的第cr行、第cc列；a?ax为取最大值。
[0008] 进一步，a。为图像序列中的某一帧图像，当f e a。时，即图像模板f为图像a。中 以U，乂）为左上点的尺寸为#ΧΛ韵区域，a。尺寸为fX//，极限情况可为Z/= a。，则： cr^ ~i0 +M
(CK。，cr。，cc。）表示图像a。与其区域A的自相关矩阵c。的最大值为CK。，CK。在相关矩 阵c。的第cr。行，第。列；依据相关原则，当ck。满足域值Z?r(·)时，图像偏移量为：
ra?为实时图像a行偏移量，ca?为实时图像a列偏移量。
[0009] 依据权利要求5所述的基于图像匹配的吊放式声纳方位修正方法，其特征在于， 修正后图像的第i行，第万!J像素值为：<
[0010] 本发明相比现有技术的优点在于： (1)本发明可提高设备的适装性，在同等条件下，不需要安装指向罗经和姿态伺服器， 减小设备研制成本。
[0011] (2)本发明具有较高的时效性和精度，对于背景缓慢变化的图像方位修正的准确 率达到100%，对于背景剧烈变化的图像方位修正的准确率优于80%。
【附图说明】
[0012] 图1是本发明的处理流程图。
[0013] 图2是本发明的原理仿真结果示意图，（a)为模版、（b)为实时图像，（c)为相关中 心，相关中心为（40, 30)，（d)为修正后的图像。
[0014] 图3是本发明的实际数据人为偏移仿真结果示意图，（a)为模版，（b)为实时图像， (C)为相关中心，相关中心为（10,0)，（d)为修正后的图像； 图4是本发明的实际数据稳态背景单一模板处理结果示意图；（a)为模版、（b)为第56 帧图像，（c)为相关中心，相关中心为（-4,0)，（d)为第56帧图像修正后的图像； 图5是本发明的实际数据非稳态背景单一模板处理结果示意图；（a)为模版、（b)为第 1412帧图像，（c)为相关中心，相关中心为（1，0)，（d)为第1412帧图像修正后的图像； 图6是本发明的实际数据非稳态背景时变模板处理结果示意图；（a)为模版、（b)为实 时图像，（c)为相关中心，相关中心为（-3,0)，（d)为修正后的图像。
【具体实施方式】
[0015] 下文中，结合附图和实施例对本发明做进一步阐述。
[0016] 参考图1，本发明所提供的基于图像匹配的吊放式声纳方位修正方法包括： 步骤1，获取二维声纳图像序列，对任一声纳图像进行图像栅格处理： 8ij = (1) Α(·)为图像分害U函数，为声纳图像以（r，c)为左上点，尺寸为成勺第i 行，第J列区域子图，为/L(r，c)经图像分割后图像的第i行、第^/列像素值。
[0017] 步骤2,对图像栅格处理后的输出进行二值分割：
为二值分割后图像的第i行，第^/列像素值，出为二值化判决域值，皿'<3/7为均值函 数，为方差函数，I为修正值d的取值大小在〇与声纳图像信噪比之间，信噪比越大，乂 的取值越大，本领域技术人员可以根据实际情况选择的取值。
[0018] 步骤3,参考模板可选择声纳图像时间序列上的任意一帧，亦可对时间序列上连续 多帧进行统计后获得：
为连续多帧的第i行，第万!J像素的统计值，i::为第二修正值。
[0019] 步骤4,根据相似性测度原则，计算当前图像和参考模板的相关矩阵：
f为模板图像，c为模板匹配输出相关矩阵，4"为相关矩阵的第4于、第λ列相关值， 为模板的第ari行、第/7~7列像素值；（ck, cc)含义为c的最大值为ck, ck在c 的第cr行、第cc列，a?ax取最大值。
[0020] 步骤5,搜索相关矩阵峰值，根据相关矩阵来判定当前图像的偏移信息： a。为图像序列中的某一帧图像，当时，即模板//为图像a。中以（iQ，j)为左上 点的尺寸为#ΧΛ韵区域，a。尺寸为fX//，极限情况可为Z/=a。，则：
(1^。，^^。，％。）为图像<3。与其区域6的自相关矩阵1?。的最大值1^。在第 1^。行，第1^。 列。依据相关原则，当CK。满足域值挪(·）时，图像偏移量为：
ra?为实时图像a行偏移量，ca?为实时图像a列偏移量。
[0021] 步骤6,图像修正：
a ^为修正后图像的第i行，第^/列像素值。
[0022] 利用本发明对某型吊放式反蛙人声纳进行图像抖动修正，图像矩阵为192个波束 (360° )，1000个距离象素（400米)，hl(·)选用中值滤波器实现，处理效果见图1~图5所 示，性能统计如表1所示。
[0023] 表1多组实际数据处理性能统计
由表1可以看出，本发明所提供的方法可以对声纳抖动引起图像变形进行修正，修正 的准确率满足应用需求。
[0024] 图2中（a)为选用的模板，（b)为实时图像，（c)为相关矩阵峰值示意图，（d)为修 正后的图像，从图2可以看出经过本发明所提供的方法相关矩阵峰值得到有效修正。
[0025] 图3中，先选择声纳图像时间序列上的任意一帧图像（a)为模板，然后将模板图 像左偏10度，得到实时图像（b)，再采用本发明所提供的方法进行修正，得到修正后的图像 (d)，比较（a)和（d)，可以看出本发明所提供的方法可以将左偏离修正。
[0026] 图4为本发明所提供的方法应用于稳态背景的示意图。
[0027] 图5和图6为本发明所提供的方法应用于非稳态背景的示意图。
[0028] 本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域 技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发 明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明 的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案 的保护范围。
【主权项】
1. 一种基于图像匹配的吊放式声纳方位修正方法，其特征在于，包括： 步骤1、对声纳图像进行图像栅格处理： Α(·)为图像分害U函数，tUc)为戸納图像以Uc)为左上点，尺寸为#乂成勺第i 行，第J列区域子图，为/^(r，c)经图像分割处理后的图像的第i行、第^/列像素值； 步骤2、对图像栅格处理的输出讲行二倌分割：为二值分割后图像的第i行，第^/列像素值，出为二值化判决域值，皿'<3/7为均值函 数，为方差函数，J为修正值，J的取值大小在0与声纳图像信噪比之间，信噪比越大，乂 的取值越大； 步骤3、选择参考模板，所述参考模板为声纳图像时间序列上的任意一帧，或者对时间 序列上连续多帧进行统计后获得； 步骤4、根据相似性测度原则，计算二值分割后图像和参考模板的相关矩阵； 步骤5、根据相关矩阵判定当前图像的偏移信息； 步骤6、根据所述偏移信息对二值分割后图像进行修正。2. 依据权利要求1所述的基于图像匹配的吊放式声纳方位修正方法，其特征在于， 所述参考模板为对声纳图像时间序列上连续多帧进行统计后获得的情况下，参考模板第i 行，第j列的像素值为％时间序列上任意帧的第i行，第 j列的像素值。3. 依据权利要求2所述的基于图像匹配的吊放式声纳方位修正方法，其特征在 于，二值分割后图像和参考模板的相关矩阵的计算公式为= 为模板图像，^为模板的第ari行、第列像素值；c为模板 匹配输出相关矩阵，4"为相关矩阵的第4于、第λ列相关值；（ck，cr，cc)含义为c的最大 值为CK，CK在c的第cr行、第cc列；a?ax为取最大值。4. 依据权利要求4所述的基于图像匹配的吊放式声纳方位修正方法，其特征在于，a。 为图像序列中的某一帧图像，当时，即图像模板//为图像 a。中以（i。，乂）为左上点 的尺寸为的区域，a。尺寸为fX"，极限情况可为f=a。，则：(ck。，cr。，cc。）表示图像a。与其区域A的自相关矩阵c。的最大值为ck。，ck。在相关矩 阵c。的第cr。行，第。列；依据相关原则，当ck。满足域值Z?r(·)时，图像偏移量为：ra?为实时图像a行偏移量，ca?为实时图像a列偏移量。5.依据权利要求5所述的基于图像匹配的吊放式声纳方位修正方法，其特征在于，修 正后图像的第i行，第J列像素值为： 。麵峰扣·:。
【文档编号】G06T5/00GK105869138SQ201510061565
【公开日】2016年8月17日
【申请日】2015年2月6日
【发明人】姚新, 毛英磊, 张磊磊, 郜扬文
【申请人】中国船舶重工集团公司第七二六研究所