发明名称 一种基于方向标志位的Bayer彩色图像插值方法
摘要 本发明提供一种基于方向标志位的Bayer彩色图像插值方法,首先对未知绿色分量进行恢复,采用边缘检测算子对每个像素点四个方向的梯度进行检测,根据检测到的梯度大小确定每个像素点处的方向标志位信息,同时计算每个像素点与相邻像素的色差系数,然后根据方向标志位信息和色差系数对未知绿色分量进行插值恢复。恢复出整幅图像的全部绿色分量后,最后基于色差法则恢复出全部的未知红色和蓝色分量,最终完成对Bayer图像的彩色插值。与传统仅仅检测水平和垂直边缘进行单通道独立插值的方法相比,本发明在保持原有方法较低复杂度和较高图像质量的基础上,进一步提高了插值精度,更适合人眼观察,为后续用于高精度图像复原提供了基础。
申请公布号 CN104537625A 申请公布日期 2015.04.22
申请号 CN201510003095.7 申请日期 2015.01.05
申请人 中国科学院光电技术研究所 发明人 刘征;王华闯;徐智勇;魏宇星;蔡敬菊
分类号 G06T5/00(2006.01)I;G06T7/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 一种基于方向标志位的Bayer彩色图像插值方法,其特征在于,包括如下步骤:步骤一、原始Bayer图像边缘梯度检测:原始Bayer图像中,每个像素点处仅存在三基色分量中的某一种,需要恢复该像素点处缺失的其它两种颜色分量,而原始图像中不同像素点处的三种颜色分量间存在相关性,为降低复杂度,仅考虑某一分量的梯度信息来确定方向标志位;所述方向标志位含义定义为:flag=00表示边缘存在像素点水平方向;flag=01表示边缘存在像素点垂直方向;flag=10表示边缘存在像素点对角方向1;flag=11表示边缘存在像素点对角方向2;所述梯度检测算子定义如下:当中心第i行,第j列(i,j)像素处于红色或蓝色分量时,其中G(i,j)代表位于第i行,j列的绿色分量像素灰度值大小,下述符号含义类似;水平梯度检测算子α:α=|G(i,j+1)‑G(i,j‑1)|  (公式1)垂直梯度检测算子β:β=|G(i‑1,j)‑G(i+1,j)|  (公式2)对角梯度检测算子1:<img file="FDA0000651880590000011.GIF" wi="1315" he="128" />(公式3)对角梯度检测算子2:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&omega;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>*</mo><mo>|</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo></mrow>]]></math><img file="FDA0000651880590000012.GIF" wi="1337" he="130" /></maths>  (公式4)当中心(i,j)像素处于绿色分量时,水平梯度检测算子α:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&alpha;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>*</mo><mo>|</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo></mrow>]]></math><img file="FDA0000651880590000013.GIF" wi="1669" he="128" /></maths>  (公式5)垂直梯度检测算子β:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&beta;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>*</mo><mo>|</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo></mrow>]]></math><img file="FDA0000651880590000014.GIF" wi="1668" he="128" /></maths>  (公式6)对角梯度检测算子1:<img file="FDA0000651880590000021.GIF" wi="1647" he="129" />(公式7)对角梯度检测算子2:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&omega;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>*</mo><mo>|</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo></mrow>]]></math><img file="FDA0000651880590000022.GIF" wi="1668" he="130" /></maths>  (公式8)根据定义的梯度检测算子计算出每个像素点处水平、垂直和两个对角方向的梯度值,选择其中最小梯度值的方向作为边缘存在的方向并用定义的标志位进行标记;同时,计算出每个绿色分量像素点与周围3×3邻域内红色分量和蓝色分量像素之间的色差值;其中像素(i,j)处的红绿色差和蓝绿色差定义分别如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>K</mi><mi>gr</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000651880590000023.GIF" wi="1031" he="133" /></maths>  (公式9)<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>K</mi><mi>gb</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mi>B</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000651880590000024.GIF" wi="1028" he="136" /></maths>  (公式10)步骤二、恢复图像未知绿色分量:根据步骤一计算得到的每个像素点标志位信息和色差系数,对图像的未知绿色分量进行恢复;由于人眼视觉特性,高效的插值算法都会沿边缘进行插值,而不是跨越边缘,否则会出现模糊和杂色;因此,针对不同边缘的插值公式如下:当flag分别等于00,01,10,11时,恢复红色分量像素处的未知绿色分量公式分别如下,其中a,b分别是相关权重系数。<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><mo>*</mo><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mi>gr</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>gr</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mo>*</mo><mrow><mo>(</mo><mfrac><mrow><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000651880590000025.GIF" wi="1674" he="149" /></maths>  (公式11)<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>g</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>a</mi><mo>*</mo><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mi>gr</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mi>gr</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mo>*</mo><mrow><mo>(</mo><mfrac><mrow><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000651880590000026.GIF" wi="1570" he="137" /></maths>  (公式12)<img file="FDA0000651880590000027.GIF" wi="1913" he="281" /><img file="FDA0000651880590000028.GIF" wi="1911" he="282" />按照上述公式进行插值,可以恢复出全部的绿色分量;步骤三、恢复未知红色和蓝色分量:对步骤二得到的马赛克图像,利用色差定律恢复出剩余的未知红色和蓝色分量;色差定律假设不同分量之间的差值在图像局部范围内,可以看作是接近于一个常数;利用步骤二恢复出的已知绿色分量和待求颜色分量之间颜色差恒定的特点,计算得到未知的红色和蓝色分量,完成对全部颜色分量的插值恢复,最终得到插值完成的彩色图像。
地址 610209 四川省成都市双流350信箱