发明名称 一种非局部均值图像去噪算法
摘要 本发明所述的属于图形图像处理领域,具体涉及一种非局部均值去噪算法。本发明的目的在于:针对现有技术存在的问题,提供一种非局部均值图像去噪算法。利用Pseudo-Zernike矩的旋转不变性特征不仅可以匹配到平移相似的patch,还可以匹配到旋转和对称而潜在相似的patch,从而为滤波提供更多的相似patch。进一步的使用二次高斯型核函数替代原算法中的指数型核函数来计算加权系数,从而有效解决图像去噪问题。本发明通过首先计算像素点i、p距离公式,结合像素点i的权重值<img file="dest_path_image002.GIF" wi="42" he="21" />,计算得到计算出i点的灰度估计值<img file="dest_path_image004.GIF" wi="136" he="33" />。
申请公布号 CN104346787A 申请公布日期 2015.02.11
申请号 CN201410685549.9 申请日期 2014.11.25
申请人 成都卫士通信息产业股份有限公司 发明人 陈磊
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 成都九鼎天元知识产权代理有限公司 51214 代理人 韩雪
主权项 一种非局部均值图像去噪算法,其特征在于包括:步骤1:在一幅数字图像中,计算以像素点i为中心的patch N<sub>i</sub>,其n阶m重的pseudo‑zernike矩Z<sub>nm</sub>(i);同理得到以像素点p为中心的patch N<sub>p</sub>,其n阶m重的pseudo‑zernike矩Z<sub>nm</sub>(p);patch N<sub>i</sub>指的是以像素点i为中心的方形区域;patch N<sub>p</sub>指的是以像素点p为中心的方形区域;步骤2:根据步骤1中,像素点i的Z<sub>nm</sub>(i)得到Z(i);像素点p的Z<sub>nm</sub>(p)得到Z(p),进而得到像素点i、p距离公式为:d(i,p)=||Z(i)‑Z(p)||<sup>2</sup>                    (1);步骤3:根据步骤2得到权重w(i,p)为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>C</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mi>exp</mi><mo>{</mo><mo>-</mo><msup><mrow><mo>[</mo><mfrac><msup><mrow><mo>|</mo><mo>|</mo><mi>Z</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mi>Z</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><msup><mi>h</mi><mn>2</mn></msup></mfrac><mo>]</mo></mrow><mn>2</mn></msup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000615652430000011.GIF" wi="1635" he="155" /></maths>其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mi>p</mi></munder><mi>exp</mi><mo>{</mo><mo>-</mo><msup><mrow><mo>[</mo><mfrac><msup><mrow><mo>|</mo><mo>|</mo><mi>Z</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mi>Z</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><msup><mi>h</mi><mn>2</mn></msup></mfrac><mo>]</mo></mrow><mn>2</mn></msup><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0000615652430000012.GIF" wi="603" he="154" /></maths>步骤4:对整幅图像I内的所有点求和,计算出i点的灰度估计值:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>NL</mi><mo>[</mo><mi>v</mi><mo>]</mo><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>p</mi><mo>&Element;</mo><mi>I</mi></mrow></munder><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mi>v</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000615652430000013.GIF" wi="1630" he="118" /></maths>其中v(i)为噪声图像I中被估计的像素点值i的灰度值。
地址 610041 四川省成都市高新区云华路333号