CN110575201B - 基于反相Golay码的超声微泡空化成像方法及装置 - Google Patents
基于反相Golay码的超声微泡空化成像方法及装置 Download PDFInfo
- Publication number
- CN110575201B CN110575201B CN201910952607.2A CN201910952607A CN110575201B CN 110575201 B CN110575201 B CN 110575201B CN 201910952607 A CN201910952607 A CN 201910952607A CN 110575201 B CN110575201 B CN 110575201B
- Authority
- CN
- China
- Prior art keywords
- cavitation
- imaging
- pulse
- pulses
- sequence
- 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
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
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- 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/488—Diagnostic techniques involving Doppler signals
-
- 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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- 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
- A61N2007/0039—Ultrasound therapy using microbubbles
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N2007/0052—Ultrasound therapy using the same transducer for therapy and imaging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
一种基于反相Golay码的超声微泡空化成像方法及装置,通过将待发射的若干空化成像脉冲采用Golay码编码后,得到一对正相及其180°反相空化脉冲,经与常规成像超声脉冲协同控制后发射;再通过对回波进行匹配滤波、分通道时延和叠加后解调,与常规成像超声脉冲得到的图像叠加后实现。本发明能够精确采集空化发生时的所有非线性的能量,检测到的空化强度更为准确。同时由于使用了脉冲压缩技术。空化成像的精度提高。
Description
技术领域
本发明涉及的是一种超声治疗领域的技术,具体是一种基于反相Golay码的超声微泡空化成像方法及装置。
背景技术
空化成像是监测微泡空化强度的一个必要手段,合适的空化强度对于医生和实验者都是非常必要参数。通常空化激励需要连续的强超声脉冲,通过回波的非线性强度来判断空化的强度,回波的非线性成分越多则表明空化强度越高;而空化成像脉冲数越多反而分辨率越低。因此空化成像的分辨率和空化强度产了一对矛盾。
发明内容
本发明针对现有技术存在的上述不足,提出一种基于反相Golay码的超声微泡空化成像方法及装置,通过编码发射和脉冲压缩技术,将空化激励超声和成像超声融合并同时发送。
本发明是通过以下技术方案实现的:
本发明涉及一种基于反相Golay码的超声微泡空化成像方法,通过将待发射的若干空化成像脉冲采用Golay码编码后,得到一对正相及其180°反相空化脉冲,经与常规成像超声脉冲协同控制后发射;再通过对回波进行匹配滤波、分通道时延和叠加后解调,与常规成像超声脉冲得到的图像叠加后实现。
所述的Golay码,具体为正交互补Golay(A,B)码,其长度与对应的空化成像脉冲相匹配。
所述的一对正相及其反相空化脉冲,包括经Golay码的A序列、B序列编码后的信号以及其反相后的信号。
所述的叠加,通过将匹配滤波后得到的一对正相及其反相空化回波信号经叠加处理后得到空化脉冲回波的基波分量和非线性分量。
所述的常规成像超声脉冲,包括正常B模式、多普勒模式或造影模式脉冲或其组合。
所述的协同控制,优选控制空化成像脉冲与正常B模式和/或造影模式脉冲的发射时序。
所述的协同控制是指:控制空化成像脉冲和常规成像超声脉冲的发射时序,使两者在时序上的完全错开从而将干扰显著降低,实现分时全双工,具体为:以满足空化脉冲发射后得到的回波强度衰减至小于等于常规成像超声脉冲的强度所需要的安全时间,即: 其中:ft为空化脉冲频率,fB为成像脉冲频率,Vt为空化脉冲电压,VB为成像脉冲电压,d为当前成像深度,α为超声在人体内的衰减系数,一般为0.5dBMHz/cm。
技术效果
相比现有技术,本发明创新地先使用匹配滤波器天然滤除二次谐波,再通过反相技术,将基波和非线性成分分离,从而精确采集空化发生时的所有非线性的能量,检测到的空化强度更为准确。同时由于使用了脉冲压缩技术。空化成像的精度提高。
附图说明
图1为本发明发射装置示意图;
图2为本发明发射流程示意图;
图3为实施例Golay码发射波形图;
图中:a、b分别为互补序列;
图4为实施例接收端匹配滤波器示意图;
图中:a、b分别用于接收AB互补序列;
图5为AB分别匹配滤波过后相加得到的脉冲压缩过后的波形图;
图6为实施例协同控制示意图。
具体实施方式
如图1所示,为本实施例涉及一种基于Golay码的超声微泡空化成像装置,包括:用于生成单个空化脉冲和成像脉冲的发射端、用于接收回波并解析出空化图像的接收端和分别与之相连发送接收开关,其中:发送接收开关与换能器相连并正对待测位置。
所述的发射端包括:双极脉冲发生器和发送序列存储器,其中:双极脉冲发生器的控制端接收发送聚焦延时以发射用于常规成像的正常B模式、多普勒模式或造影模式脉冲,发送序列存储器用于提供Golay码以编码空化成像脉冲。
所述的Golay码,具体为正交互补Golay(A,B)码,即A(a(n),n=1,...,N),B(b(n),n=1,...,N),满足a(n)*a(-n)+b(n)*b(-n)=2Nδ(n),其中:*为卷积,N为自然常数,δ(n)为脉冲序列。
所述的Golay码,其长度与空化成像脉冲相同:当空化成像脉冲为10个周期,则采用10个码字的Golay(A,B)码,例如A=1,-1,-1,1,-1,1,-1,-1,-1,1;B=1,-1,-1,-1,1,-1,1,-1,-1,1。
如图2所示,所述的发射端采用四次发射,通过分别发射A序列、B序列编码后的脉冲以及A序列、B序列编码后的反向脉冲,实现对回波中的基波成分和强非线性成分的分离。
所述的接收端包括:匹配滤波器、带有延时的接收通道、接收波束合成器、累加器以及解调器,其中:匹配滤波器与预存有滤波器系数的存储器相连并将接收到的回波信号与对应发射信号的接收码相卷积后输出至接收通道以进行对应的的时延、放大以及模数转换处理,累加器将合成后的同一发射焦点位置的匹配滤波信号叠加以实现谐波信号的脉冲压缩和基波信号的抑制,解调器经对累加器输出的信号解调并输出至后续信号处理器进行相应的包络检测、边缘增强和/或对数压缩处理。
如图4a和图4b所示,所述的匹配滤波器的导通截止设置与发送端的A序列、B序列相同。本实施例中通过发送序列存储器和滤波系数存储器实现,图中纵坐标±0.05分别对应A、B序列中的1和-1,横坐标为时间(图中20为示意性)。
所述的解调器分别对基波成分和非线性成分进行希尔伯特变换,相应得到基波成分能量EX,和非线性成分能量EY并输出至信号处理器。
优选地,所述的发射端进一步设有协同控制模块,用于控制单个空化脉冲和成像超声脉冲的发射时序。
如图6a所示,为发射间隙满足条件时的情况:理想状态下在一个周期内发射若干次空化脉冲,相邻空化脉冲的时间间隔为1/prf。在最后次空化脉冲发射后间隔安全时隙后发射若干次的常规成像脉冲Image,相邻的常规成像脉冲的时间间隔根据影像系统本身的间隔决定,通常为:其中:c为声波在人体内传播速度,d为成像深度(可通过现有技术采集获取),extra为额外等待时间;通常d与extra成反比。
如图6b所示,为发射间隙不满足条件时的情况:即根据当前成像深度信息得到的安全时隙以及上述发射端采用四次发射中的任意一次或任意四次所需时间超过了发射间隙时,协同控制模块将临时暂停双极脉冲发生器,通过强制増加空化脉冲之间的间隙(虚线框部分治疗脉冲关闭)以保证常规成像脉冲的不受干扰。
所述的强制增加,具体是指:协同控制模块以1为步长依次递增的方式缩短一个周期内空化脉冲的个数并重新尝试计算安全时隙,即扩大双极脉冲发射器的发射间隙,直至发射间隙大于安全时隙和成像脉冲所需时间之和时所对应的空化脉冲缩短的个数。
本实施例中增加的间隙由发射包数目(packetsize)决定,该参数控制治疗超声每发射n次暂停1次的比例。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (8)
1.一种基于反相Golay码的超声微泡空化成像方法,其特征在于,通过将待发射的若干空化成像脉冲采用Golay码编码后,得到一对正相及其180°反相空化脉冲,经与常规成像超声脉冲协同控制后发射;再通过对回波进行匹配滤波、分通道时延和叠加后解调,与常规成像超声脉冲得到的图像叠加后实现;
所述的Golay码,具体为正交互补Golay(A,B)码,其长度与对应的空化成像脉冲相匹配;
所述的匹配滤波,其导通截止设置与发射的设置相同;
3.根据权利要求1所述的方法,其特征是,所述的叠加,通过将匹配滤波后得到的一对正相及其反相空化回波信号经叠加处理后得到空化脉冲回波的基波分量和非线性分量。
5.一种基于Golay码的超声微泡空化成像装置,其特征在于,包括:用于生成单个空化脉冲和成像脉冲的发射端、用于接收回波并解析出空化图像的接收端和分别与之相连发送接收开关,其中:发送接收开关与换能器相连并正对待测位置;
所述的接收端包括:匹配滤波器、带有延时的接收通道、接收波束合成器、累加器以及解调器,其中:匹配滤波器与预存有滤波器系数的存储器相连并将接收到的回波信号与对应发射信号的接收码相卷积后输出至接收通道以进行对应的时延、放大以及模数转换处理,累加器将合成后的同一发射焦点位置的匹配滤波信号叠加以实现谐波信号的脉冲压缩和基波信号的抑制,解调器经对累加器输出的信号解调并输出至后续信号处理器进行相应的包络检测、边缘增强和/或对数压缩处理;
所述的匹配滤波器的导通截止设置与发送端的A序列、B序列相同;
6.根据权利要求5所述的成像装置,其特征是,所述的发射端包括:双极脉冲发生器和发送序列存储器,其中:双极脉冲发生器的控制端接收发送聚焦延时以发射用于常规成像的正常B模式、多普勒模式或造影模式脉冲,发送序列存储器用于提供Golay码以编码空化成像脉冲。
7.根据权利要求5或6所述的成像装置,其特征是,所述的发射端采用四次发射,通过分别发射A序列、B序列编码后的脉冲以及A序列、B序列编码后的反向脉冲,实现对回波中的基波成分和强非线性成分的分离。
8.根据权利要求7所述的成像装置,其特征是,所述的发射端进一步设有协同控制模块,用于控制单个空化脉冲和成像超声脉冲的发射时序;
当发射间隙满足条件时的情况,即理想状态下在一个周期内发射若干次空化脉冲,相邻空化脉冲的时间间隔为1/prf;在最后次空化脉冲发射后间隔安全时隙后发射若干次的常规成像脉冲Image,相邻的常规成像脉冲的时间间隔为其中:c为声波在人体内传播速度,d为成像深度,extra为额外等待时间;通常d与extra成反比;
当发射间隙不满足条件时的情况,即根据当前成像深度信息得到的安全时隙以及上述发射端采用四次发射中的任意一次或任意四次所需时间超过了发射间隙时,协同控制模块将临时暂停双极脉冲发生器,通过强制増加空化脉冲之间的间隙以保证常规成像脉冲的不受干扰。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910952607.2A CN110575201B (zh) | 2019-10-09 | 2019-10-09 | 基于反相Golay码的超声微泡空化成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910952607.2A CN110575201B (zh) | 2019-10-09 | 2019-10-09 | 基于反相Golay码的超声微泡空化成像方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110575201A CN110575201A (zh) | 2019-12-17 |
CN110575201B true CN110575201B (zh) | 2021-12-24 |
Family
ID=68814454
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910952607.2A Active CN110575201B (zh) | 2019-10-09 | 2019-10-09 | 基于反相Golay码的超声微泡空化成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110575201B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111134719B (zh) * | 2019-12-19 | 2021-01-19 | 西安交通大学 | 一种聚焦超声辐照相变纳米液滴的主被动超声复合成像方法及系统 |
CN113117257B (zh) * | 2019-12-30 | 2023-07-11 | 重庆融海超声医学工程研究中心有限公司 | 用于检测空化效应的方法及装置、超声治疗设备 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6042556A (en) * | 1998-09-04 | 2000-03-28 | University Of Washington | Method for determining phase advancement of transducer elements in high intensity focused ultrasound |
US6155980A (en) * | 1999-03-16 | 2000-12-05 | General Electric Company | Ultrasonic imaging system with beamforming using unipolar or bipolar coded excitation |
US6375618B1 (en) * | 2000-01-31 | 2002-04-23 | General Electric Company | Enhanced tissue-generated harmonic imaging using coded excitation |
US6491631B2 (en) * | 2001-01-11 | 2002-12-10 | General Electric Company | Harmonic golay-coded excitation with differential pulsing for diagnostic ultrasound imaging |
US20070083120A1 (en) * | 2005-09-22 | 2007-04-12 | Cain Charles A | Pulsed cavitational ultrasound therapy |
FR2912817B1 (fr) * | 2007-02-21 | 2009-05-22 | Super Sonic Imagine Sa | Procede d'optimisation de la focalisation d'ondes au travers d'un element introducteur d'aberations. |
CN101336830B (zh) * | 2007-07-03 | 2012-07-04 | 深圳迈瑞生物医疗电子股份有限公司 | 用于超声诊断成像的正交多普勒信号间隙填充方法与装置 |
CN101149433A (zh) * | 2007-11-01 | 2008-03-26 | 绵阳索尼克电子有限责任公司 | 一种使用反向Golay码译码的超声运动检测方法 |
CN102322805B (zh) * | 2011-08-26 | 2013-03-13 | 南京理工大学 | 空化泡最大泡半径的探测方法 |
CN102670258B (zh) * | 2012-04-13 | 2013-08-07 | 中国医学科学院生物医学工程研究所 | 一种眼科超声生物测量方法 |
TWI489106B (zh) * | 2013-10-28 | 2015-06-21 | Univ Nat Taiwan | 應用於超音波非線性成像之格雷編碼成像方法 |
CN104777484B (zh) * | 2015-02-13 | 2016-04-27 | 西安交通大学 | 压缩自适应波束合成的平面波超声成像和微泡成像的方法与系统 |
CN104757999B (zh) * | 2015-04-07 | 2018-03-09 | 声泰特(成都)科技有限公司 | 一种基于超声基波和谐波的非线性成像方法和系统 |
CN105548006B (zh) * | 2015-11-27 | 2019-07-16 | 西安科技大学 | 一种检测超声空化强度的系统 |
CN105548373B (zh) * | 2015-12-15 | 2018-03-27 | 华南理工大学 | 一种准单次正交互补Golay(A,B)码超声相控阵编码激励方法 |
CN106971055B (zh) * | 2016-01-12 | 2019-10-25 | 深圳迈瑞生物医疗电子股份有限公司 | 超声造影成像方法和系统 |
US10571554B2 (en) * | 2016-11-29 | 2020-02-25 | Siemens Medical Solutions Usa, Inc. | Adaptive post beamformation synthetic aperture for ultrasound imaging |
CN108882914B (zh) * | 2017-11-20 | 2021-04-30 | 深圳迈瑞生物医疗电子股份有限公司 | 超声造影成像方法及超声成像系统 |
CN107714091B (zh) * | 2017-11-27 | 2019-12-20 | 西安交通大学 | 经颅低频超声线性调频脉冲逆转微泡成像方法 |
CN109513123B (zh) * | 2018-12-28 | 2020-05-19 | 西安交通大学 | 一种基于半球阵的高分辨三维被动空化成像方法 |
CN109785271A (zh) * | 2019-02-25 | 2019-05-21 | 天津大学 | 一种基于编码激励和相干系数的内镜超声成像算法 |
CN110075430B (zh) * | 2019-04-28 | 2020-10-02 | 南京大学 | 一种基于信息熵的超声空化实时监测方法及系统 |
-
2019
- 2019-10-09 CN CN201910952607.2A patent/CN110575201B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110575201A (zh) | 2019-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100393370B1 (ko) | 직교 골레이 코드를 이용하는 초음파 영상 형성 방법 및장치 | |
CN103431873B (zh) | 超声融合谐波影像系统和方法 | |
US6241674B1 (en) | Medical ultrasound diagnostic imaging method and system with nonlinear phase modulation pulse compression | |
CN110575201B (zh) | 基于反相Golay码的超声微泡空化成像方法及装置 | |
CN102327129B (zh) | 超声波图像装置 | |
CN105120759B (zh) | 用于从旋转血管内超声设备产生图像的系统和方法 | |
US6010456A (en) | Method and apparatus for acoustic subtraction imaging using linear and nonlinear ultrasonic images | |
US20040184351A1 (en) | Ultrasonic transmitter, utrasonic transceiver and sonar apparatus | |
CN107714091B (zh) | 经颅低频超声线性调频脉冲逆转微泡成像方法 | |
WO2000057769A2 (en) | Ultrasonic imaging system using coded transmit pulses | |
ITRM20010028A1 (it) | Metodo di formazione di immagine a ultrasuoni e apparecchio basato suuna tecnica di compressione di impulsi utilizzante codici modificati d | |
CN1210003C (zh) | 基于脉冲压缩的超声波收发装置 | |
JP2003190157A (ja) | 超音波診断装置 | |
CN109254284A (zh) | 用于编码超声波回波检测的方法 | |
CN113768542B (zh) | 超声血流成像装置及超声设备 | |
Gran et al. | Coded ultrasound for blood flow estimation using subband processing | |
CN102247167A (zh) | 一种带无线探头的b超仪及其实现方法 | |
Fu et al. | Barker coded excitation using LFM carrier for improving axial resolution in ultrasound imaging | |
JP2004222824A (ja) | 超音波診断装置 | |
JP2003325508A (ja) | 超音波診断装置 | |
Gong et al. | Delay-encoded transmission in synthetic transmit aperture (DE-STA) imaging | |
Kim et al. | Multiple transmit focusing using modified orthogonal Golay codes for small scale systems | |
CN110433406B (zh) | 超声聚焦定位图像优化方法 | |
CN108318883A (zh) | 一种基于全频谱换能器自然加权的距离向计算方法 | |
JP2003310612A (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 |