专利名称:在超声成像系统中用于优化超声图像的灰度值的方法
技术领域:
本发明涉及一种在超声成像系统中用于优化超声图像的灰度值 的方法。
背景技术:
通常在运用B型成像进行超声诊断时,诊断人员需要调整TGC (时间增益补偿)、Gain (总增益)、动态范围等,使图像达到一 个比较好的水平,这些调整操作会导致完成诊断的时间加长。这一 问题往往通过"B型图像一键优化"技术来解决。B型图像(即黑白 图像) 一键优化技术的目的就是简化TGC、 Gain、动态范围的调 节,用户只需要按一个键,由超声系统自动分析出最合适的TGC、 Gain及其他和成像相关的参数并使其生效,图像因而达到比较好的 水平,可以有效加快诊断的效率和准确率。这种工作方式称为一键 优化,或者自动图像均衡。
绝大多数一键优化(或者图像均衡)技术都是基于分析图像沿 深度方向的灰度变化情况,以此为基础计算合适的TGC曲线和 Gain,例如美国专利US 6743174B2以及美国专利US 5579768、 US 6120446和US 6102859。分析灰度变化的方法各不相同,另外分析中 都考虑了对噪声的抑制,其中US 6743174B2还提供了一个动态范围 自动优化的方法。
专利US 6743174B2假定灰度沿深度是单调均匀变化的,可以用
一条直线来拟合,很多时候这和实际情况不符,另外它利用灰度均 值和中值来调整动态范围的方式,也不太容易达到一个好的效果;
US 5579768和US 6102859把B型图像沿探头径向和切向均匀地分成 若干子区域,在分析灰度变化时,主要考虑了每一子区域的平均灰 度,可能导致灰度变化分析不够准确;另外这种均匀地分隔子区域 方式,也不利于灰度变化分析。
发明内容
本发明的目的是提供一种可以克服上述现有技术的缺陷并获得
更好的图像均衡效果的超声成像系统和成像方法。
本发明提供一种在超声成像系统中用于优化超声图像的灰度值
的方法,包括以下步骤
获得响应于预定时间增益补偿参数产生的超声图像; 从该超声图像中选择不同深度下代表均匀软组织的子区域; 分析这些代表均匀软组织的子区域的灰度沿深度的变化;以及 基于该分析的结果计算该超声图像的优化灰度值。 本发明还提供一种可优化超声图像的灰度值的超声成像系统,
包括
获得响应于预定时间增益补偿参数产生的超声图像的装置;从 该超声图像中选择不同深度下代表均匀软组织的子区域的装置;分 析这些代表均匀软组织的子区域的灰度沿深度的变化的装置;以及 基于该分析的结果计算该超声图像的优化灰度值的装置。
图l是本发明的超声成像系统的实施例的结构方框图; 图2是本发明的在超声成像系统中用于优化超声图像的灰度值的 方法的实施例的流程图3是在根据本发明的实施例中分析灰度变化的流程图; 图4是在根据本发明的实施例中子区域移动方法的示意图;和
图5是在根据本发明的实施例中图像均衡参数的计算流程示意图。
具体实施例方式
图l是本发明的可实现自动图像均衡的超声成像系统的实施例 的结构方框图(其中由于发射部分和本发明无太大关系,因此未示 出发射部分)。
一个常规的B型成像流程为探头发射出脉冲,参与接收的各
阵元接收到回波,经过放大,A/D转换,在波束合成器以不同的延 时量相加得到射频数据。射频数据经过检包络和降采样,DTGC
(Digital time gain compensation),对数压縮,动态范围变化,数 据进入电影回放存储器,以便用户可以回放数据,同时进入DSC
(数字扫描转换),DSC结果显示在屏幕上,即B型图像。上述过 程中,DTGC是指给输入的数据沿深度(即时间)给予不同的增 益,目的是使显示在屏幕上的图像亮度比较均匀, 一般说来是用户 根据显示在屏幕上的图像,拨动键盘上的电位器(一般有六到八 个,对应不同深度的增益补偿值)对图像不同的深度部分调亮或者 变暗,或者旋转键盘上的Gain旋钮使图像整体变暗/亮,主控CPU读 取这些电位器和Gain旋钮的状态,换算成对应DTGC参数值(一般 会经过插值),并将该值写到DTGC模块,DTGC模块用写入的 DTGC参数值和解调后的结果相乘。
以一个具体的系统为例,假定键盘上有八个电位器,八个电位 器数值对应dtgcjo[i], i=0~7, Gain旋钮对应数值gain,检波降采样 后的一条扫描线数据为x[k],k-0 511, DTGC模块的工作方式为
1)主控CPU读入电位器的数值dtgcjo[i], i=0 7,对这8个数进 行线性插值,变成长度为128的DTGCl[i], i=0~127。
2)主控CPU读入Gain旋钮数值gain,将DTGCl和Gain旋钮对应 数值gain相加(或者相乘,不同系统略有区别),得到长度为128的 PDM[i], i=0~127,并将PDM的值提供给DTGC模块。其中,PDM是 数字增益补偿参数(Parameter for DTGC Module)。
DTGC模块将输入的降采样结果x[k], k=0 511和PDM[i〗, M) 127相乘,假定输出的数据为y[k],k-0 511 。
<formula>formula see original document page 9</formula>
上式中delta表示每delta个输入点和同一个PDM相乘, 一般的, delta可以取4。
如果键盘上的电位器或者Gain旋钮状态发生了变化,重复1 )~2) 本发明的超声成像系统就是自动寻找合适的PDM,替换由键盘 上的TGC电位器及Gain决定的PDM。参见图1,在非一键优化状态 下,主控CPU根据键盘的电位器及Gain旋钮状态按上述方法计算 PDM并送入硬件中的DTGC模块;在一键优化状态下,主控CPU从 电影回放区读取数据,通过分析这些数据,得到合适的PDM,送入 硬件中的DTGC模块。
图2为是本发明的在超声成像系统中用于优化超声图像的灰度 值的方法的实施例的流程图。当用户启动本发明的方法时,其工作 流程如下
1. 关闭发射主控CPU首先使发射停止
2. 设置默认参数然后将DTGC模块和动态范围变化模块的参 数设为默认,DTGC模块的默认参数为使输入增益为1的数,动态范 围参数使输入数据由最小到最大线性映射到0 2553. 分析噪声水平..从电影回放区采集一帧或者多帧数据,由于 此时发射关闭,无回波信号,所以电影回放区内为噪声数据。对采 集到的噪声数据进行分析,得到对应图像不同深度的噪声水平。
4. 打开发射主控CPU使发射开启
5. 采集数据从电影回放区采集一帧或者多帧数据,由于此时 发射开启,电影回放存储器内为被扫描物反射回波经过处理的数据。
6. 分析使图像自动均衡的PDM:主控CPU对采集得到的数据进 行分析,得到使显示图像整场亮度均匀,且平均亮度符合用户要求 的PDM。
7. 送参数PDM到DTGC模块将上一步分析得到的PDM写到 DTGC模块,这时显示器中显示的自然就是亮度优化后的图像。
本发明的方法实施例中的相关算法可分为三步
分析噪声水平
分析图像灰度变化曲
分析得到PDM
以下分别介绍这三步
一、噪声水平的分析
采集到 一 帧的噪声数据可以看成 一 个二维数组 Noiselmg[M][N], M为纵向的点数, 一般在360 600之间,N为一帧 数据的接收线数目,分析噪声水平按如下步骤进行
(1) 将噪声图像平均分成Ver—Num行,每行的数据量为 M/Ver—Num*N; Ver—Num的值可以根据探头,检査模式的不同设成 不同的值。
(2) 计算每一行所有数据的均值,得到噪声水平曲线 Noise[i], i-0~Ver—Num-l ,代表图象显示区域由浅到深的噪声水 平;噪声水平曲线的计算公式为
Noise[i〗=Mean(NoiseImg[i*M/Ver—Num: (i+l) *M/Ver—Num-l][O:N-l])
(3)如果采集了多帧噪声图像,可以重复l, 2,然后将对应于 每一帧噪声图像的噪声水平曲线Noise[i]平均,得到最终的噪声水平 曲线Noise[i]。
二、分析灰度变化曲线
图3为分析灰度变化的基本流程。将从电影回放区读入的图像 数据定义成一二维数组Image[M][N], M为纵向的点数, 一般在 360~600之间,N为一帧数据的接收线数目,分析过程为
1 、将图像平均分成Ver—Num行,每行的数据量为 M/Ver—Num*N; Ver—Num的值可以根据探头,检查模式的不同设成 不同的值。注意,图像划分的行数应该和噪声图像划分的行数相 同。
2、确认每一行的类型。这一步为算法中的主要部分。我们认为 主要有三种类型 一种是以均匀软组织为主,称为有效行; 一种以
组织边缘为主,称为结构行; 一种以噪声为主,称为噪声行。确认 每一行类型的方法如下
a) 设定每一行内子区域的大小为IFW (H为纵向点数,W为 横向点数),具体值可以根据探头和检查模式的不同而不同,H的 值应该小于等于M/Ver_Num ,为了方便叙述,这里使 H=M/Vcr—N亂
b) 在一行内寻找代表软组织的子区域,这些子区域的相互 位置可以不均匀排列,以第i行为例,寻找方法为
(i)从行起始处开始寻找代表均匀软组织的子区域,用 Sublmg[m][n]表示子区域的元素,m取值范围为[(i-l" M/Ver—Num i* M/Ver—Num-1] , n取值范围[startP, startP+W-1],对于行起始, startP=0。从三个方面确认是否软组织区域均值大于同行的噪声水
平A倍(A可以根据探头和检查模式的不同设为不同值, 一般可以取
1.1);方差小于一个设定的方差值.,均值小于整幅图像灰度最大值 B倍(B可以取0.8)。如果均值小于同行的噪声水平A倍,该子区域 判断为噪声子区域;如果既非均匀软组织,又非噪声子区域,该子 区域判断为结构子区域。
(ii) 如果该子区域被判断为均匀软组织或者噪声,则 StartP=StartP+W (相当于移动子区域),重复步骤(i),判断下一个 子区域;如果该子区域判断为结构子区域,则startP= startP+deltaW, deltaW是一个小于W的数, 一般取W的l/3或者一 半,可以根据探头,检査模式的不同设成不同的值,重复步骤(i), 判断该子区域的类型。这样的査找方式可以有效地避免结构区域对 灰度计算的影响(一般结构区域的平均灰度会明显比均匀组织区域 亮,引起对整行平均灰度水平的误判),又能最大程度地利用所有 代表均匀软组织的区域,而均匀划分子区域的方式难以做到这一 点。
图4为在均匀组织区域和结构区域移动子区域方法的示意图。
(iii) 按照步骤(ii)所述方法,将子区域由左向右移动到图像最右 端,设定两个计数器,分别记录均匀软组织子区域的数目和噪声子 区域的数目。需要说明的是,也可以縮小子区域的面积,让子区域 在行内基于二维最寻找均匀组织子区域,可以达到更好的效果。
(iv) 如果一行内均匀组织子区域的数目大于设定阈值T一Tis,该 行为有效行(以均匀软组织为主),将用于图像的灰度分析;如果 噪声子区域的数目大于设定阈值T—Noise,该行为噪声行;如果前两 个条件都不成立,该行为结构行。TJis和T!Noise可以根据探头和 检査模式的不同设为不同值。
c) 对整幅图像分成的Ver—Num行都重复b),得到每一行的类型。
2. 计算每一行的灰度值Amp [i], i=0 Ver—Num-l。对于有效 行,灰度值为所有代表均匀软组织的子区域的灰度均值,对于其它 行(噪声行,结构区域行),灰阶值根据相邻有效行的值进行插 值。如果无效行并未处于两有效行之间,取最邻近的有效行的灰阶 值。
3. 如果存在多帧图像数据可供分析,对其余帧的数据重复l到 3,每一帧都得到相应的灰度变化曲线,将所有的灰度变化曲线求平 均,为最终的灰度变化曲线Amp(i) , i= 0 Ver—Num-l 。
三、计算图像均衡参数PDM
图像均衡除了使全场组织的亮度基本相同,也要使亮度和理想 的平均灰度值为MeanGray (MeanGray可以设在系统预置中)接近, 在一般情况下,用户通过调节Gain旋钮达到这种效果。因此在计算 PDM时,除了考虑图像的均衡,也要考虑整体亮度和MeanGray接近。
图5为根据灰度变化曲线分析PDM的流程示意图,计算按如下 步骤进行
1、 首先计算灰度变化曲线Amp和平均灰度MeanGray的差值 deltaGray,计算公式为
cfe/toGraj(0 = Mea"G考_,' = 0 ~ Fer _ jV薩 (})
很明显,将Amp(i)和deltaGray(i)相加结果是MeanGray。设置 PDM的目的是使检波后的数据与PDM相乘后,在经过对数压縮及动 态范围变化后,产生大小为deltaGray(i)的灰度差。
2、 采集得到的电影回放数据是经过对数压缩和动态范围变换 的结果,而DTGC模块作用在对数压縮之前。利用deltaGray,进行 反动态范围变换及反对数变换,很容易得到一组PDMl(i), i=0 Ver_Num-l,将解调的数据和与其深度相对应的PDMl(i)相乘
后,进行对数压缩及动态范围变换,得到的结果必须使整个图像代
表软组织的部分灰度在MeanGray附近。
3、考虑在计算灰度变化曲线时确定的噪声行,对这些行的数据 应该考虑做出抑制,如果第j行为噪声行,和该行最接近的有效行为 第i行,则
吗)=磨%#—一 —力) (2)
Sup—K为一个大于O的数,命名为灰度抑制系数,具体数值可以 根据探头和检査模式的不同设为不同值。
如果和某一灰度抑制无效行最接近的有效行有两个,公式(2)中 的DTGC(i)可以换成两个有效行的DTGC对应的均值。
前3步得到的DTGC数据长度为Ver一Num, 一般情况下DTGC模 块需要的参数长度不会等于Ver一Num,可以将PDM1进行线性插 值,变成DTGC模块需要的长度。插值后的结果即PDM参数。
在本发明的上述实施例中提出一种灰度曲线分析方式,其中采 用了一种动态寻找均匀组织子区域的方法,并且在判断子区域类型 时考虑和噪声水平的差异、方差、和最大亮度的差异;另外按图像 划分成多行,判别每行为三类中哪一类以均匀组织为主,以结构 区域为主,以噪声为主。即,基于动态寻找均匀组织子区域,以此 为基础判别子区域所在行的类型,并在此基础上分别灰度变化曲 线。这样的方式更有利于灰度变化曲线的分析。本发明的实施例还 提出了根据灰度变化曲线分析DTGC模块使用参数PDM的方法。
与现有技术的方法相比,本发明采用了动态寻找代表软组织子 区域的方法;并且在确认子区域类型时不仅考虑平均灰度值,还考 虑与噪声水平和最大灰度的关系,以及灰度的方差;采用这样方式
可以有效地避免结构子区域、噪声子区域对灰度变化曲线分析的影 响,从而更准确地得到合适的PDM值,使图像均衡效果更好.
本发明的实施例提出了一种基于合理地分析图像灰度变化来计
算时间增益补偿(TGC)参数PDM的方法。该方法使用动态划分图 像子区域的方式,利用图像的亮度,方差特性,选出不同深度下代 表均匀软组织的子区域,(这些子区域在横向位置的间隔可能并不 均匀),分析这些代表均匀软组织子区域灰度沿深度的变化,并以 此为基础计算使图像整场达到均衡的TGC参数;同时找出主要属于 噪声的区域,在计算TGC时对这些区域考虑加以抑制。
本发明的实施例采用了一种动态寻找均匀组织子区域的方法, 并且在判断子区域类型时考虑和噪声水平的差异,方差,和最大亮 度的差异;另外按图像划分成多行,判别每行为三类中哪一类以均 匀组织为主,以结构区域为主,以噪声为主;在分析灰度变化曲线 时,只考虑均匀组织行,并对噪声行加以抑制。这样的方法相对于 以往发明,可以更准确地分析同灰度变化曲线,因而得到更好的图 像均衡结果。
在本发明的上述实施例的说明中,在一行内搜寻代表均匀软组 织的子区域是基于一维搜索,类似的概念可以扩展到二维,实现更 准确地搜索。另外,类似的概念可以扩展到三维或者四维超声成像 中,动态寻找属于均匀组织的块,用以判断空间不同位置的灰度分 布,以此为基础实现三维图像的图像均衡。
另外,本发明的超声成像系统的上述实施例主要依靠主控CPU 来实现,但是也可以用单独的数字信号处理器(DSP)或者数字式 硬件逻辑电路实现。
权利要求
1、一种在超声成像系统中用于优化超声图像的灰度值的方法,包括以下步骤获得响应于预定时间增益补偿参数产生的超声图像;从该超声图像中选择不同深度下代表均匀软组织的子区域;分析这些代表均匀软组织的子区域的灰度沿深度的变化;以及基于该分析的结果计算该超声图像的优化灰度值。
2、 根据权利要求l的方法,其中,所述选择步骤包括 将该超声图像分成若干行; 设定每一行内的子区域的大小;以及在每一行内采用动态搜寻的方式确定每一行内代表均匀软组织 的子区域。
3、 根据权利要求2的方法,其中,所述确定步骤包括 当以下三个条件满足时,判断子区域是代表均匀软组织的子区域子区域的灰度均值大于相同行的噪声水平的预定倍; 子区域的灰度方差小于一个设定的方差值;以及 子区域的灰度均值小于整幅图像的灰度最大值的预定倍。
4、 根据权利要求3的方法,其中,所述确定步骤还包括 当子区域的灰度均值小于相同行的噪声水平的预定倍时,判断该子区域为噪声子区域;以及当子区域既非均匀软组织子区域,又非噪声子区域,则判断该 子区域为结构子区域。
5、 根据权利要求3的方法,其中,所述噪声水平的预定倍数是 根据该超声成像系统的探头和检査模式的不同而设定。
6、 根据权利要求3的方法,其中,所述噪声水平的预定倍数是1.1。
7、 根据权利要求3的方法,其中,所述整幅图像的灰度最大值 的预定倍数是0.8。
8、 根据权利要求l的方法,还包括根据每一行内的灰度分布确定每一行所指示的内容是以均匀软 组织、或组织边缘、或噪声为主,从而将以均匀软组织为主的行定 义为有效行,将以组织边缘为主的行定义为结构行,将以噪声为主 的行定义为噪声行。
9、 根据权利要求8的方法,所述计算步骤包括 计算每一行的灰度值;其中,对于有效行,该灰度值为所有代表均匀软组织的子区域的灰度均值;对于噪声行和结构行,灰度值是根据相邻有效行的灰 度值进行插值而获得的;以及获得该超声图像的灰度变化曲线。
10、 根据权利要求9的方法,还包括 计算灰度变化曲线与平均灰度的差值; 利用所计算的差值,进行反动态范围变换和反对数变换; 对所确定的噪声行的图像数据执行灰度抑制处理;以及获得该超声图像的数字时间增益补偿参数。
11、 根据权利要求l的方法,还包括 由该超声成像系统采集多帧噪声图像; 将每一帧噪声图像分成若干行;计算每一行的所有灰度数据的均值,以获得噪声水平曲线; 将对应于每一帧噪声图像的噪声水平曲线求平均,得到最终的 噪声水平曲线。
12、 根据权利要求l的方法,其中,该超声图像包括多帧图像, 该方法还包括将每一帧图像分成若干行;确定每一行的类型,其中根据每一行内的灰度分布确定每一行 所指示的内容是以均匀软组织、或组织边缘、或噪声为主,从而将 以均匀软组织为主的行定义为有效行,将以组织边缘为主的行定义 为结构行,将以噪声为主的行定义为噪声行;计算每一行的灰度值;获得关于每一帧图像的灰度变化曲线;以及 将所有的灰度变化曲线求平均,获得最终的灰度变化曲线。
13、 根据权利要求l的方法,其中,所述超声图像是三维图像, 所述子区域是三维空间中的子块。
14、 根据权利要求l的方法,其中,所述超声图像是四维图像。
15、 一种可优化超声图像的灰度值的超声成像系统,包括 获得响应于预定时间增益补偿参数产生的超声图像的装置;从该超声图像中选择不同深度下代表均匀软组织的子区域的装置;分析这些代表均匀软组织的子区域的灰度沿深度的变化的装 置;以及基于该分析的结果计算该超声图像的优化灰度值的装置。
全文摘要
本发明涉及一种在超声成像系统中用于优化超声图像的灰度值的方法,包括获得响应于预定时间增益补偿参数产生的超声图像;从该超声图像中选择不同深度下代表均匀软组织的子区域;分析这些代表均匀软组织的子区域的灰度沿深度的变化;以及基于该分析的结果计算该超声图像的优化灰度值。通过本发明的实施,超声成像系统将获得更好的图像均衡效果。
文档编号A61B8/00GK101190135SQ20061016353
公开日2008年6月4日 申请日期2006年11月29日 优先权日2006年11月29日
发明者胡勤军, 剑 董, 勇 蒋 申请人:深圳迈瑞生物医疗电子股份有限公司