发明名称 基于模糊噪声图像对联合优化的图像盲复原方法
摘要 本发明公开了一种基于模糊噪声图像对联合优化的图像盲复原方法,用于解决现有图像盲复原方法复原图像效果差的技术问题。技术方案是引入噪声图像,利用噪声图像中包含未模糊的有效边缘信息,通过联合模糊图像和噪声图像的联合目标函数优化求解模糊核和清晰图像。避免了从模糊图像中估计非模糊图像边缘的不准而导致模糊核估计错误,引起复原结果差的问题,获得效果更好、更具细节的图像。提高了复原图像的效果。
申请公布号 CN103337058A 申请公布日期 2013.10.02
申请号 CN201310285234.0 申请日期 2013.07.05
申请人 西北工业大学 发明人 李海森;张艳宁;张海超;孙瑾秋
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于模糊噪声图像对联合优化的图像盲复原方法,其特征在于包括以下步骤:给定模糊图像B和噪声图像N两幅图像,通过如(1)式所示目标能量函数联合求解清晰图像L和模糊核K:<maths num="0001"><![CDATA[<math><mrow><mo>{</mo><mi>K</mi><mo>,</mo><mi>L</mi><mo>}</mo><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><mi>K</mi><mo>,</mo><mi>L</mi></mrow></munder><mo>{</mo><mi>&gamma;</mi><msubsup><mrow><mo>|</mo><mo>|</mo><mi>B</mi><mo>-</mo><mi>K</mi><mo>&CircleTimes;</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mrow><mi>&beta;</mi><mo>|</mo><mo>|</mo><mi>N</mi><mo>-</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mrow><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>K</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mi>&eta;</mi><mo>|</mo><mo>|</mo><mi>G</mi><mo>&CircleTimes;</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mi>&alpha;</mi></msup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,G为梯度提取滤波器组,<img file="FDA00003470936600012.GIF" wi="59" he="59" />为二维卷积操作,γ,β,λ,η均为正则化系数,α为范数选择参数,<img file="FDA00003470936600013.GIF" wi="205" he="77" />为计算<img file="FDA00003470936600014.GIF" wi="131" he="60" />的α范数;对于(1)式的求解,首先,对L进行初始化,令L=N,然后,循环执行以下步骤1,2直到迭代到设定的迭代次数T;1.固定L,求解图像模糊核K;此时(1)式的目标函数简化为下式:<maths num="0002"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>K</mi></munder><mo>{</mo><mi>&gamma;</mi><msubsup><mrow><mo>|</mo><mo>|</mo><mi>B</mi><mo>-</mo><mi>K</mi><mo>&CircleTimes;</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mrow><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>K</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>(2)式通过共轭梯度求解以加快求解速度,求解公式转化为:<img file="FDA00003470936600016.GIF" wi="1293" he="219" />式中,F(·)和F<sup>-1</sup>(·)分别为二维快速傅丽叶变换的正变换和逆变换,<img file="FDA00003470936600017.GIF" wi="98" he="82" />为F(·)的共轭,为矩阵的点乘操作,I为单位矩阵;求解(3)式,得到当前的图像模糊核K;2.固定K,求解复原图像L;此时(1)式的目标函数简化为(4)式:<maths num="0003"><![CDATA[<math><mrow><mi>L</mi><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>L</mi></munder><mo>{</mo><mi>&gamma;</mi><msubsup><mrow><mo>|</mo><mo>|</mo><mi>B</mi><mo>-</mo><mi>K</mi><mo>&CircleTimes;</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mrow><mi>&beta;</mi><mo>|</mo><mo>|</mo><mi>N</mi><mo>-</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msup><mrow><mi>&eta;</mi><mo>|</mo><mo>|</mo><mi>G</mi><mo>&CircleTimes;</mo><mi>L</mi><mo>|</mo><mo>|</mo></mrow><mi>&alpha;</mi></msup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,梯度滤波器组G由C个梯度滤波器[G<sub>1</sub>,G<sub>2</sub>,...,G<sub>C</sub>]组成,(4)式中包含α范数,通过迭代加权最小均方误差法进行近似求解,其步骤如下:迭代求解下面1)、2)两步,直到迭代到设定的迭代次数M;1)计算C个加权矩阵,分别为ω<sup>1</sup>,ω<sup>2</sup>,...,ω<sup>c</sup>,...,ω<sup>C</sup>,且每一个加权矩阵的维数均与L相同,对于第c个加权矩阵ω<sup>c</sup>其求法如下:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&omega;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>c</mi></msubsup><mo>=</mo><msup><mrow><mo>|</mo><msub><mrow><mo>(</mo><msub><mi>G</mi><mi>c</mi></msub><mo>&CircleTimes;</mo><mi>L</mi><mo>)</mo></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>|</mo></mrow><mrow><mi>a</mi><mo>-</mo><mn>2</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA000034709366000110.GIF" wi="224" he="78" />表示<img file="FDA000034709366000111.GIF" wi="140" he="71" />卷积后矩阵的第i行,第j列的元素;2)将式(4)求解近似为<img file="FDA00003470936600021.GIF" wi="1641" he="141" />(6)式为2范数约束求解问题,利用共轭梯度下降方法进行求解,得到当前迭代所求的L。
地址 710072 陕西省西安市友谊西路127号
您可能感兴趣的专利