发明名称 一种基于可控核的双几何非局部均值图像去噪方法
摘要 本发明公开了一种基于双几何的非局部均值图像去噪方法,其实现过程是:(1)输入含噪的自然图像;(2)在含噪图像中,分别计算以当前像素点i为中心的局部分析窗w<sub>i</sub>和搜索区域S<sub>i</sub>内的权重矩阵<img file="DDA0000468145700000011.GIF" wi="112" he="59" />和<img file="DDA0000468145700000012.GIF" wi="144" he="61" />(3)利用新的权重公式计算当前像素点i与其搜寻区域内所有像素点j的相似性,得到搜寻区域内所有像素点的权值;(4)根据计算出的搜寻区域内所有像素点的权值,对搜寻区域内所有像素点进行加权平均,得到当前像素点修正后的灰度值;(5)用修正后像素点的灰度值取代输入的含噪图像中像素点的灰度值,得到去噪后的图像。本发明的去噪效果比原始非局部方法更好,能够利用图像的几何结构信息,可以更好地在去除噪声的同时保持自然图像的边缘和纹理等细节,可用于对自然图像的去噪处理。
申请公布号 CN103839234A 申请公布日期 2014.06.04
申请号 CN201410060377.6 申请日期 2014.02.21
申请人 西安电子科技大学 发明人 杨淑媛;焦李成;赵琳芳;刘红英;刘芳;马晶晶;马文萍;侯彪;钟桦;王爽
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 西安智萃知识产权代理有限公司 61221 代理人 李东京
主权项 1.一种基于可控核的双几何非局部均值图像去噪方法,包括如下步骤:(1)输入待去噪的含噪图像Y,其中Y=X+N,X是原始干净图像,N是均值为零标准差为σ的高斯白噪声;(2)在含噪图像Y中,分别计算局部分析窗w<sub>i</sub>和搜索区域S<sub>i</sub>内的权重矩阵<img file="FDA0000468145670000011.GIF" wi="124" he="74" />和<img file="FDA0000468145670000012.GIF" wi="150" he="74" />其中w<sub>i</sub>和S<sub>i</sub>是以像素点i为中心大小分别为9×9和17×17的方形窗口;(3)利用如下公式对含噪图像Y中待估计像素点i与其搜索区域S<sub>i</sub>内像素点j之间的权值进行计算,得到搜索区域S<sub>i</sub>内所有像素点的权值;<maths num="0001"><![CDATA[<math><mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>Z</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mo>{</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow><msubsup><mi>w</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>steer</mi></msubsup><mn>2</mn></msubsup><msup><mi>h</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>&times;</mo><msubsup><mi>w</mi><mrow><mi>i</mi><mn>2</mn></mrow><mi>steer</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>Z</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></munder><mo>{</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow><msubsup><mi>w</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>steer</mi></msubsup><mn>2</mn></msubsup><msup><mi>h</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>&times;</mo><msubsup><mi>w</mi><mrow><mi>i</mi><mn>2</mn></mrow><mi>steer</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths>其中,w(i,j)表示含噪图像Y中像素点i和j之间的相似性,即像素点j在像素点i复原过程中所占的权值,w(i,j)满足条件:0≤w(i,j)≤1且<img file="FDA0000468145670000015.GIF" wi="276" he="107" />y<sub>i</sub>是含噪图像Y中局部分析窗w<sub>i</sub>内的图像块,y<sub>j</sub>是含噪图像Y中以像素点j为中心大小为9×9局部分析窗w<sub>j</sub>内的图像块,<img file="FDA0000468145670000016.GIF" wi="124" he="86" />是由权重矩阵<img file="FDA0000468145670000017.GIF" wi="130" he="74" />加权后的欧式距离的平方,平滑参数h=0.9σ,σ是高斯白噪声的标准差,exp(·)表示底数为e的指数运算,<img file="FDA0000468145670000018.GIF" wi="201" he="73" />是权重矩阵<img file="FDA0000468145670000019.GIF" wi="127" he="73" />在像素点j处的值,Z(i)表示归一化常数,S<sub>i</sub>是以像素i为中心大小为17×17的搜索区域;(4)利用(3)中得到的搜索区域S<sub>i</sub>内所有像素点的权值,对含噪图像Y中待估计像素点i的搜索区域S<sub>i</sub>内所有像素点进行加权平均,得到待估计像素点i的估计值:<maths num="0003"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><msub><mi>S</mi><mi>i</mi></msub></mrow></munder><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mi>y</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000468145670000021.GIF" wi="94" he="83" />表示对原始图像X中像素点i灰度值的估计,即为去噪后像素点i的灰度值,y(j)是含噪图像Y中像素点j的灰度值;(5)重复步骤(2)至(4),依次得到每个像素点的估计值;(6)用计算得到的所有像素点的估计值取代含噪图像中所有像素点的灰度值,得到去噪图像<img file="FDA0000468145670000027.GIF" wi="42" he="61" />。
地址 710071 陕西省西安市太白南路2号西安电子科技大学