[go: up one dir, main page]

CN110477951B - 基于宽频带声学超材料的超快复合平面波成像方法 - Google Patents

基于宽频带声学超材料的超快复合平面波成像方法 Download PDF

Info

Publication number
CN110477951B
CN110477951B CN201910813822.4A CN201910813822A CN110477951B CN 110477951 B CN110477951 B CN 110477951B CN 201910813822 A CN201910813822 A CN 201910813822A CN 110477951 B CN110477951 B CN 110477951B
Authority
CN
China
Prior art keywords
preset
signal
transmitting
receiving
frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910813822.4A
Other languages
English (en)
Other versions
CN110477951A (zh
Inventor
郑音飞
李超
蒋东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201910813822.4A priority Critical patent/CN110477951B/zh
Publication of CN110477951A publication Critical patent/CN110477951A/zh
Priority to US16/864,584 priority patent/US11829872B2/en
Application granted granted Critical
Publication of CN110477951B publication Critical patent/CN110477951B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995—Combining images from different aspect angles, e.g. spatial compounding
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/5205—Means for monitoring or calibrating
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52079—Constructional features
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Theoretical Computer Science (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了基于宽频带声学超材料的超快复合平面波成像方法。该方法通过超快复合平面波成像装置实现；该装置包括发射接收超声探头和声学超材料结构。该方法包括：控制发射接收超声探头以预设发射频率、第一预设发射角度发出超声波信号；预设发射频率等于声学超材料结构的响应频率；控制发射接收超声探头以预设接收频率，且分别以第一预设接收角度、第二预设接收角度和第三预设接收角度接收被测物体反射的回波信号；预设接收频率为预设发射频率的n倍；第一预设接收角度等于第一预设发射角度，第二预设接收角度小于第一预设发射角度，第三预设接收角度大于第一预设发射角度；采用回波信号重建被测物体的图像。本发明能够提高成像深度与成像质量。

Description

基于宽频带声学超材料的超快复合平面波成像方法
技术领域
本发明涉及超快平面波成像技术领域，特别是涉及基于宽频带声学超材料的超快复合平面波成像方法。
背景技术
目前，对医学超声图像的评判标准主要集中在成像深度及成像质量两方面。
在成像质量方面，超快平面波的出现革新了医学成像领域，成像帧频可达上千，极大地提高了常规超声成像的帧频。对比于常规超声成像的线到线的聚焦成像方式，超快平面波成像采用的是面到面的非聚焦成像方式。非聚焦的发射/接收模式是超快平面波成像的关键。2002年，Tanter等人基于超快平面波成像测量了横波速度，提出了瞬时弹性成像；2003年Bercoff团队将瞬时弹性成像首次用于乳腺癌的临床研究，也是超快平面波的首例临床应用。弹性成像的原则是对组织位移或横波速度的测量，因此牺牲图像质量(对比度和分辨率)来提高帧频的代价是可接受的。为提高图像质量，2009年Montaldo等人提出相干复合平面波成像方法，既不损失帧频的前提下，又可大幅提高成像质量。这种方法是通过增加不同发射角度的平面波数量来提高图像质量，是对成像帧频和成像质量的一种权衡。目前相干复合成像方法已成为超快成像的核心，被广泛应用于各种医学成像场景，尤其是血流成像上。2015年Tanter团队拓展了复合成像方法，提出了一种既不需要妥协帧频，同时又可以提高图像信噪比的超快多波成像方法，但是此方法计算量大、耗时长的缺陷限制了其临床应用。
在成像深度方面，低频发射序列可以带来高穿透性，但是其成像质量较差。近年来，对声学超材料的研究成为一个比较有前景的研究方向。2014年，Nicholas Fang等人利用comsol仿真设计了具有负质量密度和负体积模量的双负声学超材料，为超声成像提供了新的思路和方法。2015年，Thomas团队基于Mie共振效应实现了具有双负性质的声学超材料。但是，声学超材料的研制仍处于低频阶段，对高频超材料的研究处于摸索阶段。
综上，现有的超快平面波成像方法是基于相干复合平面波成像的，由于平面波的非聚焦特点，导致能量丢失并影响成像深度与成像质量。
发明内容
基于此，有必要提供一种基于宽频带声学超材料的超快复合平面波成像方法，以提高平面波成像的成像深度与成像质量。
为实现上述目的，本发明提供了如下方案：
基于宽频带声学超材料的超快复合平面波成像方法，所述超快复合平面波成像方法通过超快复合平面波成像装置实现；所述超快复合平面波成像装置包括发射接收超声探头和声学超材料结构；所述发射接收超声探头发出的超声波信号通过所述声学超材料结构后到达被测物体，所述发射接收超声探头还接收由所述被测物体反射的回波信号；
所述超快复合平面波成像方法，具体包括：
控制所述发射接收超声探头以预设发射频率，且以第一预设发射角度发出超声波信号，所述超声波信号通过所述声学超材料结构后到达被测物体；所述预设发射频率与所述声学超材料结构的响应频率相等；
控制所述发射接收超声探头以预设接收频率，且分别以第一预设接收角度、第二预设接收角度和第三预设接收角度接收所述被测物体反射的回波信号；所述预设接收频率为所述预设发射频率的n倍，n＞1；所述第一预设接收角度等于所述第一预设发射角度，所述第二预设接收角度小于所述第一预设发射角度，所述第三预设接收角度大于所述第一预设发射角度；
采用所述回波信号重建所述被测物体的图像。
可选的，所述采用所述回波信号重建所述被测物体的图像，具体包括：
采用深度学习算法去除所述回波信号中的伪影信号，得到回波优化信号；
采用所述回波优化信号重建所述被测物体的图像。
可选的，所述声学超材料结构的响应频率通过响应频率确定装置确定；所述响应频率确定装置包括发射探头、接收探头、与所述发射探头连接的信号发生器以及与所述接收探头连接的示波器；所述发射探头的输出端通过所述声学超材料结构与所述接收探头的输入端连接；所述信号发生器与所述示波器连接；
所述声学超材料结构的响应频率的确定方法为：
所述信号发生器在预设频率范围内以预设步长发出多个不同频率的电信号，每个频率的电信号均分为两路，第一路信号直接在示波器上显示，第二路信号依次经过所述发射探头、所述声学超材料结构和所述接收探头后，在所述示波器上显示；
比较每个频率下，在所述示波器上显示的两路电信号的波形，确定每个频率对应的第二路信号的波形的畸变程度、放大程度和延迟程度；
将畸变程度最小、放大程度最大，且延迟程度最小的第二路信号的波形对应的信号发生器的频率确定为声学超材料结构的响应频率。
可选的，所述采用深度学习算法去除所述回波信号中的伪影信号，得到回波优化信号，具体包括：
获取训练数据集；所述训练数据集包括含伪影信号的超声训练信号和去除伪影信号的超声训练信号；
构建三层卷积神经网络；
将所述含伪影信号的超声训练信号作为输入，所述去除伪影信号的超声训练信号作为输出，对所述三层卷积神经网络进行训练，得到训练好的三层卷积神经网络；
将所述回波信号输入至所述训练好的三层卷积神经网络中，得到回波优化信号。
可选的，所述采用所述回波优化信号重建所述被测物体的图像，具体包括：
根据所述回波优化信号，采用波束合成方法重建所述被测物体的图像。
可选的，1＜n＜3。
可选的，所述第一预设接收角度与所述第一预设发射角度均为A°；所述第二预设接收角度为A-δ°，所述第三预设接收角度为A+δ°，其中δ＞0。
可选的，5＜δ＜7。
可选的，所述预设频率范围为[0.5MHz,15MHz]；所述预设步长为1MHz。
与现有技术相比，本发明的有益效果是：
本发明提出了基于宽频带声学超材料的超快复合平面波成像方法。该方法通过超快复合平面波成像装置实现；该装置包括发射接收超声探头和声学超材料结构。该方法包括：控制发射接收超声探头以预设发射频率、第一预设发射角度发出超声波信号；预设发射频率等于声学超材料结构的响应频率；控制发射接收超声探头以预设接收频率，且分别以第一预设接收角度、第二预设接收角度和第三预设接收角度接收被测物体反射的回波信号；预设接收频率为预设发射频率的n倍；第一预设接收角度等于第一预设发射角度，第二预设接收角度小于第一预设发射角度，第三预设接收角度大于第一预设发射角度；采用回波信号重建被测物体的图像。本发明设置预设发射频率与声学超材料结构的响应频率相等，使得入射能量增强，提高了成像深度；提出单角度发射-三角度接收的超快复合成像方法，使得回波信息增多，并且通过使用单频发射-n倍频接收的方法，放大了接收信号能量，提高了成像质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例基于宽频带声学超材料的超快复合平面波成像方法的流程图；
图2为本发明实施例基于宽频带声学超材料的超快复合平面波成像方法的原理图；
图3为本发明实施例超快复合平面波成像装置的结构示意图；
图4为本发明实施例发射接收超声探头的结构示意图；
图5为本发明实施例单角度发射-三角度接收的示意图；
图6为本发明实施例响应频率确定装置的结构示意图；
图7为本发明实施例序列延时控制模块控制阵元发出一定发射角度的超声波的控制原理图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例基于宽频带声学超材料的超快复合平面波成像方法的流程图；图2为本发明实施例基于宽频带声学超材料的超快复合平面波成像方法的原理图。
参见图1-图2，实施例的基于宽频带声学超材料的超快复合平面波成像方法通过超快复合平面波成像装置实现。图3为本发明实施例超快复合平面波成像装置的结构示意图。参见图3，超快复合平面波成像装置包括发射接收超声探头1和声学超材料结构2；发射接收超声探头1发出的超声波信号通过声学超材料结构2后到达被测物体3，发射接收超声探头1还接收由被测物体3反射的回波信号。声学超材料结构2具有宽频带特性。
图4为本发明实施例发射接收超声探头的结构示意图。参见图4，发射接收超声探头1包括发射阵元11、接收阵元12和序列延时控制模块13，发射阵元11(发射序列)和接收阵元12(接收序列)分别与序列延时控制模块13连接，序列延时控制模块13用于控制发射阵元11的发射频率和发射角度，以及控制接收阵元12的接收频率和接收角度。
超快复合平面波成像方法，具体包括：
步骤S1：控制发射接收超声探头以预设发射频率，且以第一预设发射角度发出超声波信号，超声波信号通过声学超材料结构后到达被测物体；预设发射频率与声学超材料结构的响应频率相等。预设发射频率指发射序列的中心频率，声学超材料结构的响应频率指超材料发挥最优性能的频率。其实现原理图如图5的(a)部分所示。
声学超材料结构的响应频率是采用预实验方法，并通过响应频率确定装置确定的。图6为本发明实施例响应频率确定装置的结构示意图，参见图6，响应频率确定装置包括发射探头3、接收探头4、与发射探头3连接的信号发生器5以及与接收探头4连接的示波器6；发射探头3的输出端通过声学超材料结构2与接收探头4的输入端连接；信号发生器5与示波器6连接。发射探头3和接收探头4均为单阵元探头。
本实施例中，声学超材料结构的响应频率的确定方法为：
1)信号发生器在预设频率范围内以预设步长发出多个不同频率的电信号，每个频率的电信号均分为两路，第一路信号(通道1中的信号)直接在示波器上显示，第二路信号(通道2中的信号)依次经过发射探头、声学超材料结构和接收探头后，在示波器上显示。
2)比较每个频率下，在示波器上显示的两路电信号的波形，确定每个频率对应的第二路信号的波形的畸变程度、放大程度和延迟程度。
3)将畸变程度最小、放大程度最大，且延迟程度最小的第二路信号的波形对应的信号发生器的频率确定为声学超材料结构的响应频率。
作为一种可选的实施方式，所述预设频率范围为[0.5MHz,15MHz]；所述预设步长为1MHz。在该实施方式中，提供了一个采用预实验方法声学超材料结构的响应频率的具体过程。
将信号发射器的频率依次从0.5MHz调整至15MHz(步长为1MHz)，分别记录不同频率下通道1、通道2的信号，共15组信号。按照最小畸变、最优放大及最小延迟的原则，对比每组通道2接收信号与通道1发射信号，挑选通道2信号波形畸变最小、放大倍数最大且延迟最小的一组，此时，该组通道1信号的发射频率为与声学超材料相匹配的频率，即将此时信号发生器的频率确定为声学超材料结构的响应频率。
本实施例中，第一预设发射角度是通过序列延时控制模块13控制得到的，图7为本发明实施例序列延时控制模块控制阵元发出一定发射角度的超声波的控制原理图。参见图7，具体控制过程：序列延时控制模块13在控制发射阵元发出平面波时，当发射阵元的发射角度为零时，即当发射零角度的平面波时，序列延时控制模块13只需要控制多个发射阵元同时发射即可；当发射阵元的发射角度不为零时，即当发射带有偏角的平面波时，序列延时控制模块13控制给各个发射阵元施加延时，对于序号为x的发射阵元的延时计算公式为：
Figure BDA0002185749690000061
θ表示偏角，pitch表示阵元间隔，c表示声速，通过施加这样的延时，使不同阵元先后发出平面波，它们在波前形成一个平面(波前)，并且波前与阵元的夹角为θ。
步骤S2：控制发射接收超声探头以预设接收频率，且分别以第一预设接收角度、第二预设接收角度和第三预设接收角度接收被测物体反射的回波信号；预设接收频率为预设发射频率的n倍。即预设发射频率为BHz，预设接收频率为nBHz。
其中，n＞1；第一预设接收角度等于第一预设发射角度，第二预设接收角度小于第一预设发射角度，第三预设接收角度大于第一预设发射角度。本实施例中，所述第一预设接收角度与所述第一预设发射角度均为A°；所述第二预设接收角度为A-δ°，所述第三预设接收角度为A+δ°，其中，δ表示偏转角度的大小，δ＞0。其实现原理图如图5的(b)部分所示。
作为一种优选的实施方式，1＜n＜3。
作为一种优选的实施方式，5＜δ＜7。
步骤S3：采用回波信号重建被测物体的图像。该步骤，具体包括：
31：采用深度学习算法去除回波信号中的伪影信号，得到回波优化信号。
具体为：
获取训练数据集；训练数据集包括含伪影信号的超声训练信号和去除伪影信号的超声训练信号。
构建三层卷积神经网络。第一层特征提取，第二层非线性滤波，第三层特征组合。每层模型公式都为ReLU激活函数，即f(x)＝ReLU(W*x+b)，其中，x表示输入数据，w和b均表示模型的参数。
将含伪影信号的超声训练信号作为输入，去除伪影信号的超声训练信号作为输出，对三层卷积神经网络进行训练，得到训练好的三层卷积神经网络。
将回波信号输入至训练好的三层卷积神经网络中，得到回波优化信号。
32：采用回波优化信号重建被测物体的图像。具体为：根据回波优化信号，采用波束合成方法重建被测物体的图像。
本实施例的基于宽频带声学超材料的超快复合平面波成像方法，具有以下优点：
1、解决了平面波成像深度(穿透性差)和成像质量(分辨率及对比度差)的难题。
2、设置预设发射频率与声学超材料结构的响应频率相等，使得入射能量增强，提高了成像深度
3、目前，现有的常用的超快复合平面波成像，是对不同发射角度(…,-A°,0°,+A°,…)平面波的回波信号进行相干复合成像，本实施例步骤S1和步骤S2中，采用单角度发射-三角度接收的超快复合成像方法。当发射角度为+A°的单角度时，接收角度为A-δ°，A+0°，A+δ°的三个角度接收时，回波信号的数据量会相应增加，从而提高成像质量。
4、目前，现有的常规超声成像都是同一频率的发射/接收序列，本实施例采用单频发射-n倍频接收的超快复合平面波成像模式，且此发射/接收频率与声学超材料最优的双负特性相匹配。这样能够使接收信号的能量进一步放大，便于成像质量的提高。
5、基于接收到的回波信号，利用深度学习对回波信号进行深度学习处理，重建图像，去掉伪影，也可以进一步提高成像质量。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上，本说明书内容不应理解为对本发明的限制。

Claims (9)

1.基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，
所述超快复合平面波成像方法通过超快复合平面波成像装置实现；所述超快复合平面波成像装置包括发射接收超声探头和声学超材料结构；所述发射接收超声探头发出的超声波信号通过所述声学超材料结构后到达被测物体，所述发射接收超声探头还接收由所述被测物体反射的回波信号；
所述超快复合平面波成像方法，具体包括：
控制所述发射接收超声探头以预设发射频率，且以第一预设发射角度发出超声波信号，所述超声波信号通过所述声学超材料结构后到达被测物体；所述预设发射频率与所述声学超材料结构的响应频率相等；
控制所述发射接收超声探头以预设接收频率，且分别以第一预设接收角度、第二预设接收角度和第三预设接收角度接收所述被测物体反射的回波信号；所述预设接收频率为所述预设发射频率的n倍，n＞1；所述第一预设接收角度等于所述第一预设发射角度，所述第二预设接收角度小于所述第一预设发射角度，所述第三预设接收角度大于所述第一预设发射角度；
采用所述回波信号重建所述被测物体的图像。
2.根据权利要求1所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，所述采用所述回波信号重建所述被测物体的图像，具体包括：
采用深度学习算法去除所述回波信号中的伪影信号，得到回波优化信号；
采用所述回波优化信号重建所述被测物体的图像。
3.根据权利要求1所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，
所述声学超材料结构的响应频率通过响应频率确定装置确定；所述响应频率确定装置包括发射探头、接收探头、与所述发射探头连接的信号发生器以及与所述接收探头连接的示波器；所述发射探头的输出端通过所述声学超材料结构与所述接收探头的输入端连接；所述信号发生器与所述示波器连接；
所述声学超材料结构的响应频率的确定方法为：
所述信号发生器在预设频率范围内以预设步长发出多个不同频率的电信号，每个频率的电信号均分为两路，第一路信号直接在示波器上显示，第二路信号依次经过所述发射探头、所述声学超材料结构和所述接收探头后，在所述示波器上显示；
比较每个频率下，在所述示波器上显示的两路电信号的波形，确定每个频率对应的第二路信号的波形的畸变程度、放大程度和延迟程度；
将畸变程度最小、放大程度最大，且延迟程度最小的第二路信号的波形对应的信号发生器的频率确定为声学超材料结构的响应频率。
4.根据权利要求2所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，所述采用深度学习算法去除所述回波信号中的伪影信号，得到回波优化信号，具体包括：
获取训练数据集；所述训练数据集包括含伪影信号的超声训练信号和去除伪影信号的超声训练信号；
构建三层卷积神经网络；
将所述含伪影信号的超声训练信号作为输入，所述去除伪影信号的超声训练信号作为输出，对所述三层卷积神经网络进行训练，得到训练好的三层卷积神经网络；
将所述回波信号输入至所述训练好的三层卷积神经网络中，得到回波优化信号。
5.根据权利要求2所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，所述采用所述回波优化信号重建所述被测物体的图像，具体包括：
根据所述回波优化信号，采用波束合成方法重建所述被测物体的图像。
6.根据权利要求1所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，1＜n＜3。
7.根据权利要求1所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，所述第一预设接收角度与所述第一预设发射角度均为A°；所述第二预设接收角度为A-δ°，所述第三预设接收角度为A+δ°，其中δ＞0。
8.根据权利要求7所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，5＜δ＜7。
9.根据权利要求3所述的基于宽频带声学超材料的超快复合平面波成像方法，其特征在于，所述预设频率范围为[0.5MHz,15MHz]；所述预设步长为1MHz。
CN201910813822.4A 2019-08-30 2019-08-30 基于宽频带声学超材料的超快复合平面波成像方法 Active CN110477951B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910813822.4A CN110477951B (zh) 2019-08-30 2019-08-30 基于宽频带声学超材料的超快复合平面波成像方法
US16/864,584 US11829872B2 (en) 2019-08-30 2020-05-01 Method for ultrafast compound plane wave imaging based on broadband acoustic metamaterial

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910813822.4A CN110477951B (zh) 2019-08-30 2019-08-30 基于宽频带声学超材料的超快复合平面波成像方法

Publications (2)

Publication Number Publication Date
CN110477951A CN110477951A (zh) 2019-11-22
CN110477951B true CN110477951B (zh) 2020-08-25

Family

ID=68555414

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910813822.4A Active CN110477951B (zh) 2019-08-30 2019-08-30 基于宽频带声学超材料的超快复合平面波成像方法

Country Status (2)

Country Link
US (1) US11829872B2 (zh)
CN (1) CN110477951B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109510791B (zh) * 2017-09-15 2024-10-11 华为技术有限公司 传输方法和传输装置
US12092735B2 (en) * 2019-10-16 2024-09-17 Koninklijke Philips N.V. Method and apparatus for deep learning-based ultrasound beamforming
WO2021114071A1 (zh) * 2019-12-10 2021-06-17 深圳大学 一种经颅超声成像方法、装置及计算机可读存储介质
CN112244894B (zh) * 2020-10-19 2021-10-12 浙江大学 基于宽频带声学超材料的超声无创穿颅成像方法及系统
WO2022082627A1 (zh) * 2020-10-22 2022-04-28 深圳迈瑞生物医疗电子股份有限公司 超声成像方法和系统以及光声成像方法和系统
CN113138201B (zh) * 2021-03-24 2022-05-20 北京大学 用于无线被动式环境状态检测的超材料物联网系统及方法
CN113504307B (zh) * 2021-09-10 2021-12-21 西南石油大学 一种多频率岩心声速测量装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106481385A (zh) * 2016-11-03 2017-03-08 清华大学苏州汽车研究院（相城） 一种基于声学超材料的宽频小容积消声结构
CN109428174A (zh) * 2017-08-31 2019-03-05 深圳光启尖端技术有限责任公司 一种超材料

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10054707B2 (en) * 2016-04-15 2018-08-21 Baker Hughes, A Ge Company, Llc Bipolar acoustic hyperlens for dual-string thru-casing ultrasonic sensors
US9952343B2 (en) * 2016-07-20 2018-04-24 Baker Hughes, A Ge Company, Llc Rhodonea cell acoustic hyperlens for thru-casing ultrasonic sensors
US10783871B2 (en) * 2016-10-04 2020-09-22 Rutgers, The State University Of New Jersey Metal acoustic lens and method of manufacturing same
CN106580369A (zh) * 2016-11-30 2017-04-26 珠海威泓医疗科技有限公司 一种多角度平面波相干彩色多普勒成像方法
GB2557345B (en) * 2016-12-08 2021-10-13 Bae Systems Plc MIMO communication system and data link
DE102017209823A1 (de) * 2017-06-09 2018-12-13 Robert Bosch Gmbh Ultraschallsensor
CN107714091B (zh) * 2017-11-27 2019-12-20 西安交通大学 经颅低频超声线性调频脉冲逆转微泡成像方法
CN109975814B (zh) * 2017-12-28 2020-09-22 深圳先进技术研究院 超声成像方法、系统和设备
US11004439B2 (en) * 2018-02-26 2021-05-11 Toyota Motor Engineering & Manufacturing North America, Inc. Acoustic absorber
CN108836389B (zh) * 2018-06-01 2020-08-18 西安交通大学 平面波相关点相干自适应波束合成成像方法
CN109559359B (zh) * 2018-09-27 2023-06-02 东南大学 基于深度学习实现的稀疏角度数据重建图像的伪影去除方法
US20210393238A1 (en) * 2018-10-18 2021-12-23 The Regents Of The University Of California System and method for rendering objects transparent to ultrasound
US11067687B2 (en) * 2019-04-25 2021-07-20 Elwha, Llc Multipath acoustic holography and virtual haptics
CN110101411B (zh) * 2019-05-28 2020-11-10 飞依诺科技(苏州)有限公司 超声成像空间复合方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106481385A (zh) * 2016-11-03 2017-03-08 清华大学苏州汽车研究院（相城） 一种基于声学超材料的宽频小容积消声结构
CN109428174A (zh) * 2017-08-31 2019-03-05 深圳光启尖端技术有限责任公司 一种超材料

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Acoustic Filtering Properties of 3D Elastic Metamaterials Structured by Crack-Like Inclusions》;I. Ya. Zhbadynskyi，等;《2018 XXIIIrd International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory (DIPED)》;20180927;145-148 *
《Simulation-based conceptual design of an acoustic metamaterial with full band gap using an air-based 1-3 piezoelectric composite for ultrasonic noise control》;ShahrokhRezaei,等;《Comptes Rendus Mécanique》;20170228;第345卷(第2期);137-152 *

Also Published As

Publication number Publication date
US11829872B2 (en) 2023-11-28
CN110477951A (zh) 2019-11-22
US20210064993A1 (en) 2021-03-04

Similar Documents

Publication Publication Date Title
CN110477951B (zh) 基于宽频带声学超材料的超快复合平面波成像方法
KR101501479B1 (ko) 초음파 최적화 방법과 그를 위한 초음파 의료 장치
CN102293665B (zh) 超声波诊断装置
CN106901693B (zh) 基于环形超声换能器阵列的大景深光声-声速双模成像法
JP5851575B2 (ja) 超音波診断装置
CN107714091B (zh) 经颅低频超声线性调频脉冲逆转微泡成像方法
US20130172752A1 (en) Ultrasound transducer apparatus and ultrasound imaging system and imaging method
JPH0120899B2 (zh)
CN107280707A (zh) 用于声电成像的相控阵超声聚焦系统
CN103490754B (zh) 一种大时间带宽积的超声信号及其脉冲压缩方法及系统
CN1529567A (zh) 超声波诊断装置
JPH06209941A (ja) 超音波診断装置
CN101966088A (zh) 基于柔性相控超声阵列的口腔综合检测方法与装置
CN104622505B (zh) 一种超声颅内血流检测系统和方法
CN116077099B (zh) 基于环形阵列多子阵快速图像重建的超声ct反射成像方法
Vallet et al. Enhancement of photoacoustic imaging quality by using CMUT technology: Experimental study
CN110441398A (zh) 一种基于多层介质声速模型的合成孔径超声成像方法
CN107204021B (zh) 基于高斯函数探头响应模型和压缩感知的超声成像方法
Tarbox et al. Low-power ultrasound imaging systems using time delay spectrometry
Irie et al. Transmission of 100-MHz-range ultrasound through a fused quartz fiber
WO2006057092A1 (ja) 超音波撮像装置
Bera et al. Three-dimensional beamforming combining micro-beamformed RF datasets
CN114847983A (zh) 一种基于预扫描的乳腺超声ct反射成像参数自适应补偿方法
CN104905815B (zh) 一种线阵超声探头
KR101772098B1 (ko) 다양한 포커싱을 통해 다중 선택 가능한 초음파 프로브 및 이를 구비한 초음파 영상 진단장치

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant