CN115500865A - 用于被动声成像的快速延迟乘和自适应波束形成方法 - Google Patents
用于被动声成像的快速延迟乘和自适应波束形成方法 Download PDFInfo
- Publication number
- CN115500865A CN115500865A CN202211220827.4A CN202211220827A CN115500865A CN 115500865 A CN115500865 A CN 115500865A CN 202211220827 A CN202211220827 A CN 202211220827A CN 115500865 A CN115500865 A CN 115500865A
- Authority
- CN
- China
- Prior art keywords
- imaging
- time
- cavitation
- signal
- passive
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 8
- 239000000523 sample Substances 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 230000010354 integration Effects 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 230000006835 compression Effects 0.000 claims abstract description 5
- 238000007906 compression Methods 0.000 claims abstract description 5
- 238000012634 optical imaging Methods 0.000 claims abstract 3
- 230000001133 acceleration Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000002604 ultrasonography Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 5
- 230000003321 amplification Effects 0.000 claims description 4
- 230000003139 buffering effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000002087 whitening effect Effects 0.000 claims description 2
- 238000002601 radiography Methods 0.000 claims 2
- 230000017105 transposition Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 5
- 230000001629 suppression Effects 0.000 abstract description 2
- XREKLQOUFWBSFH-UHFFFAOYSA-N dimethyl 2-acetylbutanedioate Chemical compound COC(=O)CC(C(C)=O)C(=O)OC XREKLQOUFWBSFH-UHFFFAOYSA-N 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
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/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N7/02—Localised ultrasound hyperthermia
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一种用于被动声成像的快速延迟乘和自适应波束形成方法,首先设置位于待测物体一侧的功率探头的超声输出功率和时序,在待测物体的另一侧采集B模式信号或各阵元被动接收的空化信号后,采用改进的最小方差‑相乘延时叠加算法处理,得到聚焦信号,经能量时间积分以实现长脉冲光成像后,通过对数压缩处理和扫描转换,得到最终的显示图像。本发明将自适应波束形成器和DMAS相结合,在远场具有更好的旁瓣抑制,在被动空化成像时具有最佳聚焦和最小伪影的效果的同时,将计算复杂度从o(n2)降低至o(n)。
Description
技术领域
本发明涉及的是一种超声成像领域的技术,具体是一种用于被动声成像的快速延迟乘和(DMAS)自适应波束形成方法。
背景技术
高强度聚焦超声(HIFU)是一种通过聚焦超声脉冲来机械地摧毁组织的技术,该聚焦超声脉冲生成空化。治疗过程中需要精确定位聚焦位置,也需要检测空化的强弱和空间分布。现有的反向散射超声回波的频率的谐波的分析能够被动检测空化,但其空间分辨率较低,而传统的B型超声成像技术虽然能够实时地观察组织的改变,但是灵敏度较低,必须要聚焦超声足够强到引起组织变性才能观测到。
发明内容
本发明针对现有技术在远场会产生尾部伪影不足、计算复杂度过高的缺陷,提出一种用于被动声成像的快速延迟乘和自适应波束形成方法,将自适应波束形成器和DMAS相结合,在远场具有更好的旁瓣抑制,在被动空化成像时具有最佳聚焦和最小伪影的效果的同时,将计算复杂度从o(n2)降低至o(n)。
本发明是通过以下技术方案实现的:
本发明涉及一种用于被动声成像的快速延迟乘和自适应波束形成方法,包括:
步骤一:设置位于待测物体一侧的功率探头的超声输出功率和时序,在待测物体的另一侧采集B模式信号或各阵元被动接收的空化信号;
所述的时序是指:控制空化成像脉冲与正常B模式和/或造影模式脉冲的发射时序,使两者在时序上的完全错开从而将干扰显著降低,实现分时全双工,具体为:以满足空化脉冲发射后得到的回波强度衰减至小于等于常规成像超声脉冲的强度所需要的安全时间,即:其中:ft为空化脉冲频率,fB为成像脉冲频率,Vt为空化脉冲电压,VB为成像脉冲电压,d为当前成像深度,α为超声在人体内的衰减系数,一般为0.5dB MHz/cm。
步骤二:对步骤一所得的空化信号采用改进的最小方差-相乘延时叠加(MV-DMAS)算法处理,具体为:在xz位置进行聚焦延时,t时刻所得的每条扫描线数据si(x,z,t)=pi(t+τ(x,z)),其中:pi(t)为换能器第i个阵元收到的空化信号,z垂直于探头表面,x平行于探头表面;该空化信号对应的延时的长度对扫描线数据进行加权得到聚焦信号其中改进的权值wi,通过以下方式得到:其中:a为单位矢量;H代表转置。
步骤三:将步骤二得到的聚焦信号进行能量时间积分以实现长脉冲光成像,具体为:通过化简得到波束成型信号因此可得到xy位置的被动声场成像的声强其中:n为阵元数;为w加权过后的si;为保留符号开平方;ρ0为人体的平均密度;c0为人体的平均声速;T为积分时长;
步骤四:对步骤三得到的长脉冲光成像进行对数压缩处理和扫描转换,得到最终的显示图像。
技术效果
本发明通过在同一台设备硬件上时分复用,实现主动和被动的协同工作,通过引入可调节的高通滤波器后,通过调节不同的截止频率实现被动空化成像和被动声场映射的同时显示,即明确定位空化位置并显示声场强弱的空间分布。
附图说明
图1为本发明流程图;
图2为实施例系统示意图;
图3为实施例时序控制示意图;
图3和图4为实施例效果示意图;
图5为实施例效果示意图
图中:(a)为DAS成像算法;(b)为MV成像算法;(c)为DMAS成像算法;(d)为本方法。
具体实施方式
如图1所示,为本实施例的一种快速延迟乘和自适应波束形成方法,包括:
步骤1、调节高通滤波器参数,以同时完成被动空化成像和被动声场映射,具体为:通过设置高通滤波器的截止频率,使其仅滤去直流,从而得到全部交流的声场信号用于被动声场映射的成像;通过设置高通滤波器的截止频率为发射信号的频率时,能够在被测区域发生空化时采集到白化频谱,用于被动空化成像;当未发生空化时,则仅收到发射频率和少量的偶次倍频,故只需滤去f0和2f0后,即得到空化信号或无其他接收信号。例如发射频率1Mhz,设置高通滤波截至频率为5Mhz并对接收信号进行被动成像,即反映空化强弱的空间分辨率的被动空化成像。
步骤2、缓冲hifu单次发射期间每通道的接收数据,具体为:对于每个接收通道准备一个高速缓冲,该缓冲的零时刻是Tdelay=T0+Thifu+Ts,其中:T0为hifu发射的时间,Thifu为hifu声波到达被测物体的时间,Ts为被测物体收到hifu信号后发出的反射信号到达探头的时间。
所述的高速缓冲的位宽优选设置为与模数转换处理的位宽相同,通常为12bit或者14bit。
步骤4、根据通道延时对该通道的接收信号进行时移处理,具体为:si(x,z,t)=pi(t+τ(x,z)),其中:pi(t)为换能器第i个阵元收到的空化信号。
步骤7、化简聚焦信号q得到q*,通过时间积分得到被动声场成像的声强I,具体为:其中:化简后的聚焦信号为n为阵元数;为w加权过后的si;为保留符号开平方;ρ0为人体的平均密度;c0为人体的平均声速;T为积分时长;
步骤8、将每次发射计算得到的声强I进行时间平均作为各像素值。
步骤9、对数压缩。
步骤10、增益调节。
步骤11、扫描变换。
如图2和图3所示,为本实施例涉及的一种实现上述方法的成像系统,包括:与探头相连的波束成形单元以及图像加速单元,其中:波束成形单元包括128个接收通道、触发及时序控制单元以及发射波束赋形单元,触发及时序控制单元接收来自图像加速单元的指令并分别向HIFU单元、发射波束赋形单元和接收通道发出如图4所示的时序控制指令;接收通道缓冲HIFU单次发射期间每通道的接收数据,依次经低噪放大处理、压控放大器处理、放大处理和抗混叠滤波后经模数转换得到数字信号,该数字信号经高通滤波后缓存并通过PCIe总线批量发送至图像加速单元,图像加速单元通过计算对应的每个通道延时并计算权重后,基于改进的MV-DMAS算法得到被动声场成像的声强I,并通过时间平均得到对应被动声场成像;图像加速单元同时将收到的缓存后的数字信号进行B模式成像,最终与被动声场成像融合后输出。
经过具体实际实验,本实施例使用128单元线性探头,单元间距=0.3mm。并在40mm深度处放置空化泡。对从仿真模型生成的pi(t)进行去直流高通滤波,得到以下结果。
如图5所示,为本实施例完成的B模式和被动成像模式双工模式融合图像,图中a为DAS结果,b为最小方差(Minimum Variance)结果,c为延迟乘和(DMAS)结果,d为本发明改进的MV-DMAS结果。结果显示,本方法分辨率优于现有技术且通过GPU加速,实现了实时声场被动成像。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (7)
1.一种用于被动声成像的快速延迟乘和自适应波束形成方法,其特征在于,首先设置位于待测物体一侧的功率探头的超声输出功率和时序,在待测物体的另一侧采集B模式信号或各阵元被动接收的空化信号后,采用改进的最小方差-相乘延时叠加算法处理,得到聚焦信号,经能量时间积分以实现长脉冲光成像后,通过对数压缩处理和扫描转换,得到最终的显示图像。
2.根据权利要求1所述的用于被动声成像的快速延迟乘和自适应波束形成方法,其特征是,所述的时序是指:控制空化成像脉冲与正常B模式和/或造影模式脉冲的发射时序,使两者在时序上的完全错开从而将干扰显著降低,实现分时全双工。
6.根据权利要求1-5中任一所述的用于被动声成像的快速延迟乘和自适应波束形成方法,其特征是,具体包括:
步骤1、调节高通滤波器参数,以同时完成被动空化成像和被动声场映射,具体为:通过设置高通滤波器的截止频率,使其仅滤去直流,从而得到全部交流的声场信号用于被动声场映射的成像;通过设置高通滤波器的截止频率为发射信号的频率时,能够在被测区域发生空化时采集到白化频谱,用于被动空化成像;当未发生空化时,则仅收到发射频率和少量的偶次倍频,故只需滤去f0和2f0后,即得到空化信号或无其他接收信号;
步骤2、缓冲hifu单次发射期间每通道的接收数据,具体为:对于每个接收通道准备一个高速缓冲,该缓冲的零时刻是Tdelay=T0+Thifu+Ts,其中:T0为hifu发射的时间,Thifu为hifu声波到达被测物体的时间,Ts为被测物体收到hifu信号后发出的反射信号到达探头的时间;
步骤4、根据通道延时对该通道的接收信号进行时移处理,具体为:si(x,z,t)=pi(t+τ(x,z)),其中:pi(t)为换能器第i个阵元收到的空化信号;
步骤7、化简聚焦信号q得到q*,通过时间积分得到被动声场成像的声强I,具体为:其中:化简后的聚焦信号为n为阵元数;为w加权过后的si;为保留符号开平方;ρ0为人体的平均密度;c0为人体的平均声速;T为积分时长;
步骤8、将每次发射计算得到的声强I进行时间平均作为各像素值;
步骤9、对数压缩;
步骤10、增益调节;
步骤11、扫描变换。
7.一种实现权利要求1-6中任一所述用于被动声成像的快速延迟乘和自适应波束形成方法的成像系统,其特征在于,包括:与探头相连的波束成形单元以及图像加速单元,其中:波束成形单元包括128个接收通道、触发及时序控制单元以及发射波束赋形单元,触发及时序控制单元接收来自图像加速单元的指令并分别向HIFU单元、发射波束赋形单元和接收通道发出时序控制指令;接收通道缓冲HIFU单次发射期间每通道的接收数据,依次经低噪放大处理、压控放大器处理、放大处理和抗混叠滤波后经模数转换得到数字信号,该数字信号经高通滤波后缓存并通过PCIe总线批量发送至图像加速单元,图像加速单元通过计算对应的每个通道延时并计算权重后,基于改进的MV-DMAS算法得到被动声场成像的声强I,并通过时间平均得到对应被动声场成像;图像加速单元同时将收到的缓存后的数字信号进行B模式成像,最终与被动声场成像融合后输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211220827.4A CN115500865B (zh) | 2022-10-08 | 2022-10-08 | 用于被动声成像的快速延迟乘和自适应波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211220827.4A CN115500865B (zh) | 2022-10-08 | 2022-10-08 | 用于被动声成像的快速延迟乘和自适应波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115500865A true CN115500865A (zh) | 2022-12-23 |
CN115500865B CN115500865B (zh) | 2024-06-18 |
Family
ID=84507343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211220827.4A Active CN115500865B (zh) | 2022-10-08 | 2022-10-08 | 用于被动声成像的快速延迟乘和自适应波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115500865B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040059221A1 (en) * | 2002-09-25 | 2004-03-25 | Takashi Azuma | Ultrasonic imaging system and method |
EP2480136A1 (en) * | 2009-09-22 | 2012-08-01 | Isis Innovation Ltd | Ultrasound systems |
CN104777484A (zh) * | 2015-02-13 | 2015-07-15 | 西安交通大学 | 压缩自适应波束合成的平面波超声成像和微泡成像的方法与系统 |
CN108354627A (zh) * | 2018-04-04 | 2018-08-03 | 东北大学 | 一种提高帧频的超声波束形成方法 |
US20210212709A1 (en) * | 2018-05-25 | 2021-07-15 | Cardiawave Sa | Ultrasonic processing apparatus comprising means for imaging cavitation bubbles |
-
2022
- 2022-10-08 CN CN202211220827.4A patent/CN115500865B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040059221A1 (en) * | 2002-09-25 | 2004-03-25 | Takashi Azuma | Ultrasonic imaging system and method |
EP2480136A1 (en) * | 2009-09-22 | 2012-08-01 | Isis Innovation Ltd | Ultrasound systems |
CN104777484A (zh) * | 2015-02-13 | 2015-07-15 | 西安交通大学 | 压缩自适应波束合成的平面波超声成像和微泡成像的方法与系统 |
CN108354627A (zh) * | 2018-04-04 | 2018-08-03 | 东北大学 | 一种提高帧频的超声波束形成方法 |
US20210212709A1 (en) * | 2018-05-25 | 2021-07-15 | Cardiawave Sa | Ultrasonic processing apparatus comprising means for imaging cavitation bubbles |
Non-Patent Citations (2)
Title |
---|
ERASMIA LYKA等: "Passive Acoustic Mapping Using Data-Adaptive Beamforming Based on Higher Order Statistics", IEEE TRANSACTIONS ON MEDICAL IMAGING, vol. 37, no. 12, 31 December 2018 (2018-12-31), pages 2592 - 2592, XP011696698, DOI: 10.1109/TMI.2018.2843291 * |
周芳等: "Passive Acoustic Mapping Using Data-Adaptive Beamforming Based on Higher Order Statistics", 中国医学物理学杂志, vol. 40, no. 5, 31 May 2023 (2023-05-31), pages 568 - 575 * |
Also Published As
Publication number | Publication date |
---|---|
CN115500865B (zh) | 2024-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11717270B2 (en) | Method, apparatus and system for imaging in ultrasonic scanning | |
US5891038A (en) | Method, apparatus and applications for combining transmit wave functions to obtain synthetic waveform in ultrasonic imaging system | |
JP4717995B2 (ja) | 超音波ビーム経路の数値的最適化方式 | |
EP1832894B1 (en) | System and method for automatic gain compensation based on image processing | |
US6108572A (en) | Method and apparatus for harmonic imaging using multiple focal zones | |
JP4620261B2 (ja) | 高次の非線形成分による超音波撮像 | |
JP4424707B2 (ja) | Bモード超音波イメージングにおける自動送信波形最適化のための方法及び装置 | |
US20060094962A1 (en) | Aperture shading estimation techniques for reducing ultrasound multi-line image distortion | |
CN104614728B (zh) | 一种超声成像装置及成像方法 | |
CN107714091B (zh) | 经颅低频超声线性调频脉冲逆转微泡成像方法 | |
JPH08500023A (ja) | 自動膀胱走査装置 | |
US20050124883A1 (en) | Adaptive parallel artifact mitigation | |
EP1695110A2 (en) | Ultrasound imaging system and method having adaptive selection of image frame rate and/or number of echo samples averaged | |
EP1815795B1 (en) | Ultrasonographic device | |
US20110245676A1 (en) | Method and apparatus for ultrasound signal acquisition and processing | |
US20080030581A1 (en) | Multizone Color Doppler Beam Transmission Method | |
CN115500865B (zh) | 用于被动声成像的快速延迟乘和自适应波束形成方法 | |
Gammelmark et al. | Multi-element synthetic transmit aperture imaging using temporal encoding | |
US6231512B1 (en) | Method and apparatus for parametric harmonic imaging | |
JP5390207B2 (ja) | 超音波診断装置および超音波診断装置の制御プログラム | |
JP4580490B2 (ja) | 超音波診断装置 | |
KR101120691B1 (ko) | 공간합성을 이용한 초음파영상 합성방법 | |
JP4599408B2 (ja) | 超音波診断装置 | |
JPH03261466A (ja) | 超音波診断装置 | |
JP2015186494A (ja) | 超音波診断装置 |
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 |