CN109285113B - 一种基于梯度的改进彩色图像插值方法 - Google Patents

一种基于梯度的改进彩色图像插值方法 Download PDF

Info

Publication number
CN109285113B
CN109285113B CN201710598476.3A CN201710598476A CN109285113B CN 109285113 B CN109285113 B CN 109285113B CN 201710598476 A CN201710598476 A CN 201710598476A CN 109285113 B CN109285113 B CN 109285113B
Authority
CN
China
Prior art keywords
interpolation
color
gradient
pixel
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710598476.3A
Other languages
English (en)
Other versions
CN109285113A (zh
Inventor
叶涛
杨俊雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Remote Sensing Equipment
Original Assignee
Beijing Institute of Remote Sensing Equipment
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Remote Sensing Equipment filed Critical Beijing Institute of Remote Sensing Equipment
Priority to CN201710598476.3A priority Critical patent/CN109285113B/zh
Publication of CN109285113A publication Critical patent/CN109285113A/zh
Application granted granted Critical
Publication of CN109285113B publication Critical patent/CN109285113B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformation in the plane of the image
    • G06T3/40—Scaling the whole image or part thereof
    • G06T3/4015—Demosaicing, e.g. colour filter array [CFA], Bayer pattern

Abstract

本发明公开了一种基于梯度的改进彩色图像插值方法，本发明方法首先根据Bayer滤波阵列判断单色图像相应像素位置的颜色，然后根据相应像素位置的单色彩完成缺失颜色的插值，最后按照从上到下、从左到右扫描整幅单色图像，完成整幅图像的彩色插值。本发明方法通过卷积技术，将梯度插值方法改进为逐行扫描、实时插值的方式，每一个卷积过程，得到相应像素点的三色取值。相比于传统的彩色图像插值方法，本方法的实时性和并行性均得到很大的提升。

Description

一种基于梯度的改进彩色图像插值方法
技术领域
本发明属于彩色相机彩色图像合成的技术领域，具体涉及一种基于梯度的改进彩色图像插值方法。
背景技术
考虑降低成本和缩小数码设备的体积，绝大多数数码摄影设备只使用单个图像传感器，并在图像传感器的表面覆盖一层色彩滤波阵列，这样，每个像素只采集到红色、绿色、蓝色三基色中的一种，而缺失另外两种基色信息，缺失的两个基色需要利用周围像素点的色彩信息插值得到，进而复原一副全彩色图像，该过程称作彩色图像插值或去马赛克。目前业界存在多种彩色滤波阵列，而Bayer彩色滤波阵列应用最为普遍。
传统的彩色图像插值法通常采用先计算并存储整幅绿色通道的图像，然后通过整体插值的方式计算红色或蓝色通道的图像，实时性差，且不利于在FPGA等并行硬件平台上的实现。
发明内容
本发明目的在于提供一种基于梯度的改进彩色图像插值方法，解决传统方法实时性差和不利于并行实现的问题。
一种基于梯度的改进彩色图像插值方法，包含以下步骤：
第一步判断单色图像相应像素位置的颜色
根据Bayer滤波阵列，设像素在单色图像中的坐标位置为(i,j)，当(i+j)除以2的余数为0时，则判定像素颜色为绿色；当i除以2的余数为0且j除以2的余数不为0时，则判定像素颜色为蓝色；当i除以2的余数不为0且j除以2的余数为0时，则判定像素颜色为红色。
第二步根据像素的单色彩完成缺失颜色的插值
忽略边缘像素的影响，设像素在单色图像中的坐标位置为(i,j)，该点的像素值为d(i,j)，该点的颜色为红色、绿色或蓝色，其中G(i,j)、R(i,j)、B(i,j)分别表示图像中的坐标位置为(i,j)处像素的绿色分量、红色分量、蓝色分量的灰度值：
当单色图像坐标(i,j)处颜色为绿色时，其绿色分量为：
G(i,j)＝d(i,j) (1)
在奇数行偶数列位置的蓝色和红色分量的插值分别为：
Figure BDA0001356471830000021
Figure BDA0001356471830000022
在偶数行奇数列位置的蓝色和红色分量的插值分别为：
Figure BDA0001356471830000023
Figure BDA0001356471830000024
当单色图像坐标(i,j)处颜色为蓝色时，其蓝色分量为：
B(i,j)＝d(i,j) (6)
其绿色分量的插值由水平和竖直方向的梯度大小决定，设α为该点水平方向的梯度，β为其竖直方向的梯度：
当α＜β时，表明d(i,j)与其水平方向上的像素点关系紧密，插值将在水平方向上进行，公式如下所示：
Figure BDA0001356471830000025
当α＞β时，表明d(i,j)与其竖直方向上的像素点关系紧密，插值将在竖直方向上进行，公式如下所示：
Figure BDA0001356471830000026
当α＝β时，表明d(i,j)与其水平方向上和垂直方向上的像素点关系同等紧密，插值将在整个邻域内进行，公式如下所示：
Figure BDA0001356471830000031
其红色分量的插值为：
Figure BDA0001356471830000032
当单色图像坐标(i,j)处颜色为红色时，其红色分量为：
R(i,j)＝d(i,j) (11)
其绿色分量的插值由水平和竖直方向的梯度大小决定，设α为该点水平方向的梯度，β为其竖直方向的梯度：
当α＜β时，表明d(i,j)与其水平方向上的像素点关系紧密，插值将在水平方向上进行，公式如下所示：
Figure BDA0001356471830000033
当α＞β时，表明d(i,j)与其竖直方向上的像素点关系紧密，插值将在竖直方向上进行，公式如下所示：
Figure BDA0001356471830000034
当α＝β时，表明d(i,j)与其水平方向上和垂直方向上的像素点关系同等紧密，插值将在整个邻域内进行，公式如下所示：
Figure BDA0001356471830000035
其蓝色分量的插值为：
Figure BDA0001356471830000036
从上到下，从左到右扫描整幅单色图像，结合公式(1)-(15)，完成整幅图像的彩色插值。
更优的，其中第二步中水平方向的梯度α和竖直方向的梯度β通过二阶梯度公式获取：
Figure BDA0001356471830000041
一种基于梯度的改进彩色图像插值系统，包括：单色图像像素位置颜色判定模块和颜色插值模块；其中单色图像像素位置颜色判定模块根据Bayer滤波阵列对单色图像相应像素位置的颜色进行判定，颜色插值模块根据像素的单色彩完成缺失颜色的插值。
本发明方法通过卷积技术，将梯度插值方法改进为逐行扫描、实时插值的方式，每一个卷积过程，得到相应像素点的三色取值。相比于传统的彩色图像插值方法，本方法的实时性和并行性均得到很大的提升。
附图说明
图1现有技术中的专业数码相机的工作示意图；
图2Bayer滤波阵列示意图。
具体实例方式
一种基于梯度的改进彩色图像插值方法，包含以下步骤：
第一步判断单色图像相应像素位置的颜色
根据Bayer滤波阵列，设像素在单色图像中的坐标位置为(i,j)，当(i+j)除以2的余数为0时，则判定像素颜色为绿色；当i除以2的余数为0且j除以2的余数不为0时，则判定像素颜色为蓝色；当i除以2的余数不为0且j除以2的余数为0时，则判定像素颜色为红色。
第二步根据像素的单色彩完成缺失颜色的插值
忽略边缘像素的影响，设像素在单色图像中的坐标位置为(i,j)，该点的像素值为d(i,j)，该点的颜色为红色、绿色或蓝色，其中G(i,j)、R(i,j)、B(i,j)分别表示图像中的坐标位置为(i,j)处像素的绿色分量、红色分量、蓝色分量的灰度值：
当单色图像坐标(i,j)处颜色为绿色时，其绿色分量为：
G(i,j)＝d(i,j) (1)
在奇数行偶数列位置的蓝色和红色分量的插值分别为：
Figure BDA0001356471830000051
Figure BDA0001356471830000052
在偶数行奇数列位置的蓝色和红色分量的插值分别为：
Figure BDA0001356471830000053
Figure BDA0001356471830000054
当单色图像坐标(i,j)处颜色为蓝色时，其蓝色分量为：
B(i,j)＝d(i,j) (6)
其绿色分量的插值由水平和竖直方向的梯度大小决定，设α为该点水平方向的梯度，β为其竖直方向的梯度，二阶梯度公式为：
Figure BDA0001356471830000055
当α＜β时，表明d(i,j)与其水平方向上的像素点关系紧密，插值将在水平方向上进行，公式如下所示：
Figure BDA0001356471830000056
当α＞β时，表明d(i,j)与其竖直方向上的像素点关系紧密，插值将在竖直方向上进行，公式如下所示：
Figure BDA0001356471830000057
当α＝β时，表明d(i,j)与其水平方向上和垂直方向上的像素点关系同等紧密，插值将在整个邻域内进行，公式如下所示：
Figure BDA0001356471830000061
其红色分量的插值为：
Figure BDA0001356471830000062
当单色图像坐标(i,j)处颜色为红色时，其红色分量为：
R(i,j)＝d(i,j) (12)
其绿色分量的插值由水平和竖直方向的梯度大小决定，设α为改点水平方向的梯度，β为其竖直方向的梯度，其二阶梯度公式如下：
Figure BDA0001356471830000063
当α＜β时，表明d(i,j)与其水平方向上的像素点关系紧密，插值将在水平方向上进行，公式如下所示：
Figure BDA0001356471830000064
当α＞β时，表明d(i,j)与其竖直方向上的像素点关系紧密，插值将在竖直方向上进行，公式如下所示：
Figure BDA0001356471830000065
当α＝β时，表明d(i,j)与其水平方向上和垂直方向上的像素点关系同等紧密，插值将在整个邻域内进行，公式如下所示：
Figure BDA0001356471830000066
其蓝色分量的插值为：
Figure BDA0001356471830000071
从上到下，从左到右扫描整幅单色图像，结合公式(1)-(17)，完成整幅图像的彩色插值。
以上就是本发明一种基于梯度的改进彩色图像插值方法的技术内容和方案。需要注意的是，以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和范围之内所作的任何修改、等同替换或改进等，均包含在本发明的保护范围之内。

Claims (2)

1.一种基于梯度的改进彩色图像插值方法，其特征在于包含以下步骤：
第一步 判断单色图像相应像素位置的颜色
根据Bayer滤波阵列，设像素在单色图像中的坐标位置为(i,j)，当(i+j)除以2的余数为0时，则判定像素颜色为绿色；当i除以2的余数为0且j除以2的余数不为0时，则判定像素颜色为蓝色；当i除以2的余数不为0且j除以2的余数为0时，则判定像素颜色为红色；
第二步 根据像素的单色彩完成缺失颜色的插值
忽略边缘像素的影响，设像素在单色图像中的坐标位置为(i,j)，该点的像素值为d(i,j)，该点的颜色为红色、绿色或蓝色，其中G(i,j)、R(i,j)、B(i,j)分别表示图像中的坐标位置为(i,j)处像素的绿色分量、红色分量、蓝色分量的灰度值：
当单色图像坐标(i,j)处颜色为绿色时，其绿色分量为：
G(i,j)＝d(i,j) (1)
在奇数行偶数列位置的蓝色和红色分量的插值分别为：
Figure FDA0003809813740000011
Figure FDA0003809813740000012
在偶数行奇数列位置的蓝色和红色分量的插值分别为：
Figure FDA0003809813740000013
Figure FDA0003809813740000014
当单色图像坐标(i,j)处颜色为蓝色时，其蓝色分量为：
B(i,j)＝d(i,j) (6)
其绿色分量的插值由水平和竖直方向的梯度大小决定，设α为该点水平方向的梯度，β为其竖直方向的梯度：
当α＜β时，表明d(i,j)与其水平方向上的像素点关系紧密，插值将在水平方向上进行，公式如下所示：
Figure FDA0003809813740000021
当α＞β时，表明d(i,j)与其竖直方向上的像素点关系紧密，插值将在竖直方向上进行，公式如下所示：
Figure FDA0003809813740000022
当α＝β时，表明d(i,j)与其水平方向上和竖直方向上的像素点关系同等紧密，插值将在整个邻域内进行，公式如下所示：
Figure FDA0003809813740000023
其红色分量的插值为：
Figure FDA0003809813740000024
当单色图像坐标(i,j)处颜色为红色时，其红色分量为：
R(i,j)＝d(i,j) (11)
其绿色分量的插值由水平和竖直方向的梯度大小决定，设α为该点水平方向的梯度，β为其竖直方向的梯度：
当α＜β时，表明d(i,j)与其水平方向上的像素点关系紧密，插值将在水平方向上进行，公式如下所示：
Figure FDA0003809813740000025
当α＞β时，表明d(i,j)与其竖直方向上的像素点关系紧密，插值将在竖直方向上进行，公式如下所示：
Figure FDA0003809813740000026
当α＝β时，表明d(i,j)与其水平方向上和竖直方向上的像素点关系同等紧密，插值将在整个邻域内进行，公式如下所示：
Figure FDA0003809813740000031
其蓝色分量的插值为：
Figure FDA0003809813740000032
从上到下，从左到右扫描整幅单色图像，结合公式(1)-(15)，完成整幅图像的彩色插值。
2.如权利要求1所述的图像插值方法，其特征在于：所述第二步中水平方向的梯度α和竖直方向的梯度β通过二阶梯度公式获取：
Figure FDA0003809813740000033
CN201710598476.3A 2017-07-21 2017-07-21 一种基于梯度的改进彩色图像插值方法 Active CN109285113B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710598476.3A CN109285113B (zh) 2017-07-21 2017-07-21 一种基于梯度的改进彩色图像插值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710598476.3A CN109285113B (zh) 2017-07-21 2017-07-21 一种基于梯度的改进彩色图像插值方法

Publications (2)

Publication Number Publication Date
CN109285113A CN109285113A (zh) 2019-01-29
CN109285113B true CN109285113B (zh) 2022-11-22

Family

ID=65185582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710598476.3A Active CN109285113B (zh) 2017-07-21 2017-07-21 一种基于梯度的改进彩色图像插值方法

Country Status (1)

Country Link
CN (1) CN109285113B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113115010B (zh) * 2020-01-11 2023-04-18 晋城三赢精密电子有限公司 基于rgb-ir图像传感器的图像格式转换方法及装置
CN113824936B (zh) * 2021-09-23 2024-02-09 合肥埃科光电科技股份有限公司 色彩滤波阵列线扫描相机色彩插值方法、装置及设备
CN113824935B (zh) * 2021-09-23 2024-04-09 合肥埃科光电科技股份有限公司 伪彩色线扫描相机的时间延迟积分方法、装置及设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104537625A (zh) * 2015-01-05 2015-04-22 中国科学院光电技术研究所 一种基于方向标志位的Bayer彩色图像插值方法
CN105577981A (zh) * 2015-12-22 2016-05-11 深圳大学 一种边缘自适应的色彩还原方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7502505B2 (en) * 2004-03-15 2009-03-10 Microsoft Corporation High-quality gradient-corrected linear interpolation for demosaicing of color images

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104537625A (zh) * 2015-01-05 2015-04-22 中国科学院光电技术研究所 一种基于方向标志位的Bayer彩色图像插值方法
CN105577981A (zh) * 2015-12-22 2016-05-11 深圳大学 一种边缘自适应的色彩还原方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于彩色图像插值算法的胶囊机器人微视觉";李杰等;《计算机测量与控制》;20140225;第22卷(第2期);全文 *

Also Published As

Publication number Publication date
CN109285113A (zh) 2019-01-29

Similar Documents

Publication Publication Date Title
US20190206953A1 (en) Display apparatus, electroluminescent display panel and method of acquiring and displaying image by display apparatus
WO2015124005A1 (zh) 像素阵列的驱动方法和显示装置
CN103327220B (zh) 低照度Bayer图像上以绿色通道为导向的去噪方法
WO2017045286A1 (zh) 一种像素渲染方法、像素渲染装置和显示器
US10044952B2 (en) Adaptive shading correction
CN100366053C (zh) 一种对数字图像进行坏点补偿的方法
CN109285113B (zh) 一种基于梯度的改进彩色图像插值方法
CN102959959B (zh) 固态成像装置驱动方法、固态成像装置及成像设备
WO2017101191A1 (zh) 像素渲染方法、像素渲染装置及显示装置
TWI547169B (zh) 影像處理方法與模組
TWI449027B (zh) 適應性畫素補償方法
CN104618703A (zh) 一种白平衡调整方法
CN105141838A (zh) 去马赛克方法以及使用该方法的装置
CN108122201A (zh) 一种拜耳插值拉链效应去除方法
TW201314623A (zh) 影像處理模組及影像處理方法
CN105185352A (zh) 图像的边缘修饰方法和边缘修饰装置
US9858889B2 (en) Color compensation circuit, display apparatus, and color compensation method
CN101639932B (zh) 一种数字图像分辨率增强的方法及系统
CN104461441A (zh) 渲染方法、渲染装置及显示装置
US20150055861A1 (en) Methods and Systems for Image Demosaicing
US20160277721A1 (en) Color filtered area processing method for improving image processing
CN101150658A (zh) 一种坏点自适应的网格噪声消除装置及方法
CN103747222B (zh) 一种图像假色抑制的方法及装置
TW201616481A (zh) 對顯示器的畫面進行渲染的方法
US10140962B2 (en) Display unit, display panel and driving method thereof, and display device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant