一种非均匀声学成像方法技术

本发明专利技术涉及一种非均匀声学成像方法，生成初始成像阵列，以初始麦克风阵列进行初步采集，获取初始成像阵列中每一阵元至麦克风阵列中每一阵元的时延，根据时延生成时延矩阵，根据时延矩阵生成非均匀的优化阵列流行，根据优化阵列流行进行非均匀采样，得到声学成像。本发明专利技术的方法可以在不增加扫描阵元数量的前提下，提高声源处的扫描阵元密度，提高成像精度；且使扫描对准声源，减少了非声源区域的扫描，抑制旁瓣干扰产生，进一步提高了声源成像的精度。度。度。

【技术实现步骤摘要】
一种非均匀声学成像方法


[0001]本专利技术属于声源成像
，具体涉及一种非均匀声学成像方法。

技术介绍

[0002]在声学检测中，除了检测声源位置以外还需要检测声源物理特征。通过对声源物理特征进行分析，可以确定声源的物理状态、声源的大致形状，并且能够更加精确的定位声源位置，同时知晓该声源到底是有效声源还是反射干扰。
[0003]对于声源物理特征的分析，一般使用声学成像方法分析得到声源的图像。在现有技术中，声学成像使用由多个麦克风阵元组成的麦克风阵列扫描存在声源的一个空间区域，扫描的空间区域是由多个目标阵元组成的目标阵列，即由麦克风阵列中的多个阵元共同对组成目标区域的多个目标阵元进行扫描。
[0004]最常见的目标阵列构型为平面矩阵或者圆形矩阵，在这种目标阵列的构型下，如果想要得到声源的更详细分布，一般会将方阵的密度增大，即整个目标阵列中的目标阵元数量，从而以更密集的扫描得到清晰图像。
[0005]但上述方法将目标阵元数量提升以后不仅会大幅提升扫描及波束形成过程消耗的计算量，由于声源的主体部分仅占用目标区域中的部分空间，增加整体密度同时也会增加不重要区域的扫描密度，从而引入额外的旁瓣干扰，降低声源成像清晰度。
[0006]因此需要一种声学成像方法，不仅能够得到以较高精度对声源进行扫描，同时能够抑制旁瓣的干扰，进一步提高声源的成像精度。

技术实现思路

[0007]基于现有技术中存在的上述缺点和不足，本专利技术的目的之一是至少解决现有技术中存在的上述问题之一或多个，换言之，本专利技术的目的之一是提供满足前述需求之一或多个的一种非均匀声学成像方法。
[0008]为了达到上述专利技术目的，本专利技术采用以下技术方案：一种非均匀声学成像方法，具体包括如下步骤：S1、生成初始成像阵列；S2、以初始麦克风阵列进行初步采集，获取初始成像阵列中每一阵元至麦克风阵列中每一阵元的时延，根据时延生成时延矩阵；S3、根据时延矩阵生成非均匀的优化阵列流行；S4、根据优化阵列流行进行非均匀采样，得到声学成像。
[0009]作为一种优选方案，步骤S3具体包括如下步骤：S31、选取时延矩阵中时延最小的若干个点；S32、获取时延最小的若干个点相对应的初始成像阵列中阵元；S33、使用时延最小的若干个点相对应的初始成像阵列中阵元生成非均匀的优化阵列流行。
[0010]作为一种优选方案，步骤S4具体包括如下步骤：S41、根据优化阵列流行进行非均匀采样；S42、使用波束成形算法对非均匀采样的结果进行声学成像。
[0011]作为一种优选方案，波达方向定位算法使用波束成形算法。
[0012]作为一种优选方案，步骤S3之后、S4之前还包括步骤：S310、对优化阵列流行进行离散，使重合的点分离。
[0013]作为一种进一步优选的方案，离散使用高斯函数进行。
[0014]作为一种优选方案，步骤S4之后还包括步骤：S410、对声学成像插值运算，得到整体分辨率统一的声学成像。
[0015]本专利技术与现有技术相比，有益效果是：本专利技术的方法可以在不增加扫描阵元数量的前提下，提高声源处的扫描阵元密度，提高成像精度；且使扫描对准声源，减少了非声源区域的扫描，抑制旁瓣干扰产生，进一步提高了声源成像的精度。
附图说明
[0016]图1是本申请实施例的一种初始成像阵列的排列示意图；图2是本申请实施例的一种非均匀的优化流行阵列的排列示意图；图3是本申请实施例的离散后的非均匀的优化流行阵列的排列示意图；图4是本申请实施例的非均匀声学成像方法所成的三维声学云图；图5是本申请实施例的非均匀声学成像方法所成的二维声学云图。
具体实施方式
[0017]下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。
[0018]在下述介绍中提供了本申请的多个实施例，不同实施例之间可以替换或者合并组合，因此本申请也可认为包含所记载的相同和/或不同实施例的所有可能组合。因而，如果一个实施例包含特征A、B、C，另一个实施例包含特征B、D，那么本申请也应视为包括含有A、B、C、D的一个或多个所有其他可能的组合的实施例，尽管该实施例可能并未在以下内容中有明确的文字记载。
[0019]下面的描述提供了示例，并且不对权利要求书中阐述的范围、适用性或示例进行限制。可以在不脱离本申请内容的范围的情况下，对描述的元素的功能和布置做出改变。各个示例可以适当省略、替代或添加各种过程或组件。例如所描述的方法可以以所描述的顺序不同的顺序来执行，并且可以添加、省略或组合各种步骤。此外，可以将关于一些示例描述的特征组合到其他示例中。
[0020]一种非均匀声学成像方法，具体包括如下步骤：S1、生成初始成像阵列M，M中的每一个阵元记录了该阵元的位置。
[0021]初始成像阵列M用于初步检测声源形状及分布，由于声源的分布尚未知晓，在初始成像阵列M中必然仅有部分区域为声源区域，其余区域只有少量声音或无声源，如果直接使用初始成像阵列进行扫描，非声源区域会产生旁瓣干扰，影响最终的成像清晰度。
[0022]另外，由于声学成像是通过联合初始成像阵列中的阵元及麦克风阵列中的阵元进行扫描成像的，在相同的阵列面积下，如果想得到更高的成像清晰度需要增大阵元密度，而更高的阵元密度会极大增加两个阵列中阵元的总配对数量，从而提高算力消耗。
[0023]因此，本申请的初始成像阵列不直接用于扫描，而是先修改、使其变为更符合声源物理形态的非均匀构型。
[0024]在本申请的某些实施例中，初始成像阵列可以如图1所示，使用最简单的均匀矩形阵列；而在本申请的另一些实施例中，初始成像阵列可以使用圆形均匀阵列等多种构型，甚至可以使用非均匀的初始成像阵列，根据成像前预估的声源位置适当将阵元集中于声源所在区域附近。
[0025]生成初始成像阵列后执行步骤S2、以初始麦克风阵列进行初步采集，获取初始成像阵列中每一阵元至麦克风阵列中每一阵元的时延，根据时延生成时延矩阵。
[0026]本申请的实施例在初始成像阵列以外，需要配置一个实体的麦克风阵列，麦克风阵列用于对初始成像阵列中的阵元进行扫描以获取成像数据。
[0027]具体的，步骤S2通过获取初始成像阵列中每一个阵元的位置，同时获取用于扫描声源成像的麦克风阵列中每一个阵元的位置，从而将初始成像阵列中的阵元与麦克风阵列中的阵元依次配对，枚举两种阵元的全部配对组合。
[0028]然后，依次选择每一对初始成像阵列与麦克风阵列中阵元的配对组合，计算该对阵元间的声波传输时延，将每个配对组合的时延记录以后，整合生成时延矩阵。
[0029]时延矩阵记录了任意一对阵元间的声波传输时延，作为一种具体的实施方法，根据M矩阵计算时延矩阵的公式为：，其中ω为角频率（2*π*f / fs），f为频率范围（1~fs），fs为音频采样频率；c为声速。
[0030]得到每对阵元间的时延矩阵后，进行步骤S3、根据时延矩阵生成非均匀的优化阵列流行。
...

【技术保护点】

【技术特征摘要】
1.一种非均匀声学成像方法，其特征在于，具体包括如下步骤：S1、生成初始成像阵列；S2、获取所述初始成像阵列中每一阵元至麦克风阵列中每一阵元的时延，生成时延矩阵；S3、根据所述时延矩阵生成非均匀的优化阵列流行；S4、根据所述优化阵列流行进行非均匀采样，得到声学成像。2.如权利要求1所述的一种非均匀声学成像方法，其特征在于，所述步骤S3具体包括如下步骤：S31、选取所述时延矩阵中时延最小的若干个点；S32、获取所述时延最小的若干个点相对应的所述初始成像阵列中阵元；S33、使用所述时延最小的若干个点相对应的所述初始成像阵列中阵元生成所述非均匀的优化阵列流行。3.如权利要求1所述的一种非均匀声学成像方法，其特征在于，所述步...

【专利技术属性】
技术研发人员：曹祖杨，黄明，闫昱甫，陶慧芳，张永全，梁友贵，李佳罗，
申请(专利权)人：杭州兆华电子股份有限公司，
类型：发明
国别省市：