CN103006272B - 基于超声交织编程的速度分布测量方法 - Google Patents
基于超声交织编程的速度分布测量方法 Download PDFInfo
- Publication number
- CN103006272B CN103006272B CN201310007683.9A CN201310007683A CN103006272B CN 103006272 B CN103006272 B CN 103006272B CN 201310007683 A CN201310007683 A CN 201310007683A CN 103006272 B CN103006272 B CN 103006272B
- Authority
- CN
- China
- Prior art keywords
- interweaves
- image
- scanning line
- time quantum
- radiofrequency signal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种基于超声交织编程的速度分布测量方法,包括如下步骤1)向待测流体中添加示踪剂;2)控制器控制高频超声换能器按照交织编程方式向外发射高频超声波信号,并采集流体中示踪剂的背向散射射频信号;3)在步骤2)的基础之上,对交织编程方法采集的交织射频图像解耦,即可获得图像对;4)在步骤3)的基础之上,将上述交织编程方法采集到的两帧射频信号图像1和图像2划分成多个分析窗口即询问窗口,接着对询问窗口进行傅里叶空间内的二维互相关运算以得到该处示踪剂的局部位移,再根据两帧图像的时间间隔Δt,则可以计算得到速度向量即位移除以时间Δt。该方法可以测定较高流速状态下流体速度场。
Description
技术领域
本发明涉及一种利用交织编程原理采集超声射频信号进行速度场测量的方法。
背景技术
超声成像测速技术是一种非侵入式流场测量技术,可以实现非透明流质和流道壁速度场的测量,能够用来精确测量临床上血液流场,实现人体内血液复杂流场分布的实时测量,对各种心血管疾病准确诊断和临床研究具有重要意义。
超声成像测速过程中所使用的线性超声探头由多个单元阵列组成,在进行图像采集时首先一组单元阵列通过高频信号发生电路激发换能器在测量区域产生超声波,然后接收超声回波信号(又被称为射频RF信号),产生超声图像中的一列。然后下一组新的单元阵列按照该方式继续采样产生超声图像的第二列,按照同样方式依次扫描,当所有超声探头单元阵列扫描流场以后将产生一张完整的超声图像,所产生的超声图像是由多条经连续取样方式产生的RF信号组成。常见的线性超声探头是由128个单元构成的阵列,采集的原始图像由128条RF信号线组成。最后通过对所获得的超声粒子射频图像作粒子图像测速法(PIV)互相关分析,获取各个位置处流体的流动速度。
然而目前连续取样方式受超声设备采样频率的限制,只能获得较低流速下速度场的分布,对于流体速度较大情况下所得计算结果误差较大。例如:利用图像采集帧率f=250Hz(目前常用超声设备图像采集帧率为f=200-500Hz之间)的超声设备,从左向右扫描人体颈动脉中血液流场,设定在流场中示踪材料(造影剂)的移动速度与人体颈动脉处实际血液流动具有相同的流速U=1m/s。对于由n=128条线组成的超声图像,获取图像的第一条线至最后一条线之间的时间差为T,即获得一帧图像所花费的时间周期为T=1/f=4ms。假设所采集得到的超声粒子图像宽度为L=5cm(实际使用的探头宽度),则其横向分辨率可表示为dx=L/n=0.39mm即每条扫描线的宽度为0.39mm。示踪粒子在T时间内的位置偏移Dx=U/f=1/250=4mm,相当于10.3倍扫描线的宽度,此较大位移对于传统的光学粒子图像测速法分析是可以接受的,但在超声成像技术中,采集到的是气泡的散斑图像,噪音较大。粒子位移过大会导致PIV互相关运算时问询窗口内找不到匹配粒子图像对,从而使根据粒子图像测速法互相关运算所得的速度存在较大误差。如要达到理想状态下粒子偏移量仅是几条线的宽度,图像采集帧率f至少需增加三倍,因此需要使用更加先进的超声设备才能满足此要求,将会大大增加硬件成本。即便昂贵设备可以实现多条线同时扫描,提高图像的采集帧率,在一帧图像采集时间内使粒子偏移量减少至几条线的宽度,却存在流动过采样的缺点,例如一个典型的心动周期持续时间约为1秒,若以750-1000Hz的频率取样将产生大量的冗余数据,因为即使图像采集帧率比上述图像采集帧率降低一个数量级,也足够分辨心脏血液动力学特性。
综上所述,目前缺乏以当前超声粒子成像硬件设备为基础,可以克服较高流速条件下粒子位移偏量过大的限制,获得流体流速较大情况下的二维流场分布的测速方法。
发明内容
本发明的目的在于提供一种基于超声交织编程的速度分布测量方法,该方法可以得到较高流速状态下流体的速度场;该方法以当前超声粒子成像硬件设备为基础,克服流速较高条件下连续采样过程中粒子位移偏移量过高的限制,不增加硬件成本,同时提高测速的动态范围。
为了实现以上目的,本发明所采用的技术方案是:
基于超声交织编程的速度分布测量方法,包括以下步骤:
1)向待测流体中添加示踪剂;
2)控制器控制高频超声换能器按照交织编程方式向外发射高频超声波信号,并采集流体中示踪剂的背向散射射频信号;
设超声成像测速系统中超声探头由n个单元构成,按线扫描的方式来记录每帧图像;设一幅图像由n条扫描线构成,其中a表示扫描线所在位置,取值范围为[1,n],b表示采集该条扫描线时所处的时间单元,取值范围为[1,n];
交织编程方式采集图像时,具体方法如下:在第1时间单元采集a=1位置处的射频信号扫描线,记为L11;在第2时间单元仍采集a=1位置处的射频信号扫描线,记为L12;在第3时间单元采集a=2位置处的射频信号扫描线,记为L23;在第4时间单元仍采集a=2位置处的射频信号扫描线,记为L24;以此类推,在第2n-1时间单元采集a=n位置处的射频信号扫描线,记为Ln(2n-1);在第2n时间单元仍采集a=n位置处的射频信号扫描线,记为Ln(2n);交织编程方法采集到的图像由(L11,L12,L23,L24,…,Ln(2n-1),Ln(2n))组成,包含2n列;
3)在步骤2)的基础之上,对交织编程方法采集的交织射频图像解耦,即可获得图像对;
所述的交织编程方法采集的交织射频图像解耦具体实现方法如下:抽取步骤2)图中第1,3,5…~2n-1时间单元所采集到的图像列构成图像1(L11,L23,L35,…Ln(2n-1)),抽取第2,4,6…2n时间单元所采集到的图像列构成图像2(L12,L24,L36,…Ln(2n));
4)在步骤3)的基础之上,将上述交织编程方法采集到的两帧射频信号图像1和图像2划分成多个分析窗口即询问窗口,接着对询问窗口进行傅里叶空间内的二维互相关运算以得到该处示踪剂的局部位移,再根据两帧图像的时间间隔Δt,则可以计算得到速度向量即位移除以时间Δt。
所述的示踪剂是造影微泡或造影剂。
高频超声换能器可以发送高频信号同时接受该反射信号,其中发送高频信号由超声探头中的高频超声换能器完成,高频超声换能器是本技术领域现有技术不做详细描述。
本发明的步骤1)与步骤4)属于本技术领域公知技术,在本发明中不作详细描述。
与现有技术相比,本发明的优点在于:本发明所提出的利用交织编程原理进行超声成像测速的方法,克服流速较高条件下连续采样过程中粒子位移偏移量过高的限制,不增加硬件成本,大大提高所测流速的动态范围,使超声测速技术的应用范围进一步扩大,不仅可以用于医学上心血管内血流动力学参数的获取,而且还可用于一般尺度非透明流体成像和多相流测量领域,例如石油工业中油气的输送,化学药物溶解监测,海岸工程学研究和分析,江河入海口及海岸线的环境管理等。
附图说明
图1是传统顺序扫描方式采集第一帧图像示意图;
图2是传统顺序扫描方式采集第二帧图像示意图;
图3是交织编程成像法获取图像示意图;
图4是计算速度向量场示意图;
图5是基于超声交织编程的速度分布测量方法的流程图;
图6是测速实验系统示意图;
图7是按照传统顺序扫描方式进行超声成像测速获取的二维速度场;
图8是采用本发明交织编程原理进行超声成像测速得到的二维速度场。
具体实施方式
下面结合附图对本发明做进一步说明。
基于超声交织编程的速度分布测量方法,包括以下步骤:
1)向待测流体中添加示踪剂,示踪剂可以是造影微泡或其它声学性能较好的造影剂。
2)通过计算机控制高频超声换能器按照交织编程的方式向外发射高频超声波信号,并采集流体中示踪剂的背向散射射频信号。
传统顺序扫描方式与交织编程方式获取射频信号图像的区别之处在于:
假定超声成像测速系统中超声探头由n个单元构成,按线扫描的方式来记录每帧图像。由于超声造影剂微泡散射的信号是通过超声声束逐条扫描,再将扫描线组合在一起形成图像,从而造成不同扫描线之间时间不同。假定一幅图像由n条扫描线构成,定义图像中任意一条扫描线为Lab,其中a表示扫描线所在位置,取值范围为[1,n],b表示采集该条扫描线时所处的时间单元,取值范围为[1,n]。
按传统方法采集图像时,如图1所示,对于第一幅图像,在第1时间单元采集a=1位置处的射频信号扫描线,记为L11;在第2时间单元采集a=2位置处射频信号扫描线,记为L22;在第3时间单元采集a=3位置处射频信号扫描线,记为L33;以此类推,在第n时间单元采集a=n位置处的射频信号扫描线,记为Lnn。如图2所示,对于第二幅图像,在第n+1时间单元采集a=1位置处的射频信号扫描线,记为L1(n+1);在第n+2时间单元采集a=2位置处的射频信号扫描线,记为L2(n+2);在第n+3时间单元采集a=3位置处的射频信号扫描线,记为L3(n+3);以此类推,在第2n时间单元采集a=n位置处的射频信号扫描线,记为Ln(2n)。因此,按传统方法采集射频信号进行分析时,图像1由(L11,L22,L33,…,Lnn)组成,图像2由(L1(n+1),L2(n+2),L3(n+3),…,Ln(2n))组成。图像对的时间间隔为n。
按交织编程成像法采集图像时,如图3所示,在第1时间单元采集a=1位置处的射频信号扫描线,记为L11;在第2时间单元仍采集a=1位置处的射频信号扫描线,记为L12;由于被测流体处于运动中,因此L11,L12上记录了不同的示踪剂背向散射射频信号。在第3时间单元采集a=2位置处的射频信号扫描线,记为L23;在第4时间单元仍采集a=2位置处的射频信号扫描线,记为L24;以此类推,在第2n-1时间单元采集a=n位置处的射频信号扫描线,记为Ln(2n-1);在第2n时间单元仍采集a=n位置处的射频信号扫描线,记为Ln(2n)。交织编程法采集到的图像由(L11,L12,L23,L24,…,Ln(2n-1),Ln(2n))组成,包含2n列。
3)对利用交织编程成像法采集的交织射频图像解耦,即可获得图像对。
对交织编程成像法采集的射频信号图解耦过程如下,抽取图3中第1,3,5…~第2n-1时间单元所采集到的图像列构成图像1(L11,L23,L35,…Ln(2n-1)),抽取第2,4,6…~第2n时间单元所采集到的图像列构成图像2(L12,L24,L36,...Ln(2n))。此时,图像对的时间间隔为1。
由上述分析可见,对于高速流动的流体,若按传统方法采集射频信号进行分析时,图像对的时间间隔为n,由于示踪粒子位移过大,在第二幅图像中无法获得相应的匹配粒子图像,而按交织编程成像法采集射频信号图并解耦分析时,图像对的时间间隔仅为1,完全可以由图像对获得匹配粒子对。因此采用交织编程成像法进行分析,可以克服较高流速下粒子位移偏移量过大的限制,适用于更宽的流体测速动态范围。
4)将图像划分成多个分析窗口,所有分析窗口进行互相关运算,获得流场速度场向量图以及速度梯度分布和壁面剪切力分布。
如图4所示,将上述交织编程成像法采集得到的两帧射频信号图像划分成多个分析窗口(询问窗口),接着对询问窗口进行傅里叶空间内的二维互相关运算以得到该处示踪剂的局部位移。如果两帧图像的时间间隔Δt已知,则可以计算得到速度向量(位移除以时间Δt)。对所有的分析窗口进行互相关运算,就可以得到二维的速度向量场。速度梯度分布和壁面剪切力分布可由速度向量场经计算获得。
图5是基于超声交织编程的速度分布测量方法的流程图。基于超声交织编程的速度分布测量方法包括以下步骤。
1.向待测流体中添加示踪剂。示踪剂可以是造影微泡或其它声学性能较好的造影剂。
通过计算机控制高频超声换能器按照交织编程的方式向外发射高频超声波信号并采集图像。在第1时间单元采集a=1位置处的射频信号扫描线,记为L11;在第2时间单元仍采集a=1位置处的射频信号扫描线,记为L12;在第3时间单元采集a=2位置处的射频信号扫描线,记为L23;在第4时间单元仍采集a=2位置处的射频信号扫描线,记为L24;以此类推,在第2n-1时间单元采集a=n位置处的射频信号扫描线,记为Ln(2n-1);在第2n时间单元仍采集a=n位置处的射频信号扫描线,记为Ln(2n)。交织编程法采集到的图像(L11,L12,L23,L24,…,Ln(2n-1),Ln(2n))由2n列组成。
为获得速度场计算所需的两帧图像,需对交织编程成像法采集的交织射频信号图像解耦。抽取第1,3,5…~第2n-1时间单元所采集到的图像列构成图像1(L11,L23,L35,…Ln(2n-1)),抽取第2,4,6…~第2n时间单元所采集到的图像列构成图像2(L12,L24,L36,…Ln(2n))。图像对的时间间隔为1。
将上述交织编程成像法采集得到的两帧射频信号图像划分成多个分析窗口(询问窗口)。接着对询问窗口进行傅里叶空间内的二维互相关运算以得到该处示踪剂的局部位移。如果两帧图像的时间间隔Δt已知,则可以计算得到速度向量。对所有的分析窗口进行互相关运算,就可以得到二维的速度向量场。速度梯度分布和壁面剪切力分布可由速度向量场经计算获得。
实验系统如图6所示,乳胶管被水平放置在水槽中且在水槽底部放置吸声材料以减少超声波的反射,乳胶管的内径为5mm壁厚为1mm。通过一台离心泵驱动流体在管内以恒定流量循环流动。利用超声设备获取超声射频(RF)数据,该超声设备带有一个频率带宽为5-14MHz的线性探头,探头包含128个单元。调整流量值,使管道中心处最大速度约为1.8m/s。
图7是采用传统顺序扫描方式进行超声成像测速获取的二维速度向量场。其中箭头表示速度的方向,灰度表示速度的大小。由于高流速条件下,粒子偏移量过大,所得速度场计算结果误差较大,总体流速分布杂乱无章,管道中心流场处的部分速度向量分布甚至出现了与实际流动速度方向相反的情况。
图8是采用本发明交织编程进行超声成像测速得到的二维速度向量场,准确完整的获取了全场速度矢量分布。
Claims (2)
1.基于超声交织编程的速度分布测量方法,其特征在于,包括以下步骤:
1)向待测流体中添加示踪剂,所述待测流体为一般尺度非透明流体和多相流;
2)控制器控制高频超声换能器按照交织编程方式向外发射高频超声波信号,并采集流体中示踪剂的背向散射射频信号;
设超声成像测速系统中超声探头由n个单元构成,按线扫描的方式来记录每帧图像;设一幅图像由n条扫描线构成,其中a表示扫描线所在位置,取值范围为[1,n],b表示采集该条扫描线时所处的时间单元,取值范围为[1,n];
交织编程方式采集图像时,具体方法如下:在第1时间单元采集a=1位置处的射频信号扫描线,记为L11;在第2时间单元仍采集a=1位置处的射频信号扫描线,记为L12;在第3时间单元采集a=2位置处的射频信号扫描线,记为L23;在第4时间单元仍采集a=2位置处的射频信号扫描线,记为L24;以此类推,在第2n-1时间单元采集a=n位置处的射频信号扫描线,记为Ln(2n-1);在第2n时间单元仍采集a=n位置处的射频信号扫描线,记为Ln(2n);交织编程方法采集到的图像由(L11,L12,L23,L24,…,Ln(2n-1),Ln(2n))组成,包含2n列;
3)在步骤2)的基础之上,对交织编程方法采集的交织射频图像解耦,即可获得图像对;
所述的交织编程方法采集的交织射频图像解耦具体实现方法如下:抽取步骤2)图中第1,3,5…~2n-1时间单元所采集到的图像列构成图像1(L11,L23,L35,…Ln(2n-1)),抽取第2,4,6…2n时间单元所采集到的图像列构成图像2(L12,L24,L36,…Ln(2n));
4)在步骤3)的基础之上,将上述交织编程方法采集到的两帧射频信号图像1和图像2划分成多个分析窗口即询问窗口,接着对询问窗口进行傅里叶空间内的二维互相关运算以得到询问窗口示踪剂的局部位移,再根据两帧图像的时间间隔Δt,则可以计算得到速度向量即询问窗口示踪剂的位移除以时间Δt。
2.根据权利要求1所述的基于超声交织编程的速度分布测量方法,其特征在于,所述的示踪剂是造影微泡或造影剂。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310007683.9A CN103006272B (zh) | 2013-01-09 | 2013-01-09 | 基于超声交织编程的速度分布测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310007683.9A CN103006272B (zh) | 2013-01-09 | 2013-01-09 | 基于超声交织编程的速度分布测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103006272A CN103006272A (zh) | 2013-04-03 |
CN103006272B true CN103006272B (zh) | 2014-12-03 |
Family
ID=47955756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310007683.9A Expired - Fee Related CN103006272B (zh) | 2013-01-09 | 2013-01-09 | 基于超声交织编程的速度分布测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103006272B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104546012B (zh) * | 2014-12-31 | 2018-07-24 | 中国科学院深圳先进技术研究院 | 心脏功能评估方法和装置 |
CN105411624B (zh) * | 2015-12-25 | 2018-08-10 | 中国科学院深圳先进技术研究院 | 一种超声三维流体成像与测速方法 |
CN107970041B (zh) * | 2017-12-29 | 2021-05-14 | 深圳开立生物医疗科技股份有限公司 | 对造影剂的流动方向标识的方法、装置及设备 |
CN111588410B (zh) * | 2020-05-18 | 2023-05-16 | 南京超维景生物科技有限公司 | 基于超声造影图像的图像处理方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101002690A (zh) * | 2006-01-19 | 2007-07-25 | 株式会社东芝 | 超声波检查装置 |
EP1977692A2 (en) * | 2007-04-05 | 2008-10-08 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus, breast imaging system, and breast imaging method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6709395B2 (en) * | 2002-06-25 | 2004-03-23 | Koninklijke Philips Electronics N.V. | System and method for electronically altering ultrasound scan line origin for a three-dimensional ultrasound system |
JP5738507B2 (ja) * | 2006-01-19 | 2015-06-24 | 東芝メディカルシステムズ株式会社 | 超音波プローブの軌跡表現装置及び超音波診断装置 |
JP2012139489A (ja) * | 2010-12-16 | 2012-07-26 | Toshiba Corp | 超音波診断装置及びその制御方法 |
-
2013
- 2013-01-09 CN CN201310007683.9A patent/CN103006272B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101002690A (zh) * | 2006-01-19 | 2007-07-25 | 株式会社东芝 | 超声波检查装置 |
EP1977692A2 (en) * | 2007-04-05 | 2008-10-08 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus, breast imaging system, and breast imaging method |
Non-Patent Citations (3)
Title |
---|
《一种扇形扫描三维超声成像系统三维重建方法》;杨金耀,李德来;《中国医疗器械信息》;20120731(第7期);47-50 * |
JP特开2008-86742A 2008.04.17 * |
杨金耀,李德来.《一种扇形扫描三维超声成像系统三维重建方法》.《中国医疗器械信息》.2012,(第7期),47-50. * |
Also Published As
Publication number | Publication date |
---|---|
CN103006272A (zh) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108490068B (zh) | 超声平面波扫描式多相流可视化测量装置 | |
Tan et al. | Ultrasonic Doppler technique for application to multiphase flows: A review | |
Al-Jumaili et al. | Acoustic emission source location in complex structures using full automatic delta T mapping technique | |
CN104101687B (zh) | 基于超声多普勒与电学多传感器的多相流可视化测试方法 | |
CN108490069A (zh) | 两相流分布参数双模式超声可视化测量装置 | |
CN103006272B (zh) | 基于超声交织编程的速度分布测量方法 | |
CN107102166B (zh) | 超声多普勒多相流流速分布检测设备 | |
CN102589627B (zh) | 一种用于超声波流量计的绝对传播时间测量方法 | |
CN104155358A (zh) | 超声多普勒与电学传感器组合多相流可视化测试装置 | |
CN112067691B (zh) | 油气水三相介质衰减谱融合多频超声层析成像方法 | |
Carvalho et al. | Application of the ultrasonic technique and high-speed filming for the study of the structure of air–water bubbly flows | |
CN110470744A (zh) | 多模式曲面相控阵超声层析成像装置 | |
CN109541690B (zh) | 一种浅层介质结构面松散程度评价方法 | |
KR102326149B1 (ko) | 모델-기반 이미지 재구성 방법 | |
CN105411624A (zh) | 一种超声三维流体成像与测速方法 | |
CN104665877A (zh) | 颈动脉血管局部脉搏波传播速度测量方法 | |
CN102860842A (zh) | 一种实时准静态超声弹性成像方法 | |
CN103954962B (zh) | 一种基于压缩感知的isar成像脉冲估计算法 | |
CN109884183B (zh) | 透射反射模态融合的超声层析成像方法 | |
CN205352653U (zh) | 一种低速风洞内的超声精确测距系统 | |
CN101738489B (zh) | 一种测量散射性流体横向流速的方法 | |
CN110097608A (zh) | 修正路径追踪描述的连续波超声层析成像重建方法 | |
CN109188016A (zh) | 油气水三相流分相流速声电双模态测量方法 | |
CN110706298B (zh) | 正则化加权最小二乘的透射反射双模式超声成像重建方法 | |
Hefft et al. | A step towards the live identification of pipe obstructions with the use of passive acoustic emission and supervised machine learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent for invention or patent application | ||
CB02 | Change of applicant information |
Address after: Center branch No. 3 ancient Tan Avenue in Gaochun County of Nanjing City, Jiangsu province 211300 Room 405 Applicant after: Southeast University Address before: 210096 Jiangsu city Nanjing Province four pailou No. 2 Applicant before: Southeast University |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141203 Termination date: 20180109 |