基于超声交织编程的速度分布测量方法技术

本发明专利技术公开了一种基于超声交织编程的速度分布测量方法，包括如下步骤1）向待测流体中添加示踪剂；2）控制器控制高频超声换能器按照交织编程方式向外发射高频超声波信号，并采集流体中示踪剂的背向散射射频信号；3）在步骤2）的基础之上，对交织编程方法采集的交织射频图像解耦，即可获得图像对；4）在步骤3）的基础之上，将上述交织编程方法采集到的两帧射频信号图像1和图像2划分成多个分析窗口即询问窗口，接着对询问窗口进行傅里叶空间内的二维互相关运算以得到该处示踪剂的局部位移，再根据两帧图像的时间间隔Δt，则可以计算得到速度向量即位移除以时间Δt。该方法可以测定较高流速状态下流体速度场。

【技术实现步骤摘要】

本专利技术涉及一种利用交织编程原理采集超声射频信号进行速度场测量的方法。
技术介绍
超声成像测速技术是一种非侵入式流场测量技术，可以实现非透明流质和流道壁速度场的测量，能够用来精确测量临床上血液流场，实现人体内血液复杂流场分布的实时测量，对各种心血管疾病准确诊断和临床研究具有重要意义。超声成像测速过程中所使用的线性超声探头由多个单元阵列组成，在进行图像采集时首先一组单元阵列通过高频信号发生电路激发换能器在测量区域产生超声波，然后接收超声回波信号(又被称为射频RF信号)，产生超声图像中的一列。然后下一组新的单元阵 列按照该方式继续采样产生超声图像的第二列，按照同样方式依次扫描，当所有超声探头单元阵列扫描流场以后将产生一张完整的超声图像，所产生的超声图像是由多条经连续取样方式产生的RF信号组成。常见的线性超声探头是由128个单元构成的阵列，采集的原始图像由128条RF信号线组成。最后通过对所获得的超声粒子射频图像作粒子图像测速法(PIV)互相关分析，获取各个位置处流体的流动速度。然而目前连续取样方式受超声设备采样频率的限制，只能获得较低流速下速度场的分布，对于流体速度较大情况下所得计算结果误差较大。例如利用图像采集帧率f=250Hz(目前常用超声设备图像采集帧率为f=200-500Hz之间)的超声设备，从左向右扫描人体颈动脉中血液流场，设定在流场中示踪材料(造影剂)的移动速度与人体颈动脉处实际血液流动具有相同的流速U=lm/s。对于由n=128条线组成的超声图像,获取图像的第一条线至最后一条线之间的时间差为T，即获得一帧图像所花费的时间周期为T=l/f=4ms。假设所采集得到的超声粒子图像宽度为L=5cm(实际使用的探头宽度)，则其横向分辨率可表示为dx=L/n=0. 39mm即每条扫描线的宽度为O. 39mm。示踪粒子在T时间内的位置偏移Dx=U/f=l/250=4mm，相当于10. 3倍扫描线的宽度，此较大位移对于传统的光学粒子图像测速法分析是可以接受的，但在超声成像技术中，采集到的是气泡的散斑图像，噪音较大。粒子位移过大会导致PIV互相关运算时问询窗口内找不到匹配粒子图像对，从而使根据粒子图像测速法互相关运算所得的速度存在较大误差。如要达到理想状态下粒子偏移量仅是几条线的宽度，图像采集帧率f至少需增加三倍，因此需要使用更加先进的超声设备才能满足此要求，将会大大增加硬件成本。即便昂贵设备可以实现多条线同时扫描，提高图像的采集帧率，在一帧图像采集时间内使粒子偏移量减少至几条线的宽度，却存在流动过采样的缺点，例如一个典型的心动周期持续时间约为I秒，若以750-1000HZ的频率取样将产生大量的冗余数据，因为即使图像采集帧率比上述图像采集帧率降低一个数量级，也足够分辨心脏血液动力学特性。综上所述，目前缺乏以当前超声粒子成像硬件设备为基础，可以克服较高流速条件下粒子位移偏量过大的限制，获得流体流速较大情况下的二维流场分布的测速方法。
技术实现思路
本专利技术的目的在于提供一种，该方法可以得到较高流速状态下流体的速度场；该方法以当前超声粒子成像硬件设备为基础，克服流速较高条件下连续采样过程中粒子位移偏移量过高的限制，不增加硬件成本，同时提高测速的动态范围。为了实现以上目的，本专利技术所采用的技术方案是，包括以下步骤I)向待测流体中添加示踪剂；2)控制器控制高频超声换能器按照交织编程方式向外发射高频超声波信号，并采集流体中示踪剂的背向散射射频信号；设超声成像测速系统中超声探头由η个单元构成，按线扫描的方式来记录每帧图像；设一幅图像由η条扫描线构成，其中a表示扫描线所在位置，取值范围为[1，n]，b表示采集该条扫描线时所处的时间单元，取值范围为[l，n];交织编程方式采集图像时，具体方法如下在第I时间单元采集a=l位置处的射频信号扫描线，记为Lll ;在第2时间单元仍采集a=l位置处的射频信号扫描线，记为L12 ；在第3时间单元采集a=2位置处的射频信号扫描线，记为L23 ;在第4时间单元仍采集a=2位置处的射频信号扫描线，记为L24 ;以此类推，在第2n-l时间单元采集a=n位置处的射频信号扫描线，记为Ln (2n-l);在第2n时间单元仍采集a=n位置处的射频信号扫描线，记为Ln (2n);交织编程方法采集到的图像由(Lll，L12, L23, L24,…，Ln (2n_l)，Ln (2n))组成，包含2n列；3)在步骤2)的基础之上，对交织编程方法采集的交织射频图像解耦，即可获得图像对；所述的交织编程方法采集的交织射频图像解耦具体实现方法如下抽取步骤2) 图中第1，3，5·:2η-1时间单元所采集到的图像列构成图像1(L11，L23，L35，…Ln(2n_l))，抽取第2，4，6吣211时间单元所采集到的图像列构成图像2 (L12，L24，L36，…Ln(2n));4)在步骤3)的基础之上，将上述交织编程方法采集到的两帧射频信号图像I和图像2划分成多个分析窗口即询问窗口，接着对询问窗口进行傅里叶空间内的二维互相关运算以得到该处示踪剂的局部位移，再根据两帧图像的时间间隔At，则可以计算得到速度向量即位移除以时间At。所述的示踪剂是造影微泡或造影剂。高频超声换能器可以发送高频信号同时接受该反射信号，其中发送高频信号由超声探头中的闻频超声换能器完成，闻频超声换能器是本
现有技术不做详细描述。本专利技术的步骤I)与步骤4)属于本
公知技术，在本专利技术中不作详细描述。与现有技术相比，本专利技术的优点在于本专利技术所提出的利用交织编程原理进行超声成像测速的方法，克服流速较高条件下连续采样过程中粒子位移偏移量过高的限制，不增加硬件成本，大大提高所测流速的动态范围，使超声测速技术的应用范围进一步扩大，不仅可以用于医学上心血管内血流动力学参数的获取，而且还可用于一般尺度非透明流体成像和多相流测量领域，例如石油工业中油气的输送，化学药物溶解监测，海岸工程学研究和分析，江河入海口及海岸线的环境管理等。附图说明图1是传统顺序扫描方式采集第一帧图像示意图；图2是传统顺序扫描方式采集第二帧图像示意图；图3是交织编程成像法获取图像示意图；图4是计算速度向量场示意图； 图5是的流程图；图6是测速实验系统示意图；图7是按照传统顺序扫描方式进行超声成像测速获取的二维速度场；图8是采用本专利技术交织编程原理进行超声成像测速得到的二维速度场。具体实施例方式下面结合附图对本专利技术做进一步说明。，包括以下步骤I)向待测流体中添加示踪剂，示踪剂可以是造影微泡或其它声学性能较好的造影剂。2)通过计算机控制高频超声换能器按照交织编程的方式向外发射高频超声波信号，并采集流体中示踪剂的背向散射射频信号。传统顺序扫描方式与交织编程方式获取射频信号图像的区别之处在于假定超声成像测速系统中超声探头由η个单元构成，按线扫描的方式来记录每帧图像。由于超声造影剂微泡散射的信号是通过超声声束逐条扫描，再将扫描线组合在一起形成图像，从而造成不同扫描线之间时间不同。假定一幅图像由η条扫描线构成，定义图像中任意一条扫描线为Lab，其中a表示扫描线所在位置，取值范围为[1，n]，b表示采集该条扫描线时所处的时间单元，取值范围为[l，n]。按传统方法采集图像时，如图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。...

【技术特征摘要】
1.基于超声交织编程的速度分布测量方法，其特征在于，包括以下步骤 1)向待测流体中添加示踪剂； 2)控制器控制高频超声换能器按照交织编程方式向外发射高频超声波信号，并采集流体中示踪剂的背向散射射频信号； 设超声成像测速系统中超声探头由η个单元构成，按线扫描的方式来记录每帧图像；设一幅图像由η条扫描线构成，其中a表示扫描线所在位置，取值范围为[l，n]，b表示采集该条扫描线时所处的时间单元，取值范围为[l，n]; 交织编程方式采集图像时，具体方法如下在第I时间单元采集a=l位置处的射频信号扫描线，记为Lll ;在第2时间单元仍采集a=l位置处的射频信号扫描线，记为L12 ;在第3时间单元采集a=2位置处的射频信号扫描线，记为L23 ;在第4时间单元仍采集a=2位置处的射频信号扫描线，记为L24 ;以此类推，在第2n-l时间单元采集a=n位置处的射频信号扫描线，记为Ln(2n-1);在第2n时间单元仍采集a=n位置处的射频信号扫描线，记为Ln (2n)；交...

【专利技术属性】
技术研发人员：周宾，姚琛，许传龙，张龙，
申请(专利权)人：东南大学，
类型：发明
国别省市：