发明名称 结合非局部先验的复合正则化图像去噪方法
摘要 本发明公开了一种基于图像多特征融合的非局部均值去噪方法,该方法改进了非局部均值去噪方法中图像块间相似度的计算方法,属于图像处理中的图像去噪的研究领域;传统的非局部均值方法采用欧氏距离或高斯加权的欧式距离来衡量块间相似度,这种方式存在一定的缺陷,容易引入不相似图像块,累计造成误差,尤其是在纹理细节区域,无法较好地保持图像的结构信息,降低了去噪性能。针对这样的问题,本方法根据LBP纹理特征来计算图像块之间的相似度,并同时结合LBP特征和灰度特征,进行一种混合相似度计算,最终确定相似图像块的权值分配。基于图像多特征融合的非局部均值去噪方法在去噪效果上得到了显著提升。
申请公布号 CN106204461A 申请公布日期 2016.12.07
申请号 CN201510220372.X 申请日期 2015.05.04
申请人 南京邮电大学 发明人 周宁宁;曹璟
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 汪旭东
主权项 结合非局部先验的复合正则化图像去噪方法,其特征在于,以全变差变分去噪模型为基础,加入非局部结构自相似性的先验信息,建立非局部自相似性先验正则和全变差正则相结合的复合正则化图像去噪模型,并将一种混合相似度计算权值的方法引入非局部先验正则项,包括如下步骤:步骤一:建立图像去噪的全变差模型:假设图像的去噪模型为:f=u+n,其中,f为观测到的含噪图像,u为原始真实图像,n为噪声,图像去噪的全变差模型为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>u</mi><mo>^</mo></mover><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>u</mi></munder><mo>{</mo><mfrac><mi>&mu;</mi><mn>2</mn></mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><mi>u</mi><mo>-</mo><mi>f</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msub><mrow><mo>|</mo><mo>|</mo><mo>&dtri;</mo><mi>u</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000710456910000011.GIF" wi="666" he="133" /></maths>其中,第一项为数据保真项,要求f和u在L<sup>2</sup>范数意义下最相似即差异最小,第二项为全变差正则项,是对图像分片光滑的先验,用以保持边缘结构特征,μ为正则化参数;步骤二:根据图像中的结构自相似性先验信息,建立非局部结构自相似性正则项:假设一幅大小为N×N的图像u经向量化后,表示成u=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>]∈R<sup>N</sup>,P<sub>i</sub>和P<sub>j</sub>分别表示以x<sub>i</sub>和x<sub>j</sub>为中心的图像块,对于像素点x<sub>i</sub>,在整幅图像中找出所有满足给定条件的相似点x<sub>j</sub>,即当认为x<sub>i</sub>和x<sub>j</sub>是相似的,只要P<sub>i</sub>和P<sub>j</sub>的相似度大于某一阈值T,即图像中大量的结构自相似性表示为<img file="FDA0000710456910000012.GIF" wi="571" he="141" />其中,w(x<sub>i</sub>,x<sub>j</sub>)为权函数,满足<img file="FDA0000710456910000013.GIF" wi="401" he="141" />δ(x<sub>i</sub>)是x<sub>i</sub>的非局部邻域,是与x<sub>i</sub>相似度大于阈值T的所有像素点x<sub>j</sub>的集合,满足δ(x<sub>i</sub>)={x<sub>j</sub>|w(x<sub>i</sub>,x<sub>j</sub>)>T},此处,T=0,为非局部邻域内的像素均分配一个大于0的权值,n<sub>i</sub>是残差,<img file="FDA0000710456910000014.GIF" wi="512" he="134" />写成矩阵—向量的形式u=Au+N,其中,向量N∈R<sup>N</sup>的第i个分量为n<sub>i</sub>,矩阵A∈R<sup>N×N</sup>由非局部权值构成,表示为<img file="FDA0000710456910000021.GIF" wi="687" he="167" />u=Au+N作为正则化中的先验信息来刻画图像的非局部结构自相似特征,非局部结构自相似正则化项表示为<img file="FDA0000710456910000022.GIF" wi="407" he="124" />其中,I为N×N大小的单位矩阵;步骤三:在图像去噪的全变差模型中加入非局部结构自相似性正则项,建立复合正则化模型:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>u</mi><mo>^</mo></mover><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>u</mi></munder><mo>{</mo><mfrac><mi>&mu;</mi><mn>2</mn></mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><mi>u</mi><mo>-</mo><mi>f</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msub><mrow><mo>|</mo><mo>|</mo><mo>&dtri;</mo><mi>u</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>+</mo><mfrac><mi>&alpha;</mi><mn>2</mn></mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><mi>A</mi><mo>)</mo></mrow><mi>u</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000710456910000023.GIF" wi="972" he="118" /></maths>α为正则化参数;步骤四:采用分裂Bregman迭代法,求解复合正则化模型;步骤五:输出去噪后图像。
地址 210003 江苏省南京市鼓楼区新模范马路66号