CN101846693A - 超声粒子图像测速系统和超声粒子图像测速方法 - Google Patents
超声粒子图像测速系统和超声粒子图像测速方法 Download PDFInfo
- Publication number
- CN101846693A CN101846693A CN200910106237A CN200910106237A CN101846693A CN 101846693 A CN101846693 A CN 101846693A CN 200910106237 A CN200910106237 A CN 200910106237A CN 200910106237 A CN200910106237 A CN 200910106237A CN 101846693 A CN101846693 A CN 101846693A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- frequency
- particle image
- fluid
- image velocimetry
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
本发明涉及一种超声粒子图像测速系统和超声粒子图像测速方法。所述超声粒子图像测速系统包括注射装置、高频超声换能器、高频超声电子系统和计算机,所述注射装置用于向流体中添加示踪剂,所述高频超声电子系统激励高频超声换能器发射高频超声波信号,所述高频超声波信号被流体中的示踪剂背向散射,所述高频超声换能器接收背向散射射频信号,所述高频超声电子系统对背向散射射频信号进行处理,得到亮度模式超声图像,所述计算机对相邻两帧亮度模式超声图像进行粒子图像测速分析,得到流场速度向量图和流体动力学信息。所述超声粒子图像测速系统及超声粒子图像测速方法可以用于测量微尺度非透明流体。
Description
技术领域
本发明涉及一种超声粒子图像测速系统和超声粒子图像测速方法。
背景技术
目前比较成熟的微尺度流场成像和测量技术是光学显微粒子成像技术。光学显微粒子成像技术是利用激光照射在播散于流体并跟随流体运动的微米级荧光示踪粒子,荧光示踪粒子被入射光激发后发射出的荧光经显微物镜、三棱镜、滤光镜和中继镜后进入CCD相机进行成像,由计算机进行图像分析和处理,从而获得速度场分布。但是光学显微粒子成像技术仅适用于测量光学透明的微尺度流场,而对于非透明的流道壁或流质,该技术无法适用。
发明内容
有鉴于此,有必要针对光学显微粒子成像技术不能测量微尺度非透明流体的问题,提供一种可以测量微尺度非透明流体的超声粒子图像测速系统。
一种超声粒子图像测速系统,包括注射装置、高频超声换能器、高频超声电子系统和计算机,所述注射装置用于向流体中添加示踪剂,所述高频超声电子系统激励高频超声换能器发射高频超声波信号,所述高频超声波信号被流体中的示踪剂背向散射,所述高频超声换能器接收背向散射射频信号,所述高频超声电子系统对背向散射射频信号进行处理,得到亮度模式超声图像,所述计算机对相邻两帧亮度模式超声图像进行粒子图像测速分析,得到流场速度向量图和流体动力学信息。
优选的,所述流体动力学信息包括流场剪切力分布和速度梯度分布。
优选的,所述高频超声换能器是线阵高频超声换能器、扇扫高频超声换能器和线扫高频超声换能器中的一种。
优选的,所述高频超声电子系统包括信号发生器、功率放大器、保护电路、衰减器、带通滤波器、放大器和数据采集卡,所述信号发生器在计算机的控制下产生激励信号,激励信号由功率放大器放大后经保护电路传送至高频超声换能器,背向散射射频信号经过衰减器、带通滤波器和放大器处理后由数据采集卡传送至计算机。
优选的,所述流体为非透明流体。
还提供一种超声粒子图像测速方法。
一种超声粒子图像测速方法,包括:向待测流体中添加示踪剂;发射高频超声波信号;采集背向散射射频信号;获得亮度模式超声粒子图像;获得流场速度向量图和流体动力学信息。
优选的,所述流体动力学信息包括流场剪切力分布和速度梯度分布。
优选的,所述示踪剂是造影微泡或造影微粒。
优选的,所述高频超声波信号的频率范围是20MHz~100MHz。
优选的,所述流体为非透明流体。
上述超声粒子图像测速系统采用声学的手段得到亮度模式超声图像,并最终得到速度向量图,获得流场剪切力分布等流体信息。因为不是采用光学的手段,所以上述超声粒子图像测速系统及超声粒子图像测速方法可以用于测量微尺度非透明流体。
附图说明
图1是超声粒子图像测速系统的示意图。
图2是超声粒子图像测速方法的流程图。
图3是连续的两帧亮度模式超声粒子图像的示意图。
图4是二维的速度向量场的示意图。
图5是直径为600微米的导管内层流流体的速度向量图。
图6是弯曲流体的速度向量图。
图7是导管内局部狭窄部位的流体速度向量图。
具体实施方式
图1是超声粒子图像测速系统的示意图。超声粒子图像测速系统100包括注射装置101、高频超声换能器103、高频超声电子系统105和计算机107。注射装置101用于向流体中添加示踪剂。高频超声电子系统105激励高频超声换能器103发射高频超声波信号,高频超声波信号被流体中的示踪剂背向散射,高频超声换能器103接收背向散射射频信号,高频超声电子系统105对背向散射射频信号进行处理,得到亮度模式超声图像,计算机107对相邻两帧亮度模式超声图像进行粒子图像测速分析,最终得到速度向量图和流体动力学信息。
高频超声换能器103可以是线阵高频超声换能器、扇扫高频超声换能器、线扫高频超声换能器或者其他类型的高频超声换能器。待测微流场的尺寸决定了超声粒子图像测速系统100所需的纵向分辨率。高频超声换能器103的信号发射特性(比如脉冲宽度、带宽)也与空间分辨率的需求有关。高的横向分辨率意味着较高的图像质量,而低的横向分辨率将限制较小速度的识别能力。视场的设置基于感兴趣的流道几何尺寸的大小。对于弯曲或者有分支的流道,一般需要较大的视场,目的是为了更好的表示由于局部几何特征的变化引起的速度向量的变化。采用较大的视场将会影响帧频,同样也会影响可测量的速度动态范围。
高频超声电子系统105主要实现驱动超声探头和接受超声背向散射射频信号两大功能。高频超声电子系统105主要包括信号发生器、功率放大器、保护电路、衰减器、带通滤波器、放大器、和数据采集卡。高频超声电子系统105的工作流程是:计算机107控制信号发生器产生激励信号,由功率放大器放大后经保护电路激励高频超声换能器103发射高频超声波信号。高频超声波信号被流体中的示踪剂背向散射,即产生背向散射射频信号。背向散射射频信号经过衰减器、带通滤波器、放大器等硬件电路处理,最终由集成在计算机107上的数据采集卡采集到计算机107内存中,再经过一系列的信号处理过程得到所需的超声粒子图像。
图2是超声粒子图像测速方法的流程图。超声粒子图像测速方法包括以下步骤:
S201:向待测流体中添加示踪剂。
示踪剂可以是造影微泡或其他声学性能较好的造影剂。为了得到质量最优的超声粒子图像测速数据,必须选择合适的系统参数,比如窗口宽度、成像深度、聚焦深度以及互相关分析窗口的大小和重叠率等等。而示踪剂的浓度在超声粒子图像测速技术中也扮演着重要的角色,合适的示踪剂浓度下得到的亮度模式超声图像对于超声粒子图像测速分析是非常重要的。利用从互相关函数(用来标示两个分析窗口之间模式匹配的程度)得到的互相关指数来评价示踪剂浓度对于超声粒子图像测速结果的影响。对于质量较好的超声粒子图像测速结果,互相关指数一般介于0.2~0.8之间。互相关指数可以用来实时地指示成像时最佳的示踪剂浓度:当互相关指数曲线指示已经达到最佳的浓度时,超声粒子图像测速系统100就可以开始采集数据了。
S202:发射高频超声波信号。
计算机107控制信号发生器产生激励信号,由功率放大器放大后经保护电路激励高频超声换能器103发射高频超声波信号。高频超声波信号的频率范围是20MHz~100MHz。
S203:采集背向散射射频信号。
高频超声波信号被流体中的示踪剂背向散射,即产生背向散射射频信号。背向散射射频信号由高频超声换能器103所接收。
S204:获得亮度模式超声粒子图像。
对超声背向散射射频信号进行整合以产生该点时间和空间上所需的信号强度。分析处理射频信号并提取基波成分及其他谐波成分(例如分频谐波、超频谐波或二次谐波等),最终获得亮度模式超声粒子图像,进行超声基波成像或超声谐波成像。利用造影微泡或其他造影剂产生的超声谐波射频数据,目的就是将造影微泡或其他造影剂的背向散射射频信号从组织的背向散射信号中分离出来,提高超声成像的敏锐度和信噪比。
S205:获得流场速度向量图和流体动力学信息。
首先,连续的两帧亮度模式超声粒子图像被划分成多个分析窗口(次窗口),如图3所示。接着对次窗口进行傅立叶空间的二维互相关运算以得到该处微粒的局部位移。在两帧图像的时间间隔Δt已知的情况下,可以计算得到速度向量。对所有的分析窗口进行互相关运算,就可以得到二维的速度向量场,如图4所示。图5是直径为600微米的导管内层流流体的速度向量图。从图5可以看出,显微超声粒子图像测速得到的导管纵向速度分布满足层流流体的速度分布,这也证实了超声粒子图像测速技术的可行性和可靠性。图6是弯曲流体的速度向量图。图7是导管内局部狭窄部位的流体速度向量图。在获得速度矢量场分布的基础上,可以通过流体力学公式计算出流场剪切力分布和速度梯度分布等流体动力学信息。
上述超声粒子图像测速方法和超声粒子图像测速系统100同时具有较好的空间和时间分辨率,可以捕获瞬变流场信息;不依赖声束和速度之间的夹角,可以提供二维速度向量场,并能准确得到剪切应力的分布;操作简单,安全便捷;是一种非侵入、高分辨率(包括时间分辨率和空间分辨率)、多维、精确的微尺度流体实时成像和测量技术。
上述超声粒子图像测速技术和超声粒子图像测速系统100可以提供实时、非侵入的微细血管内多维血流速度成像和测量,可根据需要提供定量的血流动力学信息,如血流速度、剪切率、剪切应力及流线图等,这在心血管疾病的发生发展、颈动脉易损斑块的评价、动脉瘤破裂危险性的预测等方面都有重要意义。
超声粒子图像测速技术和超声粒子图像测速系统100不仅可用于微细血管内血液成像和测量领域,在微尺度非透明流体流场无干扰测量的非临床应用成像方面也非常有用。例如生物微流体系统、微结构设备(如微机电系统)中的流体、微导管流体等。而且还可用于一般微尺度非透明流体成像和测量领域,例如石油的加工,饮料(啤酒、可乐、葡萄酒、香槟、果汁、牛奶、豆浆等)、香水、墨水、自来水、染料、胶水及某些塑料的生产和加工,化学药物溶解监测,海岸工程学研究和分析,江河入海口及海岸线的环境管理等。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种超声粒子图像测速系统,其特征在于:包括注射装置、高频超声换能器、高频超声电子系统、计算机;所述注射装置用于向流体中添加示踪剂;所述高频超声电子系统激励高频超声换能器发射高频超声波信号;所述高频超声波信号被流体中的示踪剂背向散射;所述高频超声换能器接收背向散射射频信号;所述高频超声电子系统对背向散射射频信号进行处理,得到亮度模式超声图像;所述计算机对相邻两帧亮度模式超声图像进行粒子图像测速分析,得到流场速度向量图和流体动力学信息。
2.根据权利要求1所述的超声粒子图像测速系统,其特征在于:所述流体动力学信息包括流场剪切力分布和速度梯度分布。
3.根据权利要求1所述的超声粒子图像测速系统,其特征在于:所述高频超声换能器是线阵高频超声换能器、扇扫高频超声换能器和线扫高频超声换能器中的一种。
4.根据权利要求1所述的超声粒子图像测速系统,其特征在于:所述高频超声电子系统包括信号发生器、功率放大器、保护电路、衰减器、带通滤波器、放大器、和数据采集卡,所述信号发生器在计算机的控制下产生激励信号,激励信号由功率放大器放大后经保护电路传送至高频超声换能器,背向散射射频信号经过衰减器、带通滤波器和放大器处理后由数据采集卡传送至计算机。
5.根据权利要求1所述的超声粒子图像测速系统,其特征在于:所述流体为非透明流体。
6.一种超声粒子图像测速方法,包括:
向待测流体中添加示踪剂;
发射高频超声波信号;
采集背向散射射频信号;
获得亮度模式超声粒子图像;
获得流场速度向量图和流体动力学信息。
7.根据权利要求6所述的超声粒子图像测速方法,其特征在于:所述流体动力学信息包括流场剪切力分布和速度梯度分布。
8.根据权利要求6所述的超声粒子图像测速方法,其特征在于:所述示踪剂是造影微泡或造影微粒。
9.根据权利要求6所述的超声粒子图像测速方法,其特征在于:所述高频超声波信号的频率范围是20MHz~100MHz。
10.根据权利要求6所述的超声粒子图像测速方法,其特征在于:所述流体为非透明流体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910106237 CN101846693B (zh) | 2009-03-26 | 2009-03-26 | 超声粒子图像测速系统和超声粒子图像测速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910106237 CN101846693B (zh) | 2009-03-26 | 2009-03-26 | 超声粒子图像测速系统和超声粒子图像测速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101846693A true CN101846693A (zh) | 2010-09-29 |
CN101846693B CN101846693B (zh) | 2013-08-21 |
Family
ID=42771378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910106237 Active CN101846693B (zh) | 2009-03-26 | 2009-03-26 | 超声粒子图像测速系统和超声粒子图像测速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101846693B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102151152A (zh) * | 2011-03-01 | 2011-08-17 | 深圳市一体医疗科技股份有限公司 | 一种用于测量粘弹性介质弹性的测量探头、系统及方法 |
CN102331511A (zh) * | 2011-06-17 | 2012-01-25 | 河海大学 | Piv图像高频采集方法 |
CN102818690A (zh) * | 2012-01-15 | 2012-12-12 | 洛阳理工学院 | 一种对消声器内部流场及压力损失分析的试验装置 |
CN103698554A (zh) * | 2013-12-17 | 2014-04-02 | 华中科技大学 | 一种流场实时精确测量系统及方法 |
WO2014190541A1 (zh) * | 2013-05-31 | 2014-12-04 | 中国科学院深圳先进技术研究院 | 一种管壁应力相位角的测量方法和系统 |
CN105411624A (zh) * | 2015-12-25 | 2016-03-23 | 中国科学院深圳先进技术研究院 | 一种超声三维流体成像与测速方法 |
CN106919111A (zh) * | 2017-04-11 | 2017-07-04 | 上海瑞影医疗科技有限公司 | 一种高频超声设备及射频信号处理装置 |
CN108122207A (zh) * | 2016-11-30 | 2018-06-05 | 展讯通信(上海)有限公司 | 图像分频方法、装置及电子设备 |
CN109163769A (zh) * | 2018-09-03 | 2019-01-08 | 山西省地质矿产研究院(山西省煤层气测试技术研究院) | 一种管道流量电磁阵列传感器及其检测方法 |
CN110811687A (zh) * | 2015-06-05 | 2020-02-21 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体成像方法及超声流体成像系统 |
CN111413516A (zh) * | 2020-03-11 | 2020-07-14 | 中南大学 | 一种基于极化特征的铁水流速检测装置 |
CN114563077A (zh) * | 2020-11-27 | 2022-05-31 | 嘉善飞阔医疗科技有限公司 | 一种基于piv的超声手术刀输出声功率测量系统 |
CN115128299A (zh) * | 2022-08-31 | 2022-09-30 | 之江实验室 | 一种测量非透明流场的光声粒子图像测速系统和方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5855557A (en) * | 1996-12-30 | 1999-01-05 | Siemens Medical Systems, Inc. | Ultrasonic imaging system and method for generating and displaying velocity field information |
JP4422421B2 (ja) * | 2003-03-17 | 2010-02-24 | 株式会社日立メディコ | 超音波撮像装置 |
-
2009
- 2009-03-26 CN CN 200910106237 patent/CN101846693B/zh active Active
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102151152A (zh) * | 2011-03-01 | 2011-08-17 | 深圳市一体医疗科技股份有限公司 | 一种用于测量粘弹性介质弹性的测量探头、系统及方法 |
CN102331511A (zh) * | 2011-06-17 | 2012-01-25 | 河海大学 | Piv图像高频采集方法 |
CN102331511B (zh) * | 2011-06-17 | 2014-05-07 | 河海大学 | Piv图像高频采集方法 |
CN102818690A (zh) * | 2012-01-15 | 2012-12-12 | 洛阳理工学院 | 一种对消声器内部流场及压力损失分析的试验装置 |
CN102818690B (zh) * | 2012-01-15 | 2015-04-29 | 洛阳理工学院 | 一种对消声器内部流场及压力损失分析的试验装置 |
WO2014190541A1 (zh) * | 2013-05-31 | 2014-12-04 | 中国科学院深圳先进技术研究院 | 一种管壁应力相位角的测量方法和系统 |
CN103698554A (zh) * | 2013-12-17 | 2014-04-02 | 华中科技大学 | 一种流场实时精确测量系统及方法 |
CN103698554B (zh) * | 2013-12-17 | 2017-04-19 | 华中科技大学 | 一种流场实时精确测量系统及方法 |
CN110811687A (zh) * | 2015-06-05 | 2020-02-21 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体成像方法及超声流体成像系统 |
CN110811687B (zh) * | 2015-06-05 | 2022-04-22 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体成像方法及超声流体成像系统 |
CN105411624A (zh) * | 2015-12-25 | 2016-03-23 | 中国科学院深圳先进技术研究院 | 一种超声三维流体成像与测速方法 |
CN105411624B (zh) * | 2015-12-25 | 2018-08-10 | 中国科学院深圳先进技术研究院 | 一种超声三维流体成像与测速方法 |
CN108122207A (zh) * | 2016-11-30 | 2018-06-05 | 展讯通信(上海)有限公司 | 图像分频方法、装置及电子设备 |
CN106919111A (zh) * | 2017-04-11 | 2017-07-04 | 上海瑞影医疗科技有限公司 | 一种高频超声设备及射频信号处理装置 |
CN109163769B (zh) * | 2018-09-03 | 2020-02-07 | 山西省地质矿产研究院(山西省煤层气测试技术研究院) | 一种管道流量电磁阵列传感器的检测方法 |
CN109163769A (zh) * | 2018-09-03 | 2019-01-08 | 山西省地质矿产研究院(山西省煤层气测试技术研究院) | 一种管道流量电磁阵列传感器及其检测方法 |
CN111413516A (zh) * | 2020-03-11 | 2020-07-14 | 中南大学 | 一种基于极化特征的铁水流速检测装置 |
CN114563077A (zh) * | 2020-11-27 | 2022-05-31 | 嘉善飞阔医疗科技有限公司 | 一种基于piv的超声手术刀输出声功率测量系统 |
CN114563077B (zh) * | 2020-11-27 | 2023-10-13 | 嘉善飞阔医疗科技有限公司 | 一种基于piv的超声手术刀输出声功率测量系统 |
CN115128299A (zh) * | 2022-08-31 | 2022-09-30 | 之江实验室 | 一种测量非透明流场的光声粒子图像测速系统和方法 |
CN115128299B (zh) * | 2022-08-31 | 2022-12-27 | 之江实验室 | 一种测量非透明流场的光声粒子图像测速系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101846693B (zh) | 2013-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101846693B (zh) | 超声粒子图像测速系统和超声粒子图像测速方法 | |
Poelma et al. | 3D Flow reconstruction using ultrasound PIV | |
CN205374308U (zh) | 多功能肿瘤细胞或其它病理细胞检测装置 | |
Shelton et al. | Volumetric imaging of erythrocytes using label‐free multiphoton photoacoustic microscopy | |
CN110179446A (zh) | 一种联合光声与激光散斑的多模态成像设备 | |
CN105222833A (zh) | 气液两相流中气泡大小、数目和运动速度的测量方法 | |
CN103054610B (zh) | 无超声换能器频带限制的光声成像装置及其检测方法 | |
CN102854142A (zh) | 基于光束扫描的光学分辨式光声显微镜 | |
CN115128299B (zh) | 一种测量非透明流场的光声粒子图像测速系统和方法 | |
CN105572049A (zh) | 光声定量弹性成像方法及装置 | |
CN102322805A (zh) | 空化泡最大泡半径的探测装置及其探测方法 | |
CN109946236A (zh) | 一种偏振光声显微成像装置及成像方法 | |
CN115901178B (zh) | 多体海工结构间波浪共振流场特性的测量系统和分析方法 | |
He et al. | In vivo cell characteristic extraction and identification by photoacoustic flow cytography | |
McDicken et al. | Physics: principles, practice and artefacts | |
Li et al. | Separation of free and bound harmonics in waves | |
KR20110118723A (ko) | 형광 검출 방법, 형광 검출 장치 및 프로그램을 기록한 컴퓨터로 읽을 수 있는 매체 | |
Yoshida et al. | Comparable analysis of bubble translation due to acoustic radiation force based on simultaneous acoustical and optical observation | |
CN111938579A (zh) | 一种基于光纤延时的提升光声图像质量的多脉冲光声显微成像方法及其装置 | |
Casey et al. | Single bubble acoustic characterization and stability measurement of adherent microbubbles | |
Brunker et al. | Acoustic resolution photoacoustic Doppler flowmetry: practical considerations for obtaining accurate measurements of blood flow | |
CN110353624A (zh) | 一种基于声子晶体共振技术放大眼角膜散射信号的方法 | |
CN108051369A (zh) | 一种非接触式全光学光声成像装置及其方法 | |
Coutinho et al. | A new method for ultrasound detection of interfacial position in gas-liquid two-phase flow | |
CN107449379A (zh) | 用于sim卡芯片的无损测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160106 Address after: 518000 Guangdong city of Shenzhen province Nanshan District Guangdong streets Nanshan Science Park Road 16 No. 16 Taiho Technology Building 2203 Patentee after: Shenzhen Zhongke Medical Technology Co., Ltd. Address before: 518067 Nanshan medical equipment Industrial Park, No. 1019 Nanhai Road, Shenzhen, Guangdong, Nanshan District 3A Patentee before: Shenzhen Advanced Technology Research Inst. |