基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法.pdf

1、10申请公布号CN102073994A43申请公布日20110525CN102073994ACN102073994A21申请号201010617508822申请日20101231G06T5/00200601A61B8/0020060171申请人哈尔滨工业大学地址150001黑龙江省哈尔滨市南岗区西大直街92号72发明人谷延锋于潇崔兆宇74专利代理机构哈尔滨市松花江专利商标事务所23109代理人张宏威54发明名称基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法57摘要基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，它涉及一种医学超声图像噪声抑制方法。它为解决将现有各项异性扩散技术应用于

2、超声图像时，用梯度作为边缘检测算子难以准确区分图像边缘和噪声的问题而提出。一对图像进行小波分解，计算标准模值；二若为首次迭代，估计标准模值的分布参数，否则跳到步骤四；三分类得到各尺度下的噪声平稳区域；四用噪声区域内标准模值的均值确定扩散阈值；并代入扩散函数得到扩散参数；五用扩散参数对小波系数进行加权后，进行小波逆变换；六比较相邻两次迭代中的扩散阈值，判断是否终止。本发明在多尺度下对图像进行处理，结果更加精确，可广泛适用于各种需要对超声图像进行处理的场合。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书5页附图2页CN102074000A1/2页21基于多尺度各

3、项异性扩散的超声医学图像散斑噪声抑制方法，其特征在于所述方法由如下步骤完成步骤一使用MALLAT定义的二次样条小波滤波器对图像进行非抽样小波分解，得到三个尺度下的近似分量，水平方向的小波系数和竖直方向的小波系数，计算小波系数标准模值MJF；步骤二判断是否为首次迭代；是，则用混合瑞利模型来模拟图像扇形扫描区内小波系数标准模值MJF的分布，并设定出每个尺度下模型参数的值，执行步骤三；否，则跳到步骤四；步骤三使用极大似然分类的方法得出各尺度下的噪声平稳区域；步骤四通过计算每个尺度下噪声区域内MJF的均值确定扩散阈值；并将该扩散阈值代入扩散函数求得扩散参数PJ；步骤五用步骤四求得的扩散参数PJ对步骤一

4、分解得到的小波系数进行加权计算，然后对图像进行小波逆变换，完成一次多尺度各向异性扩散过程；步骤六比较相邻两次迭代过程中的扩散阈值，判断它们之间相差的百分比是否小于阈值TH；是，则抑制完成；否，则返回步骤一进行下次迭代。2根据权利要求1所述的基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，其特征在于步骤一中归一化小波系数标准模值MJF采用如下算法得到首先计算小波系数模公式1式中参数和分别代表J尺度下水平方向的小波系数和竖直方向的小波系数；然后计算标准模值公式2式中S是图像中一活动窗口内的灰度均值，各尺度下窗的大小定义为DJ2J1D011，参数D03为J1时窗的大小。3根据权利要求2所述的基

5、于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，其特征在于步骤二中用混合瑞利模型来模拟图像扇形扫描区内小波系数标准模值MJF分布的方法为边缘和噪声处的标准模值MJF分别服从方差为和的瑞利分布，概率密度函数分别由PEX和PNX表示则图像扇形扫描区内的MJF服从混合瑞利分布PXNPNX1NPEX公式3式中参数和N的值通过期望最大化方法估计得出。4根据权利要求3所述的基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，其特征在于步骤三使用极大似然分类的方法得出各尺度下的噪声平稳区域的方法如下利用极大似然分类，错误概率最小时得到分类阈值权利要求书CN102073994ACN102074000A2

6、/2页3公式4然后，按照从大尺度到小尺度的方法对图像进行分类公式5其中，参数U1的区域表示噪声平稳区域。5根据权利要求4所述的基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，其特征在于步骤四中扩散参数PJ的求得方法如下选择非线性扩散函数公式6令MJF；扩散阈值J取为对应尺度下噪声平稳区域内标准模值MJF的均值与系数的加权公式7扩散参数PJ由下式确定PJ1CJ公式8。6根据权利要求5所述的基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，其特征在于步骤五中通过扩散参数PJ对小波系数进行加权计算的方法如下将小波系数WJF乘以扩散参数PJ，得到新的小波系数WJF；公式9。权利要求书CN10

7、2073994ACN102074000A1/5页4基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法技术领域0001本发明涉及一种超声图像散斑噪声抑制方法。背景技术0002超声成像技术现广泛用于对诸如心脏的体内结构进行成像。由于成像的物理因素，相干波的叠加形成散斑噪声，严重影响了图像质量，尤其是使边缘和细节变得模糊，降低了图像的可观性和分割精度，限制了超声医学图像在临床诊断中的进一步应用。对超声图像进行降噪处理，可以为图像分割、斑点追踪等后续处理提供更有利的条件。因此，研究在保持图像边缘和细节信息的基础上有效抑制散斑噪声，有助于提高临床诊断的准确性，具有重要意义。0003从目前在该领域内对超

8、声医学图像进行噪声抑制的方法来看，主要有空域滤波、小波变换和各向异性扩散算法等。空域滤波主要指邻域平均法和中值滤波等技术，其本质是根据像素点邻域窗口内的信息确定一个新值来取代原值，在抑制噪声的同时往往会引起图像的模糊和重要细节的丢失，限制了处理的效果。基于小波变换的方法将超声图像变换到小波域，通过阈值方法去除某些小波系数，再逆变换以去除噪声，该方法主要用于加性噪声的去除，对于如散斑噪声等乘性噪声阈值难以确定。基于各向异性扩散方程的降噪技术使用选择性扩散方式，在图像的噪声处有较大的扩散系数，有利于图像的平滑，而在边缘处有较小的扩散系数，保持了图像的细节，但超声图像中，由于乘性噪声的影响，难以从灰

9、度变化上准确的区分图像的边缘和噪声，加大了算法的难度。0004各项异性扩散的过程可通过小波变换的方式来实现，过程如图7所示，将分解得到的小波系数逐一加权后再进行逆变换输出结果，相当于在多个尺度下进行各向异性扩散。该方法能够在多尺度下对图像进行处理，结果更加精确，同时继承了各向异性扩散能够保持图像细节的特点，不失为超声图像处理的好办法。该方法中，如何针对超声医学图像的特点选用合适的边缘检测器，以及如何调整扩散参数是关系到降噪质量好坏的最关键问题。发明内容0005本发明为了解决现有基于各向异性扩散方程的降噪技术使用梯度检测边缘的方法易受乘性噪声影响，难以从灰度变化上准确区分图像边缘和噪声的问题，而

10、提出的基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法。0006基于多尺度各项异性扩散的超声医学图像散斑噪声抑制方法，所述方法由如下步骤完成0007步骤一使用MALLAT定义的二次样条小波滤波器对图像进行非抽样小波分解，得到三个尺度下的近似分量，水平方向的小波系数和竖直方向的小波系数，计算小波系数标准模值MJF；说明书CN102073994ACN102074000A2/5页50008步骤二判断是否为首次迭代；是，则用混合瑞利模型来模拟图像扇形扫描区内小波系数标准模值MJF的分布，并设定出每个尺度下模型参数的值，执行步骤三；否，则跳到步骤四；0009步骤三使用极大似然分类的方法得出各尺度下的

11、噪声平稳区域；0010步骤四通过计算每个尺度下噪声区域内MJF的均值确定扩散阈值；并将该扩散阈值代入扩散函数求得扩散参数PJ；0011步骤五用步骤四求得的扩散参数PJ对步骤一分解得到的小波系数进行加权计算，然后对图像进行小波逆变换，完成一次多尺度各向异性扩散过程；0012步骤六比较相邻两次迭代过程中的扩散阈值，判断它们之间相差的百分比是否小于阈值TH；是，则抑制完成；否，则返回步骤一进行下次迭代。0013本发明排除了乘性噪声影响，可以从图像灰度变化上更加准确地区分图像边缘和噪声。本发明所述方法在多尺度下对图像进行处理，结果更加精确，同时继承了各向异性扩散能够保持图像细节的特点。本发明可广泛适用

12、于各种需要对超声图像进行处理的场合。附图说明0014图1为本发明所述方法的流程图；图2为实际图像中相同尺度下标准模值的分布直方图与混合瑞利分布模型曲线的比较示意图，图中实心点为标准模值的分布直方图，曲线为混合瑞利分布概率密度曲线；图3为扩散阈值随迭代次数的变化情况；图4为原始未处理的超声图像；图5为采用本发明所述方法当K05时输出的图像；图6为采用本发明所述方法当K2时输出的图像；图7为用小波变换实现一次多尺度各向异性扩散的示意框图。具体实施方式0015具体实施方式一结合图1至图7至说明本实施方式，所述方法由如下步骤完成0016步骤一使用MALLAT定义的二次样条小波滤波器对图像进行非抽样小波

13、分解，得到三个尺度下的近似分量，水平方向的小波系数和竖直方向的小波系数，计算小波系数标准模值MJF；0017步骤二判断是否为首次迭代；是，则用混合瑞利模型来模拟图像扇形扫描区内小波系数标准模值MJF的分布，并设定出每个尺度下模型参数的值，执行步骤三；否，则跳到步骤四；0018步骤三使用极大似然分类的方法得出各尺度下的噪声平稳区域；0019步骤四通过计算每个尺度下噪声区域内MJF的均值确定扩散阈值；并将该扩散阈值代入扩散函数求得扩散参数PJ；0020步骤五用步骤四求得的扩散参数PJ对步骤一分解得到的小波系数进行加权计算，然后对图像进行小波逆变换，完成一次多尺度各向异性扩散过程；0021步骤六比较

14、相邻两次迭代过程中的扩散阈值，判断它们之间相差的百分比是否小于阈值TH；是，则抑制完成；否，则返回步骤一进行下次迭代。说明书CN102073994ACN102074000A3/5页60022步骤一中小波系数标准模值MJF采用如下算法得到0023首先计算小波系数模0024公式10025式中参数和分别代表J尺度下水平方向的小波系数和竖直方向的小波系数；然后计算标准模值0026公式20027式中S是图像中一活动窗口内的灰度均值，各尺度下窗的大小定义为DJ2J1D011，参数D03为J1时窗的大小。0028步骤一采用的小波滤波器的系数参加下表002900300031步骤二中步骤二中用混合瑞利模型来模拟

15、图像扇形扫描区内小波系数标准模值MJF分布的方法为边缘和噪声处的标准模值MJF分别服从方差为和的瑞利分布，概率密度函数分别由PEX和PNX表示00320033则图像扇形扫描区内的MJF服从混合瑞利分布0034PXNPNX1NPEX公式30035式中参数和N的值通过期望最大化方法估计得出。实际分布曲线与标准混合瑞利模型曲线如图2所示。0036步骤三使用极大似然分类的方法得出各尺度下的噪声平稳区域的方法如下图像扫描区可分为两部分细节区域和噪声区域；利用极大似然分类，错误概率最小时得到分类阈值0037公式40038按照从大尺度到小尺度的方法对图像进行分类说明书CN102073994ACN102074

16、000A4/5页70039公式50040当参数U1的区域表示噪声平稳区域，且在最大尺度JJ下认为UJ1；K为可调节参数，K较小时如K05可用于增强图像可读性，K较大时如K2可用于后续分割。0041步骤四中扩散参数PJ的计算方法如下选择非线性扩散函数0042公式60043令MJF；扩散阈值J取为对应尺度下噪声平稳区域内标准模值MJF的均值与系数的加权0044公式70045扩散参数PJ由下式确定0046PJ1CJ公式80047步骤五中通过扩散参数PJ对小波系数进行加权计算的方法如下0048将小波系数WJF乘以扩散参数PJ，得到新的小波系数WJF；0049公式9。0050随着迭代次数的增加，扩散阈值

17、逐渐趋于稳定，比较前后两次迭代过程中的扩散阈值，相差百分比小于指定阈值TH时停止迭代；设N和N1分别为由公式7算得的第N次和第N1次迭代中的扩散阈值，当所有尺度下的JJ1，J均满足以下条件时迭代停止0051公式100052阈值TH一般取05，在各尺度下随迭代次数的变化情况如图3所示。0053本发明采用小波系数标准模值MJF作为超声图像边缘检测器的原理0054原始的散斑噪声是一种乘性噪声，幅值具有瑞利分布的特性。在医学超声成像系统中，为改善质量图像在输出前往往要经过同态滤波等预处理，输出图像数学模型可表示如下0055公式110056其中，参数GX为原始无噪图像。设图像中一移动窗口内原始灰度均值为

18、，则该窗口内原始图像可定义为GXRX，则公式11可写为0057公式120058对上式进行小波变换0059公式130060窗口很小时可认为RX1，由小波变换性质易知则上式中等号说明书CN102073994ACN102074000A5/5页8右侧第一项为零，即故将小波系数模值除以局部窗内灰度均值可去除原始图像自身灰度的影响。同理，小波系数模值除以局部窗内灰度均值得到的标准模值，去除了原始图像自身灰度的影响，使它相对于现有的使用梯度检测边缘的方法更准确更有效，成为对超声医学图像更准确更有效的边缘检测器。0061以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明所提交的权利要求书确定的专利保护范围。说明书CN102073994ACN102074000A1/2页9图1说明书附图CN102073994ACN102074000A2/2页10图2图3图4图5图6图7说明书附图CN102073994A