具体实施方式
为了在成像帧率变化很大时,也能够保证图像的帧相关处理结果非常接近,本发明给出了一种帧相关系数的控制方法,其精确计算了当前帧率下的帧相关系数。如图2所示,本发明所给出的帧相关系数控制方法的主要过程如下:
首先,比较当前成像帧率Frt下的当前滤波器幅频响应|H(f)|与一基准成像帧率Fr0下的基准滤波器幅频响应|H0(f)|,获得当前滤波器幅频响应|H(f)|与基准滤波器幅频响应|H0(f)|在重叠频带内的差异E,即E=F(|H(f)|,|H0(f)|),F表示差异函数;
然后,将所述差异E为最小时的当前滤波器幅频响应|H(f)|所对应的帧相关系数,判定为当前成像帧率Frt所对应的最佳帧相关系数。
可见,本发明在现有技术的基础上,给出了一种新的帧相关系数计算方法,该方法定义了一个最小差异原则,其计算获得的帧相关系数使得当前滤波器幅频响应|H(f)|与基准滤波器幅频响应|H0(f)|的差异最小,也就是说此时在重叠频带内,滤波器在当前成像帧率Frt和基准成像帧率Fr0下具有相近的幅频响应特性。只有在图像处理时保证不同成像帧率下的滤波器幅频响应都接近于基准幅频响应,才能保证不同成像帧率下的滤波器幅频响应之间在重叠频带内具有相近的幅频响应特性,从而有效保证了图像处理的平滑处理效果以及噪声的抑制作用。
如图2所示,本发明的帧相关系数控制方法需要设定一个基准的频率Fr0(如实时成像帧率,即50HZ或60H)下的帧相关系数a(Fr0)。此帧相关系数a(Fr0)可以根据前述公式(3)、通过选择帧相关控制系数k确定来计算,帧相关控制系数k对应超声成像***的帧相关档位,其表征图像的帧相关效果。以帧相关系数a(Fr0)对应的幅频响应|H0(f)|作为参考函数,当当前成像帧率Frt大于Fr0的时候,平滑处理要求当前成像帧率Frt下帧相关系数a(Frt)对应的滤波器幅频响应|H(f)|在0~Fr0/2的范围内具有相近的响应;而,当当前成像帧率Frt小于Fr0的时候,要求当前成像帧率Frt下帧相关系数a(Frt)对应的滤波器幅频响应|H(f)|在0~Frt/2范围内具有相近的响应。可见,两个滤波器幅频响应的重叠频带可以表示为[0,min(Fr0/2,Frt/2)],通过设定当前成像帧率Frt下的帧相关系数a(Frt)使得对应的滤波器幅频响应在频带[0,min(Fr0/2,Frt/2)]内相近,实现图像的平滑处理,从而保证在成像帧率变化较大的情况下,也能够获得基本一致的帧相关处理效果。
图2所示的帧相关系数的控制方法,按照差异函数E=F(|H(f)|,|H0(f)|)的建立方式,可以存在多种实现方式,以下针对常用的几种方式进行详细说明。
实施例1:如图3所示,本实施例中所述差异E为:基准滤波器幅频响应|H0(f)|与当前滤波器幅频响应|H(f)|之差在重叠频带[0,min(Fr0/2,Frt/2)]内的平方和,即
如图3所示,本实施例的具体过程如下:
步骤200,预先设定的基准成像帧率Fr0,这里的基准成像帧率Fr0为实时成像帧率,如50HZ或60HZ。
步骤210,将预先设定的基准成像帧率Fr0,代入前述公式(3)中计算其帧相关系数a(Fr0),结果如下所示:
a(Fr0)=exp{-1/k×Fr0}
其中,k为帧相关控制系数,其由经验来确定,此k值决定了图像的帧相关效果,并与超声探测仪器的帧相关档位相对应,其取值范围为0~1。
步骤220,将上述计算获得帧相关系数a(Fr0)输入到数字滤波器的传递函数模型中,计算滤波器的基准幅频响应|H0(f)|,其对应的频带范围为0~Fr0/2。
在这里,以数字滤波器的传递函数模型属于一阶无限脉冲响应(InfiniteImpulse Response)数字滤波器Z域的传递函数模型为例,详细说明基准幅频响应|H0(f)|的计算过程。
根据前述公式(2)
将z=exp{j2πf'},f'=f/Fr,f∈[0,Fr/2],代入上述H(z)中即获得滤波器的频率响应H(f),其幅度对应的幅频响应|H(f)|表达式为:
其中,f'为归一化数字频率,f为实际频率,Fr为成像帧率。
根据上述公式,将帧相关系数a(Fr0)、基准成像帧率Fr0代入到上述公式(5)中,即可获得关于基准成像帧率Fr0下的基准幅频响应|H0(f)|的函数表达式。
步骤230,设定当前成像帧率Frt下的N个候选帧相关系数a`(Frt),其中N为大于2的自然数;然后再分别将此N个候选帧相关系数a`(Frt)输入到数字滤波器的传递函数模型(如上述公式(5)所示的一阶无限脉冲响应数字滤波器Z域的传递函数模型)中,计算各个候选帧相关系数a`(Frt)对应的滤波器当前幅频响应|H(f)|,其对应的频带范围为0~Frt/2。
这里的N个候选帧相关系数a`(Frt)可以根据经验来确定,也可以按照以下方式来确定:
第一步,按照常规方法计算当前成像帧率Fr
t下的帧相关系数估计值
例如,将当前成像帧率Fr
t代入前述公式(3)所示的帧相关系数与成像帧率的指数函数中,计算其帧相关系数估计值
结果如下所示:
从上述公式可以看出,不同的帧相关控制系数k对应不同帧相关效果的候选帧相关系数a`(Fr
t),在计算该估计值
时,其k值可以参照步骤210中所选的值。
第二步,在所述估计值
的附近选择N个值作为候选帧相关系数a`(Fr
t),该N为大于等于2的自然数。例如,在以
值为中心的区间内取100个数值,作为候选帧相关系数a`(Fr
t)。
步骤240,在基准幅频响应|H0(f)|和当前幅频响应|H(f)|的重叠频带内,分别比较所述基准幅频响应|H0(f)|与各个候选帧相关系数a`(Frt)所对应的滤波器当前幅频响应|H(f)|之间的差异E,即按照上述公式(4)计算在[0,min(Fr0/2,Frt/2)]内各个a`(Frt)所对应的E值,根据差异最小原则,从各个候选帧相关系数a`(Frt)中选择一个最优的作为当前成像帧率Frt下的最佳帧相关系数a(Frt),即将E值最小时的|H(f)|所对应的候选帧相关系数a`(Frt)作为当前成像帧率Frt下的最佳帧相关系数a(Frt)。
本实施例滤波器系数是根据当前滤波器幅频响应|H(f)|在频带[0,min(Fr0/2,Frt/2)]内与a(Fr0)对应的滤波器幅频响应|H0(f)|具有最小的均方误差原则来确定的。可见,幅频响应函数为帧相关系数a和频率f的非线性函数,当前成像帧率Frt下对应的帧相关系数a(Frt)可以通过数值法求解。利用本实施例的实现方式获得的不同成像帧率下的帧相关系数对应的幅频响应如图4所示。由图4可见,在重叠频带内、不同帧率下的幅频响应也非常相近,确保了不同帧率情况下组织运动的模糊程度相近,且对噪声也有较好地抑制效果。
当然上述差异的计算方式并不只一种,请见下述实施例。
实施例2:本实施例与实施例1的不同点在于,在步骤240中,本实施例中所述差异E为:基准滤波器幅频响应|H0(f)|与当前滤波器幅频响应|H(f)|之差在重叠频带[0,min(Fr0/2,Frt/2)]内的绝对值之和,即
通过比较该差异值E的大小,将该差异值E为最小时的当前滤波器幅频响应|H(f)|所对应的候选帧相关系数a`(Frt)确定为当前成像帧率Frt下的最佳帧相关系数a(Frt)。本实施例的其余步骤与实施例1相同。本实施例滤波器系数是根据当前滤波器幅频响应|H(f)|在频带[0,min(Fr0/2,Frt/2)]内与a(Fr0)对应的滤波器幅频响应|H0(f)|具有最小的误差绝对值之和的原则来确定的。同实施例1,当前帧率Frt下对应的帧相关系数a(Frt)可以通过数值法来求解。利用本实施例的实现方式获得的不同成像帧率下的帧相关系数对应的幅频响应与图4类似。
实施例3:本实施例与实施例1的不同点在于,在步骤240中,本实施例中所述差异E为:在重叠频带[0,min(Fr0/2,Frt/2)]内,在一设定频率f1下基准滤波器幅频响应|H0(f)|与当前滤波器幅频响应|H(f)|之差的绝对值,即
E=‖H(f1)|-|H0(f1)‖, (7)
其中,设定频率f1为上述重叠频带(0,min(Fr0/2,Frt/2)]内除零外的任一值。通过设定一个频率f1,然后根据上述公式(7)计算并比较差异值E的大小,将该差异值E为最小时的当前滤波器幅频响应|H(f)|所对应的候选帧相关系数a`(Frt)确定为当前成像帧率Frt下的最佳帧相关系数a(Frt)。本实施例的其余步骤与实施例1相同。本实施例滤波器系数是根据在频带[0,min(Fr0/2,Frt/2)]内当前滤波器幅频响应|H(f)|与a(Fr0)对应的滤波器幅频响应|H0(f)|之差的绝对值最小原则来确定的。同实施例1,当前帧率Frt下对应的帧相关系数a(Frt)可以通过数值法来求解。
实施例4:本实施例与实施例3的不同点在于,在步骤240中,本实施例中所述差异为最小的情况是指:在一设定频率f1下基准滤波器幅频响应|H0(f)|与当前滤波器幅频响应|H(f)|之差为零,即
E=|H(f1)|-|H0(f1)|=0, (8)
也就是,在设定频率f1上使|H0(f)|=|H(f)|,此时当前滤波器幅频响应|H(f)|所对应的候选帧相关系数a`(Frt)被判定为当前成像帧率Frt下的最佳帧相关系数a(Frt),其中f1=(0,min(Fr0/2,Frt/2)]。本实施例的其余步骤与实施例3相同,该其余步骤是指除步骤230设定候选帧相关系数、及与候选帧相关系数相关的步骤外。
本实例是通过设定a(Frt)使得当前滤波器幅频响应|H(f)|在某预先设定的频率f1(f1不等于0)上具有与基准滤波器幅频响应|H0(f)|相同的数值,即令|H0(f)|=|H(f)|,其中,频率f1最优可选择为min(Fr0/2,Frt/2),即Fr0/2和Frt/2两者中的较小者。并,通过将a(Fr0)、f1代入幅频响应函数|H0(f)|,并将f1代入幅频响应函数|H(f)|,通过解析求解上述等式(8)即获得当前成像帧率Frt下的最佳帧相关系数a(Frt)。
本实施例属于实施例3中的一种特殊情况,采用本实施例的解析求解后,可以简化上述如图3所示的计算规则,如可以减少设定候选帧相关系数的步骤等等,因为上述公式(8)是一个确定的等式,其可以直接计算出最佳帧相关系数a(Frt)。本实施例相对于实施例1至3而言,本实施例的实现方式由于计算量较小,可以由成像***软件根据当前成像帧率Frt和帧相关档位对应的帧相关控制系数k实时计算获得当前成像帧率Frt的帧相关系数a(Frt)。利用本实施例的实现方式获得的不同成像帧率下的帧相关系数对应的幅频响应如图5所示。由图5可见,在重叠频带内幅频响应也非常相近,确保了不同帧率情况下组织运动的模糊程度相近,且对噪声也有较好的抑制效果。
实施例5:为了简化计算步骤,本实施例在上述实施例1、2、3或4的基础上,还增加以下的建表步骤:
通过实施例1、2、3或4确定了当前成像帧率Frt下的最佳帧相关系数a(Frt)后,将此最佳帧相关系数a(Frt)保存于一关系表中,该表记录当前成像帧率Frt、当前成像帧率Frt下的最佳帧相关系数a(Frt)、及表征帧相关效果的帧相关控制系数k三者之间的对应关系;
然后,在实时图像处理时,可以根据当前成像帧率Frt及帧相关控制系数k,查找所述关系表,实时确定当前成像帧率Frt下的最佳帧相关系数a(Frt)。利用本实施例建立的关于最佳帧相关系数的关系表可以存储于一计算机可读存储介质中。
由于上述实施例1至3的数值法求解运算量较大,因此本实施例通过预先计算好不同成像帧率下的最佳帧相关系数,将这些帧相关系数存储为一个系数表。当需要设置不同的帧相关效果的时候,可以通过更改基准帧率Fr0下的帧相关控制系数k,获得不同的|H0(f)|,然后数值计算出相应的不同成像帧率Frt下的帧相关系数a,从而组成了不同帧相关效果和不同成像帧率的一个二维系数表,然后将利用实施例1、2或3求解出的当前成像帧率下的最佳帧相关系数a(Frt)加入其中,即可获得上述关系表。超声成像***根据不同帧相关效果的档位和成像帧率对上述关系表进行查表,即可获得当前成像帧率下的帧相关系数。
上述分析只是列举了常见的几种差异的计算方式,本发明对差异如何计算并无限制。比如,根据上述各个实施例的差异数学模型可以概括的说,差异E的数学模型可以表示为基准滤波器幅频响应|H0(f)|与当前滤波器幅频响应|H(f)|之差的绝对值在重叠频带[0,min(Fr0/2,Frt/2)]内的n次方之和,即
其中,n为大于零的自然数。当上述公式(9)中n=1时,即上述实施例2所示的情况;当n=2时,即上述实施例1所示的情况;当n=3时也是可以采用的,只是计算相对复杂一些。另外,根据上述各个差异E的数学模型以及公式(9),其实还存在一种等效的情况,即差异E的数学模型可以表示为基准滤波器幅频响应|H0(f)|与当前滤波器幅频响应|H(f)|之差在重叠频带[0,min(Fr0/2,Frt/2)]内的2n次方之和,即
其中,n为大于零的自然数。当上述公式(10)中n=1时,即上述实施例1所示的情况;当n=2时,同样也是可取的计算方式。
由上可见,本发明提出的帧相关系数的控制方法不局限于组织灰度成像处理,还可以适用于彩色血流成像、组织多普勒成像等各类成像模式,其能精确的确定最佳帧相关系数用以帧相关处理中,保证图像平滑处理的效果。
通常有效的图像信息在相邻帧之间变化不大,而无效的噪声信息则通常变化较大,利用帧相关所采用的低通滤波器特性,可以较好保持帧间变化不大的低频有效图像信息,并抑制帧间变化较大的高频无效噪声。采用现有技术获得的图1所示的帧相关滤波器幅频响应曲线中,可以明显地看出各个曲线间相差很大,存在较多的噪声信息;而采用本发明的方法所获得的图4和图5所示的帧相关滤波器幅频响应曲线中,各个相邻曲线之间基本具有一致的帧相关效果,并且各个曲线在相互重叠的频带内具有相似的幅频特性。可见,本发明通过对不同的成像帧率,优化选择相应的帧相关系数从而获得不同成像帧率下的基本一致的帧相关效果,有效地保证了图像处理的平滑处理效果以及噪声的抑制作用。
依据上述帧相关系数的控制方法,本发明还相应提供了一种帧相关系数控制器,其能够直接获得当前成像帧率下的最佳帧相关系数,保证不同成像帧率下的滤波器幅频响应在重叠频带内具有相近的频率响应特性。如图6所示,所述控制器包括:
第一单元300,用于根据预设的基准成像帧率Fr0计算基准成像帧率Fr0下的帧相关系数a(Fr0);
第二单元310,用于根据第一单元300获得的帧相关系数a(Fr0)计算基准滤波器幅频响应|H0(f)|;及
第三单元320,用于比较当前成像帧率Frt下的当前滤波器幅频响应|H(f)|与所述基准滤波器幅频响应|H0(f)|在重叠频带内的差异E,输出所述差异E为最小时的当前滤波器幅频响应|H(f)|所对应的帧相关系数a(Frt)。
这里的参数a(Fr0)、Fr0、|H0(f)|、Frt、E、|H(f)|、a(Frt)、重叠频带的计算或设定原则,可以参见上述帧相关系数的控制方法的实施过程。
如图7所示,为了计算差异E,可以在控制器的第三单元320中建立以下单元:
差异模型单元321,用于存储所述当前滤波器幅频响应|H(f)|与所述基准滤波器幅频响应|H0(f)|的差异数学模型E=F(|H(f)|,|H0(f)|),并利用该数学模式计算所述当前滤波器幅频响应|H(f)|与所述基准滤波器幅频响应|H0(f)|的差异E;及
比较单元322,用于比较所述差异E,输出所述差异E为最小时的当前滤波器幅频响应|H(f)|所对应的帧相关系数。采用此结构,并利用上述实施例4的话,即可以实现实时计算当前成像帧率下的最佳帧相关系数;但是如果采用上述实施例1、2或3时,涉及到需要从多个候选帧相关系数中进行选择时,所述第三单元320就需要在差异模型单元321和比较单元322的基础上再增加以下几个单元,如图7所示,所述第三单元320还包括:
候选组计算模块323,用于生成当前成像帧率Frt下的候选帧相关系数a`(Frt);
滤波器幅频响应计算模块324,用于将所述候选帧相关系数a`(Frt)输入到数字滤波器的传递函数模型325中,计算候选帧相关系数a`(Frt)对应的当前滤波器幅频响应|H(f)|,将此当前滤波器幅频响应|H(f)|输入到所述差异模型单元321中进行差异计算,并通过比较单元322比较后获得相应的帧相关系数。
基于数据结果需要保存、以及建立快捷用关系表的需求,本发明的帧相关系数控制器可以在图6所示结构的基础上,增加一存储单元330,其用于记录当前成像帧率Frt下的最佳帧相关系数a(Frt),该最佳帧相关系数a(Frt)为所述差异E为最小时的当前滤波器幅频响应|H(f)|所对应的帧相关系数,即从选帧相关系数a`(Frt)中按照差异最小原则选出的帧相关系数。这里的存储单元330可以设置一关系表,并将记录的最佳帧相关系数a(Frt)添加其中,利用该关系表记录当前成像帧率Frt、当前成像帧率Frt下的最佳帧相关系数a(Frt)、及表征帧相关效果的帧相关控制系数k三者之间的对应关系。
从上述帧相关系数控制器的结构描述可以看出,其主要是利用上述帧相关系数的控制方法而建立一种产品,该产品采用上述帧相关系数的控制方法能够获得当前帧率下的最佳帧相关系数,利用该帧相关系数进行超声图像处理,可以保证在重叠频带内不同成像帧率下的滤波器幅频响应具有相近的频率响应特性,并且有效保证了图像处理的平滑效果及噪声抑制效果。
基于上述帧相关系数的控制方法在帧相关处理中的应用,本发明还提供了一种超声图像帧相关处理方法,该方法包括以下步骤:
步骤S1,采用上述帧相关系数的控制方法获得当前成像帧率Frt所对应的最佳帧相关系数a(Frt);
步骤S2,将获得的最佳帧相关系数a(Frt)输入到帧相关处理的数学模型中,进行帧相关运算,获得处理后的当前帧图像。这里的帧相关处理数学模型可以采用前述公式(1),其具体实施例可参见图3。由图3可见,在上述实施例1、2或3计算出最佳帧相关系数a(Frt)后,只要利用该帧相关系数进行图像的帧相关处理即可。
本发明所提出的超声图像帧相关处理方法中,关于最佳帧相关系数a(Frt)的计算过程以及其拓展的技术方案参见上述关于帧相关系数控制方法的说明,在此不再重复说明。
在上述帧相关系数控制方法及其帧相关处理方法的基础上,本发明还提供了一种新的超声成像***,如图8所示,所述超声成像***包括:超声图像处理模块和用于获取超声图像的探头;所述超声图像处理模块中设置有一帧相关系数控制器400和一帧相关处理单元410,这里的帧相关系数控制器400采用如图6所示的结构,具体可参见前述说明;所述帧相关处理单元410用于接收帧相关系数控制器400的第三单元320输出的帧相关系数,并将该帧相关系数输入到帧相关处理的数学模型中进行帧相关运算,获得平滑处理后的当前帧图像。这一超声成像***中,关于帧相关系数控制器400的内容以及其拓展的技术方案参见上述关于图6所示帧相关系数控制器的说明,在此不再重复说明。
基于上述最佳帧相关系数可以采用查表方式进行的实现方式,本发明还提供了一种新的超声成像***,如图9所示,该***中包括:超声图像处理模块和用于获取超声图像的探头;并且,在超声图像处理模块中设置有一存储介质430、查找单元420和一图6所示的帧相关处理单元410,所述存储介质430保存有一关系表,该表记录有当前成像帧率、当前成像帧率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数三者之间的对应关系,所述最佳帧相关系数对应的当前滤波器幅频响应与一基准成像帧率下的基准滤波器幅频响应在重叠频带内的差异最小;所述帧相关处理单元410通过所述查找单元420与所述存储介质430相连,所述查找单元420用于查找所述关系表中与当前成像帧率及当前帧相关控制系数相对应的最佳帧相关系数,并将该最佳帧相关系数送入到所述帧相关处理单元410中;所述帧相关处理单元410用于将所述最佳帧相关系数输入到帧相关处理的数学模型中进行帧相关运算,获得处理后的当前帧图像。
图9所示的超声成像***在于:采用上述帧相关系数控制方法所获得的最佳帧相关系数来建立一关系表,用于在图像处理时实时查表获得相应的帧相关系数,这样做可以减少计算量,提高图像处理的速度。建表的关键问题在于如何计算最佳帧相关系数,其可以参见上述关于帧相关系数控制方法的具体说明。比如:如图采用上述实施例1、2或3的方式,最佳帧相关系数可以采用如下方法获得:
首先,设定当前成像帧率下的多个候选帧相关系数,该候选帧相关系数的选择方法如下:
按照帧相关系数与成像帧率的指数函数,计算当前成像帧率下的帧相关系数估计值;在所述估计值的附近选择N个值作为候选帧相关系数,该N为大于等于2的自然数;
然后,分别将所述N个候选帧相关系数输入到数字滤波器的传递函数模型中,计算各个候选帧相关系数对应的当前滤波器幅频响应;
最后,分别将所述N个候选帧相关系数对应的当前滤波器幅频响应与所述基准滤波器幅频响应(其计算过程可以参见上述实施例1的相关说明)进行比较,根据当前滤波器幅频响应与所述基准滤波器幅频响应在重叠频带内的差异为最小的原则,从所述N个候选帧相关系数中选出一个作为当前成像帧率所对应的最佳帧相关系数,即所述差异为最小时的当前滤波器幅频响应所对应的候选帧相关系数。但如果采用实施4或5的话,则可以参见上述具体说明。
综上所述,本发明通过设定不同帧率下的帧相关系数,使得对应的滤波器幅频响应在重叠的频带内具有相近的频率响应特性的要求,提供了一种帧相关系数的计算方法,基于此计算方法,本发明又衍生出了一种可以制成独立产品并用于获得最佳帧相关系数的帧相关控制器、一种利用上述帧相关系数计算方法的帧相关处理方法及利用该方法的超声成像***、以及一种采用含有最佳帧相关系数关系表的存储介质的超声成像***。对于最后一个超声成像***,其存储介质中含有的最佳帧相关系数关系表采用上述帧相关系数计算方法获得。本发明所公开的各种方法和***,均可以根据最小差异原则计算获得当前成像帧率下的最佳帧相关系数,并用于后续的图像处理中,使得本发明在成像帧率发生较大改变时,超声图像帧相关的处理效果基本保持不变。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。