基于贝叶斯模型的超声图像斑点噪声滤波方法技术

本发明专利技术公开了一种基于贝叶斯模型的超声图像斑点噪声滤波方法，具体包括：S1、建立基于贝叶斯模型的非局部滤波模型；S2、采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声，进而得出基于Pearson统计距离的权重函数；S3、求解非局部滤波模型中的变量；S4、采用求解后的非局部滤波模型对超声图像中斑点噪声进行滤波处理。本发明专利技术能够有效的滤除超声图像中的斑点噪声，提高超声图像的成像清晰度。

【技术实现步骤摘要】
基于贝叶斯模型的超声图像斑点噪声滤波方法
本专利技术涉及超声图像滤波
，尤其涉及一种基于贝叶斯模型的超声图像斑点噪声滤波方法。
技术介绍
在超声信号的采集过程中，由于斑点噪声的存在，使得重建的超声图像质量明显比CT、MRI等成像模态的差。为了提高超声的图像质量，通常需要对斑点噪声进行滤波去噪。但由于斑点噪声对人体组织的依赖性，导致对斑点噪声的建模和滤除非常困难。对超声图像的斑点噪声滤波问题，虽然已有大量的研究报道，但经典的滤波方法在抑制噪声的同时会丢失图像中的边缘等细节信息。近些年来，基于偏微分方程的各向异性扩散滤波模型开始应用于图像的平滑去噪。各向异性扩散模型的优点是在对图像进行平滑去噪的同时会很好的保留图像中的边缘等细节信息。最简单的扩散方程是把图像平滑看作是各向同性的热传导方程的解。为了获得对加性噪声模型的平滑滤波结果，一个直观的想法就是最小化由估计的图像变化值。1963年Tikhonov提出了对应的变分问题，即在物理学领域著名的热扩散方程。热扩散方程的解等价于观察图像与高斯核函数的卷积，因为卷积算子是线性的，所以热扩散方程的解也是线性的。由于一幅图像与高斯核地卷积等价于它的Fourier变换与另一个高斯核的乘积，所以各向同性的扩散过程本质上是一个低通滤波的过程，它将抑制图像中的高频信号。不幸的是，图像的许多特征，如边缘、纹理等，通常表现为高频信号。因此边缘、纹理等有用信息也将与噪声一起被各向同性的滤波器滤除掉。为了克服上述线性滤波方法导致的各向同性平滑，Perona和Malik对热扩散方程进行了扩展，提出了著名的P-M非线性滤波方法。P-M方程通过把热扩散方程改写成散度的形式，并在散度算子内添加递减函数作为图像特征检测算子，以允许对图象滤波过程进行更加精确的控制。但P-M方程在一些图像点上产生逆扩散。而我们知道，逆扩散虽然会增强图像边缘特征，但噪声同样也有可能得到增强，因此是一个不稳定的扩散过程。但经典的扩散方程只是基于图像局部边缘信息的，没有从全局图像去考虑。Gilboa和Osher引入一种非局部平滑各向异性扩散方程，该方法基本思想是利用图像的冗余信息以及自相似性对图像进行去噪。非局部滤波方法将在全局图像中搜索与当前像素块相似的像素块，最后将所有这些相似像素块的加权平均值作为噪声点处的像素值。但非局部扩散方程只适用于高斯白噪声模型的平滑去噪，但并不适用于服从Gamma统计分布的超声图像斑点噪声滤除。因此，针对上述技术问题，有必要提供一种基于贝叶斯模型的超声图像斑点噪声滤波方法。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种基于贝叶斯模型的超声图像斑点噪声滤波方法。为了达到上述目的，本专利技术实施例提供的技术方案如下：一种基于贝叶斯模型的超声图像斑点噪声滤波方法，所述方法包括：S1、建立基于贝叶斯模型的非局部滤波模型其中，第一项λ为平滑项，第二项(u-u0)2为数据保真项，λ是一个正常数，用于控制平滑项与数据保真项之间的平衡；数据保真项用于使求解的值u不会偏移原始观测值u0太远；S2、采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声，进而得出基于Pearson统计距离的权重函数；S3、求解非局部滤波模型中的变量u和d，其中S4、采用求解后的非局部滤波模型对超声图像中斑点噪声进行滤波处理。作为本专利技术的进一步改进，所述步骤S2包括：采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声：u(x)＝v(x)+vγ(x)η(x)，其中，γ是一个依赖于超声设备并与后续成像过程相关的参数。作为本专利技术的进一步改进，所述步骤S2中γ为0.5。作为本专利技术的进一步改进，所述步骤S2中：两两像素块之间相似度的Pearson统计距离为：基于Pearson统计距离的权重函数为：其中，u(x+·)和u(y+·)为图像中的两个像素块。作为本专利技术的进一步改进，所述步骤S3中非局部滤波模型的求解采用Split-Bregman分解法、最陡梯度下降法、或对偶投影法。作为本专利技术的进一步改进，所述步骤S3中的非局部滤波模型为：作为本专利技术的进一步改进，所述步骤S3中非局部滤波模型的求解通过求解如下以离散形式给出的迭代方程获得：其中，b是与Bregman迭代算法相关的一个变量，也称作Bregman变量，β是一个大于零的常数。作为本专利技术的进一步改进，所述步骤S3中非局部滤波模型中变量u的求解具体包括：固定变量d，则非局部滤波模型简化为求关于变量u的最小化子问题：通过Eular-Lagrange方程，非局部滤波模型的最小解满足：使用高效的Gauss-Seidel迭代法进行求解，即对n≥0有：通过至少两次Gauss-Seidel迭代收敛到方程的最小解处。作为本专利技术的进一步改进，所述步骤S3中非局部滤波模型中变量d的求解具体包括：固定变量u，则非局部滤波模型简化为求关于变量d的最小化子问题：采用软阈值法，则非局部滤波模型的Eular-Lagrange方程为：根据上述公式求解变量d；最后，变量b可根据如下公式进行更新：本专利技术具有以下有益效果：本专利技术根据经对数变换后的超声图像中斑点噪声数理统计特性，获得了基于Pearson统计距离的权重计算函数，是一种非参数统计的方法，避免了传统扩散方程中对参数的复杂设计，对斑点噪声取得了较好的滤波效果。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术基于贝叶斯模型的超声图像斑点噪声滤波方法的流程图。具体实施方式为了使本
的人员更好地理解本专利技术中的技术方案，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本专利技术保护的范围。参图1所示，本专利技术公开了一种基于贝叶斯模型的超声图像斑点噪声滤波方法，具体包括：S1、建立基于贝叶斯模型的非局部滤波模型；S2、采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声，进而得出基于Pearson统计距离的权重函数；S3、求解非局部滤波模型中的变量；S4、采用求解后的非局部滤波模型对超声图像中斑点噪声进行滤波处理。以下结合具体实施方式对本专利技术作进一步说明。Bresson的研究表明，在现有的非局部扩散滤波模型中非局部全变分扩散模型具有更好的平滑滤波效果。根据贝叶斯定理，给非局部扩散项添加数据保真项(fidelityterm)的限制：其中，λ是一个正的常数，用于控制平滑项与数据保真项之间的平衡；第二项(u-u0)2为数据保真项，它的作用是使求解的值u不会偏移原始观测值u0太远，否则也容易导致对图像的过度平滑滤波。从数理统计的角度看，公式(1)中的平滑项对应贝叶斯定理中有关噪声模型的先验项(priorterm)，而数据保真项则对应其中的似然项(likelihoodterm)，因此该扩散方程的求解过程也就是使后验概率(posterior)最大化的本文档来自技高网...

【技术保护点】
一种基于贝叶斯模型的超声图像斑点噪声滤波方法，其特征在于，所述方法包括：S1、建立基于贝叶斯模型的非局部滤波模型JBayesian-NL-TV(u)=∫Ω||▿NLu||+λ2(u-u0)2dx]]>=∫Ω∫Ω(u(y)-u(x))2w(x,y)dy+λ2(u-u0)2dx,]]>其中，第一项λ为平滑项，第二项(u‑u0)2为数据保真项，λ是一个正常数，用于控制平滑项与数据保真项之间的平衡；数据保真项用于使求解的值u不会偏移原始观测值u0太远；S2、采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声，进而得出基于Pearson统计距离的权重函数；S3、求解非局部滤波模型中的变量u和d，其中S4、采用求解后的非局部滤波模型对超声图像中斑点噪声进行滤波处理。

【技术特征摘要】
1.一种基于贝叶斯模型的超声图像斑点噪声滤波方法，其特征在于，所述方法包括：S1、建立基于贝叶斯模型的非局部滤波模型其中，第一项λ为平滑项，第二项(u-u0)2为数据保真项，λ是一个正常数，用于控制平滑项与数据保真项之间的平衡；数据保真项用于使求解的值u不会偏移原始观测值u0太远；S2、采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声，进而得出基于Pearson统计距离的权重函数；S3、求解非局部滤波模型中的变量u和d，其中S4、采用求解后的非局部滤波模型对超声图像中斑点噪声进行滤波处理；在所述步骤S2中，采用Gamma分布拟合经对数压缩后超声图像中的斑点噪声：u(x)＝v(x)+vγ(x)η(x)，其中，γ是一个依赖于超声设备并与后续成像过程相关的参数；两两像素块之间相似度的Pearson统计距离为：基于Pearson统计距离的权重函数为：其中，u(x+·)和u(y+·)为图像中的两个像素块。2.根据权利要求1所述的方法，其特征在于，所述步骤S2中γ为0.5。3.根据权利要求1所述的方法，其特征在于，所述步骤S3中非局部滤波模型的求解采用Split-Bregman分解法、最陡梯度下降法、或对偶投影法。4.根据权利要求1所述的方法，其特征在于，所述步骤S3中的非局部滤波模型为：5.根据权利要求4所述的方法，其特征在于，所述步骤S3中非局部滤波模型的求解通过求解如下以离散形式给出的迭代方程获得：

【专利技术属性】
技术研发人员：温铁祥，李凌，辜嘉，
申请(专利权)人：深圳先进技术研究院，
类型：发明
国别省市：广东;44