一种超声彩色血流成像的脉冲重复频率调整方法及装置技术领域
本发明实施例涉及超声信号处理技术领域,尤其涉及一种超声彩色血流成像的脉
冲重复频率调整方法及装置。
背景技术
超声成像因为其无创性、实时性、操作方便、价格便宜等诸多优势,使其成为临床
上应用最为广泛的诊断工具之一。超声彩色血流成像是超声成像中的一种重要成像模式,
在临床上的应用极其普遍,是多种疾病诊断的重要依据,因此提高超声彩色血流图像的质
量对于临床医生的诊断有重要意义。
超声彩色血流成像通常是在二维组织成像的基础上,由用户确定一个感兴趣区域
(ROI框),然后对该区域进行血流成像,再将血流图像叠加在二维组织图像上一起显示。
超声彩色血流成像的基本原理是超声多普勒效应,即发射一定频率的超声波(超
声波是一种频率高于20000赫兹的声波,它的方向性好,穿透能力强,易于获得较集中的声
能)进入人体组织后,在血流运动作用下接收的超声回波信号会产生频偏,频偏的大小与血
流速度成正比,用公式表示如下:
Δf=2*f*Vcos(θ)/C
其中f表示发射超声波信号的中心频率,Δf表示接收超声波信号的频偏,V表示血
流速度,θ表示声束与血流的夹角,C表示超声波传播速度。上述公式也可以表示成:
V=C*Δf/(2*f*cos(θ))
通常超声彩色血流成像有一定的速度范围,在发射超声波信号的中心频率固定的
情况下,该速度范围由脉冲重复频率(PRF)决定,脉冲重复频率实际上也是慢时间方向(即
扫描次数方向)的采样率。对于超声彩色血流成像而言,信号实际上有两个维度:一个维度
是扫描深度方向,也称为快时间方向;另一个维度是扫描次数方向,也称为慢时间方向。按
照奈奎斯特采样定理,接收超声波信号中能检测到的最大频偏是脉冲重复频率的一半(即
PRF/2),由于超声彩色血流成像通常不考虑或难以考虑角度的影响,因此超声彩色血流成
像能检测到的最大速度为:
Vmax=C*PRF/(4*f)
考虑血流运动的方向,超声彩色血流成像的速度范围就是-Vmax~Vmax。
实际计算时通常以相位来反应血流速度大小和方向,自相关算法是最早也是目前
最常用的超声彩色血流成像算法,其计算公式如下所示:
其中,I和Q分别是超声多普勒血流信号的同相分量和正交分量,脉冲重复频率实
际上就是I/Q信号在慢时间方向的采样率,为了去除组织低频杂波信号,I和Q必须先经过滤
波处理,PS是对同一扫描线位置按照脉冲重复频率重复扫描的次数,为保证一定的帧频和
计算精度,PS的取值范围通常在8~16之间。上述自相关运算得到的N、D可通过以下公式计
算相位和频偏:
其中,tan-1表示反正切运算,T是脉冲重复频率PRF的倒数。反正切运算得到的相位
范围始终是-π~π,对应的频偏范围是对应的速度范围是-Vmax~Vmax。当实际
血流速度超过这个范围时,会发生混叠现象,此时计算的血流速度和方向不能反映真实的
血流情况。
理想的超声彩色血流成像应该尽可能避免混叠现象,同时尽可能体现出血流的速
度梯度。如上所述,速度范围主要由脉冲重复频率控制,较大的脉冲重复频率可以确保不引
起混叠,但灵敏度降低,同时丢失速度梯度;较小的脉冲重复频率灵敏度高,但容易引起混
叠。因此,超声彩色血流成像需要设置合适的速度范围,以获得更好的血流成像表现力。由
于不同的人的各个应用部位之间血流速度有显著差异,因此脉冲重复频率需要随时调整,
但目前是用户根据应用部位和实际图像,依据经验手动调节脉冲重复频率实现控制速度范
围。
发明内容
本发明实施例提供了一种超声彩色血流成像的脉冲重复频率调整方法及装置,以
解决现有技术中PRF依靠使用经验进行调整,准确性差且效率低的技术问题。
第一方面,本发明实施例提供了一种超声彩色血流成像的脉冲重复频率调整方
法,该方法包括:
按照当前脉冲重复频率PRF进行扫查,根据扫查回波信号判断是否发生混叠;
在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心动周期调整后的
所有PRF;
选取所有PRF中的极值PRF作为最终优化的脉冲重复频率
进一步地,所述根据扫查回波信号判断是否发生混叠包括:
获取当前采样率对应的第一信号相位,所述当前采样率实为当前脉冲重复频率
PRF;
获取二分之一当前采样率对应的第二信号相位;
根据所述第一信号相位和第二信号相位判断是否发生混叠。
进一步地,所述根据所述第一信号相位和第二信号相位判断是否发生混叠,包括:
在第一信号相位和第二信号相位方向符号相同且第二信号相位不近似等于二倍
第一信号相位时,确定第二信号相位发生混叠;或者
在第一信号相位和第二信号相位方向符号相反时,确定第二信号相位发生混叠。
进一步地,所述根据混叠状态调整PRF,包括:
计算按照当前PRF扫查帧中第二信号相位中发生混叠的血流像素比例;
根据所述血流像素比例调整PRF。
进一步地,所述根据所述血流像素比例调整PRF,包括:
在所述血流像素比例低于比例阈值时,减小PRF;
在所述血流像素比例高于比例阈值时,增大PRF。
第二方面,本发明实施例还提供了一种超声彩色血流成像的脉冲重复频率调整装
置,该装置包括:
判断模块,用于按照当前脉冲重复频率PRF进行扫查,根据扫查回波信号判断是否
发生混叠;
调整模块,用于在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心
动周期调整后的所有PRF;
选取模块,用于选取所有PRF中的极值PRF作为最终优化的脉冲重复频率。
进一步地,所述判断模块包括:
第一相位获取单元,用于获取当前采样率对应的第一信号相位,所述当前采样率
实为当前脉冲重复频率PRF;
第二相位获取单元,用于获取二分之一当前采样率对应的第二信号相位;
判断单元,用于根据所述第一信号相位和第二信号相位判断是否发生混叠。
进一步地,所述判断单元具体用于:
在第一信号相位和第二信号相位方向符号相同且第二信号相位不近似等于二倍
第一信号相位时,确定第二信号相位发生混叠;或者
在第一信号相位和第二信号相位方向符号相反时,确定第二信号相位发生混叠。
进一步地,所述调整模块包括:
比例计算单元,用于计算按照当前PRF扫查帧中第二信号相位中发生混叠的血流
像素比例;
调整单元,用于根据所述血流像素比例调整PRF。
进一步地,所述调整单元具体用于:
在所述血流像素比例低于比例阈值时,减小PRF;
在所述血流像素比例高于比例阈值时,增大PRF。
本发明实施例通过按照当前脉冲重复频率PRF进行扫查,根据扫查回波信号判断
是否发生混叠;在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心动周期调
整后的所有PRF;选取所有PRF中的极值PRF作为最终优化的脉冲重复频率,解决了用户只能
根据应用部位和实际图像手动调节脉冲重复频率来控制血流速度范围的问题,使得超声彩
色血流成像功能可以根据当前扫查部位和扫查切面的实时血流图像对脉冲重复频率进行
自动调整和优化,使得当前状态下的血流图像尽可能避免混叠现象,又保持较好的血流速
度梯度。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本发明的其它
特征、目的和优点将会变得更明显:
图1是本发明实施例一提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图;
图2是本发明实施例二提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图;
图3是本发明实施例三提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图;
图4是本发明实施例四提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图;
图5是本发明实施例五提供的一种超声彩色血流成像的脉冲重复频率调整装置的
结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描
述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便
于描述,附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图,该方法可以由超声彩色血流成像的脉冲重复频率调整装置来执行,该装置可以由
软件和/或硬件的方式实现,该装置可以集成在具有超声彩色血流成像功能的电子设备中,
例如可以是彩超机。该方法具体包括如下步骤:
步骤110、按照当前脉冲重复频率PRF进行扫查,根据扫查回波信号判断是否发生
混叠。
其中,脉冲重复频率为上一次发射超声波信号与下一次发射超声波信号之间的时
间间隔的倒数。当前扫查部位可以包括多条扫描线位置。按照当前脉冲重复频率PRF进行扫
查,可以按照当前脉冲重复频率PRF向同一扫描线位置多次发射超声波信号。扫查回波信号
为接收的超声回波信号,每发射一次超声波信号将接收一次超声回波信号。示例性的,可以
根据扫查回波信号判断是否发生混叠,即将连续信号经过离散采样后,得到的离散信号的
傅氏谱为原信号傅氏谱SF倍的周期延拓,如果原信号中包含的最高频率成分fmax>SF/2,则
在离散信号谱中相应周期的谱会出现频率混叠。反之,如果fmax<SF/2,即采样频率大于分
析信号中最高频率成分的两倍,则采样后离散信号频谱中不会出现频率混叠。
步骤120、在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心动周期
调整后的所有PRF。
其中,当前扫查部位包括多条扫描线位置,对每条扫描线位置对应的扫查回波信
号进行分析判断是否发生混叠,在发生混叠时,根据混叠状态调整PRF。优选的可以是每生
成一帧超声彩色血流图像调整一次PRF,作为下一帧超声彩色血流图像的脉冲重复频率。该
完整心动周期指从一次心跳的起始到下一次心跳的起始,心血管系统所经历的过程,该心
动周期因人而异。由于每个心动周期内的血流速度一直在某一范围内变化。因此,需要根据
血流速度不断调整PRF。故在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心
动周期调整后的所有PRF。每生成一帧超声彩色血流图像所用时间小于一个完整心动周期,
在一个完整心动周期内将生成多帧超声彩色血流图像,对应的会多次调整PRF。在发生混叠
时,说明实际血流速度超过当前脉冲重复频率所对应的血流速度范围,需增大PRF。
步骤130、选取所有PRF中的极值PRF作为最终优化的脉冲重复频率。
其中,选取所有PRF中的极值PRF作为最终优化的脉冲重复频率,即无需再调整
PRF。该极值可以为极大值。由于心脏的周期性收缩和舒张,人体中的血流速度呈周期性变
化,要尽可能避免超声彩色血流成像发生混叠应该以血流速度最大的那一帧作为参考标
准,因此直至取得至少一个完整心动周期调整后的所有PRF,选取所有PRF中的极值PRF作为
最终优化的脉冲重复频率。
本实施例提供的技术方案,通过按照当前脉冲重复频率PRF进行扫查,根据扫查回
波信号判断是否发生混叠;在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整
心动周期调整后的所有PRF;选取所有PRF中的极值PRF作为最终优化的脉冲重复频率,解决
了用户只能根据应用部位和实际图像手动调节脉冲重复频率来控制血流速度范围的问题,
使得超声彩色血流成像功能可以根据当前扫查部位和扫查切面的实时血流图像对脉冲重
复频率进行自动调整和优化,使得当前状态下的血流图像尽可能避免混叠现象,又保持较
好的血流速度梯度。
实施例二
图2是本发明实施例二提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图,本实施例以上述实施例为基础进行优化,提供了一种混叠判断的具体方法,具体是
所述根据扫查回波信号判断是否发生混叠包括:获取当前采样率对应的第一信号相位,所
述当前采样率实为当前脉冲重复频率PRF;获取二分之一当前采样率对应的第二信号相位;
根据所述第一信号相位和第二信号相位判断是否发生混叠。相应的,本实施例的方法包括:
步骤210、按照当前脉冲重复频率PRF进行扫查。
步骤220、获取当前采样率对应的第一信号相位,所述当前采样率实为当前脉冲重
复频率PRF。
其中,自相关算法是目前最常用的超声彩色血流成像算法,根据超声多普勒血流
信号的同相分量和正交分量,来获取第一信号相位,用相位可以反应血流速度大小和方向。
脉冲重复频率就是超声多普勒血流信号的同相分量和正交分量信号在慢时间方向的采样
率,即当前采样率实为当前脉冲重复频率PRF。当前采样率与当前脉冲重复频率PRF相等,随
着PRF的调整,采样率会随着PRF变化,始终与PRF相等。可以根据按照当前脉冲重复频率PRF
向当前扫查部位发射预定次数的超声波信号获得的超声回波信号,获得当前采样率对应的
第一信号相位。
步骤230、获取二分之一当前采样率对应的第二信号相位。
按照奈奎斯特采样定理,当采样频率大于信号中最高频率的2倍时,采样之后的数
字信号完整地保留了原始信号中的信息。因此,需要将采样率降低至原有采样频率的二分
之一,以方便确定是否会发生混叠。其中,获取二分之一当前采样率对应的第二信号相位,
可以是从按照当前脉冲重复频率PRF向当前扫查部位发射预定次数的超声波信号获得的超
声回波信号中,按照二分之一当前采样率在慢时间方向进行抽取获得的超声回波信号,来
求取对应的第二信号相位,还可以是按照二分之一当前脉冲重复频率PRF向当前扫查部位
发射预定次数的超声波信号获得的超声回波信号,来求取对应的第二信号相位。
步骤240、根据所述第一信号相位和第二信号相位判断是否发生混叠。
其中,由于根据频偏与脉冲重复频率和超声多普勒血流信号的同相分量和正交分
量的相位的关系,以及频偏与血流速度的对应关系,可知若血流速度为定值,则频偏也为定
值,故脉冲重复频率和超声多普勒血流信号的同相分量和正交分量的相位的乘积为定值。
若当前采样率(即当前脉冲重复频率)和对应的超声多普勒血流信号的同相分量和正交分
量的第一信号相位的乘积,与二分之一当前采样率和对应的超声多普勒血流信号的同相分
量和正交分量的第二信号相位的乘积相等,则实际血流速度未超过当前脉冲重复频率所对
应的血流速度范围,故未发生混叠;否则,实际血流速度超过当前脉冲重复频率所对应的血
流速度范围,故发生混叠。可以继续重复获取2倍降采样率(即四分之一当前采样率、八分之
一当前采样率等)的超声回波信号,来计算相应的信号相位来判断是否发生混叠,来提高混
叠判断的准确率。
步骤250、在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心动周期
调整后的所有PRF。
步骤260、选取所有PRF中的极值PRF作为最终优化的脉冲重复频率。
本实施例的技术方案提供了一种超声彩色血流成像的脉冲重复频率调整方法,该
实施例在上述实施例的基础上,通过获取当前采样率对应的第一信号相位,所述当前采样
率实为当前脉冲重复频率PRF;获取二分之一当前采样率对应的第二信号相位;根据所述第
一信号相位和第二信号相位判断是否发生混叠,使得超声彩色血流成像功能可以根据当前
扫查部位和扫查切面的实时血流图像对脉冲重复频率进行自动调整和优化,使得当前状态
下的血流图像尽可能避免混叠现象,又保持较好的血流速度梯度。
实施例三
图3是本发明实施例三提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图,本实施例以上述实施例为基础进行优化,具体是所述根据所述第一信号相位和第
二信号相位判断是否发生混叠,包括:在第一信号相位和第二信号相位方向符号相同且第
二信号相位不近似等于二倍第一信号相位时,确定第二信号相位发生混叠;或者在第一信
号相位和第二信号相位方向符号相反时,确定第二信号相位发生混叠。相应的,本实施例的
方法包括:
步骤310、按照当前脉冲重复频率PRF进行扫查。
步骤320、获取当前采样率对应的第一信号相位,所述当前采样率实为当前脉冲重
复频率PRF。
步骤330、获取二分之一当前采样率对应的第二信号相位。
步骤340、在第一信号相位和第二信号相位方向符号相同且第二信号相位不近似
等于二倍第一信号相位时,确定第二信号相位发生混叠;或者在第一信号相位和第二信号
相位方向符号相反时,确定第二信号相位发生混叠。
其中,在第一信号相位和第二信号相位方向符号相同且第二信号相位不近似等于
二倍第一信号相位,即第二信号相位与二倍第一信号相位的差值不满足预设阈值范围时,
确定第二信号相位发生混叠;或者在第一信号相位和第二信号相位方向符号相反时,确定
第二信号相位发生混叠。而在第一信号相位和第二信号相位方向符号相同且第二信号相位
近似等于二倍第一信号相位(即第二信号相位与二倍第一信号相位的差值满足预设阈值范
围)时,确定第二信号相位未发生混叠。
步骤350、在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整心动周期
调整后的所有PRF。
步骤360、选取所有PRF中的极值PRF作为最终优化的脉冲重复频率。
本实施例的技术方案提供了一种超声彩色血流成像的脉冲重复频率调整方法,该
实施例在上述实施例的基础上,通过在第一信号相位和第二信号相位方向符号相同且第二
信号相位不近似等于二倍第一信号相位时,确定第二信号相位发生混叠;或者在第一信号
相位和第二信号相位方向符号相反时,确定第二信号相位发生混叠,使得超声彩色血流成
像功能可以根据当前扫查部位和扫查切面的实时血流图像对脉冲重复频率进行自动调整
和优化,使得当前状态下的血流图像尽可能避免混叠现象,又保持较好的血流速度梯度。
实施例四
图4是本发明实施例四提供的一种超声彩色血流成像的脉冲重复频率调整方法的
流程图,本实施例以上述实施例为基础进行优化,提供了一种调整PRF的具体方法,具体是
所述根据混叠状态调整PRF,包括:计算按照当前PRF扫查帧中第二信号相位中发生混叠的
血流像素比例;根据所述血流像素比例调整PRF。相应的,本实施例的方法包括:
步骤410、按照当前脉冲重复频率PRF进行扫查。
步骤420、获取当前采样率对应的第一信号相位,所述当前采样率实为当前脉冲重
复频率PRF。
步骤430、获取二分之一当前采样率对应的第二信号相位。
步骤440、根据所述第一信号相位和第二信号相位判断是否发生混叠。
步骤450、在发生混叠时,计算按照当前PRF扫查帧中第二信号相位中发生混叠的
血流像素比例。
其中,扫查帧包括多个血流像素点,每个血流像素点的像素值有与其对应第一信
号相位和第二信号相位。在发生混叠时,计算按照当前PRF扫查帧中第二信号相位中发生混
叠的血流像素比例。示例性的,可以采用发生混叠的血流像素点个数与总的血流像素点个
数的商来表示。
步骤460、根据所述血流像素比例调整PRF,直至取得至少一个完整心动周期调整
后的所有PRF。
其中,根据血流像素比例调整PRF,直至取得至少一个完整心动周期调整后的所有
PRF。
优选的,所述根据所述血流像素比例调整PRF,包括:在所述血流像素比例低于比
例阈值时,减小PRF;在所述血流像素比例高于比例阈值时,增大PRF。
其中,在血流像素比例低于比例阈值时,减小PRF;在血流像素比例高于比例阈值
时,增大PRF。如果没有发生混叠或发生混叠的比例较小,说明当前脉冲重复频率过高,应适
当减小脉冲重复频率;如果发生混叠的比例较高,说明当前脉冲重复频率过低,应适当提高
脉冲重复频率。PRF可以根据当前脉冲重复频率、比例阈值和血流像素比例来确定,对应关
系可以是PRF为血流像素比例和比例阈值的差值加1的和,与当前脉冲重复频率的乘积。
示例性的,可以采用如下公式计算调整PRF:
PRF_Opt=PRF*(1+(alias_ratio–threshold));其中,PRF为当前帧脉冲重复频
率,PRF_Opt即为优化后的脉冲重复频率,threshold为预设的阈值,alias_ratio为发生混
叠的血流像素比例。
步骤470、选取所有PRF中的极值PRF作为最终优化的脉冲重复频率。
本实施例的技术方案提供了一种超声彩色血流成像的脉冲重复频率调整方法,该
实施例在上述实施例的基础上,通过计算按照当前PRF扫查帧中第二信号相位中发生混叠
的血流像素比例;根据血流像素比例调整PRF,使得超声彩色血流成像功能可以根据当前扫
查部位和扫查切面的实时血流图像对脉冲重复频率进行自动调整和优化,使得当前状态下
的血流图像尽可能避免混叠现象,又保持较好的血流速度梯度。
实施例五
图5是本发明实施例五提供的一种超声彩色血流成像的脉冲重复频率调整装置的
结构示意图,该装置适用于执行本发明实施例提供的超声彩色血流成像的脉冲重复频率调
整方法,如图5所示,该装置包括判断模块510、调整模块520和选取模块530。
其中,判断模块510,用于按照当前脉冲重复频率PRF进行扫查,根据扫查回波信号
判断是否发生混叠;调整模块520,用于在发生混叠时,根据混叠状态调整PRF,直至取得至
少一个完整心动周期调整后的所有PRF;选取模块530,用于选取所有PRF中的极值PRF作为
最终优化的脉冲重复频率。
本实施例提供的技术方案,通过按照当前脉冲重复频率PRF进行扫查,根据扫查回
波信号判断是否发生混叠;在发生混叠时,根据混叠状态调整PRF,直至取得至少一个完整
心动周期调整后的所有PRF;选取所有PRF中的极值PRF作为最终优化的脉冲重复频率,解决
了用户只能根据应用部位和实际图像手动调节脉冲重复频率来控制血流速度范围的问题,
使得超声彩色血流成像功能可以根据当前扫查部位和扫查切面的实时血流图像对脉冲重
复频率进行自动调整和优化,使得当前状态下的血流图像尽可能避免混叠现象,又保持较
好的血流速度梯度。
优选的,所述判断模块510包括:第一相位获取单元511,用于获取当前采样率对应
的第一信号相位,所述当前采样率实为当前脉冲重复频率PRF;第二相位获取单元512,用于
获取二分之一当前采样率对应的第二信号相位;判断单元513,用于根据所述第一信号相位
和第二信号相位判断是否发生混叠。
优选的,所述判断单元513具体用于:在第一信号相位和第二信号相位方向符号相
同且第二信号相位不近似等于二倍第一信号相位时,确定第二信号相位发生混叠;或者,在
第一信号相位和第二信号相位方向符号相反时,确定第二信号相位发生混叠。
优选的,所述调整模块520包括:比例计算单元521,用于计算按照当前PRF扫查帧
中第二信号相位中发生混叠的血流像素比例;调整单元522,用于根据所述血流像素比例调
整PRF。
优选的,所述调整单元522具体用于:在所述血流像素比例低于比例阈值时,减小
PRF;在所述血流像素比例高于比例阈值时,增大PRF。
上述超声彩色血流成像的脉冲重复频率调整装置可执行本发明任意实施例所提
供的超声彩色血流成像的脉冲重复频率调整方法,具备执行方法相应的功能模块和有益效
果。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,
本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、
重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行
了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还
可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。