专利名称：一种提升超声图像质量的方法
技术领域：
本发明涉及超声成像技术领域，特别是超声成像技术中的数字图像处理技术，提升超声图像质量的方法。
背景技术：
一个典型的超声成像系统如图1所示，它包括一个超声波成像装置，主控制器控制整个系统的正常运转。在主控制器的控制下，探头发射出超声波，延时一定时间后接收从组织反射回来的超声波。波束形成器将各通道反射回波信号进行聚焦延迟、加权、求和，形成一条或多条扫描线，再经过滤波、检波等处理获得扫描回波的包络信号。回波包络送到数字扫描变换模块完成坐标变换，形成图像数据。所形成的图像数据再送到显示模块，叠加上字符后输出到显示器显示。
由于超声束的空间分辨力有限，组织的反射面不光滑，不同的反射源产生的回波信号相位也不一样，在波束形成的时候可能相叠加也可能相消，及电子电路的噪声等因素导致超声图像中存在许多斑点噪声。斑点噪声的存在将会掩盖超声图像中的部分有用信息，对医生的诊断造成干扰，影响医生的视觉感受甚至影响到医生的诊断信心。
有很多的方法抑制斑点噪声，例如多帧图像平均，扫描线之间的平滑，图像像素的邻域平均等。另外为了突出有用的组织结构特征，有的超声成像设备采用锐化图像的方式，使得超声图像中的组织结构、边缘信息得到加强。
上述常用的办法虽然抑制了噪声或增强了组织结构特征，但在抑制了斑点噪声的同时，却模糊了组织结构或损失了动态特性，或在增强了组织结构信息的同时，却将噪声也一起加强了，无法真正提高超声图像质量。
与解决上述问题相关的是，GE公司在美国专利US 6,208,763和US 6,592,523中公开了一种基于梯度分割的图像增强方法。该方法依据梯度信息，将图像分割为结构区域与非结构区域。对结构区域进行基于灰度加权方向性二阶倒数的各向异性锐化处理，以增强图像边缘的对比度；对非结构区域进行各向同性的平滑处理以消除斑点噪声。
深圳迈瑞生物医疗电子股份有限公司在中国专利200510035913.8和200510036841.9中公开了一种基于梯度分割的图像增强方法。该方法依据梯度信息，将图像分割为结构区域与非结构区域。对结构区域进行基于中间像素的加权平均各向异性增强处理；对非结构区域进行各向同性的平滑处理以消除斑点噪声。
发明内容
上述方法对特征区域处理时，不能建立有效的数学模型增强图像的边缘信息和结构信息。本发明要解决的就是上述方法的不足，即利用热扩散方程模型，达到在抑制斑点噪声的同时，不损失反而增强超声图像中组织结构信息的目的，从而提升超声诊断图像质量。
本发明的目的在于克服现有技术的缺点，提供一种提升超声图像质量的方法。
为了在抑制斑点噪声的同时增强组织结构特征，本发明的基本思想在现有超声图像基础上，通过后处理的方式，利用由结构张量表示的图像像素的局部特征，将图像划分成组织结构区域，不含组织结构的平坦区域以及组织结构的边缘区域。对三类区域分别采用不同的处理方式，从而达到抑制斑点噪声的同时增强组织结构信息。为了抑制各种大小范围的斑点噪声，增强不同尺度的组织结构特征，所述的图像处理分别在多个尺度范围内进行。
为实现上述目的，本发明提出了一种提升超声图像质量方法，该方法增强超声图像结构信息与抑制斑点噪声，用于超声成像系统对扫描图像显示数据的优化处理。其的步骤如下 A.利用高斯-拉普拉斯金字塔分解算法分解超声显示图像； B.对分解产生的高斯层图像数据进行各向异性扩散滤波； C.对扩散结果进行锐化处理； D.利用高斯-拉普拉斯金字塔重构算法将各向异性扩散滤波结果与拉普拉斯层数据重构回去； E.A、B、C、D三个步骤迭代若干次。
金字塔分解可以获得不同分辨尺度下原始图像的表达，金字塔重构可以把不同分辨尺度下的图像几乎没有任何信息丢失重构回去。本发明正是基于这种思想，将图像分解到不同的尺度，然后在不同的尺度下对图像做扩散滤波处理，将不同尺度下处理过的图像重构回去，从而获得最终的结果。
各向异性扩散是基于热传导模型，实现对图像的平滑的同时可以保持甚至增强图像的边界或结构信息。采用高斯核与输入图像做卷积的图像滤波方法，等同于扩散系数为常值的热传导方程的解，即各向同性扩散。由于这种基于高斯核的图像滤波方法没有考虑图像特征在空间上的分布，因此在平滑图像的同时也滤除了图像的一些重要的特征信号。本发明考虑到各向同性扩散滤波的缺点，采用各向异性扩散滤波，根据图像的局部特征计算扩散系数。图像的局部特征用结构张量的特征向量和特征值来表征。特征值间的差值越大，表明这个区域的结构特征越强，大特征值方向采用较小的扩散系数，而小特征值方向采用较大的扩散系数；反之如果特征值间的差值越小，表明这个区域有结构信息的可能性就越小，即越平坦，两个特征方向的扩散系数基本相同。如此一来，扩散结果就是，有结构信息的地方沿着结构边缘扩散，而垂直于结构的边缘方向则几乎没有扩散，从而实现了结构信息的保留甚至增强。对于扩散方程的求解，本发明采用托马斯求解方法，以最高的效率完成扩散滤波算法，以便于实现实时处理超声扫描显示图像。
上述方案中，所述步骤B，对由输入图像分解出来的高斯层进行4个方向的各向异性扩散滤波，由以下几个步骤组成 a.计算4个方向扩散系数，4个方向分别是直角坐标系中的-45°、0°、45°和90°方向； b.在上述4个方向进行扩散滤波。
所述过程a，由以下几个步骤组成 1.高斯平滑； 2.计算梯度； 3.计算结构张量； 4.计算扩散张量； 5.裂解为4个方向的扩散系数。
有益效果 本发明达到在抑制斑点噪声的同时，不损失反而增强超声图像中组织结构信息的目的，从而提升超声诊断图像质量。
图1是典型的超声成像系统图。
图2是超声成像系统中本发明的图像处理模块图。
图3是提升超声图像质量的方法总体流程图。
图4是高斯-拉普拉斯金字塔分解流程图。
图5是高斯-拉普拉斯金字塔重构流程图。
图6是各向异性扩散滤波算法流程图。
图7是进行扩散的4个方向示意图。
具体实施例方式 图1是典型的超声成像系统图。图中由主控制器、探头、波束形成器、数字扫描变换(DSC，Digital Scan Conversion)、显示器组成，在主控制器的控制下，探头发射出超声波，延时一定时间后接收从组织反射回来的超声波。波束形成器将各通道反射回波信号进行聚焦延迟、加权、求和，形成一条或多条扫描线，再经过滤波、检波等处理获得扫描回波的包络信号。回波包络送到数字扫描变换模块完成坐标变换，形成图像数据。所形成的图像数据再送到显示模块，叠加上字符后输出到显示器显示。
图2是超声成像系统中本发明的图像处理模块图。
本发明提升超声图像质量的方法，抑制超声图像中的斑点噪声，增强超声图像中组织结构信息，用于在超声成像系统中对超声扫描图像显示数据的优化处理，其在系统中的位置如图2所示。是在DSC与显示器之间加入一个图像处理模块，图中虚线所示。该图像处理模块是提升超声图像质量方法的核心，主要包括图像数据输入接口，处理参数输入接口，图像显示数据输出接口，高斯-拉普拉斯金字塔分解重构算法，结构张量求解算法，扩散张量求解算法，各向异性扩散滤波算法，加性裂解算法，数字图像锐化算法。
图3是提升超声图像质量的方法总体流程图。
本发明提升超声图像质量的方法，其总体实现如图3所示。首先把将要送出显示的图像数据利用高斯-拉普拉斯金字塔分解算法，将其分解成若干层。接着对数据量最小的高斯层进行各向异性扩散滤波处理，再对处理结果做锐化处理。将处理过的高斯层数据与未处理的拉普拉斯层数据利用高斯-拉普拉斯重构算法进行重构。以上四个步骤经过若干次循环即可实现提升超声图像质量的目标，图像斑点噪声得到有效抑制，同时组织结构信息得到增强。
图4、图5是高斯-拉普拉斯金字塔分解、重构流程图。
本发明为了减少整个算法的计算量，以及实现在多尺度下对图像进行处理，对输入图像进行了高斯-拉普拉斯金字塔分解，将图像中的能量集中于数据量较少的高斯层数据。所述高斯-拉普拉斯金字塔分解与重构算法，为PETER J.BURT等在1983年提出来的。分解算法的流程如图4所示，重构算法的流程如图5所示。
高斯金字塔或拉普拉斯金字塔的分解和重构均由两个重要的算子组成，分别为Reduce和Expand，Reduce是数据经过数字滤波器后再下采样的过程，Expand与Reduce相反，是对数据上采样再滤波的过程。
分解过程中，高斯金字塔的求法 Gauss0＝x； Gauss1＝Reduce(Gauss1-1)。
而拉普拉斯金字塔的求法 Laplace1＝Gauss1-Expand(Gauss1+1) 重构算法 Gauss1＝Expand(Gauss1+1)+Laplace1 其中x为输入阵列，Gauss1为第1层高斯层，Laplace1为第L层拉普拉斯层。
本发明在金字塔分解算法中的滤波器，采用若干阶FIR滤波器。下采样前的滤波器和上采样后的滤波器具有特定的关系。
图6、图7是各向异性扩散滤波算法、进行扩散流程图。
本发明通过高斯-拉普拉斯金字塔分解算法对输入图像分解后，对最底层的高斯层进行各向异性扩散滤波，目的之一是减少计算量，目的之二是不同层数的金字塔分解，其最底层高斯层尺度也不同。因而通过多次循环可以实现在不同尺度下处理图像。本发明采用四个方向扩散滤波，分别是直角坐标系中的-45°、0°、45°和90°方向。而要进行各向异性扩散滤波，必须求出相应每个像素的4个方向的扩散系数。
本发明提升超声图像质量的方法，为了抑制超声图像中的噪声，增强组织结构信息，在计算扩散系数的时候将图像的局部特征考虑进去。本发明采用结构张量来表示局部的特征，某个像素的结构张量为其中，Ix为x方向的梯度，Iy为y方向的梯度。结构张量可以分解为两个正交的向量。利用这两个向量的方向和特征值(u1，u2)来表征局部特性。两个特征值的差值越大，说明这个区域的结构特征越强，大特征值方向的扩散小，小特征值方向的扩散大，反之如果两个特征值的差值越小，说明这个地方越没有结构信息，即越平坦，两个特征方向的扩散可基本相同。如此一来，扩散结果就是，有结构信息的地方沿着结构边缘有扩散，而垂直于结构的边缘方向则几乎没有扩散，从而实现了结构信息的保留甚至增强。
结构张量的计算要考虑到噪声的影响，因此在计算梯度的时候，首先采用若干阶FIR滤波器平滑图像，平滑之后的图像噪声得到抑制。基于平滑之后的图像求出每个像素的x，y方向的梯度，获得结构张量的组成部分，即
IxIy，并分别对这三个数据组成的阵列再次作平滑处理，进一步减小噪声对整个算法的影响，得到最终需要的结构张量 结构张量的特征值和特征向量为 扩散张量的计算由结构张量的特征值，计算出扩散张量的特征值λ1和λ2，再由新的特征值结合结构张量的特征向量计算出扩散张量。
dd＝(λ2-λ1)(j11-j22)/(μ1-μ2) d11＝0.5*[(λ2+λ1)+(λ2-λ1)(j11-j22)/(μ1-μ2)] d12＝(λ1-λ2)j12/(μ1-μ2) d22＝0.5*[(λ2+λ1)-(λ2-λ1)(j11-j22)/(μ1-μ2)] 4个扩散系数的获得。再由扩散张量转化为4个方向的扩散系数 α0＝d11-|d12|， α2＝d22-|d12|。
扩散滤波采用托马斯算法，以达到处理的高效率。
本发明的一次处理循环包括两次各项异性扩散滤波，第一次的扩散系数计算方法为 λ2＝α 初步实现在结构信息较少的区域做平滑，而在结构信息较多的地方保留其结构信息。
第二次扩散系数的计算方法为 if(μ1-μ2)2<ss_low λ1＝α e/se λ1＝0 if(μ1-μ2)2<ss_low or(μ1-μ2)2>ss_high λ2＝α e/se λ2＝0 ss_low，ss_high为预设的两个阈值，(μ1-μ2)2值低于ss_low认为是平坦区域，高于ss_high认为是强结构区域。第二次各向异性扩散滤波进一步平滑平坦区域，而在组织结构的边界上沿着边缘方向扩散，使组织结构信息增强。
两次扩散滤波之后，对结果进行锐化处理，以避免扩散过程可能带来的图像模糊，同时增强组织结构信息。
本发明在完成一次处理循环之后将继续在下一个尺度下对图像做类似的处理。不同尺度下判定区域是否平坦或是包含丰富结构信息，所用的阈值与所在的尺度相关，从而保证处理后的图像质量大幅提升。
权利要求
1、一种提升超声图像质量方法，用于超声成像系统对扫描图像显示数据的优化处理，其的步骤如下
A.利用高斯-拉普拉斯金字塔分解算法分解超声显示图像；
B.对分解产生的高斯层图像数据进行各向异性扩散滤波；
C.对扩散结果进行锐化处理；
D.利用高斯-拉普拉斯金字塔重构算法将各向异性扩散滤波结果与拉普拉斯层数据重构回去；
E.A、B、C、D三个步骤迭代若干次。
2、根据权利要求1所述的提升超声图像质量方法，其特征在于，所述步骤B各向异性扩散滤波过程中，利用图像的局部结构特征来控制扩散滤波的强度。
3、根据权利要求2所述的提升超声图像质量方法，其特征在于，所述图像的局部结构特征是利用结构张量来表征的，利用结构张量，可以表征图像的局部特征，判别是否为图像的结构区域。
4、根据权利要求2所述的提升超声图像质量方法，其特征在于，所述各向异性扩散滤波算法中，由扩散张量表征的图像局部特征，并由此控制扩散滤波强度，扩散滤波的强度由扩散系数确定，在图像的结构区域，扩散系数很小，图像的结构信息得以保留甚至增强；在图像的非结构区域，扩散系数大，噪声得到有效抑制；在图像或组织结构的边缘区域，沿边缘等高线方向扩散。
5、根据权利要求1所述的提升超声图像质量方法，其特征在于，所述超声图像局部特征由结构张量对应的两个特征向量和两个特征值表示，两个特征值差的平方越大，所在区域结构信息越丰富，相反地，两个特征值差的平方越小，所在区域越平坦，设置两个阈值，特征值差的平方小于小阈值的区域认为是平坦区域，大于大阈值的区域认为是结构边缘，介于两者之间的区域认为是结构组织，在这三种区域采用不同的扩散系数计算算法。
6、根据权利要求1所述的提升超声图像质量方法，其特征在于所述步骤B，对由输入图像分解出来的高斯层进行4个方向的各向异性扩散滤波，由以下几个步骤组成
a.计算4个方向扩散系数，4个方向分别是直角坐标系中的-45°、0°、45°和90°方向；
b.在上述4个方向进行扩散滤波。
7、根据权利要求6所述的提升超声图像质量方法，其特征在于，
所述过程a，由以下几个步骤组成
1.高斯平滑；
2.计算梯度；
3.计算结构张量；
4.计算扩散张量；
5.裂解为4个方向的扩散系数。
8、根据权利要求7所述的提升超声图像质量方法，其特征在于，所述1.高斯平滑；2.计算梯度；3.计算结构张量；4.计算扩散张量；5.裂解为4个方向的扩散系数，其步骤具体如下
采用结构张量来表示局部的特征，某个像素的结构张量为其中，Ix为x方向的梯度，Iy为y方向的梯度，结构张量可以分解为两个正交的向量，利用这两个向量的方向和特征值(u1，u2)来表征局部特性；
结构张量的计算要考虑到噪声的影响，计算梯度的时候，首先采用若干阶FIR滤波器平滑图像，平滑之后的图像噪声得到抑制，基于平滑之后的图像求出每个像素的x，y方向的梯度，获得结构张量的组成部分，即
IxIy，并分别对这三个数据组成的阵列再次作平滑处理，进一步减小噪声对整个算法的影响，得到最终需要的结构张量结构张量的特征值和特征向量为
扩散张量的计算由结构张量的特征值，计算出扩散张量的特征值λ1和λ2，再由新的特征值结合结构张量的特征向量计算出扩散张量
dd＝(λ2-λ1)(j11-j22)/(μ1-μ2)
d11＝0.5*[(λ2+λ1)+(λ2-λ1)(j11-j22)/(μ1-μ2)]
d12＝(λ1-λ2)j12/(μ1-μ2)
d22＝0.5*[(λ2+λ1)-(λ2-λ1)(j11-j22)/(μ1-μ2)]
4个扩散系数的获得，再由扩散张量转化为4个方向的扩散系数
α0＝d11-|d12|，
α2＝d22-|d12|。
9、根据权利要求1所述的提升超声图像质量方法，其特征在于，所述步骤D中包含两次各向异性扩散滤波，一次主要是增强组织结构信息，另一次主要是平滑非结构区域同时增强结构边缘信息。
10、根据权利要求9所述的提升超声图像质量方法，其特征在于，所述一次处理循环包括两次各项异性扩散滤波，第一次的扩散系数计算方法为
λ2＝α
初步实现在结构信息较少的区域做平滑，而在结构信息较多的地方保留其结构信息；
第二次扩散系数的计算方法为
if(μ1-μ2)2<ss_low
λ1＝α
else
λ1＝0
if(μ1-μ2)2<ss_low or(μ1-μ2)2>ss_high
λ2＝α
else
λ2＝0
ss_low，ss_high为预设的两个阈值，(μ1-μ2)2值低于ss_low认为是平坦区域，高于ss_high认为是强结构区域，第二次各向异性扩散滤波进一步平滑平坦区域，而在组织结构的边界上沿着边缘方向扩散，使组织结构信息增强。
全文摘要
本发明提供一种提升超声图像质量的方法，用于超声成像系统中对超声扫描图像显示数据的优化处理，包括步骤采用高斯-拉普拉斯金字塔分解方法将超声图像数据分解若干层，对分解结果中的高斯层数据进行各向异性扩散滤波处理，然后对处理后的超声图像数据逆向重构，经过若干次的迭代处理，提升了超声图像质量。其中，各向异性扩散滤波算法根据超声图像的局部结构信息计算出不同方向上的扩散系数，对图像数据进行滤波操作。结果是图像的组织结构区域信息得到增强，而非组织结构区域噪声得到有效抑制。
文档编号G06T5/00GK101452574SQ20081018821
公开日2009年6月10日 申请日期2008年12月19日 优先权日2008年12月19日
发明者刘林泉, 周传涛 申请人:深圳市恩普电子技术有限公司