帧相关系数控制及帧相关处理方法、装置及超声成像系统的制作方法

专利名称：帧相关系数控制及帧相关处理方法、装置及超声成像系统的制作方法
技术领域：
本发明涉及超声系统中超声图像时域平滑滤波系数的优化技术领域， 尤其涉及的是， 一种帧相关系数的控制方法及其装置、以及利用该方法的 帧相关处理方法及超声成像系统。
背景纟支术
超声图像由于存在较大的电路噪声和斑点(Speckle)噪声(其中，斑 点噪声是由于大量尺寸小于超声波波长的粒子的散射信号相干叠加而形成 的，这是超声显像系统中一个固有的问题)，所以需要经过一些时域和空间 域的平滑处理，以改善图像的信噪比。其中，时域的平滑处理通常也称帧 相关处理。常用的帧相关处理采用上一帧输出图像和当前帧成像图像做加 权求和的方式实现，即为一种简单的递归滤波方法，其可以表示为
:KO，(卜l) + (l-"WO (1) 其中，少(/)为当前帧帧相关处理后输出的图像，；c(O为当前帧成像图像，
少G-1)为上一帧的输出图像，fl为帧相关系数。
其实，帧相关处理属于一种滤波方法，其在信号处理领域可以描述为一
个一阶IIR (Infinite Impulse Response )数字滤波器，该数字滤波器Z域的传
递函lt/Z(z)可以表示为将z-exp[/2;^，代入上述/f(z)即获得滤波器的频率响应H(/')。通常情况下， 帧相关处理时仅考虑该滤波器的幅频响应，即上述频率响应的幅度I/Z(f)1, 其中，/'为归一化数字频率。
才艮据上述滤波器的频率响应可知，该滤波器为一个低通滤波器，由于不 同帧的电路噪声通常是独立，而不同帧的斑点噪声由于人体组织或者探头 的移动，导致不同帧的斑点噪声也具有较弱的相关性，而且不同帧的组织 或者血流的图像具有较强的相关性，因此通过该滤波器可以有效抑制图像 中的噪声，而增强有用信号，从而改善图像的对比分辨率。通过设置不同 的帧相关系数fl可以实现不同程度的噪声抑制，随着帧相关系数a的增大， 噪声抑制效果越显著，但是也容易引起运动组织图像的模糊。同样的帧相 关系数a在不同的图像成像帧率下表现出不同的噪声抑制和运动组织模糊 效果，如在成像帧率较高的时候，设置较大的帧相关系数a则可以有效地降 低图像噪声，且不会引起运动模糊，但是当成像帧率较低的时候，同样的 帧相关系数"则会引起明显的运动模糊。
现有的超声诊断系统中，常将帧相关系数"设置为成像帧率的指数函数， 其具体表示为 -
其中，A:为帧相关控制系数，Fr为成像帧率。可见，随着成像帧率Fr的降低, 帧相关系数。也随之降低， 一定程度上改善了成像帧率&对图像时域平滑的
将公式(3)所表示的帧相关系数fl代入到前述公式(2)所表示的数字滤波 器Z域的传递函数中，
为例子，不同成像帧率下，相同帧相关控制系数A所对应的滤波器幅频响 应如图1所示，其分别给出了帧相关控制系数A:为0.02,帧率为15Hz,30Hz，
其中，
<formula>formula see original document page 9</formula>
，即可获得对应的滤波器幅频响应|//(/)|。作45Hz和60Hz的情况下，通过上述方法获得的四条幅频响应曲线。实系数 数字滤波器的幅频响应对应的频率范围是零到当前成像帧率的二分之一， 即戶0 i^/2，由图1可见，不同帧率下的幅频响应曲线随着帧率的降低所 获得滤波器对高频噪声的抑制能力也随之降低，帧率变化越大的时候噪声 抑制效果差别也越大，因此在帧率变化较大时，帧相关处理的高频衰减明 显下降，从而降低了帧相关对斑点噪声和电路噪声的抑制作用。
依据上述分析，在成像帧率/^变化不大的时候，采用上述公式(3)筒单 求解帧相关系数a进行帧相关处理时，其帧相关效果相近，但是，在成像 帧率Fr变化较大时，就可能出现噪声抑制效果变差的情况，比如超声诊断 系统从组织灰度(B-mode)成像模式切换到彩色血流成像模式，通常成像 帧率i^从几十帧下降到几帧，此时组织灰度图像的噪声会突然增大，如果 这时还采用上述公式(3)简单求解帧相关系数a进行帧相关处理的话，就无 法保证图像的帧相关效果，会出现图像失真，降低图像噪声的滤除效果。
可见，现有技术中获得帧相关系数"的方法还存在一定的弊端，需要 进一步地改进。

发明内容
本发明的目的在于提供一种帧相关系数的控制方法及其装置、以及利 用该方法的帧相关处理方法及超声成像系统，其能够获得帧相关系数，并 且在成像帧率变化很大时，使用本发明获得的帧相关系数能够使图像的帧 相关处理结果非常接近，保证超声图像噪声滤除效果，以及有效改善图像 的信p桑比。
为了实现上述目的，本发明釆用如下技术方案
本发明提供了一种帧相关系数的控制方法，所述方法包括A、比较当 前成像帧率下的当前滤波器幅频响应与一基准成像帧率下的基准滤波器幅 频响应，获得所述当前滤波器幅频响应与基准滤波器幅频响应在重叠频带内的差异；B、将所述差异为最小时的当前滤波器幅频响应所对应的帧相关 系数，判定为当前成像帧率所对应的最佳帧相关系数。
在上述方法的基础上，本发明还提供了一种帧相关系数控制器，所述 控制器包括第一单元，用于根据预设的基准成像帧率计算基准成像帧率 下的帧相关系数；第二单元，用于根据第一单元获得的帧相关系数计算基 准滤波器幅频响应；及第三单元，用于比较当前成像帧率下的当前滤波器 幅频响应与所述基准滤波器幅频响应在重叠频带内的差异，输出所述差异 为最小时的当前滤波器幅频响应所对应的帧相关系数。
在上述方法的基础上，本发明还提供了一种超声图像帧相关处理方法， 所述方法包括以下步骤
51、 比较当前成像帧率下的当前滤波器幅频响应与一基准成像帧率下 的基准滤波器幅频响应，获得所述当前滤波器幅频响应与基准滤波器幅频 响应在重叠频带内的差异；
52、 将所述差异为最小时的当前滤波器幅频响应所对应的帧相关系数， 判定为当前成像帧率所对应的最佳帧相关系数；
53、 利用获得的最佳帧相关系数进行帧相关运算，获得处理后的当前 帧图像。
针对上述方法，本发明还提供了一种超声成像系统，所述系统包括 超声图像处理模块和用于获取超声图像的探头；所述超声图像处理模块中 设置有一帧相关系数控制器和一帧相关处理单元，该控制器包括第一单 元，用于根据预设的基准成像帧率计算基准成像帧率下的帧相关系数；第 二单元，用于根据第一单元获得的帧相关系数计算基准滤波器幅频响应； 及第三单元，用于比较当前成^f象帧率下的当前滤波器幅频响应与所述基准 滤波器幅频响应在重叠频带内的差异，输出所述差异为最小时的当前滤波 器幅频响应所对应的帧相关系数；所述帧相关处理单元用于接收所述第三 单元输出的帧相关系数，并利用该帧相关系数进行帧相关运算，获得处理后的当前帧图像。
另外，本发明还提供了一种超声成像系统，所述系统包括超声图像 处理模块和用于获取超声图像的探头；所述超声图像处理模块中设置有一 存储介质、查找单元和一帧相关处理单元，所述存储介质保存有一关系表， 该表记录有当前成像帧率、当前成像帧率下的最佳帧相关系数、及表征帧 相关效果的帧相关控制系数三者之间的对应关系，所述最佳帧相关系数对 应的当前滤波器幅频响应与一基准成像帧率下的基准滤波器幅频响应在重 叠频带内的差异最小；所述帧相关处理单元通过所述查找单元与所述存储 介质相连，所述查找单元用于查找所述关系表中与当前成像帧率及当前帧 相关控制系数相对应的最佳帧相关系数，并将该最佳帧相关系数送入到所 述帧相关处理单元中；所述帧相关处理单元用于利用所述最佳帧相关系数 进行帧相关运算，获得处理后的当前帧图像。
从上述技术方案可知，本发明通过建立当前滤波器幅频响应与一基准 成像帧率下的基准滤波器幅频响应在重叠频带内的最小差异原则，选择最 佳的帧相关系数，用以超声图像的帧相关处理，相对于现有技术而言，本 发明采用逆向思维从所要获得的效果出发计算最佳的帧相关系数，从而能 获得比现有技术更加合适的帧相关系数，优化超声图像时域平滑滤波系数, 使得利用该方法的帧相关处理方法及其超声成像系统，具有更好的滤波效 果、并且在成像帧率变化很大时，也能够保证图像的显示效果、以及使图 像的巾贞相关处理结果非常接近。


图l是现有技术获得的滤波器幅频响应特性图； 图2是本发明帧相关系数控制方法的流程图； 图3是本发明实施例1至3的方法流程图4是本发明利用实施例1的方法所获得的帧相关控制系数A为0.02,帧率依次为15Hz、 30Hz、 45Hz和60Hz所对应的幅频响应特性图5是本发明利用实施例4的方法所获得的帧相关控制系数A为0.02,
帧率依次为15Hz、 30Hz、 45Hz和60Hz所对应的幅频响应特性图； 图6是本发明帧相关系数控制器的结构示意图； 图7是图6中第三单元的内部结构示意图8是本发明提供的一种超声成像系统中涉及到图像处理部分的结构 示意图9是本发明提供的另一种超声成像系统中涉及到图像处理部分的结 构示意图。
具体实施例方式
为了在成像帧率变化很大时，也能够保证图像的帧相关处理结果非常 接近，本发明给出了一种帧相关系数的控制方法，其精确计算了当前帧率 下的帧相关系数。如图2所示，本发明所给出的帧相关系数控制方法的主 要过程如下
首先，比较当前成像帧率i^>t下的当前滤波器幅频响应|/7(/")|与一基准成 像帧率i^o下的基准滤波器幅频响应l//。01,获得当前滤波器幅频响应|//(01 与基准滤波器幅频响应li/o(/" )l在重叠频带内的差异五，即￡=F(IW(/)|， |//0(/)|), F表示差异函数；
然后，将所述差异￡为最小时的当前滤波器幅频响应|//(/)|所对应的帧 相关系数，判定为当前成像帧率所对应的最佳帧相关系数。
可见，本发明在现有技术的基础上，给出了一种新的帧相关系数计算 方法，该方法定义了一个最小差异原则，其计算获得的帧相关系数使得当 前滤波器幅频响应|//(/)|与基准滤波器幅频响应li/o(/)i的差异最小，也就是 说此时在重叠频带内，滤波器在当前成像帧率i^和基准成像帧率fh下具 有相近的幅频响应特性。只有在图像处理时保证不同成像帧率下的滤波器幅频响应都接近于基准幅频响应，才能保证不同成像帧率下的滤波器幅频 响应之间在重叠频带内具有相近的幅频响应特性，从而有效保证了图像处 理的平滑处理效果以及噪声的抑制作用。
如图2所示，本发明的帧相关系数控制方法需要设定一个基准的频率 /70(如实时成像帧率，即50HZ或60H)下的帧相关系数a(F")。此帧相关 系数a(F")可以根据前述公式(3)、通过选择帧相关控制系数A:确定来计算， 帧相关控制系数A:对应超声成像系统的帧相关档位，其表征图像的帧相关 效果。以帧相关系数aCFVo)对应的幅频响应l^(/)l作为参考函数，当当前成 像帧率i^大于i^。的时候，平滑处理要求当前成像帧率fVt下帧相关系数 a(i^t)对应的滤波器幅频响应l/7(/)i在0 ~7^)/2的范围内具有相近的响应；而， 当当前成像帧率^Vv卜于fh)的时候，要求当前成像帧率Frt下帧相关系数 a(i^t)对应的滤波器幅频响应l/Z(/)l在范围内具有相近的响应。可见， 两个滤波器幅频响应的重叠频带可以表示为[O, min(/^/2， /^>t/2)]，通过设 定当前成像帧率i^下的帧相关系数"(i^)使得对应的滤波器幅频响应在频 带[O, min(^V2， i^/2)]内相近，实现图像的平滑处理，从而保证在成像帧 率变化较大的情况下，也能够获得基本一致的帧相关处理效杲。
图2所示的帧相关系数的控制方法，按照差异函数￡=尸(1//(/)|， 的建立方式，可以存在多种实现方式，以下针对常用的几种方式进行详细 说明。
实施例1:如图3所示，本实施例中所述差异五为基准滤 波器幅频响应li/o(/)l与当前滤波器幅频响应|//(/)|之差在重叠频带
内的平方和，即
如图3所示，本实施例的具体过程如下
步骤200，预先设定的基准成像帧率i^，这里的基准成像帧率&0为实时成像帧率，如50HZ或60HZ。
步骤210，将预先设定的基准成像帧率i^，代入前述公式(3)中计算其 帧相关系数a(F"),结果如下所示
其中，A为帧相关控制系数，其由经验来确定，此A:值决定了图像的帧 相关效果，并与超声探测仪器的帧相关档位相对应，其取值范围为0 1。
步骤220,将上述计算获得帧相关系数a(F")输入到数字滤波器的传递 函数模型中，计算滤波器的基准幅频响应|//0(/)|,其对应的频带范围为0~ f>0/2。
在这里，以数字滤波器的传递函数模型属于一阶无限脉冲响应(Infinite Impulse Response ) ^t字滤波器Z域的传递函数才莫型为例，详细说明基准幅 频响应IZ/o(/)i的计算过程。
根据前述公式(2)<formula>formula see original document page 15</formula>
将z-exp(/2,)， /'= //_F>，/e
，代入上述H(z)中即获得滤波器的频率 响应H(/)，其幅度对应的幅频响应IH(/)I表达式为
其中，/'为归一化数字频率，/为实际频率，」FV为成像巾贞率。
根据上述公式，将帧相关系数a(fh))、基准成像帧率i^V。代入到上述公式 (5)中，即可获得关于基准成像帧率片o下的基准幅频响应li/()(/)l的函数表达 式。
步骤230，设定当前成像帧率Frt下的N个候选帧相关系数a、(Frt)，其中N 为大于2的自然数；然后再分别将此N个候选帧相关系数"XPrt)输入到数字 滤波器的传递函数模型(如上述公式(5)所示的一阶无限脉冲响应数字滤波器Z域的传递函数冲莫型)中，计算各个候选帧相关系数fl、(尸^t)对应的滤波器 当前幅频响应l^(/)卜其对应的频带范围为0 Frt/2。
这里的N个候选帧相关系数a 、(fh)可以根据经验来确定，也可以按照以 下方式来确定
第 一步，按照常规方法计算当前成像帧率i^下的帧相关系数估计值 例如，将当前成像帧率i^t代入前述公式(3)所示的帧相关系数与成像 帧率的指数函数中，计算其帧相关系数估计值^，结果如下所示
从上述公式可以看出，不同的帧相关控制系凄U对应不同帧相关效果的 候选帧相关系数a、(/、)，在计算该估计值^。时，其Hi可以参照步骤210 中所选的值。
第二步，在所述估计值^刁的附近选择N个值作为候选帧相关系数 a、(^Vt),该N为大于等于2的自然数。例如，在以^^值为中心的区间内取 IOO个数值，作为候选帧相关系数o、(F^)。
步骤240,在基准幅频响应|//0(/)|和当前幅频响应|//(/)|的重叠频带内，分 别比较所述基准幅频响应1/7。(/)l与各个候选帧相关系数"V^)所对应的滤波 器当前幅频响应|//(/)|之间的差异￡,即按照上述公式(4)计算在
内各个^(Frt)所对应的五值，根据差异最小原则，从各 个候选帧相关系数"XFrt)中选择一个最优的作为当前成像帧率i^h下的最佳 帧相关系数fl(Fn),即将￡值最小时的|7/(/)|所对应的候选帧相关系数 a V^t)作为当前成像帧率i^下的最佳帧相关系数a(i^t)。
本实施例滤波器系数是根据当前滤波器幅频响应|//(/)|在频带
内与a(i^。)对应的滤波器幅频响应l/fo(/)l具有最小的均 方误差原则来确定的。可见，幅频响应函数为帧相关系数a和频率/的非线 性函数，当前成像帧率^V,下对应的帧相关系数a(F^)可以通过数值法求解。利用本实施例的实现方式获得的不同成像帧率下的帧相关系数对应的幅频
响应如图4所示。由图4可见，在重叠频带内、不同帧率下的幅频响应也 非常相近，确保了不同帧率情况下组织运动的模糊程度相近，且对噪声也 有较好地抑制效果。
当然上述差异的计算方式并不只一种，请见下述实施例。
实施例2:本实施例与实施例1的不同点在于，在步骤240中，本实施
例中所述差异￡为基准滤波器幅频响应IZ/o(/)l与当前滤波器幅频响 应|//(/)|之差在重叠频带
内的绝对值之和，即
￡= SH/)卜。 (6)
/e0，min(Fro/2，F/v/2)
通过比较该差异值￡的大小，将该差异值￡为最小时的当前滤波器幅 频响应|//(/)|所对应的候选帧相关系数a'CFVt)确定为当前成像帧率Fn下的最 佳帧相关系数a(Frt)。本实施例的其余步骤与实施例1相同。本实施例滤 波器系数是根据当前滤波器幅频响应|7/(/)|在频带
内与《fy。)对应的滤波器幅频响应i7/o(/)l具有最小的误差绝对值之和的原则 来确定的。同实施例1，当前帧率i^下对应的帧相关系数a(i^)可以通过 数值法来求解。利用本实施例的实现方式获得的不同成像帧率下的帧相关 系数对应的幅频响应与图4类似。'
实施例3:本实施例与实施例1的不同点在于，在步骤240中，本实施 例中所述差异￡为在重叠频带
内，在一设定频率 力下基准滤波器幅频响应l/Zo(/)l与当前滤波器幅频响应IZ/(/)i之差的绝对值， 即
E卞(/i)l-l仏(川卜 (7) 其中，设定频率力为上述重叠频带(O, minCF)V2， &/2)]内除零外的任 一值。通过设定一个频率/1,然后根据上述公式(7)计算并比较差异值五的 大小，将该差异值五为最小时的当前滤波器幅频响应|//(/)|所对应的候选帕相关系数a、(fVt)确定为当前成像帧率Frt下的最佳帧相关系数"CFVt)。本实 施例的其余步骤与实施例l相同。本实施例滤波器系数是根据在频带
内当前滤波器幅频响应与"(尸")对应的滤波器 幅频响应IFo(/)l之差的绝对值最小原则来确定的。同实施例l,当前帧率i^ 下对应的帧相关系数a(/^t)可以通过数值法来求解。
实施例4:本实施例与实施例3的不同点在于，在步骤240中，本实施 例中所述差异为最小的情况是指在一设定频率/i下基准滤波器幅频响应 li/o(/)l与当前滤波器幅频响应l/7(/)l之差为零，即
^1//(/梓。(川0， (8)
也就是，在设定频率力上使|//0(/)|=|//(/)1，此时当前滤波器幅频响应|//(/)| 所对应的候选帧相关系数aX^h)被判定为当前成像帧率下的最佳帧相关 系数fl(i^t)，其中/产(0， min(i^)/2,&t/2)]。本实施例的其余步骤与实施例3 相同，该其余步骤是指除步骤230设定候选帧相关系数、及与候选帧相关 系数相关的步骤外。
本实例是通过设定"(PVt)使得当前滤波器幅频响应li/(/)l在某预先设定 的频率力(Zi不等于0)上具有与基准滤波器幅频响应|//0(/)|相同的数值，即 令|^(/)|=|//(/)|，其中，频率力最优可选择为min(F /2， &/2),即Fr0/2和f>t/2 两者中的较小者。并，通过将"(i^j)、 yi代入幅频响应函数|//()(/)|，并将/i 代入幅频响应函数|//(/)|,通过解析求解上述等式(8 )即获得当前成像帧率 下的最佳帧相关系数a(F、)。
本实施例属于实施例3中的一种特殊情况，采用本实施例的解析求解 后，可以简化上述如图3所示的计算规则，如可以减少设定候选帧相关系 数的步骤等等，因为上述公式(8)是一个确定的等式，其可以直接计算出最 佳帧相关系数flCFVt)。本实施例相对于实施例1至3而言，本实施例的实现 方式由于计算量较小，可以由成像系统软件根据当前成像帧率Frt和帧相关 档位对应的帧相关控制系数A:实时计算获得当前成像帧率^PVt的帧相关系数aCF>t)。利用本实施例的实现方式获得的不同成像帧率下的帧相关系数对应 的幅频响应如图5所示。由图5可见，在重叠频带内幅频响应也非常相近， 确保了不同帧率情况下组织运动的模糊程度相近，且对噪声也有较好的抑 制效果。
实施例5:为了简化计算步骤，本实施例在上述实施例l、 2、 3或4的 基础上，还增加以下的建表步骤
通过实施例1、 2、 3或4确定了当前成像帧率^Vt下的最佳帧相关系数 a(^Vt)后，将此最佳帧相关系数a(F^t)保存于一关系表中，该表记录当前成 像帧率F V当前成像帧率Frt下的最佳帧相关系数fl(Frt)、及表征帧相关效 果的帧相关控制系数A三者之间的对应关系；
然后，在实时图像处理时，可以根据当前成像帧率i^及帧相关控制系 数A;，查找所述关系表，实时确定当前成像帧率i^下的最佳帧相关系数 fl(^>t)。利用本实施例建立的关于最佳帧相关系数的关系表可以存储于一计 算机可读存储介质中。
由于上述实施例1至3的数值法求解运算量较大，因此本实施例通过 预先计算好不同成像帧率下的最佳帧相关系数，将这些帧相关系数存储为 一个系数表。当需要设置不同的帧相关效果的时候，可以通过更改基准帧 率i^o下的帧相关控制系数A:,获得不同的IWo(/)1，然后数值计算出相应的不 同成像帧率Frt下的帧相关系数a，从而组成了不同帧相关效果和不同成像 帧率的一个二维系数表，然后将利用实施例l、 2或3求解出的当前成像帧 率下的最佳帧相关系数"(^Vt)加入其中，即可获得上述关系表。超声成像系 统根据不同帧相关效果的档位和成像帧率对上述关系表进行查表，即可获 得当前成像帧率下的帧相关系数。
上述分析只是列举了常见的几种差异的计算方式，本发明对差异如何 计算并无限制。比如，根据上述各个实施例的差异数学模型可以概括的说， 差异￡的数学模型可以表示为基准滤波器幅频响应|//。(/)|与当前滤波器幅频响应|7/(/)|之差的绝对值在重叠频带
内 的w次方之和，即
五=z^(/)i一i压wr, ("
./"e
内的2w次方之和，即
/e[0,min — /2，fA/2)
其中，w为大于零的自然数。当上述公式(10)中"=1时，即上述实施例 l所示的情况；当"=2时，同样也是可取的计算方式。
由上可见，本发明提出的帧相关系数的控制方法不局限于组织灰度成 像处理，还可以适用于彩色血流成像、组织多普勒成像等各类成像模式， 其能精确的确定最佳帧相关系数用以帧相关处理中，保证图像平滑处理的 效果。
通常有效的图像信息在相邻帧之间变化不大，而无效的噪声信息则通 常变化较大，利用帧相关所采用的低通滤波器特性，可以较好保持帧间变 化不大的低频有效图像信息，并抑制帧间变化较大的高频无效噪声。采用 现有技术获得的图l所示的帧相关滤波器幅频响应曲线中，可以明显地看
出各个曲线间相差很大，存在较多的噪声信息；而采用本发明的方法所获 得的图4和图5所示的帧相关滤波器幅频响应曲线中，各个相邻曲线之间 基本具有一致的帧相关效果，并且各个曲线在相互重叠的频带内具有相似 的幅频特性。可见，本发明通过对不同的成像帧率，优化选择相应的帧相关系数从而获得不同成像帧率下的基本一致的帧相关效果，有效地保证了 图像处理的平滑处理效果以及噪声的抑制作用。
依据上述帧相关系数的控制方法，本发明还相应提供了一种帧相关系 数控制器，其能够直接获得当前成像帧率下的最佳帧相关系数，保证不同 成^(象帧率下的滤波器幅频响应在重叠频带内具有相近的频率响应特性。如
图6所示，所述控制器包括
第一单元300,用于根据预设的基准成像帧率^V。计算基准成像帧率Fr0 下的帧相关系数"(i^o);
第二单元310，用于根据第一单元300获得的帧相关系数a(i^)计算基 准滤波器幅频响应|7/()(/)|;及
第三单元320，用于比较当前成像帧率Frt下的当前滤波器幅频响应 |//(/)|与所述基准滤波器幅频响应|//()(/)|在重叠频带内的差异￡，输出所述差 异五为最小时的当前滤波器幅频响应|//(/)|所对应的帧相关系数"(^Vt)。
这里的参数a(i^)、 Fr0、 |//0(/)|、 /^、五、|7/(/)1、 "(&0、重叠频带的计 算或设定原则，可以参见上述帧相关系数的控制方法的实施过程。
如图7所示，为了计算差异五，可以在控制器的第三单元320中建立以 下单元
差异模型单元321,用于存储所述当前滤波器幅频响应|//(/)|与所述基准 滤波器幅频响应l7/o(/)l的差异数学模型￡=F(l//(/)|， li/o(/)卩，并利用该数学模
式计算所述当前滤波器幅频响应|//(/)|与所述基准滤波器幅频响应|//0(/)1的 差异￡;及
比较单元322，用于比较所述差异五，输出所述差异五为最小时的当前 滤波器幅频响应I/Z(/)I所对应的帧相关系数。采用此结构，并利用上述实施 例4的话，即可以实现实时计算当前成像帧率下的最佳帧相关系数；但是 如果采用上述实施例l、 2或3时，涉及到需要从多个候选帧相关系数中进 行选择时，所述第三单元320就需要在差异模型单元321和比较单元322的基础上再增加以下几个单元，如图7所示，所述第三单元320还包括 候选组计算模块323,用于生成当前成像帧率i^t下的候选帧相关系数
滤波器幅频响应计算模块324，用于将所述候选帧相关系数fl、(fh)输入 到数字滤波器的传递函数模型325中，计算候选帧相关系数a、(/^t)对应的 当前滤波器幅频响应l/f(/)1,将此当前滤波器幅频响应|//(/)|输入到所述差异 模型单元321中进行差异计算，并通过比较单元322比较后获得相应的帧 相关系数。
基于数据结果需要保存、以及建立快捷用关系表的需求，本发明的帧 相关系数控制器可以在图6所示结构的基础上，增加一存储单元330,其用 于记录当前成像帧率F^下的最佳帧相关系数"(i^Vt)，该最佳帧相关系数 fl(FrO为所述差异￡为最小时的当前滤波器幅频响应li/(/)l所对应的帧相关系 数，即从选帧相关系数a、(Frt)中按照差异最小原则选出的帧相关系数。这 里的存储单元330可以设置一关系表，并将记录的最佳帧相关系数"(Frt) 添加其中，利用该关系表记录当前成^f象帧率Frt、当前成^^帧率Fn下的最 佳帧相关系数flCF>t)、及表征帧相关效果的帧相关控制系数A:三者之间的对 应关系。
从上述帧相关系数控制器的结构描述可以看出，其主要是利用上述帧 相关系数的控制方法而建立一种产品，该产品采用上述帧相关系数的控制 方法能够获得当前帧率下的最佳帧相关系数，利用该帧相关系数进行超声 图像处理，可以保证在重叠频带内不同成像帧率下的滤波器幅频响应具有 相近的频率响应特性，并且有效保证了图像处理的平滑效果及噪声抑制效 果。
基于上述帧相关系数的控制方法在帧相关处理中的应用，本发明还提 供了一种超声图像帧相关处理方法，该方法包括以下步骤
步骤Sl，采用上述帧相关系数的控制方法获得当前成像帧率片t所对应的最佳帧相关系数a(Frt);
步骤S2，将获得的最佳帧相关系数a(/^)输入到帧相关处理的数学模 型中，进行帧相关运算，获得处理后的当前帧图像。这里的帧相关处理数 学模型可以采用前述公式(l),其具体实施例可参见图3。由图3可见，在 上述实施例1、 2或3计算出最佳帧相关系数a(i^t)后，只要利用该帧相关 系数进行图像的帧相关处理即可。
本发明所提出的超声图像帧相关处理方法中，关于最佳帧相关系数 "(/^t)的计算过程以及其拓展的技术方案参见上述关于帧相关系数控制方法 的说明，在此不再重复说明。
在上述帧相关系数控制方法及其帧相关处理方法的基础上，本发明还 提供了一种新的超声成像系统，如图8所示，所述超声成像系统包括超 声图像处理模块和用于获取超声图像的探头；所述超声图像处理模块中设 置有一帧相关系数控制器400和一帧相关处理单元410，这里的帧相关系数 控制器400采用如图6所示的结构，具体可参见前述说明；所述帧相关处 理单元410用于接收帧相关系数控制器400的第三单元320输出的帧相关 系数，并将该帧相关系数输入到帧相关处理的数学模型中进行帧相关运算， 获得平滑处理后的当前帧图像。这一超声成像系统中，关于帧相关系数控 制器400的内容以及其拓展的技术方案参见上述关于图6所示帧相关系数 控制器的说明，在此不再重复说明。
基于上述最佳帧相关系数可以采用查表方式进行的实现方式，本发明 还提供了一种新的超声成像系统，如图9所示，该系统中包括超声图像 处理模块和用于获取超声图像的探头；并且，在超声图像处理模块中设置 有一存储介质430、查找单元420和一图6所示的帧相关处理单元410,所 述存储介质430保存有一关系表，该表记录有当前成像帧率、当前成像帧 率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数三者之间的 对应关系，所述最佳帧相关系数对应的当前滤波器幅频响应与一基准成像帧率下的基准滤波器幅频响应在重叠频带内的差异最小；所述帧相关处理 单元410通过所述查找单元420与所述存储介质430相连，所述查找单元 420用于查找所述关系表中与当前成像帧率及当前帧相关控制系数相对应 的最佳帧相关系数，并将该最佳帧相关系数送入到所述帧相关处理单元410 中；所述帧相关处理单元410用于将所述最佳帧相关系数输入到帧相关处 理的数学模型中进行帧相关运算，获得处理后的当前帧图像。
图9所示的超声成像系统在于釆用上述帧相关系数控制方法所获得 的最佳帧相关系数来建立一关系表，用于在图像处理时实时查表获得相应 的帧相关系数，这样做可以减少计算量，提高图像处理的速度。建表的关 键问题在于如何计算最佳帧相关系数，其可以参见上述关于帧相关系数控 制方法的具体说明。比如如图采用上述实施例l、 2或3的方式，最佳帧 相关系数可以采用如下方法获得
首先，设定当前成像帧率下的多个候选帧相关系数，该候选帧相关系 数的选择方法如下
按照帧相关系数与成像帧率的指数函数，计算当前成像帧率下的帧相 关系数估计值；在所述估计值的附近选择N个值作为候选帧相关系数，该 N为大于等于2的自然数；
然后，分别将所述N个候选帧相关系数输入到数字滤波器的传递函数 模型中，计算各个候选帧相关系数对应的当前滤波器幅频响应；
最后，分别将所述N个候选帧相关系数对应的当前滤波器幅频响应与 所述基准滤波器幅频响应(其计算过程可以参见上述实施例1的相关说明) 进行比较，根据当前滤波器幅频响应与所述基准滤波器幅频响应在重叠频 带内的差异为最小的原则，从所述N个候选帧相关系数中选出 一个作为当 前成像帧率所对应的最佳帧相关系数，即所述差异为最小时的当前滤波器 幅频响应所对应的候选帧相关系数。但如果采用实施4或5的话，则可以 参见上述具体说明。综上所述，本发明通过设定不同帧率下的帧相关系数，使得对应的滤 波器幅频响应在重叠的频带内具有相近的频率响应特性的要求，提供了一 种帧相关系数的计算方法，基于此计算方法，本发明又衍生出了一种可以 制成独立产品并用于获得最佳帧相关系数的帧相关控制器、 一种利用上述 帧相关系数计算方法的帧相关处理方法及利用该方法的超声成像系统、以 及一种采用含有最佳帧相关系数关系表的存储介质的超声成像系统。对于 最后一个超声成像系统，其存储介质中含有的最佳帧相关系数关系表采用 上述帧相关系数计算方法获得。本发明所公开的各种方法和系统，均可以 根据最小差异原则计算获得当前成像帧率下的最佳帧相关系数，并用于后 续的图像处理中，使得本发明在成像帧率发生较大改变时，超声图像帧相 关的处理效果基本保持不变。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以 改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护 范围。
权利要求
1、一种帧相关系数的控制方法，其特征在于，所述方法包括A、比较当前成像帧率下的当前滤波器幅频响应与一基准成像帧率下的基准滤波器幅频响应，获得所述当前滤波器幅频响应与基准滤波器幅频响应在重叠频带内的差异；B、将所述差异为最小时的当前滤波器幅频响应所对应的帧相关系数，判定为当前成像帧率所对应的最佳帧相关系数。
2、 根据权利要求1所述的方法，其特征在于，所述重叠频带的范围是 零到T值的闭区间，该T值为&0/2和Fn/2两者中的较小者，其中i^o为 基准成像帧率，i^为当前成像帧率。
3、 根据权利要求1所述的方法，其特征在于，所述步骤A包括以下步骤Al、设定当前成像帧率下的多个候选帧相关系数；A2、依据数字滤波器的传递函数，计算各个候选帧相关系数对应的当 前滤波器幅频响应；A3、分别将所述各个候选帧相关系数对应的当前滤波器幅频响应与所 述基准滤波器幅频响应进行比较，获得各个当前滤波器幅频响应与基准滤 波器幅频响应在重叠频带内的差异；并且，所述步骤B中，通过比较各个候选帧相关系数所对应的所述差 异，获取最小的所述差异。
4、 根据权利要求l、 2或3所述的方法，其特征在于，所述差异为 所述基准滤波器幅频响应与所述当前滤波器幅频响应之差的绝对值在重叠 频带内的w次方之和，其中，w为大于零的自然数。
5、 根据权利要求1、 2或3所述的方法，其特征在于，所述差异为 在一设定频率下的所述基准滤波器幅频响应与所述当前滤波器幅频响应之 差。
6、 根据权利要求5所述的方法，其特征在于，所述设定频率为所述重 叠频带内除零外的任一值。
7、 根据权利要求5所述的方法，其特征在于，所述差异为最小的情况 是指在一设定频率下的所述基准滤波器幅频响应与所述当前滤波器幅频 响应之差为零。
8、 根据权利要求7所迷的方法，其特征在于，所述设定频率为F"/2 和f>t/2两者中的较小者。
9、 根据权利要求3所述的方法，其特征在于，所述候选帧相关系数的 选择方法包括以下步骤按照帧相关系数与成像帧率的指数函数，计算当前成像帧率下的帧相 关系数估计值；在所述估计值的附近选择N个值作为候选帧相关系数，该N为大于等 于2的自然数。
10、 根据权利要求l、 2或3所述的方法，其特征在于，所述方法还包 括以下步骤C、 将所述最佳帧相关系数保存于一关系表中，该表记录当前成像帧率、 当前成像帧率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数 三者之间的对应关系；D、 根据当前成像帧率及帧相关控制系数，查找所述关系表，确定当前 成像帧率下的最佳帧相关系数。
11、 一种帧相关系数控制器，其特征在于，所述控制器包括 第一单元，用于根据预设的基准成像帧率计算基准成像帧率下的帧相关系数；第二单元，用于根据第 一单元获得的帧相关系数计算基准滤波器幅频 响应；及第三单元，用于比较当前成像帧率下的当前滤波器幅频响应与所述基 准滤波器幅频响应在重叠频带内的差异，输出所述差异为最小时的当前滤 波器幅频响应所对应的帧相关系数。
12、 根据权利要求11所述的控制器，其特征在于，所述第三单元包括 差异模型单元，用于计算所述当前滤波器幅频响应与所述基准滤波器幅频响应的差异；及比较单元，用于比较所述差异，输出所述差异为最小时的当前滤波器 幅频响应所对应的帧相关系凄t。
13、 根据权利要求12所述的控制器，其特征在于，所述第三单元还包 括 候选组计算模块，用于生成当前成像帧率下的候选帧相关系数；及 滤波器幅频响应计算模块，用于将所述候选帧相关系数输入到数字滤波器的传递函数中，计算候选帧相关系数对应的当前滤波器幅频响应，将此当前滤波器幅频响应输入到所述差异^t型单元中。
14、 根据权利要求11或12所述的控制器，其特征在于，所述控制器 还包括设置有一关系表的存储单元，所述关系表用于记录当前成像帧率、当 前成像帧率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数三 者之间的对应关系；该最佳帧相关系数为所述差异为最小时的当前滤波器幅频响应所对应的帧相关系凄史。
15、 一种超声图^f象帧相关处理方法，其特征在于，所述方法包括以下 步骤S1 、比较当前成像帧率下的当前滤波器幅频响应与 一基准成像帧率下 的基准滤波器幅频响应，获得所述当前滤波器幅频响应与基准滤波器幅频 响应在重叠频带内的差异； S2将所述差异为最小时的当前滤波器幅频响应所对应的帧相关系数， 判定为当前成像帧率所对应的最佳帧相关系数；S3、利用获得的最佳帧相关系数进行帧相关运算，获得处理后的当前 帧图像。
16、 根据权利要求15所述的方法，其特征在于，所述步骤S2与步骤 S3之间包括以下步骤将所述最佳帧相关系数保存于 一关系表中，该表记录当前成像帧率、 当前成像帧率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数 三者之间的对应关系；根据当前成像帧率及帧相关控制系数，查找所述关系表，确定当前成 像帧率下的最佳帧相关系数。
17、 一种超声成像系统，所述系统包括超声图像处理模块和用于获 取超声图像的探头；其特征在于，所述超声图像处理模块中设置有一帧相 关系数控制器和一帧相关处理单元，该控制器包括第一单元，用于^f艮据预设的基准成像帧率计算基准成像帧率下的帧相 关系数；第二单元，用于才艮据第 一单元获得的帧相关系数计算基准滤波器幅频 响应.，及第三单元，用于比较当前成像帧率下的当前滤波器幅频响应与所述基 准滤波器幅频响应在重叠频带内的差异，输出所述差异为最小时的当前滤波器幅频响应所对应的帧相关系数；所述帧相关处理单元用于接收所述第三单元输出的帧相关系数，并利 用该帧相关系数进行帧相关运算，获得处理后的当前帧图像。
18、 根据权利要求17所述的超声成像系统，其特征在于，所述第三单 元包括差异模型单元，用于计算所述当前滤波器幅频响应与所述基准滤波器 幅频响应的差异；及比较单元，用于比较所述差异，输出所述差异为最小时的当前滤波器 幅频响应所对应的4美选帧相关系数。
19、 根据权利要求18所述的超声成像系统，其特征在于，所述第三单 元还包括候选组计算模块，用于生成当前成像帧率下的候选帧相关系数；及 滤波器幅频响应计算模块，用于将所述候选帧相关系数输入到数字滤波器的传递函数中，计算候选帧相关系数对应的当前滤波器幅频响应，将此当前滤波器幅频响应输入到所述差异模型单元中。
20、 根据权利要求17或18所述的超声成像系统，其特征在于，所述 超声成像系统还包括设置有一关系表的存储单元，所述关系表用于记录当前成像帧率、当 前成像帧率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数三 者之间的对应关系；该最佳帧相关系数为所述差异为最小时的当前滤波器 幅频响应所;寸应的帧相关系#t。
21、 一种超声成像系统，所述系统包括超声图像处理模块和用于获取超声图像的探头；其特征在于，所述超声图像处理模块中设置有一存储 介质、查找单元和一帧相关处理单元，所述存储介质保存有一关系表，该表记录有当前成像帧率、当前成像 帧率下的最佳帧相关系数、及表征帧相关效果的帧相关控制系数三者之间 的对应关系，所述最佳帧相关系数对应的当前滤波器幅频响应与一基准成 像帧率下的基准滤波器幅频响应在重叠频带内的差异最小；所述帧相关处理单元通过所述查找单元与所述存储介质相连，所述查 找单元用于查找所述关系表中与当前成像帧率及当前帧相关控制系数相对 应的最佳帧相关系数，并将该最佳帧相关系数送入到所述帧相关处理单元 中；所述帧相关处理单元用于利用所述最佳帧相关系数进行帧相关运算， 获得处理后的当前帧图像。
22、根据权利要求21所述的超声成像系统，其特征在于，所述最佳帧 相关系数的计算过程包括以下步骤Cl、设定当前成像帧率下的多个候选帧相关系数，该候选帧相关系数 的选择方法如下按照帧相关系数与成像帧率的指数函数，计算当前成像帧率下的帧相 关系数估计值；在所述估计值的附近选择N个值作为候选帧相关系数，该 N为大于等于2的自然数；C2、依据数字滤波器的传递函数，计算各个候选帧相关系数对应的当 前滤波器幅频响应；C3、分别将所述N个候选帧相关系数对应的当前滤波器幅频响应与所 述基准滤波器幅频响应进行比较，将所述差异为最小时的当前滤波器幅频 响应所对应的候选帧相关系数，判定为所述最佳帧相关系数。
全文摘要
本发明公开了一种帧相关系数控制及帧相关处理方法、装置及超声成像系统，其帧相关系数的控制方法包括比较当前成像帧率下的当前滤波器幅频响应与一基准成像帧率下的基准滤波器幅频响应，获得所述当前滤波器幅频响应与基准滤波器幅频响应在重叠频带内的差异；将所述差异为最小时的当前滤波器幅频响应所对应的帧相关系数，判定为当前成像帧率所对应的最佳帧相关系数。本发明能够精确地获得帧相关系数，并且在成像帧率变化很大时，使得图像的帧相关处理结果非常接近，保证超声图像噪声滤除效果，以及有效改善图像的信噪比。
文档编号G06T5/00GK101645167SQ20081014278
公开日2010年2月10日 申请日期2008年8月6日 优先权日2008年8月6日
发明者羽 张 申请人:深圳迈瑞生物医疗电子股份有限公司