发明名称 基于非局部低秩的色彩滤波阵列图像去马赛克方法
摘要 本发明公开了一种基于非局部低秩的色彩滤波阵列图像去马赛克方法,主要解决现有技术易出现边缘模糊、拉链效应和块效应的问题。其实现步骤为:1.输入一幅色彩滤波阵列图像;2.对绿色通道作方向插值;3.取绿色插值图的待修正块;4.取待修正块的图像块集合并从集合中找出相似块,组成相似块矩阵;5.定义相似块矩阵的拉格朗日优化函数并求解,得到修正的相似块矩阵;6.将修正的相似块矩阵的中间列向量重组,得到修正的块;7.对红色和蓝色通道作方向插值;8.取红色和蓝色插值图的待修正块,重复执行步骤4-6,修正红色和蓝色插值图的所有块;9.输出全彩色图像。本发明能避免边缘模糊,抑制拉链效应和块效应,用于对色彩滤波阵列图像的恢复。
申请公布号 CN103595981A 申请公布日期 2014.02.19
申请号 CN201310512909.0 申请日期 2013.10.25
申请人 西安电子科技大学 发明人 张小华;焦李成;张平;马文萍;马晶晶;田小林;钟桦;白婷
分类号 H04N9/04(2006.01)I 主分类号 H04N9/04(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于非局部低秩的色彩滤波阵列图像去马赛克方法,包括如下步骤:(1)输入一幅色彩滤波阵列图像I;(2)通过方向插值法对色彩滤波阵列图像I中绿色通道缺失的像素进行估计,得到绿色通道图像的插值图像<img file="FDA00004025398400000113.GIF" wi="40" he="57" />;(3)在绿色通道图像的插值图像<img file="FDA00004025398400000114.GIF" wi="38" he="57" />中逐像素取一个5×5大小的图像块,作为当前待修正图像块X;(4)在当前待修正图像块X的中心像素的32×32大小的邻域中取所有5×5的块,组成当前待修正图像块X的图像块集合Ω;(5)计算当前待修正块X与其图像块集合Ω中每个块之间的像素欧式距离d,将欧式距离从小到大排序,取前100个块作为当前待修正块的相似块;(6)将当前待修正块X和其相似块拉成列向量,组成相似块矩阵I<sub>X</sub>,并使待修正块X拉成的待修正列向量Y作为I<sub>X</sub>的中间列向量;(7)定义相似块矩阵I<sub>X</sub>的拉格朗日优化函数F(I′<sub>X</sub>):F(I′<sub>X</sub>)=||I<sub>X</sub>-I′<sub>X</sub>||<sub>F</sub>+λ<sub>1</sub>||I<sub>X</sub>-I′<sub>X</sub>||<sub>1</sub>+λ<sub>2</sub>||I′<sub>X</sub>||<sub>*</sub>,其中,I′<sub>X</sub>为函数F(I′<sub>X</sub>)的自变量,<img file="FDA00004025398400000115.GIF" wi="478" he="138" />c为常数,<img file="FDA00004025398400000116.GIF" wi="79" he="55" />为相似块矩阵I<sub>X</sub>中插值误差的方差,<img file="FDA00004025398400000117.GIF" wi="64" he="61" />为相似块矩阵I<sub>X</sub>的奇异值;(8)最小化约束求解相似块矩阵I<sub>X</sub>的拉格朗日优化函数F(I′<sub>X</sub>),得到修正后的相似块矩阵<img file="FDA0000402539840000013.GIF" wi="84" he="75" /><maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>I</mi><mo>^</mo></mover><mi>X</mi></msub><mo>=</mo><munder><mrow><mi>min</mi><mi>F</mi><mrow><mo>(</mo><msub><mrow><mi>I</mi><mo>'</mo></mrow><mi>X</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow><msub><mrow><mi>I</mi><mo>'</mo></mrow><mi>X</mi></msub></munder></mrow></math>]]></maths>其中,I′<sub>X</sub>为函数F(I′<sub>X</sub>)的自变量;(9)根据修正后的相似块矩阵<img file="FDA0000402539840000015.GIF" wi="82" he="75" />求修正后的列向量<img file="FDA0000402539840000016.GIF" wi="64" he="64" /><maths num="0002"><![CDATA[<math><mrow><mover><mi>Y</mi><mo>^</mo></mover><mo>=</mo><mi>mid</mi><mrow><mo>(</mo><msub><mover><mi>I</mi><mo>^</mo></mover><mi>X</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,<img file="FDA0000402539840000018.GIF" wi="164" he="77" />为取修正后的相似块矩阵<img file="FDA0000402539840000019.GIF" wi="54" he="75" />的中间列向量;(10)将修正后的列向量<img file="FDA00004025398400000110.GIF" wi="38" he="63" />重组为5×5的矩阵,得到修正后的图像块<img file="FDA00004025398400000111.GIF" wi="74" he="65" />(11)对绿色通道图像的插值图像<img file="FDA00004025398400000112.GIF" wi="47" he="66" />的所有图像块重复执行步骤(3)-(10),完成对绿色通道图像的最终估计;(12)通过方向插值法对色彩滤波阵列图像I中红色和蓝色通道缺失的像素进行估计,得到红色通道图像的插值图像<img file="FDA00004025398400000215.GIF" wi="39" he="56" />和得到蓝色通道图像的插值图像<img file="FDA00004025398400000216.GIF" wi="37" he="60" />;(13)分别在红色通道图像的插值图像<img file="FDA00004025398400000217.GIF" wi="34" he="52" />和蓝色通道图像的插值图像<img file="FDA00004025398400000218.GIF" wi="35" he="52" />中逐像素取一个5×5大小的图像块,作为当前待修正的图像块X,重复执行步骤(4)-(10),完成对红色和蓝色通道图像的最终估计;(14)输出含有绿色、红色和蓝色的全彩色图像。
地址 710071 陕西省西安市太白南路2号