专利名称：一种超声图像噪声抑制方法
技术领域：
本发明涉及超声图像处理领域，具体涉及到一种基于非降采样轮廓波域的超声图像抑制方法。
背景技术：
超声成像是利用超声波照射检测体，通过接收和处理载有检测体内部组织或结构性质特征信息的回波，获得检测体组织性质与结构的可见图像的方法和技木。超声成像具 有价格低廉、成像实时性好、对患者无侵害和便于携带使用等优点，因此，在医学成像领域获得了广泛的应用。当前，依据超声图像的定量分析从事辅助治疗的研究逐渐兴起，例如血管内超声图像的血管壁边缘提取和血管腔面积测量、HIFU治疗中的基于超声图像监控的病灶区域提取和治疗区温度实时测量等，这对超声成像的质量有了更高的要求。由于超声成像的相干特性，图像中固有的斑点噪声(Speckle noise)呈现出密集可见的颗粒状或波纹状物体，在图像中产生众多伪边缘，严重影响到图像中细微特征的分辨和对比度，并给后续的定量分析造成困难。因此，通常需要对超声图像的噪声加以抑制。现有的针对超声图像的噪声抑制方法，主要分为空域滤波和变换域滤波两大类。空域滤波器主要包括Lee、Kun, Frost等人提出的一系列基于局部统计特性的空域滤波器，此类滤波器利用了图像的局部统计參数，在抑制图像中噪声时，可在一定程度上保留图像的边缘。Yu提出的斑点噪声各向异性扩散模型SRAD(Speckle reducinganisotropic diffusion)算法,在图像中平坦区域进行平滑滤波,而在图像的边缘细节区域算法对细节信息进行保留。在基于变换域滤波的超声图像降噪方法中，最为广泛应用的方法是基于小波变换的降噪滤波技木。应用小波变换对图像进行多分辨率分析的方法，改善了在空域和频域中图像信号与噪声提取和分离的难题。如小波变换的软阈值收缩方法考虑到信号与噪声的系数在幅值上的差异，利用软阈值处理对噪声的系数进行滤除；PizuriCa提出ー种基于联合估计的广义似然比降噪GenLik算法，在对小波系数进行降噪处理时，并不需要利用先验知识对噪声的分布形式进行假设，而是选择从系数自身中提取信号和噪声的分布信息；而基于贝耶斯估计的降噪方法针对对数变换后超声图像的小波系数，利用最小均方误差检测(,minimum mean squared error, MMSE)或最大后验概率(maximum a posteriori, MAP)检测理论结合对信号和噪声的分布概率模型(包括高斯分布、瑞利分布、Nakagami分布和广义高斯分布等)作參数化近似，可在降噪滤波中尽可能好的保留图像信号。但是，现有技术在采用以上的方法对超声图像中斑点噪声进行抑制时，很难保证在有效抑制噪声的同时能保证图像中的关键有用信息不会被破坏。为了更加有效的抑制超声图像中的斑点噪声，各类提升小波变换域的降噪算法完善了小波域方法的各种不足，成为了当前的ー个重要研究方向。如轮廓波变换(Contourlet)是ー种功能强大且用途广泛的ニ维信号变换工具，与小波变换相比有着更好的多分辨率分析、多方向性表示性能，可以准确的捕捉图像中内在的几何结构。非降采样轮廓波变换(Nonsubsampling contourlet, NSCT)在Contourlet变换的基础上，利用非降采样塔式分解和方向滤波器组构造分解结构，避免了降采样过程，使得变换具有平移不变性，并可由多孔算法实现，进ー步提高了 Contourlet变换在图像降噪领域的性能。
发明内容
有鉴于此，本发明的目的是提供一种基于非降采样轮廓波变换的超声图像噪声抑制方法，以在降噪过程中同时实现对噪声的充分抑制和对细节信息的有效保留，以满足医用超声成像仪器在成像时进行快速、有效降噪处理的需要。为实现上述目的，本发明提供了一种超声图像噪声抑制方法，该方法包括步骤一、对超声图像进行对数变换，以获取对数图像；步骤ニ、提取对数图像的边缘信号区；步骤三、对对数图像进行非降采样轮廓波变换，以得到多个高频子带；步骤四、提取出各高频子带的边缘信号区和信号平稳区；步骤五、对各高频子带在边缘信号区和信号平缓区分别进行降噪处理；步骤六、对降噪后的高频子带进行非降采样轮廓波反变换；步骤七、对步骤六得到的图像进行指数变换得到最終降噪图像。较佳地，在步骤ニ中提取到的对数图像的边缘信号区之外的点组成其信号平稳区，以mi表示边缘信号区，以Hl2表示信号平稳区。较佳地,在步骤三中所述的多个高频子带具有尺度和方向，以bj, i表示高频子带，其中j表示子带的尺度，I表示子带的方向。较佳地,在步骤三中对各尺度下各方向的高频子带by,首先以=mi se/提取各高频子带内的边缘信号区味,〃，以=、、,〃提取信号平稳区Hl2J1 ;其中，@代表形态学膨胀运算，Sej表示j尺度下膨胀使用的结构元素，信号平稳区m2,M为Hi^1的补集。较佳地，在步骤五中进ー步包括步骤5. 1，在边缘信号区Hi1,い，保持其中各点的系数值；以及步骤5. 2，在信号平稳区以软阈值收缩方法进行降噪处理。较佳地，在所述步骤5. 2中，以软阈值收缩方法对信号平稳区m2, j, x进行降噪处理的方法进ー步包括计算所述信号平稳区m2, j, x中各非降采样轮廓波变换系数的标准差Sb ;选择t · Sb作为阈值，利用软阈值收缩方程对信号平稳区进行降噪处理，得到降噪后的NSCT子帯，软域值收缩方程为
权利要求
1.一种超声图像噪声抑制方法，其特征在于，包括 步骤一、对超声图像进行对数变换，以获取对数图像； 步骤ニ、提取对数图像的边缘信号区； 步骤三、对对数图像进行非降采样轮廓波变换，以得到多个高频子带； 步骤四、提取出各高频子带的边缘信号区和信号平稳区； 步骤五、对各高频子带在边缘信号区和信号平缓区分别进行降噪处理； 步骤六、对降噪后的高频子带进行非降采样轮廓波反变换； 步骤七、对步骤六得到的图像进行指数变换得到最終降噪图像。
2.如权利要求I所述的方法，其特征在于，在步骤ニ中提取到的对数图像的边缘信号区之外的点组成其信号平稳区，以Hl1表示边缘信号区，以Hl2表示信号平稳区。
3.如权利要求2所述的方法，其特征在干，在步骤三中所述的多个高频子带具有尺度和方向，以by表示高频子带,其中j表示子带的尺度，I表示子带的方向。
4.如权利要求3所述的方法，其特征在于，对各尺度下各方向的高频子带\1；首先以mU,/ =mi se/提取各高频子带内的边缘信号区mm,以m2,j,i = 提取信号平稳区m2,j,!;其中，@代表形态学膨胀运算，Sち表示j尺度下膨胀使用的结构元素，信号平稳区I为mI, j,I的补集。
5.如权利要求4所述的方法，其特征在于，在步骤五中进ー步包括 步骤5. I，在边缘信号区保持其中各点的系数值；以及 步骤5. 2，在信号平稳区以软阈值收缩方法进行降噪处理。
6.如权利要求5所述的方法，其特征在于，在所述步骤5.2中，以软阈值收缩方法对信号平稳区m2, 口进行降噪处理的方法进ー步包括 计算所述信号平稳区%ja中各非降采样轮廓波变换系数的标准差δ b ； 选择t · Sb作为阈值，利用软阈值收缩方程对信号平稳区进行降噪处理，得到降噪后的NSCT子带，软域值收缩方程为:
7.如权利要求6所述的方法，其特征在于，所述參数t取值为2 3。
8.如权利要求3所述的方法，其特征在于，所述尺度最大值为5层，每尺度下分解的方向数为16。
全文摘要
本发明提供了一种超声图像噪声抑制方法，该方法包括步骤一、对超声图像进行对数变换，以获取对数图像；步骤二、提取对数图像的边缘信号区；步骤三、对对数图像进行非降采样轮廓波变换，以得到多个高频子带；步骤四、提取出各高频子带的边缘信号区和信号平稳区；步骤五、对各高频子带在边缘信号区和信号平缓区分别进行降噪处理；步骤六、对降噪后的高频子带进行非降采样轮廓波反变换；步骤七、对步骤六得到的图像进行指数变换得到最终降噪图像。采用本发明提供的噪声抑制方法，能够有效地对斑点噪声进行抑制，并且较大程度上保留超声图像中的细节信息，达到较佳的降噪效果。
文档编号G06T5/00GK102646270SQ20121006734
公开日2012年8月22日 申请日期2012年3月14日 优先权日2012年3月14日
发明者侯朝焕, 原建平, 闫晟 申请人:中国科学院声学研究所