基于优化的非局部均值的医学图像去噪方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN107025640A(43)申请公布日2017.08.08(21)申请号201710294540.9(22)申请日2017.04.28(71)申请人南京觅踪电子科技有限公司地址210000江苏省南京市马台街141号3号楼403室(72)发明人王晓芳(74)专利代理机构北京思创大成知识产权代理有限公司11614代理人王尧(51)Int.Cl.G06T5/00(2006.01)G06T7/00(2017.01)权利要求书2页说明书4页附图1页(54)发明名称基于优化的非局部均值的医学图像去噪方法(57)摘要一种基于优化的非局部均值医学图像去噪方法，它包括以下步骤：步骤一、针对待处理图像，确定邻域搜索范围和搜索参数；步骤二、计算当前像素点的邻域范围内所有像素与当前像素的相似度，得到邻域高斯权重矩阵；步骤三、根据邻域高斯权重矩阵计算当前像素点的滤波结果，限定滤波结果在有效地图像范围内；步骤四、重复步骤二和步骤三，完成对整幅图像的遍历，得到去噪结果图并输出。本发明的图像去噪方法，建立有效的计算边界条件，提高参数的有效性，去噪核心算法有效的平滑了医疗图像由于软硬件产生的各种伪影和噪声，同时有效且良好的保留了细节对比度。CN107025640ACN107025640A权利要求书1/2页1.一种基于优化的非局部均值医学图像去噪方法，其特征是它包括以下步骤：步骤一、针对待处理图像，确定邻域搜索范围和搜索参数；步骤二、计算当前像素点的邻域范围内所有像素与当前像素的相似度，得到邻域高斯权重矩阵；步骤三、根据邻域高斯权重矩阵计算当前像素点的滤波结果，限定滤波结果在有效地图像范围内；步骤四、重复步骤二和步骤三，完成对整幅图像的遍历，得到去噪结果图并输出。2.根据权利要求1所述的基于优化的非局部均值医学图像去噪方法，其特征是步骤一具体为：步骤101：获取待处理图像的直方图，调节控制系统a，计算无效像素区域和获取待处理图像直方图中的最小有效像素值k1和最大有效像素值k2；其中：i是直方图灰阶序号，k1是直方图有效像素计算起始位置即最小有效像素值k1，k2是直方图有效像素计算结束位置即最大有效像素值k2，h是直方图数组，a是判断系数，0＜a＜1，n是直方图灰阶的最大值；步骤102：根据待处理图像的直方图，去除无效像素区域，获取邻域搜索的范围；步骤103：设定邻域搜索窗宽的宽度S和计算邻域窗宽宽度C。3.根据权利要求1基于优化的非局部均值医学图像去噪方法，其特征是所述的步骤二具体为：步骤201：遍历图像中所有像素点，在当前像素p上，按照所述的邻域搜索窗宽S遍历以p为中心的SXS邻域；步骤202：对前述SXS邻域中的每个搜索点q，取其CXC大小的邻域Nq，同时取以p为中心CXC大小的邻域Np，计算Np与Nq的相似度，得到搜索点q的高斯权重矩阵，计算公式为：其中：d代表距离运算，d(NP,Nq)＝|f(Np)-f(Nq)|，f(Np)是邻域Np内各点的像素值，f(Nq)是邻域Nq内各点的像素值，h是滤波层度系数；步骤203：对W进行归一化处理。4.根据权利要求3基于优化的非局部均值医学图像去噪方法，其特征是所述的归一化处理计算公式为：其中，2CN107025640A权利要求书2/2页5.根据权利要求1基于优化的非局部均值医学图像去噪方法，其特征是所述的步骤三具体为：步骤301：根据邻域高斯权重矩阵，计算当前像素p的滤波结果NLM(f(p))，：NLM(f(p))＝∑qW(NP,Nq)f(q)公式4其中：NLM为非局部均值函数，f(q)代表当前像素p的SXS邻域中各搜索点q的像素值；步骤302：对滤波结果进行归一化处理，防止像素值溢出。3CN107025640A说明书1/4页基于优化的非局部均值的医学图像去噪方法技术领域[0001]本发明是一种图像处理方法。具体地，是一种基于优化的非局部均值的医学图像去噪方法。背景技术[0002]图像去噪是图像处理的链条的第一步，目的是改善给定的图像，解决实际图像由于噪声干扰而导致图像质量下降的问题。通过去噪技术可以有效地提高图像质量，增大信噪比，更好的体现原来图像所携带的信息，是重要的预处理手段。目前，根据噪声的产生原因、频谱、信噪关系等分类，广泛使用的去噪方法有几十种。[0003]其中，非局部均值方法克服了局部均值模糊了细节信息的缺点，而有效的保留了细节。我们在非局部均值方法的基础上添加了医学人体器官模型，增加了计算边界条件和参数的调控，使得计算精度更进一步的提升。但是，由于非局部均值的原始思路是在对一个像素进行计算时搜索范围是整幅图像，计算量较大，算法效率低。发明内容[0004]本发明的目的是提供一种结合器官模型，建立有效的计算边界条件，提高参数的有