CN107753058B - 一种剪切波动态滤波方法 - Google Patents
一种剪切波动态滤波方法 Download PDFInfo
- Publication number
- CN107753058B CN107753058B CN201711176353.7A CN201711176353A CN107753058B CN 107753058 B CN107753058 B CN 107753058B CN 201711176353 A CN201711176353 A CN 201711176353A CN 107753058 B CN107753058 B CN 107753058B
- Authority
- CN
- China
- Prior art keywords
- deformation
- time
- shear wave
- detection position
- curve
- 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/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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开一种剪切波动态滤波方法,该方法中,将动态滤波技术应用于剪切波波速估计之中,通过对不同检测位置的“形变‑时间”曲线进行频谱分析,计算出该检测位置剪切波对位移贡献最大的频率范围,并根据此范围对“形变‑时间”曲线进行滤波,利用滤波后的“形变‑时间”结果进行剪切波波速估计。本发明改善了剪切波在组织传播中衰减大等因素对波速估计的影响,使得剪切波波速估计结果更加精确,提高了结果的有效性。
Description
技术领域
本发明涉及剪切波弹性成像领域,特别涉及一种剪切波动态滤波方法。
背景技术
剪切波是传播方向与介质质点的振动方向垂直的波,又称横波,S波。
剪切波弹性成像技术可实现生物组织实时硬度定量检测,为临床判断组织病变情况提供依据。它的基本原理如下:声辐射力聚焦冲击能在组织内部产生剪切波,由于剪切波在不同硬度的组织传播速度有差异,通过检测剪切波在不同位置的传播速度可以间接反映该位置的软硬情况。根据剪切波波速大小进行伪彩色映射即为剪切波弹性成像。因此,剪切波波速估计的准确性对组织弹性情况的判断至关重要。
剪切波在组织传播过程中会引起组织的位移形变,利用超快帧频成像技术观察和检测目标区域的位移形变,通过位移形变重建出观察位置的剪切波震动曲线。传统的剪切波波速估计方法为Time-of-Filght(TOF)方法,即对观察位置不同位置所重建的剪切波震动曲线进行波峰匹配,找出波峰相差时间间隔,此间隔可以理解为剪切波从该两个检测位置之间传播所用时间。而两检测位置之间的距离为已知,通过距离除以时间得出剪切波经过两个检测位置之间的平均速度。
声辐射力冲击的剪切波在组织传播时衰减非常快,距离冲击源越远,位移估计结果的信噪比就越低。而且与振动弹性成像中机械振动产生的剪切波不一样,声辐射力聚焦冲击所产生的剪切波不是单一频率的剪切波,而是含有一定频宽的信号。这些因素都会直接影响到剪切波波速估计的准确性。
发明内容
为了使剪切波波速估值更加精确,本发明提出一种剪切波动态滤波方法。
本发明为实现其目的所采用的技术方案是:一种剪切波动态滤波方法,在生物组织弹性测量过程中,选取弹性观察区域,通过声辐射力冲击激发剪切波经过该区域传播,并利用超高帧频成像装置采集该区域回波信号进行位移估计运算,得到形变估计结果后;对位移估计结果进行如下步骤:
步骤1、确定待观察深度的位置,选取至少两个横向检测位置,根据形变估计结果矩阵作出所选位置对应的“形变-时间”曲线;
步骤1、对所选深度,各横向检测位置对应的“形变-时间”曲线进行傅里叶变换,得到对应的频谱曲线;
步骤3、分别对各横向检测位置的“形变-时间”频谱曲线进行频谱分析,确定各个对应检测位置的频谱曲线的截止频率范围;
步骤4、根据各检测位置对应截止频率生成滤波器系数;
步骤5、根据各位置对应滤波器对“形变-时间”曲线进行延时间方向进行滤波;
步骤6、搜索各检测位置滤波后的“形变-时间”曲线的峰值,并记录所对应的时间;
步骤7、根据各横向检测位置和波峰值对应时间点进行直线拟合得出“时间-距离”直线;
步骤8、求出“时间-距离”直线的斜率,得到剪切波速度值。
本发明中,将动态滤波技术应用于剪切波波速估计之中,通过对不同检测位置的“形变-时间”曲线进行频谱分析,计算出该检测位置剪切波对位移贡献最大的频率范围,并根据此范围对“形变-时间”曲线进行滤波,利用滤波后的“形变-时间”结果进行剪切波波速估计。本发明改善了剪切波在组织传播中衰减大等因素对波速估计的影响,使得剪切波波速估计结果更加精确,提高了结果的有效性。
进一步的,上述的剪切波动态滤波方法中:在生成滤波器系数时,滤波器系数生成方法选择FIR滤波器窗函数法。
下面结合具体实施例对本发明作较为详细的描述。
附图说明
图1是本发明实施例1流程图。
图2剪切波位移估计结果。
图3剪切波位移信噪比。
图4各位置“形变-时间”曲线。
图5拟合“时间-距离”直线。
图6弹性质量评估显示。
具体实施方式
本实施例中,在进行剪切波波速估计之前,先对不同检测位置的剪切波位移估计结果进行动态频谱分析,根据实际中心频率和频宽确定滤波器截止频率范围,这样能有效选出该检测位置对组织位移贡献最大的剪切波频率范围。根据此滤波器对剪切波位移估计结果进行滤波,并用滤波后的结果进行剪切波波速估计。经过动态滤波后的剪切波位移估计结果能够使波峰信息更加明显,提高了剪切波波速估计结果的准确性,具体流程如图1所示。
在生物组织弹性测量过程中,选取弹性观察区域,通过声辐射力冲击激发剪切波经过该区域传播,并利用超高帧频成像技术采集该区域回波信号进行位移估计运算得到形变估计结果后,对位移估计结果进行如下步骤:
步骤1:确定待观察深度的位置,选取两个或多个横向检测位置,根据形变估计结果矩阵作出所选位置对应的“形变-时间”曲线,如图2所示。图中显示了五个不同的横向检测位置的“形变-时间”曲线,虽然位置不同但波形基本相似,频率一致,只是波峰到来的时刻不同,另外,可以看出,声辐射力冲击是在50ms。
步骤2:对各横向检测位置对应的“形变-时间”曲线进行傅里叶变换,得到对应的频谱曲线。如附图3所示,该图展示了检测位置一的曲线频谱。
步骤3:分别对各横向检测位置的“形变-时间”频谱曲线进行频谱分析,确定各个对应检测位置的频谱曲线的截止频率范围,如附图4所示。截止频率方法可根据频谱幅值确定,可选择与频谱峰值相差30dB~50dB以上的位置作为截止频率范围。
步骤4:根据各检测位置对应截止频率生成滤波器系数。滤波器系数生成方法可选择例如FIR滤波器窗函数法。
步骤5:根据各位置对应滤波器对“形变-时间”曲线进行延时间方向进行滤波。
步骤6:搜索各检测位置滤波后的“形变-时间”曲线的峰值,并记录所对应的时间,如图5所示,是位置1的曲线经过滤波后的曲线,与图2中相应的曲线相比较基本形状没有变化,说明所携带的信息没有丢失。
步骤7:根据各横向检测位置和波峰值对应时间点进行直线拟合得出“时间-距离”直线,如图6所示。
步骤8:求出“时间-距离”直线的斜率,得到剪切波速度值。
Claims (2)
1.一种剪切波动态滤波方法,其特征在于:在生物组织弹性测量过程中,选取弹性观察区域,通过声辐射力冲击激发剪切波经过该区域传播,并利用超高帧频成像装置采集该区域回波信号进行位移估计运算,得到形变估计结果后;对位移估计结果进行如下步骤:
步骤1、确定待观察深度的位置,选取至少两个横向检测位置,根据形变估计结果矩阵作出所选位置对应的“形变-时间”曲线;
步骤2、对所选深度,各横向检测位置对应的“形变-时间”曲线进行傅里叶变换,得到对应的频谱曲线;
步骤3、分别对各横向检测位置的“形变-时间”频谱曲线进行频谱分析,确定各个对应检测位置的频谱曲线的截止频率范围;该步骤中,选择与频谱峰值相差30dB~50dB以上的位置作为截止频率范围;
步骤4、根据各检测位置对应截止频率生成滤波器系数;
步骤5、根据各位置对应滤波器对“形变-时间”曲线进行延时间方向进行滤波;
步骤6、搜索各检测位置滤波后的“形变-时间”曲线的峰值,并记录所对应的时间;
步骤7、根据各横向检测位置和波峰值对应时间点进行直线拟合得出“时间-距离”直线;
步骤8、求出“时间-距离”直线的斜率,得到剪切波速度值。
2.根据权利要求1所述的剪切波动态滤波方法,其特征在于:在生成滤波器系数时,滤波器系数生成方法选择FIR滤波器窗函数法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711176353.7A CN107753058B (zh) | 2017-11-22 | 2017-11-22 | 一种剪切波动态滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711176353.7A CN107753058B (zh) | 2017-11-22 | 2017-11-22 | 一种剪切波动态滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107753058A CN107753058A (zh) | 2018-03-06 |
CN107753058B true CN107753058B (zh) | 2021-03-02 |
Family
ID=61278696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711176353.7A Active CN107753058B (zh) | 2017-11-22 | 2017-11-22 | 一种剪切波动态滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107753058B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020019254A1 (zh) * | 2018-07-26 | 2020-01-30 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波成像方法及*** |
CN109589138A (zh) * | 2018-11-26 | 2019-04-09 | 深圳中科乐普医疗技术有限公司 | 一种剪切波波速计算方法及弹性成像设备 |
CN112799078B (zh) * | 2021-04-15 | 2021-07-02 | 深圳中科乐普医疗技术有限公司 | 一种剪切波传播速度的检测方法、***和超声成像设备 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5431169A (en) * | 1992-08-03 | 1995-07-11 | Olympus Optical Co., Ltd. | Ultrasonic diagnosing apparatus |
US6850168B2 (en) * | 2000-11-13 | 2005-02-01 | Baker Hughes Incorporated | Method and apparatus for LWD shear velocity measurement |
FI20065217L (fi) * | 2006-04-03 | 2007-10-04 | Metso Automation Oy | Menetelmä synkronoidun aikakeskiarvon laskemiseksi |
CN1976224A (zh) * | 2006-12-15 | 2007-06-06 | 天津大学 | 低通连续可调滤波器和滤波方法 |
CN101179545B (zh) * | 2007-12-20 | 2010-06-09 | 清华大学 | 基于多普勒频率对消的全数字主载波跟踪方法 |
US8456952B2 (en) * | 2008-10-03 | 2013-06-04 | Baker Hughes Incorporated | Curve-fitting technique for determining dispersion characteristics of guided elastic waves |
CN101567676B (zh) * | 2009-05-08 | 2012-05-23 | 江苏大学 | 一种滤波方法 |
US8734352B2 (en) * | 2009-07-17 | 2014-05-27 | Koninklijke Philips N.V. | Spatially-fine shear wave dispersion ultrasound vibrometry sampling |
US10172527B2 (en) * | 2009-07-31 | 2019-01-08 | Supersonic Imagine | Method and apparatus for measuring a physical parameter in mammal soft tissues by propagating shear waves |
CN201682468U (zh) * | 2010-02-11 | 2010-12-22 | 北京国能子金电气技术有限公司 | 变频率信号的自适应低通数字滤波器 |
US8944183B2 (en) * | 2010-08-11 | 2015-02-03 | Baker Hughes Incorporated | Low frequency formation shear slowness from drilling noise derived quadrupole array data |
US8469891B2 (en) * | 2011-02-17 | 2013-06-25 | Siemens Medical Solutions Usa, Inc. | Viscoelasticity measurement using amplitude-phase modulated ultrasound wave |
CN103399348B (zh) * | 2013-08-15 | 2015-09-30 | 电子科技大学 | 基于Shearlet变换的地震信号去噪方法 |
US10278671B2 (en) * | 2013-09-30 | 2019-05-07 | Siemens Medical Solutions Usa, Inc. | Shear wave detection in medical ultrasound imaging |
CN105266843B (zh) * | 2014-07-23 | 2018-06-08 | 通用电气公司 | 超声***及使用超声***获得目标组织的二维弹性图的方法 |
CN104546014A (zh) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | 一种用于生物组织弹性测量的剪切波速度估计方法 |
CN105212961B (zh) * | 2015-08-20 | 2018-08-31 | 深圳市红源资产管理有限公司 | 一种声辐射剪切波波速检测方法及*** |
CN106618638B (zh) * | 2016-11-04 | 2019-02-26 | 声泰特(成都)科技有限公司 | 一种定量剪切波弹性成像*** |
-
2017
- 2017-11-22 CN CN201711176353.7A patent/CN107753058B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107753058A (zh) | 2018-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2688299C1 (ru) | Способ и устройство для определения параметра вязкоупругости вязкоупругой среды | |
KR101398948B1 (ko) | 진폭―위상 변조된 초음파를 이용한 점탄성 측정 | |
CN107753058B (zh) | 一种剪切波动态滤波方法 | |
CN107510474B (zh) | 剪切波弹性成像方法及*** | |
CN103269639B (zh) | 利用质心估计剪切波速度 | |
CN103293515B (zh) | 一种舰船线谱噪声源纵向分布特性测量方法 | |
US8864671B2 (en) | Methods and systems for color flow imaging | |
CN103380386A (zh) | 超声波测量 | |
EP2903530B1 (en) | Shear wave attenuation from k-space analysis system | |
CN101784234A (zh) | 用于测量感兴趣区域的平均粘弹值的方法和设备 | |
CN111337171B (zh) | 一种应用于临界折射纵波应力检测的声时差测量方法 | |
JPS59174152A (ja) | 超音波媒体特性値測定方式 | |
Kiefer et al. | Simultaneous ultrasonic measurement of thickness and speed of sound in elastic plates using coded excitation signals | |
CN109512465B (zh) | 一种声辐射力双向剪切波复合成像方法及装置 | |
Malo et al. | Wave mode discrimination of coded ultrasonic guided waves using two-dimensional compressed pulse analysis | |
CN104897779B (zh) | 利用线性调频信号测量超声波传播时间的方法 | |
JP6144038B2 (ja) | 非接触音響検査装置および非接触音響検査方法 | |
CN104897777A (zh) | 基于Burg算法的自回归谱外推技术提高TOFD检测纵向分辨率的方法 | |
NO338715B1 (no) | Fremgangsmåte for karakterisering av en suspensjon | |
CN111388012B (zh) | 用于检测组织硬度的方法、设备及*** | |
JP4534309B2 (ja) | 金属薄板の厚み共振スペクトル測定方法及び金属薄板の電磁超音波計測方法 | |
CN112799078B (zh) | 一种剪切波传播速度的检测方法、***和超声成像设备 | |
CN103033154B (zh) | 一种用于提高超声检测时间分辨率的自回归谱分析方法 | |
US20170156701A1 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
CN110108797A (zh) | 利用声阻抗变化信息的介质分界面超声检测方法 |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 405, building 9, baiwangxin Industrial Park, No. 1002, Songbai Road, Yangguang community, Xili street, Nanshan District, Shenzhen City, Guangdong Province Applicant after: Shenzhen Zhongke Medical Technology Co., Ltd. Address before: Baoan District Shiyan street Shenzhen city Guangdong province 518000 White Pine Road in Yuntai technology industrial factory building 9 building seven floor East Applicant before: Shenzhen Zhongke Medical Technology Co., Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |