CN101511031B - 一种用于拜尔数据插值的方法 - Google Patents
一种用于拜尔数据插值的方法 Download PDFInfo
- Publication number
- CN101511031B CN101511031B CN2008101478636A CN200810147863A CN101511031B CN 101511031 B CN101511031 B CN 101511031B CN 2008101478636 A CN2008101478636 A CN 2008101478636A CN 200810147863 A CN200810147863 A CN 200810147863A CN 101511031 B CN101511031 B CN 101511031B
- Authority
- CN
- China
- Prior art keywords
- component
- pixel point
- current pixel
- disappearance
- edge
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Color Television Image Signal Generators (AREA)
Abstract
本发明公开了一种用于Bayer Data插值的方法,通过当前像素点所包含的R或B分量,在三个方向计算梯度值,在边缘方向两侧选择边缘差值最小的一侧包含G分量的像素点恢复当前像素点的G分量,解决了双线性插值所带来的边缘过渡不平滑等问题;根据边缘方向及当前像素点包含的R或B分量,对当前像素点所缺失的B或R分量进行恢复,可以尽量避免梯度颜色双线性插值由于判断误差所带来的颜色失真。
Description
技术领域
本发明涉及电子图像尤其是彩色图像处理技术领域,具体来讲,涉及一种用于拜尔数据插值的方法。
背景技术
信号传感相机可以检测到与照片中像素点对应位置的光线密度变化。在彩色信号传感相机系统中,色彩滤波镜阵列(CFA)仅能有规则的测量任何给定像素单元的其中一种颜色,这种规则颜色分布的图像称为拜尔数据,导致信号传感相机只能有规则获取图像中每个像素点的一种颜色,如图1所示。因此,需要用颜色插值技术来恢复图像中每个像素点的其他两种颜色。
目前相机系统可以产生红、绿、蓝三种颜色。用颜色插值技术将每个像素点只有一种颜色分量的图像恢复为三种颜色的像素点。大多数色彩滤波镜阵列(CFA)只有一种颜色信号有较高的频率,相较其他两种颜色信号采样也较多。在RGB图像中,这种信号为绿色(G)。绿色由于其较高的细节频率并且是人类视觉集中的信号通常被作为亮度信号。因此目前的颜色插值技术一般先将图像每个像素点的绿色(G)分量恢复后再恢复其他两种颜色分量。目前较为典型的插值方法有以下三种:
双线性插值技术:该技术首先恢复图像每个像素点的绿色(G)分量,其中缺失G分量的像素点的G分量通过旁边像素点双线性插值得到。然后根据就近原则恢复每个像素点的其它两种分量。
边缘检测插值技术:在图像边缘采用不同于平滑区域的插值方法。并且该插值方法根据边缘与预先存储的模板进行匹配确定。
梯度插值技术:1、恢复图像所有像素点的G分量。缺失G分量像素点,分别计算垂直以及水平方向梯度,用最小值方向上的拥有G分量的相邻像素点进行双线性插值恢复。2、在G分量恢复的基础上使用双线性插值恢复像素点的其他两个颜色分量。
最初的双线性插值技术由于只是简单的通过行相邻像素点进行恢复,因此在图像边缘部分存在过渡不平滑,出现块效应。针对该不足,方向性检测的插值技术增加了边缘检测的技术,并且通过模板比较在边缘部分使用不同的插值技术,但是该方法计算复杂。梯度插值技术针对以上不足,计算垂直以及水平方向的梯度,使用最小梯度方向的像素点恢复当前点。该技术虽然解决了计算复杂度以及对边缘方向问题有所处理,但是该技术只有G分量进行垂直与水平方向的梯度计算,然后在最小方向上进行双线性插值,由于G分量集中了图像的大部分信息,因此梯度误差插值导致颜色失真,而且其他两个分量在G分量的基础之上恢复,因此造成误差扩散。
发明内容
本发明的目的在于克服现有拜尔数据插值方法的不足,提供一种边缘过渡不平滑、颜色失真小的用于拜尔数据插值的方法。
为达到上述发明目的,本发明用于拜尔数据插值的方法,包括缺失G分量像素点的恢复和G分量像素点缺失的R和B分量的恢复,其特征在于,缺失G分量像素点的恢复包括以下步骤:
(1)、使用当前像素点所包含的R或B分量,在垂直方向,水平方向,两个对角线方向,计算梯度值,选择最小梯度值方向作为边缘方向;
(2)、在边缘方向两侧选择边缘差值最小的一侧包含G分量的像素点恢复当前像素点的G分量;
(3)、对当前像素点缺失B或R分量的恢复。
在本发明的一具体实施方式中,步骤(3)中,当前像素点缺失B或R分量的恢复,在步骤(1)的基础上,根据边缘方向及当前像素点包含的R或B分量,如果边缘方向为垂直或水平方向,则计算边缘方向两侧差值,选择最小的一侧,离当前像素点最近且包含颜色分量为当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复;如果边缘方向为对角线方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复。
在本发明的另一具体实施方式中,G分量像素点缺失的R和B分量的恢复为:G分量像素点缺失的R分量通过垂直或水平的两个R分量求平均得到,而G分量像素点缺失的B分量则通过水平或垂直的两个R分量求平均得到。
在本发明中,通过当前像素点所包含的R或B分量,在三个方向计算梯度值,在边缘方向两侧选择边缘差值最小的一侧包含G分量的像素点恢复当前像素点的G分量,解决了双线性插值所带来的边缘过渡不平滑等问题,并且计算简单、可行;同时,在判断纹理情况的前提下,当前像素点缺失B或R分量的恢复,在步骤(1)的基础上,根据边缘方向及当前像素点包含的R或B分量,如果边缘方向为垂直或水平方向,则计算边缘方向两侧差值,选择最小的一侧,离当前像素点最近且包含颜色分量为当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复;如果边缘方向为对角线方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复,可以尽量避免梯度颜色双线性插值由于判断误差所带来的颜色失真。
附图说明
图1是拜尔数据的一种排列方式图。
具体实施方式
为更好地理解本发明,下面结合附图和具体实施方式对本发明进行更为详细描述。在以下的描述中,当已有的现有技术的详细描述也许会淡化本发明的主题内容时,这些描述在这儿将被忽略。
图1是拜尔数据的一种排列方式图。本实施例结合图1给出的拜尔数据的排列方式进行说明。
一、缺失G分量像素点的恢复
(1)、使用当前像素点所包含的R或B分量,在垂直方向,水平方向,两个对角线方向,计算梯度值,选择最小梯度值方向作为边缘方向。
在本实施例中,用X表示当前像素点包含的颜色分量值,即R或B分量的值,计算多方向梯度值:
diff1=|(X(i-2),j+X(i+2),j)/2-Xi,j|(1)
diff2=|(Xi,(j-2)+Xi,(j+2))/2-Xi,j|(2)
diff3=|(X(i-2),(j-2)+X(i+2),(j+2))/2-Xi,j|(3)
diff4=|(X(i-2),(j+2)+X(i+2),(j-2))/2-Xi,j|(4)
其中i,j表示当前像素点坐标位置,X表示当前像素点包含的颜色分量值,diff1、diff2、diff3、diff4分别为垂直方向,水平方向,两个对角线方向的梯度值。
(2)、在边缘方向两侧选择边缘差值最小的一侧包含G分量的像素点恢复当前像素点的G分量。
在本实施例中:
1)、如果diff1最小,垂直方向作为边缘方向进行G分量恢复。
a)、如果|Xi,j-Xi,j-2|>|Xi,j-Xi,j+2|,即右侧边缘差值最小,当前像素点的G分量Gij
Gi,j=Gi,j+1(5)
b)、如果|Xi,j-Xi,j-2|≤|Xi,j-Xi,j+2|,即左侧边缘差值最小,当前像素点的G分量Gij
Gi,j=Gi,j-1(6)
2)、如果diff2最小,水平方向作为边缘方向进行G分量恢复。
a)、|Xi,j-Xi-2,j|>|Xi,j-Xi+2,j|,即下侧边缘差值最小,当前像素点的G分量G′i,j
Gi,j=Gi+1,j (7)
b)、|Xi,j-Xi-2,j|≤|Xi,j-Xi+2,j|,即上侧边缘差值最小,当前像素点的G分量Gi,j
Gi,j=Gi-1,j (8)
3)、如果diff3最小,斜下方向作为边缘方向进行G分量恢复,使用第2)的a)、b)步骤进行G分量恢复。
4)、如果diff4最小,斜上方向作为边缘方向进行G分量恢复,使用第1)的a)、b)步骤进行G分量恢复。
(3)、当前像素点缺失B或R分量的恢复
根据边缘方向及当前像素点包含的R或B分量,计算边缘方向两侧差值,选择最小的一侧,离当前像素点最近且包含颜色分量为当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复。
在本实施例,根据计算多方向梯度值,进行判断:
1)、如果diff1最小,垂直方向作为边缘方向,进行缺失B或R分量的恢复。
a)、|Xi,j-Xi,j-2|>|Xi,j-Xi,j+2|,即右侧边缘差值最小,缺失Y分量为:
Yi,j=(|Yi-1,j+1-Gi-1,j+1|+|Yi+1,j+1-Gi+1,j+1|)/2+Gi,j (9)
b)、|Xi,j-Xi,j-2|≤|Xi,j-Xi,j+2|,即左侧边缘差值最小,缺失Y分量为:
Yi,j=(|Yi-1,j-1-Gi-1,j-1|+|Yi+1,j-1-Gi+1,j-1|)/2+Gi,j (10)
2)、如果diff2最小,水平方向作为边缘方向,进行缺失B或R分量的恢复。
a)、|Xi,j-Xi-2,j|>|Xi,j-Xi+2,j|,即下侧边缘差值最小,缺失Y分量为:
Yi,j=(|Yi+1,j-1-Gi+1,j-1|+|Yi+1,j+1-Gi+1,j+1|)/2+Gi,j (11)
b)、|Xi,j-Xi-2,j|≤|Xi,j-Xi+2,j|,即上侧边缘差值最小,缺失Y分量为:
Yi,j=(|Yi-1,j-1-Gi-1,j-1|+|Yi-1,j+1-Gi-1,j+1|)/2+Gi,j (12)
3)、如果diff3最小,边缘方向为对角线方向的斜下方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复,即
Yi,j=(|Yi-1,j-1-Gi-1,j-1|+|Yi+1,j+1-Gi+1,j+1|)/2+Gi,j (13)
4)、如果diff4最小,边缘方向为对角线方向的斜上方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复,即
Yi,j=(|Yi+1,j-1-Gi+1,j-1|+|Yi-1,j+1-Gi-1,j+1|)/2+Gi,j (14)
其中X表示当前像素点的颜色分量,Y表示当前像素点缺失B或R颜色分量,i,j表示当前像素点的坐标位置。
二、G分量像素点缺失的R和B分量的恢复
G分量像素点缺失的R分量通过垂直或水平的两个R分量求平均得到,而G分量像素点缺失的B分量则通过水平或垂直的两个R分量求平均得到,即:
Yi,j=(Yi-1,j+Yi+1,j)/2或Yi,j=(Yi,j-1+Yi,j+1)/2(15)
步骤(1)实例
下面结合图1的具体排列方式进一步对本发明进行说明和描述,以便理解本发明的方法。
一、缺失G分量像素点的恢复
(1)、使用当前像素点所包含的R或B分量,在垂直方向,水平方向,两个对角线方向,计算梯度值,选择最小梯度值方向作为边缘方向。
在本实例中,以缺失G分量像素点R3,4进行说明。计算多方向梯度值:
diff1=|(R1,4+R5,4)/2-R3,4|
diff2=|(R3,2+X3,6)/2-R3,4|
diff3=|(R1,2+R5,6)/2-R3,4|
diff4=|(R1,6+R5,2)/2-R3,4|
(2)、在边缘方向两侧选择边缘差值最小的一侧包含G分量的像素点恢复当前像素点的G分量。
1)、如果diff1最小,垂直方向作为边缘方向进行G分量恢复。
a)、如果|R3,4-R3,2|>|R3,4-R3,6|,即右侧边缘差值最小,当前像素点的G分量G′3,4
G3,4=G3,5
b)、如果|R3,4-R3,2|≤|R3,4-R3,6|,即左侧边缘差值最小,当前像素点的G分量G3,4
G3,4=G3,3
2)、如果diff2最小,水平方向作为边缘方向进行G分量恢复。
a)、|R3,4-R1,4|>|R3,4-R5,4|,即下侧边缘差值最小,当前像素点的G分量Gi,jG3,4=G4,4
b)、|R3,4-R1,4|≤|R3,4-R5,4|,即上侧边缘差值最小,当前像素点的G分量Gi,jG3,4=G2,4
3)、如果diff3最小,斜下方向作为边缘方向进行G分量恢复,使用第2)的a、b步骤进行G分量恢复。
4)、如果diff4最小,斜上方向作为边缘方向进行G分量恢复,使用第1)的a、b步骤进行G分量恢复。
(3)、当前像素点缺失B分量的恢复。
根据计算多方向梯度值,进行判断:
1)、如果diff1最小,垂直方向作为边缘方向,进行缺失B分量的恢复。
a)、如果|R3,4-R3,2|>|R3,4-R3,6|,即右侧边缘差值最小,缺失B分量为:
B3,4=(|B2,5-G2,5|+|B4,5-G4,5|)/2+G3,4
b)、如果|R3,4-R3,2|≤|R3,4-R3,6|,即左侧边缘差值最小,缺失B分量为:
B3,4=(|B2,3-G2,3|+|B4,3-G4,3|)/2+G3,4
2)、如果diff2最小,水平方向作为边缘方向,进行缺失B分量的恢复。
a)、|R3,4-R1,4|>|R3,4-R5,4|,即下侧边缘差值最小,缺失B分量为:
B3,4=(|B4,3-G4,3|+|B4,5-G4,5|)/2+G3,4
b)、|R3,4-R1,4|≤|R3,4-R5,4|,即上侧边缘差值最小,缺失B分量为:
B3,4=(|B2,3-G2,3|+|B2,5-G2,5|)/2+G3,4
3)、如果diff3最小,边缘方向为对角线方向的斜下方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复,即
B3,4=(|B2,3-G2,3|+|B4,5-G4,5|)/2+G3,4
4)、如果diff4最小,边缘方向为对角线方向的斜上方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复,即
B3,4=(|B4,3-G4,3|+|B2,5-G2,5|)/2+G3,4
其中X表示当前像素点的颜色分量,Y表示当前像素点缺失颜色分量,i,j表示当前像素点的坐标位置。
二、G分量像素点缺失的R和B分量的恢复
G分量像素点缺失的R分量通过垂直或水平的两个R分量求平均得到,而G分量像素点缺失的B分量则通过水平或垂直的两个R分量求平均得到,在本实例中,以G3,3为例:
B3,3=(B2,3+B4,3)/2或R3,3=(R3,2+R3,4)/2
尽管上面对本发明说明性的具体实施方式进行了描述,但应当清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (3)
1.一种用于拜尔数据插值的方法,包括缺失G分量像素点的恢复和G分量像素点缺失的R和B分量的恢复,其特征在于,缺失G分量像素点的恢复包括以下步骤:
(1)使用当前像素点所包含的R或B分量,在垂直方向,水平方向,两个对角线方向,计算梯度值,选择最小梯度值方向作为边缘方向;
(2)在边缘方向两侧选择边缘差值最小的一侧包含G分量的像素点恢复当前像素点的G分量;
(3)对当前像素点缺失B或R分量的恢复。
2.根据权利要求1所述的用于拜尔数据插值的方法,其特征在于,步骤(3)所述的对当前像素点缺失B或R分量的恢复为:
在步骤(1)的基础上,根据边缘方向及当前像素点包含的R或B分量,如果边缘方向为垂直或水平方向,则计算边缘方向两侧差值,选择最小的一侧,离当前像素点最近且包含颜色分量为当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复;如果边缘方向为对角线方向,直接用该边缘方向上的当前像素点所缺失的B或R分量的两个像素点的B、G或R、G分量之差的绝对值平均后,加上当前像素点的G分量进行恢复。
3.根据权利要求1所述的用于拜尔数据插值的方法,其特征在于,所述的G分量像素点缺失的R和B分量的恢复为:
G分量像素点缺失的R分量通过垂直或水平的两个R分量求平均得到,而G分量像素点缺失的B分量则通过水平或垂直的两个R分量求平均得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101478636A CN101511031B (zh) | 2008-12-15 | 2008-12-15 | 一种用于拜尔数据插值的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101478636A CN101511031B (zh) | 2008-12-15 | 2008-12-15 | 一种用于拜尔数据插值的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101511031A CN101511031A (zh) | 2009-08-19 |
CN101511031B true CN101511031B (zh) | 2010-07-14 |
Family
ID=41003261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101478636A Expired - Fee Related CN101511031B (zh) | 2008-12-15 | 2008-12-15 | 一种用于拜尔数据插值的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101511031B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102713974B (zh) * | 2010-01-06 | 2015-09-30 | 日本电气株式会社 | 学习装置、识别装置、学习识别系统和学习识别装置 |
CN101917629B (zh) * | 2010-08-10 | 2012-03-07 | 浙江大学 | 一种基于绿色分量和色差空间的Bayer格式颜色插值方法 |
CN102625133B (zh) * | 2011-01-31 | 2014-11-05 | 英属开曼群岛商恒景科技股份有限公司 | 影像感测坏像素点的检测方法 |
CN102630019B (zh) * | 2012-03-27 | 2014-09-10 | 上海算芯微电子有限公司 | 去马赛克的方法和装置 |
CN102938843B (zh) * | 2012-11-22 | 2014-12-03 | 华为技术有限公司 | 图像处理方法、装置以及成像设备 |
CN103248907B (zh) * | 2013-04-22 | 2015-07-22 | 雷欧尼斯(北京)信息技术有限公司 | 一种三维图像还原方法 |
CN103313066B (zh) * | 2013-05-08 | 2015-07-08 | 华为技术有限公司 | 插值方法及装置 |
CN103716606B (zh) * | 2013-12-30 | 2015-10-14 | 上海富瀚微电子股份有限公司 | 一种Bayer域图像下采样的方法、装置及摄像设备 |
CN105577981A (zh) * | 2015-12-22 | 2016-05-11 | 深圳大学 | 一种边缘自适应的色彩还原方法及系统 |
CN105578160A (zh) * | 2015-12-23 | 2016-05-11 | 天津天地伟业数码科技有限公司 | 一种基于fpga平台的高清晰度去马赛克插值方法 |
CN105959657A (zh) * | 2016-05-27 | 2016-09-21 | 乐视控股(北京)有限公司 | 图像处理方法和系统 |
CN108024100A (zh) * | 2017-12-15 | 2018-05-11 | 上海交通大学 | 基于改进边缘导向的Bayer格式图像插值方法 |
CN111340693B (zh) * | 2019-08-29 | 2023-12-01 | 杭州海康慧影科技有限公司 | 一种Bayer图像的插值方法、装置及电子设备 |
CN112788310B (zh) * | 2019-11-05 | 2022-11-01 | 瑞昱半导体股份有限公司 | 影像处理方法及相关的影像处理电路 |
-
2008
- 2008-12-15 CN CN2008101478636A patent/CN101511031B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101511031A (zh) | 2009-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101511031B (zh) | 一种用于拜尔数据插值的方法 | |
KR100782812B1 (ko) | 에지 적응적 컬러 보간 방법 및 장치 | |
US7825965B2 (en) | Method and apparatus for interpolating missing colors in a color filter array | |
JP4054184B2 (ja) | 欠陥画素補正装置 | |
KR100818447B1 (ko) | 컬러 필터를 사용하여 검출된 영상의 컬러 보간 방법 | |
US10565681B2 (en) | System and method for image demosaicing | |
US7995840B2 (en) | Image processing apparatus, image processing method, and program for attaining image processing | |
US9111365B2 (en) | Edge-adaptive interpolation and noise filtering method, computer-readable recording medium, and portable terminal | |
US7653240B1 (en) | Color filter array and method | |
CN103347190B (zh) | 一种边相关的联合彩色去马赛克和放大方法 | |
US9030579B2 (en) | Image processing apparatus and control method that corrects a signal level of a defective pixel | |
US20130077858A1 (en) | Image processing module and image processing method | |
JP4253655B2 (ja) | デジタル撮影装置の色補間法 | |
CN100589581C (zh) | 一种色彩滤波阵列插值方法和装置 | |
US8606016B2 (en) | Edge detection apparatus and computing circuit employed in edge detection apparatus | |
US8045826B2 (en) | Detecting edges in a digital images | |
WO1997048231A1 (en) | Method and system for reconstructing missing chrominance values with interpolation for a single-sensor color imaging systems | |
CN115471420A (zh) | 图像处理装置、成像设备、方法、电子设备及存储介质 | |
CN103384334A (zh) | 图像处理装置、图像处理方法及程序 | |
KR101327790B1 (ko) | 영상 보간 방법 및 장치 | |
CN103313066A (zh) | 插值方法及装置 | |
CN116416126A (zh) | 一种Bayer图像插值算法中解决伪彩色问题的方法 | |
US11128842B2 (en) | Image processing method and associated image processing circuit | |
CN110544210B (zh) | 一种基于模糊边界插值的Bayer CFA图像恢复方法 | |
KR100776385B1 (ko) | 에지 향상과 노이즈 제거를 위한 컬러 보간 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100714 Termination date: 20151215 |
|
EXPY | Termination of patent right or utility model |