CN101416885B - 用于超声成像系统的帧平均电路 - Google Patents

用于超声成像系统的帧平均电路 Download PDF

Info

Publication number
CN101416885B
CN101416885B CN200710166607A CN200710166607A CN101416885B CN 101416885 B CN101416885 B CN 101416885B CN 200710166607 A CN200710166607 A CN 200710166607A CN 200710166607 A CN200710166607 A CN 200710166607A CN 101416885 B CN101416885 B CN 101416885B
Authority
CN
China
Prior art keywords
vpv
frame
present frame
comparator
average
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.)
Active
Application number
CN200710166607A
Other languages
English (en)
Other versions
CN101416885A (zh
Inventor
熊伟
李勇
王建永
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Mindray Scientific Co Ltd
Original Assignee
Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Mindray Bio Medical Electronics Co Ltd filed Critical Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority to CN200710166607A priority Critical patent/CN101416885B/zh
Priority to US12/255,570 priority patent/US8282560B2/en
Publication of CN101416885A publication Critical patent/CN101416885A/zh
Priority to US13/647,186 priority patent/US8721549B2/en
Application granted granted Critical
Publication of CN101416885B publication Critical patent/CN101416885B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Hematology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明提供一种用于在医学超声成像系统中实现彩色血流帧的帧平均处理的帧平均电路。该帧平均电路由两部分组成,迭代计算部分和非迭代部分。迭代计算部分采用IIR滤波器和查找表相结合,主要处理速度方向没有发生变化的情况。非迭代部分对彩色血流回波信号的速度、能量和壁滤波前的能量进行阈值判断,实现对速度方向发生变化的处理。通过本发明,能够有效地减弱回波数据中的闪烁点现象,减弱彩色血流图像中的噪声影响;能够判断血流速度的变化,较好地实时表现血流速度的变化;能够判断血流速度的混叠,表现真实的血流图像。

Description

用于超声成像系统的帧平均电路
技术领域
本发明涉及超声成像领域,并且更具体地,涉及一种用于超声成像系统的帧平均电路。
背景技术
对于医学超声系统成像,在B型成像时,二维剖面图中像素点的灰度代表超声回波数据的幅度;在彩色血流成像时,不同的颜色代表血管中红细胞的运动速度或者组织的运动速度。心脏和血管中的血流速度根据红细胞运动产生的多普勒效应来检测。超声成像系统通过计算超声回波信号的频率偏移来推算出血流的运动速度。超声回波信号的频率偏移正比于血管内红细胞的移动速度。频率的正向偏移表示血管内红细胞沿着探头扫描方向运动;频率的负向偏移表示血管内红细胞远离探头扫描方向运动。频率的偏移采用颜色编码后,形成二维彩色血流图。彩色血流图像通过颜色来表示血流速度的大小和方向。彩色血流图像一般是叠加在B型图像上一起显示,在血管外显示B型图像,在血管内显示血流图像,形成二维彩色图像。
如图1所示,医学超声成像系统一般由以下七部分组成:发射和接收模块100、波束合成模块200、信号处理模块300、后处理模块400、扫描显示控制模块500、显示监视器600和系统控制模块700。系统控制模块700是软硬件接口,包含软件控制器和硬件控制器。硬件控制器接收由软件写入的参数,进行相关处理,转换成以下模块的控制信号:波束合成模块200、信号处理模块300、后处理模块400和扫描显示控制模块500。另外,系统控制模块700通过硬件控制和扫描控制总线,对整个系统的时序和运行进行总体控制。
数据流处理如下。发射和接收模块100、如探头接收超声回波数据,接收的超声回波数据经过A/D转换后,进入到波束合成模块200中。波束合成模块200对超声回波数据进行处理。根据扫描模式的不同,分别生成两种不同类型的数据。一种和B型成像相关,另一种和彩色血流成像相关。与B型成像相关的超声回波数据进入到黑白信号处理模块310中;与彩色血流成像相关的超声回波数据进入到彩色信号处理模块320中。黑白信号处理模块310和彩色信号处理模块320根据系统参数和控制信号对输入的超声回波数据进行处理,然后将处理结果输出到后处理模块400。后处理模块400的输出结果输入到扫描显示控制模块500。扫描显示控制模块500将超声成像结果转换成为标准视频数据,然后输出到显示监视器600进行显示。
黑白信号处理模块310对波束合成数据进行对数压缩、包络信号检测和降采样等处理,并将和B型成像相关的超声回波数据处理成8比特的灰阶数据流。数据流中的每一个数据点表示扫描线上的采样点,数据大小表示该点对应的灰阶。
彩色信号处理模块320可以计算采样点上的速度、方差、能量和壁滤波前的能量。彩色信号处理模块320利用红细胞运动产生的多普勒效应,采用自相关算法计算出扫描剖面上血流的速度、方差、能量和壁滤波前的能量。这些数据分量分别采用8比特数据进行表示。
后处理模块400对来自黑白信号处理模块310的数据和来自彩色信号处理模块320的数据进行电影回放保存、帧平均和优先权编码等处理。处理完成后的结果输出到扫描显示控制模块500。本发明的帧平均处理在后处理模块400中完成。
在超声成像系统中,传统上需要对图像进行帧平均处理,以减弱随机噪声对图像的影响。一般是采用IIR滤波器来实现帧平均处理功能。这种滤波器采用当前输入的图像数据和上一帧图像数据作为输入数据来进行计算。滤波器函数如下:
Y(n)=A×Y(n-1)+(1-A)×X(n)
公式中,Y(n)是帧平均处理的结果,A为常数值。Y(n-1)为上一帧帧平均处理的结果,X(n)是当前新的一帧图像。在彩色血流成像时,该方法能够减少噪声和斑点,但是同时也对有效的速度信号在时间上进行了平均处理,从而不能检测出速度的微小变化和慢速血流。这种方法尤其是对速度方向的变化不敏感,而速度方向的变化对实际临床诊断有重要的帮助。同时,由于常数A使得速度衰减比较慢,所以不能够非常快地反映运动的变化。当探头已经离开移动物体后,由于帧平均的原因,仍然能够在血流图像中显示血流。
在美国专利“Adaptive Persistence Processing”US5595179中,J.Nelson Wright等人提出了一种帧平均结构。这种帧平均结构由两部分组成,一部分是迭代计算部分,另一部分是非迭代计算部分。迭代计算部分中IIR滤波器的系数为常数,非迭代计算部分采用FIR滤波器来实现,其系数通过查找表获得,实现相对简单。这种结构与传统结构相比有了很大改进,但是仍然有显而易见的缺点,不能处理速度方向快速变化的情况,特别是血流方向发生改变时的小速度血流信号。在血流速度出现混叠时,计算出的速度值和真实速度相差较大。
在美国专利“Persistence for Ultrasonic Flow Imaging”US5897502中,Tommy King-Yuen Wong等人提出了一种帧平均结构。在该结构中,IIR滤波器采用简单固定的系数,没有根据不同情况来变化系数。该发明对回波速度进行阈值判断,根据对速度的阈值判断,决定是输出当前速度值,还是输出上一帧通过处理过后的值。但是,在该发明结构中,没有对速度混叠的情况进行任何处理。该方法对速度方向发生快速变化有很好的处理,但是在速度方向没有发生变化的情况下,处理结果不够理想,不能有效地消除回波数据中的噪声;同时,在速度快速减小的情况下,该结构不能反映出这种变化。
根据目前已有的专利,可以将帧平均处理结构划分为两个部分,迭代计算部分和非迭代计算部分,大部分专利都是在对迭代计算部分的查找表和非迭代计算部分的算法进行改进,以减少噪声、减少血流图像中的闪烁点、提高对速度方向变化的反应速度并减弱混叠对彩色血流图像的影响。目前,在超声血流图形处理中,存在如下几个问题:
1.血流图像中存在闪烁点,影响图像质量;
2.混叠现象对血流图像的影响;以及
3.不能实时反映出实际的速度方向变化及血流形态。
因此,在超声成像系统中需要一种解决上述问题的帧平均电路。
发明内容
本发明主要针对上面的问题来改进帧平均处理。本发明提供一种用于超声成像系统的帧平均电路,包括:迭代计算部分,用于对当前帧的血流速度和上一帧帧平均后的血流速度进行滤波;和非迭代部分,用于进行关于当前帧的血流速度的判断并根据判断结果来输出当前帧的血流速度或迭代计算部分的输出值作为当前帧帧平均后的血流速度。
在本发明的一个实施例中,迭代计算部分包括:滤波器;和系数提供单元,用于向滤波器提供系数。
在本发明的一个实施例中,系数提供单元包括:存储器,用于存储查找表,查找表用于存储滤波器的系数;和地址生成器,用于通过组合当前帧的血流速度和上一帧帧平均后的血流速度来生成查找表地址。根据查找表地址对存储器进行查表来获得滤波器的系数。
在本发明的一个实施例中,滤波器是IIR滤波器。
在本发明的一个实施例中,非迭代部分包括:第一阈值比较器,用于将当前帧的各个血流分量与对应的阈值进行比较,如果当前帧的各个血流分量均大于对应的阈值,则当前帧的血流速度有效,否则,当前帧的血流速度是噪声;第二阈值比较器,用于将当前帧的血流速度和上一帧帧平均后的血流速度进行比较,以确定是否是以下情况之一:上一帧帧平均后的血流速度发生混叠;当前帧的血流速度是闪烁噪声;第一组合器,用于将第一阈值比较器的输出结果和第二阈值比较器的输出结果进行组合;第一选择器,用于根据第一组合器的输出值来选择是输出当前帧的血流速度还是输出迭代计算部分的输出值,如果当前帧的血流速度有效,并且是情况之一,或者如果当前帧的血流速度是噪声,则输出迭代计算部分的输出值,否则,输出当前帧的血流速度;速度比较器,用于将当前帧的血流速度的绝对值和迭代计算部分的输出值的绝对值进行比较,如果当前帧的血流速度的绝对值大于迭代计算部分的输出值的绝对值,则当前帧的血流速度大于上一帧帧平均后的血流速度,否则,当前帧的血流速度小于上一帧帧平均后的血流速度;第二选择器,用于根据速度比较器的输出值来选择是输出当前帧的血流速度还是输出迭代计算部分的输出值,如果当前帧的血流速度大于上一帧帧平均后的血流速度,则输出迭代计算部分的输出值,否则,输出当前帧的血流速度;速度方向判断器,用于判断当前帧的血流速度的方向是否和上一帧帧平均后的血流速度的方向相同;和第三选择器,用于根据速度方向判断器的输出值来选择是输出第一选择器的输出值还是输出第二选择器的输出值,如果方向相同,则输出第二选择器的输出值,否则,输出第一选择器的输出值。
在本发明的一个实施例中,第一组合器是与门。
在本发明的一个实施例中,第一阈值比较器包括:速度阈值比较器,用于将当前帧的血流速度与速度阈值进行比较;第一能量阈值比较器,用于将当前帧的能量与第一能量阈值进行比较;第二能量阈值比较器,用于将当前帧的壁滤波前的能量与第二能量阈值进行比较;和第二组合器,用于将速度阈值比较器的输出值、第一能量阈值比较器的输出值和第二能量阈值比较器的输出值进行组合。
在本发明的一个实施例中,第二组合器是与门。
在本发明的一个实施例中,第二阈值比较器包括:速度差阈值比较器,用于将当前帧的血流速度和上一帧帧平均后的血流速度的差的绝对值与速度差阈值进行比较,如果当前帧的血流速度和上一帧帧平均后的血流速度的差的绝对值大于速度差阈值,则上一帧帧平均后的血流速度发生混叠;速度绝对值比较器,用于将当前帧的血流速度的绝对值和上一帧帧平均后的血流速度的绝对值进行比较,如果当前帧的血流速度的绝对值小于上一帧帧平均后的血流速度的绝对值,则当前帧的血流速度是闪烁噪声;和第三组合器,用于将速度差阈值比较器的输出值和速度绝对值比较器的输出值进行组合。
在本发明的一个实施例中,第三组合器是异或门。
与现有技术相比,本发明具有以下优点:
1.有效地消除血流图像中的闪烁点;
2.有效地减少回波数据中的噪声影响;
3.可以减少回波数据的混叠现象对图像的影响;
4.能够实时反映出血流的变化形态;以及
5.能够及时反映出速度方向改变。
附图说明
通过结合以下附图,并且参考以下对具体实施方式的详细说明,可以对本发明有更透彻的理解:
图1是医学超声成像系统的结构框图;
图2是根据本发明的帧平均电路的结构框图;
图3是根据本发明的迭代计算部分的电路图;
图4是根据本发明的非迭代部分的结构框图;
图5是根据本发明的第一阈值比较器的电路图;
图6是根据本发明的第二阈值比较器的电路图;以及
图7是根据本发明的速度比较器的电路图。
具体实施方式
根据本发明的帧平均电路主要由两部分组成。一部分为迭代计算部分,采用查找表方式,对当前帧的血流速度和上一帧帧平均后的血流速度进行IIR滤波。滤波系数由查找表决定,查找表的输入为当前帧的血流速度和上一帧帧平均后的血流速度。该部分能够很好地处理同方向血流速度发生变化的情况。另一部分为非迭代部分,采用了能量、壁滤波前的能量和速度阈值来对当前帧的血流速度进行判断。该部分主要在血流速度方向发生变化后对血流速度的输出进行控制。
采用根据本发明的帧平均电路对血流速度进行处理后,输出的血流速度能够正确反映真实血流的变化。在根据本发明的帧平均电路中,对血流速度的同方向变化和反向变化都进行了处理。在同方向变化中,采用IIR滤波器能够减弱血流信号中噪声的影响。当速度方向发生变化时,通过非迭代部分的阈值判断处理,不但能够实时反映出血流速度的快速变化,而且能够减弱由于其他原因而导致的速度变化对图像的影响。
以下将结合附图对本发明进行详细说明。
本发明主要针对在医学超声成像系统的彩色血流信号处理中的帧平均处理。如图2所示,根据本发明的帧平均电路主要由以下两部分组成:迭代计算部分420和非迭代部分440。当前帧的血流速度信号X(n)410分别进入迭代计算部分420的滤波器和非迭代部分440的选择器中。迭代计算部分420的滤波器有四个输入,其中两个是计算的数据,另外两个是系数。两个计算的数据分别为当前帧的血流速度值X(n)410和上一帧帧平均后的血流速度值460。上一帧帧平均后的血流速度数据存储在RAM中。在当前帧的血流速度数据X(n)410进入迭代计算部分420的滤波器时,从RAM中读出上一帧帧平均后的血流速度数据进行滤波计算。两个系数由迭代计算部分420的查找表RAM430提供,它们是当前帧的血流速度值和上一帧帧平均后的血流速度值的函数。非迭代部分440接收迭代计算部分的输出值480并进行相关阈值判断。当前帧的血流计算结果中的各个血流分量Z(n)470(包括能量、壁滤波前的能量和速度)进入到非迭代部分440的阈值比较器中。阈值比较器将输入的数据进行比较,输出比较结果。非迭代部分440的选择器根据比较结果并结合速度方向变化来控制帧平均电路的输出结果,即决定帧平均电路是输出当前帧的血流速度值还是输出迭代计算部分420的滤波器处理后的速度值(即迭代计算部分420的输出值)。延迟单元450通过插入一帧的延迟来允许使用上一帧帧平均后的血流速度。
如图3所示,迭代计算部分由IIR滤波器和查找表组成。当新的一帧数据到来后,系统控制模块中的硬件控制器根据新来帧的数据产生电影回放RAM的地址信号A(n),从电影回放RAM中读出上一帧帧平均后的血流速度值421。该血流速度值421经过阈值判断,送出处理后的血流速度值423。血流速度值423和经过寄存后的当前帧的血流速度值425组合,生成查找表地址424。根据查找表地址424从查找表RAM中找出相应的滤波器系数1-K和K,分别和经过寄存后的当前帧的血流速度数据426和上一帧帧平均后的血流速度数据422进行乘法运算,计算结果427和428再进行加法运算,得到迭代计算部分的输出值V1(n)。V1(n)经过寄存后进入非迭代部分,再进行相应的处理。通过该滤波器处理,可以减弱彩色血流数据中噪声和其他干扰信号的影响,保留血流数据中的有效成分。
非迭代部分原理如图4所示。硬件控制器通过节点444将各个阈值数据写入到第一阈值比较器441和第二阈值比较器442中。第一阈值比较器441主要是将当前输入的血流分量、如能量、速度和壁滤波前的能量等与对应的阈值进行比较,将输出结果445送入到与门逻辑中。第二阈值比较器442主要是将连续两帧回波数据中的血流速度进行比较,比较结果446也送入到与门逻辑中。与门逻辑输出控制信号来控制第一选择器451。速度比较器443将当前回波速度V(n)和滤波器的输出V1(n)进行比较,比较结果447控制第二选择器452。速度方向判断器454将当前回波速度的方向和上一帧帧平均后血流速度的方向进行比较,输出比较结果来控制第三选择器453。第三选择器最终输出帧平均处理后的结果Y(n)。以下将详细描述第一阈值比较器、第二阈值比较器和速度比较器。
第一阈值比较器的电路如图5所示。在速度方向发生改变后,第一阈值比较器对当前回波数据中的速度、能量和壁滤波前的能量进行阈值判断,以判断出当前回波数据的各个分量是否同时满足要求。如果满足要求,即如果当前回波数据中的速度V(n)、能量P(n)和壁滤波前的能量BW(n)分别大于速度阈值Vth、第一能量阈值Pth和第二能量阈值BWth,则说明发生的速度改变是有效的改变,否则认为是噪声的影响。在噪声影响的情况下,需要减弱噪声对图像的影响。
第二阈值比较器的电路如图6所示。在速度方向发生改变后,第二阈值比较器利用当前帧的血流速度值V(n)和上一帧帧平均后的血流速度值V(n-1)进行比较。第二阈值比较器包括速度差阈值比较器、速度绝对值比较器和与或门。速度差阈值比较器将当前帧的血流速度V(n)和上一帧帧平均后的血流速度V(n-1)的差的绝对值与速度差阈值Vrevth进行比较。如果当前帧的血流速度V(n)和上一帧帧平均后的血流速度V(n-1)的差的绝对值大于速度差阈值Vrevth,则上一帧帧平均后的血流速度发生混叠。速度绝对值比较器将当前帧的血流速度V(n)的绝对值和上一帧帧平均后的血流速度V(n-1)的绝对值进行比较。如果当前帧的血流速度V(n)的绝对值小于上一帧帧平均后的血流速度V(n-1)的绝对值,则当前帧的血流速度是闪烁噪声。与或门将速度差阈值比较器的输出值和速度绝对值比较器的输出值进行组合。通过速度差阈值比较器将连续两帧的血流速度值的差的绝对值与速度差阈值进行比较,可以判断出是上一帧帧平均后的血流速度发生混叠还是速度方向发生了改变。通过速度绝对值比较器,可以区分速度是混叠还是闪烁噪声。闪烁噪声的血流速度值与正常的血流速度值相比相对较小。通过将当前帧的血流速度与上一帧帧平均后对应点的血流速度进行比较可以减弱噪声对图像的影响。
速度比较器的电路如图7所示。速度比较器将当前帧的血流速度V(n)的绝对值和迭代计算部分的输出值V1(n)的绝对值进行比较。如果当前帧的血流速度V(n)的绝对值大于迭代计算部分的输出值V1(n)的绝对值,则当前帧的血流速度V(n)大于上一帧帧平均后的血流速度V(n-1),否则,当前帧的血流速度V(n)小于上一帧帧平均后的血流速度V(n-1)。显然,该电路主要是用在速度方向没有发生变化的情况下对当前帧的血流速度值进行判断,从而实现对低速血流的检测。如果当前回波速度V(n)大于滤波器的输出V1(n),则检测到的是加速血流,否则,检测到的是减速血流。
现在回到图4,如果当前帧的血流速度有效,并且如果当前帧的血流速度是闪烁噪声或者如果上一帧帧平均后的血流速度发生混叠,或者如果当前帧的血流速度是噪声,则第一选择器451输出迭代计算部分的输出值,否则,第一选择器451输出当前帧的血流速度;如果当前帧的血流速度大于上一帧帧平均后的血流速度,则第二选择器452输出迭代计算部分的输出值,否则,第二选择器452输出当前帧的血流速度;如果当前帧的血流速度的方向与上一帧帧平均后的血流速度的方向相同,则第三选择器453输出第二选择器的输出值,否则,第三选择器453输出第一选择器的输出值。因此,本发明的帧平均电路能够在保证减弱噪声影响的同时,较好地实时表现血流速度的变化。本发明的帧平均电路对血流速度方向的变化和不变化分开进行处理。血流速度方向没有发生变化时,采用一套专门的处理电路进行输出处理,以减弱噪声的影响。当速度方向发生变化后,采用另一套专门的处理电路,以判断速度方向变化是否是有效的变化,从而较好地实时表现血流速度的变化。
综上所述,本发明提供一种用于在医学超声成像系统中实现彩色血流帧的帧平均处理的帧平均电路。该帧平均电路由两部分组成,迭代计算部分和非迭代部分。迭代计算部分采用IIR滤波器和查找表相结合,主要处理速度方向没有发生变化的情况。非迭代部分对彩色血流回波信号的速度、能量和壁滤波前的能量进行阈值判断,实现对速度方向发生变化的处理。通过本发明,能够有效地减弱回波数据中的闪烁点现象,减弱彩色血流图像中的噪声影响;能够判断血流速度的变化,较好地实时表现血流速度的变化;能够判断血流速度的混叠,表现真实的血流图像。
本发明已经应用于商用彩超中,通过对彩色血流回波数据进行后处理,使显示的彩色血流图像能够更加逼真地反映血管中血流的变化,并且能够保留有效的彩色血流回波速度信号,消除血流图像中的闪烁点,减少血流图像中的自由噪声,同时还能在血流速度出现混叠时,仍然能够正确表现血流速度的形态,从而提高了血流图像的质量。
以上通过特定的实施例对本发明进行了详细的描述,但本发明并不限于上述实施例。在不脱离本发明范围的前提下,可以对本发明进行各种修改和变更。本发明的范围由所附权利要求书限定。在本发明的说明书和权利要求书中所用的诸如“第一”、“第二”和“第三”等词语仅仅是为了便于描述,而不具有任何限制意义。

Claims (10)

1. 一种用于超声成像系统的帧平均电路,包括:
迭代计算部分,用于对当前帧的血流速度和上一帧帧平均后的血流速度进行滤波;和
非迭代部分,用于进行关于所述当前帧的血流速度的判断并根据判断结果来输出所述当前帧的血流速度或所述迭代计算部分的输出值作为当前帧帧平均后的血流速度。
2. 如权利要求1所述的电路,其特征在于,所述迭代计算部分包括:
滤波器;和
系数提供单元,用于向所述滤波器提供系数。
3. 如权利要求2所述的电路,其特征在于,所述系数提供单元包括:
存储器,用于存储查找表,所述查找表用于存储所述滤波器的系数;和
地址生成器,用于通过组合所述当前帧的血流速度和所述上一帧帧平均后的血流速度来生成查找表地址,
根据所述查找表地址对所述存储器进行查表来获得所述滤波器的系数。
4. 如权利要求2或3所述的电路,其特征在于,所述滤波器是IIR滤波器。
5. 如权利要求1至3中任一项所述的电路,其特征在于,所述非迭代部分包括:
第一阈值比较器,用于将当前帧的各个血流分量与对应的阈值进行比较,如果所述当前帧的各个血流分量均大于所述对应的阈值,则所述当前帧的血流速度有效,否则,所述当前帧的血流速度是噪声;
第二阈值比较器,用于将所述当前帧的血流速度和所述上一帧帧平均后的血流速度进行比较,以确定是否是以下情况之一:
所述上一帧帧平均后的血流速度发生混叠;
所述当前帧的血流速度是闪烁噪声;
第一组合器,用于将所述第一阈值比较器的输出结果和所述第二阈值比较器的输出结果进行组合;
第一选择器,用于根据所述第一组合器的输出值来选择是输出所述当前帧的血流速度还是输出所述迭代计算部分的输出值,如果所述当前帧的血流速度有效,并且是所述情况之一,或者如果所述当前帧的血流速度是噪声,则输出所述迭代计算部分的输出值,否则,输出所述当前帧的血流速度;
速度比较器,用于将所述当前帧的血流速度的绝对值和所述迭代计算部分的输出值的绝对值进行比较,如果所述当前帧的血流速度的绝对值大于所述迭代计算部分的输出值的绝对值,则所述当前帧的血流速度大于所述上一帧帧平均后的血流速度,否则,所述当前帧的血流速度小于所述上一帧帧平均后的血流速度;
第二选择器,用于根据所述速度比较器的输出值来选择是输出所述当前帧的血流速度还是输出所述迭代计算部分的输出值,如果所述当前帧的血流速度大于所述上一帧帧平均后的血流速度,则输出所述迭代计算部分的输出值,否则,输出所述当前帧的血流速度;
速度方向判断器,用于判断所述当前帧的血流速度的方向是否和所述上一帧帧平均后的血流速度的方向相同;和
第三选择器,用于根据所述速度方向判断器的输出值来选择是输出所述第一选择器的输出值还是输出所述第二选择器的输出值,如果方向相同,则输出所述第二选择器的输出值,否则,输出所述第一选择器的输出值。
6. 如权利要求5所述的电路,其特征在于,所述第一组合器是与门。
7. 如权利要求5所述的电路,其特征在于,所述第一阈值比较器包括:
速度阈值比较器,用于将所述当前帧的血流速度与速度阈值进行比较;
第一能量阈值比较器,用于将当前帧的能量与第一能量阈值进行比较;
第二能量阈值比较器,用于将当前帧的壁滤波前的能量与第二能量阈值进行比较;和
第二组合器,用于将所述速度阈值比较器的输出值、所述第一能量阈值比较器的输出值和所述第二能量阈值比较器的输出值进行组合。
8. 如权利要求7所述的电路,其特征在于,所述第二组合器是与门。
9. 如权利要求5所述的电路,其特征在于,所述第二阈值比较器包括:
速度差阈值比较器,用于将所述当前帧的血流速度和所述上一帧帧平均后的血流速度的差的绝对值与速度差阈值进行比较,如果所述当前帧的血流速度和所述上一帧帧平均后的血流速度的差的绝对值大于所述速度差阈值,则所述上一帧帧平均后的血流速度发生混叠;
速度绝对值比较器,用于将所述当前帧的血流速度的绝对值和所述上一帧帧平均后的血流速度的绝对值进行比较,如果所述当前帧的血流速度的绝对值小于所述上一帧帧平均后的血流速度的绝对值,则所述当前帧的血流速度是闪烁噪声;和
第三组合器,用于将所述速度差阈值比较器的输出值和所述速度绝对值比较器的输出值进行组合。
10. 如权利要求9所述的电路,其特征在于,所述第三组合器是异或门。
CN200710166607A 2007-10-22 2007-10-22 用于超声成像系统的帧平均电路 Active CN101416885B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN200710166607A CN101416885B (zh) 2007-10-22 2007-10-22 用于超声成像系统的帧平均电路
US12/255,570 US8282560B2 (en) 2007-10-22 2008-10-21 Frame averaging circuit for use in an ultrasound imaging system
US13/647,186 US8721549B2 (en) 2007-10-22 2012-10-08 Frame averaging circuit for use in an ultrasound imaging system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200710166607A CN101416885B (zh) 2007-10-22 2007-10-22 用于超声成像系统的帧平均电路

Publications (2)

Publication Number Publication Date
CN101416885A CN101416885A (zh) 2009-04-29
CN101416885B true CN101416885B (zh) 2012-10-17

Family

ID=40564155

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200710166607A Active CN101416885B (zh) 2007-10-22 2007-10-22 用于超声成像系统的帧平均电路

Country Status (2)

Country Link
US (2) US8282560B2 (zh)
CN (1) CN101416885B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101416885B (zh) 2007-10-22 2012-10-17 深圳迈瑞生物医疗电子股份有限公司 用于超声成像系统的帧平均电路
CN102113900B (zh) * 2010-01-05 2015-07-15 深圳迈瑞生物医疗电子股份有限公司 彩色血流动态帧相关方法和装置
CN105686851B (zh) * 2016-01-14 2018-10-09 深圳开立生物医疗科技股份有限公司 一种血流3d成像方法、装置及其超声设备
CN105640592B (zh) * 2016-03-18 2019-01-01 深圳安盛生物医疗技术有限公司 一种彩色血流帧相关计算方法
CN106037793B (zh) * 2016-05-11 2020-06-02 飞依诺科技(苏州)有限公司 超声成像系统及其帧平均处理方法
CN109363722B (zh) * 2018-11-29 2021-08-06 武汉中旗生物医疗电子有限公司 彩色血流成像中运动伪像的抑制方法及设备
US11506771B2 (en) * 2019-09-24 2022-11-22 GE Precision Healthcare LLC System and methods for flash suppression in ultrasound imaging
CN111265246B (zh) * 2020-02-17 2023-02-24 东软医疗系统股份有限公司 超声彩色成像处理方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5595179A (en) * 1995-05-02 1997-01-21 Acuson Corporation Adaptive persistence processing
CN1208599A (zh) * 1997-05-07 1999-02-24 通用电气公司 彩色流速超声波成象中提高分辨率和灵敏度的方法和装置
US5897502A (en) * 1996-11-26 1999-04-27 Siemens Medical Systems, Inc. Persistence for ultrasonic flow imaging
EP1067400A2 (en) * 1999-07-09 2001-01-10 General Electric Company Method and apparatus for adaptive wall filtering in spectral doppler ultrasound imaging
CN1518669A (zh) * 2001-06-22 2004-08-04 �ʼҷ����ֵ������޹�˾ 超声杂波滤波器

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5467770A (en) 1994-11-25 1995-11-21 General Electric Company Color adaptive frame averaging
JPH10511588A (ja) * 1994-12-30 1998-11-10 アキュソン コーポレイション 体液の流れ又は組織の運動の画像を増強するための適応式時間フィルタリング
US6086539A (en) * 1996-12-04 2000-07-11 Acuson Corporation Methods and apparatus for ultrasound image quantification
US5823964A (en) * 1996-12-30 1998-10-20 Siemens Medical Systems, Inc. Ultrasound doppler wall filter
US5895358A (en) 1997-05-07 1999-04-20 General Electric Company Method and apparatus for mapping color flow velocity data into display intensities
US6123672A (en) * 1998-12-31 2000-09-26 General Electric Company Color flow imaging for enhancing segmentation and flow dynamics
US6322505B1 (en) * 1999-06-08 2001-11-27 Acuson Corporation Medical diagnostic ultrasound system and method for post processing
US6870885B2 (en) * 2001-05-16 2005-03-22 Qualcomm Incorporated Apparatus and method for decoding and computing a discrete cosine transform using a butterfly processor
US7139611B1 (en) * 2003-05-05 2006-11-21 Pacesetter, Inc. System and method for rejecting far-field signals using an implantable cardiac stimulation device
CN101416885B (zh) 2007-10-22 2012-10-17 深圳迈瑞生物医疗电子股份有限公司 用于超声成像系统的帧平均电路

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5595179A (en) * 1995-05-02 1997-01-21 Acuson Corporation Adaptive persistence processing
US5897502A (en) * 1996-11-26 1999-04-27 Siemens Medical Systems, Inc. Persistence for ultrasonic flow imaging
CN1208599A (zh) * 1997-05-07 1999-02-24 通用电气公司 彩色流速超声波成象中提高分辨率和灵敏度的方法和装置
EP1067400A2 (en) * 1999-07-09 2001-01-10 General Electric Company Method and apparatus for adaptive wall filtering in spectral doppler ultrasound imaging
CN1518669A (zh) * 2001-06-22 2004-08-04 �ʼҷ����ֵ������޹�˾ 超声杂波滤波器

Also Published As

Publication number Publication date
US20090105593A1 (en) 2009-04-23
US8282560B2 (en) 2012-10-09
CN101416885A (zh) 2009-04-29
US8721549B2 (en) 2014-05-13
US20130046180A1 (en) 2013-02-21

Similar Documents

Publication Publication Date Title
CN101416885B (zh) 用于超声成像系统的帧平均电路
US5409007A (en) Filter to reduce speckle artifact in ultrasound imaging
US5282471A (en) Ultrasonic imaging system capable of displaying 3-dimensional angiogram in real time mode
US6406430B1 (en) Ultrasound image display by combining enhanced flow imaging in B-mode and color flow mode
US5860930A (en) Energy weighted parameter spatial/temporal filter
CN100544677C (zh) 处理多普勒信号间隙的方法
DE10058452B4 (de) Verfahren und Gerät zur Bewegungsdarstellung bei Ultraschall-Fluss-Bilddarstellung unter Verwendung von Paketdatenerfassung
CN105590315B (zh) 一种超声彩色多普勒图像后处理方法
EP2245988B1 (en) Clutter signal filtering using eigenvectors in an ultrasound system
US4386528A (en) Method and means for generating pixel data in an ultrasonic scanning system
JPH04317639A (ja) 超音波診断装置
CN101866480A (zh) 一种空间复合成像中的帧相关处理方法及系统
CN110832343A (zh) 具有改善的杂波抑制的能量多普勒成像系统和方法
US5897502A (en) Persistence for ultrasonic flow imaging
US5865752A (en) Method and apparatus for ultrasound imaging using normalized difference between successive frames
US5976086A (en) Ultrasonic diagnostic apparatus and filter to obtain sequential frames of image data
US20040225221A1 (en) Diagnostic ultrasound imaging system with adaptive persistence
US6893399B2 (en) Method and apparatus for B-mode image banding suppression
JP2763126B2 (ja) カラー超音波診断装置
JPH04297246A (ja) 超音波診断装置
JP4426472B2 (ja) 超音波診断装置
JP4749540B2 (ja) 超音波カラーフローの適応的スキャン方法及びその装置
US11234671B2 (en) Ultrasonic diagnostic apparatus
JP2950600B2 (ja) カラードプラ超音波診断装置
JPH03277357A (ja) 超音波診断装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20180515

Address after: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong

Co-patentee after: SHENZHEN MINDRAY SCIENTIFIC Co.,Ltd.

Patentee after: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Address before: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong

Patentee before: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

TR01 Transfer of patent right
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20090429

Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd.

Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Contract record no.: X2022440020009

Denomination of invention: Frame Averaging Circuits for Ultrasound Imaging Systems

Granted publication date: 20121017

License type: Common License

Record date: 20220804

EE01 Entry into force of recordation of patent licensing contract