CN106204481A - 一种新型图像色彩恒常性的光照估计方法 - Google Patents

一种新型图像色彩恒常性的光照估计方法 Download PDF

Info

Publication number
CN106204481A
CN106204481A CN201610537416.6A CN201610537416A CN106204481A CN 106204481 A CN106204481 A CN 106204481A CN 201610537416 A CN201610537416 A CN 201610537416A CN 106204481 A CN106204481 A CN 106204481A
Authority
CN
China
Prior art keywords
image
derivative
reflectance
point
light source
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.)
Granted
Application number
CN201610537416.6A
Other languages
English (en)
Other versions
CN106204481B (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.)
Hefei Kuxin Microelectronics Co ltd
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201610537416.6A priority Critical patent/CN106204481B/zh
Publication of CN106204481A publication Critical patent/CN106204481A/zh
Application granted granted Critical
Publication of CN106204481B publication Critical patent/CN106204481B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10024Color image

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明涉及一种新型图像色彩恒常性的光照估计方法。它是基于图像表面各点反射率的四阶导数计算,将图像等分成20个模块,对每个模块中反射率四阶导数值最大的点进行记录,利用每个模块中的最值,分别对相应模块各点四阶导数进行局部归一化处理,大量实验表明原图像各点反射率四阶导数总和与经过归一化处理过后的各点值的总和之比为常数,然后利用该线性规律推测出未知光源,结合Von Kries对角模型将图像颜色校正到标准白光下。本发明为图像色彩恒常算法提供新的方向,适用于计算机视觉和数字图像处理技术领域。

Description

一种新型图像色彩恒常性的光照估计方法
技术领域
本发明涉及计算机视觉和数字图像处理技术领域,具体涉及一种新型图像色彩恒常性的光照估计方法,用于解决色彩还原和色调处理的一系列问题。
背景技术
在人类的视觉构成中,视网膜上物体的影像会随着物体位置、大小、颜色的变化而随之改变。但是,人类不会因为某物体这些外在因素的变化就认为这是不同的物体。人类的眼睛能够稳定地感知到这个物体的大小、颜色、形状,这种现象被称为知觉恒常性,其中色彩恒常性为其中的一种。
研究色彩恒常性的目的就是要消除光照强度、光照颜色等一些外因素对物体图像的颜色、亮度的影响,获得一种与外界条件无关、只与物体本身元素有关的稳定的色彩描述。这对一个稳定的计算机视觉系统来说意义重大。
发明内容
本发明的目的在于针对现有技术的不足,提出一种新型图像色彩恒常性的光照估计方法,对未知光照估计具有较好的准确率,能够实现较好的色彩还原效果。
为达到上述目的,本发明的构思是:
本发明方法基于朗伯特反射模型和图像底层数据,对于给定的图像和光源,计算出图像表面反射率的四阶导数,将图像进行分割,然后将各部分四阶导数进行局部归一化处理,最后将未经局部归一化处理的图像四阶导数值的总和与归一化处理后图像各点值的总和进行对比,由两者之间的线性规律对未知光源进行预测,测试表明此方法具有良好的色彩恒常效果。
上述的方法中,对给定图像和光源的反射率四阶导数进行计算,同时对该图像进行等区间分割处理,求取反射率四阶导数最大值,用该最大值对本区间其他各点进行归一化处理,其他区域相同方式处理。
上述的方法中,利用图像反射率四阶导数总和与归一化后反射率四阶导数总和之比的线性规律,结合现有的色彩恒常性算法,做出目标光源的准确估计计算,并将图像矫正到标准白光下。
根据上述发明构思,本发明采用下述技术方案:
一种新型图像色彩恒常性的光照估计方法,其特征在于:具体操作步骤如下:
1)采集图像:获取给定光源照射下相机采集图像fc(x)。
2)计算图像反射率:计算图像中每个像素点的反射率,并计算反射率对该点坐标的四阶导数。
3)归一化处理:原图像等分成20块,每一块求出反射率四阶导数值最大的点,分别对各自模块局部归一化处理。
4)获得线性规律:利用图像各点反射率四阶导数总和与归一化处理后各点的值总和之比为常数k。
5)求出未知光源:结合上述线性规律,求出未知光源,经过对角模型转换矫正图像。
根据上述新型图像色彩恒常性的光照估计方法,所述步骤1)采集图像:根据给定特殊光源对某个场景进行照射,此光源一般不同于自然光的光源,利用相机对此场景进行拍摄获取所需待处理的原图像fc(x)。
根据上述新型图像色彩恒常性的光照估计方法,所述步骤2)计算图像反射率:计算出图像表面空间坐标为x点的反射率R(x),并进一步将反射率对x求四阶导,计算出x点反射率的四阶导数R4 x(x),同理,对整个图像每个坐标点都可以进行反射率四阶导数的计算。
根据上述新型图像色彩恒常性的光照估计方法,所述步骤3)归一化处理:将整个原图像分成20个互不重叠且尺寸相同的模块,分别计算每个模块内各坐标点的表面反射率四阶导数,其中在第j个模块内,反射率四阶导数值最大的像素点的空间坐标值记为Xj,max,用点Xj,max的四阶导数值对第j个模块内的其他点反射率四阶导数进行局部归一化处理,即把第j个模块各个点反射率四阶导数除以Xj,max的四阶导数值,把处理完成后像素点x的值记为LNRDC(x),同理,对每个模块都进行类似j模块的归一化处理,把原图像各像素点表面反射率四阶导数的全局总和记为GSRDC,再把图像经上述归一化处理后各点的值LNRDC(x)的全局总和记为GSLNRDC
根据上述新型图像色彩恒常性的光照估计方法,所述步骤4)获得线性规律:大量实验表明GSRDC和GSLNRDC两者之间存在线性关系,分别针对R,G,B通道得出两者线性关系的比值为k。
根据上述新型图像色彩恒常性的光照估计方法,所述步骤5)求出未知光源:利用步骤4)中的比值k,预测未知光源:
∫ | f 4 x ( x ) | d x ∫ 1 J ∫ ζ | f 4 x ( x ) | | f 4 x ( x j , max ) | d x d j = e c ∫ Ω | R 4 x ( x ) | d x ∫ 1 J ∫ ζ e c · | R 4 x , j ( x ) | e c · | R 4 x , j ( x j , max ) | d x d j = ke c
式中f(x)为相机成像模型下采集的原图像,f4 x(x)表示f(x)对x求四阶导数,ec为未知光源,得到图像的光源信息后,利用Von Kries对角模型将图像颜色校正到标准白光下。
本发明与现有技术相对比,具有如下显而易见的突出实质性特点及显著技术进步:
1)上述方法处理过程中,既考虑了图像的全局反射特性又结合了图像的局部信息,将整个原图像中的有效信息充分融合利用;
2)上述方法中,对光照估计的过程中,充分利用了原图像的边缘结构特征,同时结合亮度信息,可以更加有效地对光照进行估计;
3)上述方法中,通过对原图像多信息的融合处理,光照估计结果的误差小,比以往的方法更加准确。
附图说明
图1:本发明方法的流程框图
图2:本发明方法的处理照片图
图3:本发明方法的处理效果照片图
具体实施方式
本发明的优选实例结合附图详述如下:
实施例一:
参见图1,本新型图像色彩恒常性的光照估计方法,其特征在于:具体操作步骤如下:
1)采集图像:获取给定光源照射下相机采集图像fc(x)。
2)计算图像反射率:计算图像中每个像素点的反射率,并计算反射率对该点坐标的四阶导数。
3)归一化处理:原图像等分成20块,每一块求出反射率四阶导数值最大的点,分别对各自模块局部归一化处理。
4)获得线性规律:利用图像各点反射率四阶导数总和与归一化处理后各点的值总和之比为常数k。
5)求出未知光源:结合上述线性规律,求出未知光源,经过对角模型转换矫正图像。
实施例二:
本实施例与实施例一基本相同,特别之处如下:
所述步骤1)采集图像:根据给定特殊光源对某个场景进行照射,此光源一般不同于自然光的光源,利用相机对此场景进行拍摄获取所需待处理的原图像fc(x)。
所述步骤2)计算图像反射率:计算出图像表面空间坐标为x点的反射率R(x),并进一步将反射率对x求四阶导,计算出x点反射率的四阶导数R4 x(x),同理,对整个图像每个坐标点都可以进行反射率四阶导数的计算。
所述步骤3)归一化处理:将整个原图像分成20个互不重叠且尺寸相同的模块,分别计算每个模块内各坐标点的表面反射率四阶导数,其中在第j个模块内,反射率四阶导数值最大的像素点的空间坐标值记为Xj,max,用点Xj,max的四阶导数值对第j个模块内的其他点反射率四阶导数进行局部归一化处理,即把第j个模块各个点反射率四阶导数除以Xj,max的四阶导数值,把处理完成后像素点x的值记为LNRDC(x),同理,对每个模块都进行类似j模块的归一化处理,把原图像各像素点表面反射率四阶导数的全局总和记为GSRDC,再把图像经上述归一化处理后各点的值LNRDC(x)的全局总和记为GSLNRDC
所述步骤4)获得线性规律:大量实验表明GSRDC和GSLNRDC两者之间存在线性关系,分别针对R,G,B通道得出两者线性关系的比值为k。
所述步骤5)求出未知光源:利用步骤4)中的比值k,预测未知光源:
∫ | f 4 x ( x ) | d x ∫ 1 J ∫ ζ | f 4 x ( x ) | | f 4 x ( x j , max ) | d x d j = e c ∫ Ω | R 4 x ( x ) | d x ∫ 1 J ∫ ζ e c · | R 4 x , j ( x ) | e c · | R 4 x , j ( x j , max ) | d x d j = ke c
式中f(x)为相机成像模型下采集的原图像,f4 x(x)表示f(x)对x求四阶导数,ec为未知光源,得到图像的光源信息后,利用Von Kries对角模型将图像颜色校正到标准白光下。
实施例三:
新型图像色彩恒常性的光照估计方法如下:
首先根据给定光源对图像的照射情况,得到所需处理的原图像,如图2(a)所示,计算出图像表面空间坐标为x点的反射率R(x),并进一步将反射率对x求四阶导,计算出x点反射率的四阶导数R4 x(x),同理,对整个图像每个坐标点都可以进行反射率四阶导数的计算。
然后把整个原图像分成20个互不重叠且尺寸相同的模块,如图2(b)所示,分别计算每个模块内各坐标点的表面反射率四阶导数,其中在第j个模块内,反射率四阶导数值最大的像素点的空间坐标值记为Xj,max
基于上一步的求解,用点Xj,max的四阶导数值对第j个模块内的其他点反射率四阶导数进行局部归一化处理,即把第j个模块各个点反射率四阶导数除以Xj,max的四阶导数值,把处理完成后像素点x的值记为LNRDC(x),同理,对每个模块都进行类似j模块的归一化处理。
紧接着,把原图像各像素点表面反射率四阶导数的全局总和记为GSRDC,再把图像经上述归一化处理后各点的值LNRDC(x)的全局总和记为GSLNRDC。大量实验表明GSRDC和GSLNRDC两者之间存在线性关系,分别针对R,G,B通道得出两者线性关系的比值为k。
利用上述线性比值k,对未知光源进行预测:
∫ | f 4 x ( x ) | d x ∫ 1 J ∫ ζ | f 4 x ( x ) | | f 4 x ( x j , max ) | d x d j = e c ∫ Ω | R 4 x ( x ) | d x ∫ 1 J ∫ ζ e c · | R 4 x , j ( x ) | e c · | R 4 x , j ( x j , max ) | d x d j = ke c
式中f(x)为相机成像模型下采集的原图像,f4 x(x)表示f(x)对x求四阶导数,ec为未知光源,得到图像的光源信息后,利用Von Kries对角模型将图像颜色校正到标准白光下,达到良好的实践效果,如图3所示。

Claims (6)

1.一种新型图像色彩恒常性的光照估计方法,其特征在于:具体操作步骤如下:
1)采集图像:获取给定光源照射下相机采集图像fc(x);
2)计算图像反射率:计算图像中每个像素点的反射率,并计算反射率对该点坐标的四阶导数;
3)归一化处理:原图像等分成20块,每一块求出反射率四阶导数值最大的点,分别对各自模块局部归一化处理;
4)获得线性规律:利用图像各点反射率四阶导数总和与归一化处理后各点的值总和之比为常数k;
5)求出未知光源:结合上述线性规律,求出未知光源,经过对角模型转换矫正图像。
2.根据权利要求1所述的新型图像色彩恒常性的光照估计方法,其特征在于:所述步骤1)采集图像:根据给定特殊光源对某个场景进行照射,此光源一般不同于自然光的光源,利用相机对此场景进行拍摄获取所需待处理的原图像fc(x)。
3.根据权利要求1所述的新型图像色彩恒常性的光照估计方法,其特征在于:所述步骤2)计算图像反射率:计算出图像表面空间坐标为x点的反射率R(x),并进一步将反射率对x求四阶导,计算出x点反射率的四阶导数R4 x(x),同理,对整个图像每个坐标点都可以进行反射率四阶导数的计算。
4.根据权利要求1所述的新型图像色彩恒常性的光照估计方法,其特征在于:所述步骤3)归一化处理:将整个原图像分成20个互不重叠且尺寸相同的模块,分别计算每个模块内各坐标点的表面反射率四阶导数,其中在第j个模块内,反射率四阶导数值最大的像素点的空间坐标值记为Xj,max,用点Xj,max的四阶导数值对第j个模块内的其他点反射率四阶导数进行局部归一化处理,即把第j个模块各个点反射率四阶导数除以Xj,max的四阶导数值,把处理完成后像素点x的值记为LNRDC(x),同理,对每个模块都进行类似j模块的归一化处理,把原图像各像素点表面反射率四阶导数的全局总和记为GSRDC,再把图像经上述归一化处理后各点的值LNRDC(x)的全局总和记为GSLNRDC
5.根据权利要求4所述的新型图像色彩恒常性的光照估计方法,其特征在于:所述步骤4)获得线性规律:大量实验表明GSRDC和GSLNRDC两者之间存在线性关系,分别针对R,G,B通道得出两者线性关系的比值为k。
6.根据权利要求5所述的新型图像色彩恒常性的光照估计方法,其特征在于:所述步骤5)求出未知光源:利用比值k,预测未知光源:
∫ | f 4 x ( x ) | d x ∫ 1 J ∫ ζ | f 4 x ( x ) | | f 4 x ( x j , max ) | d x d j = e c ∫ Ω | R 4 x ( x ) | d x ∫ 1 J ∫ ζ e c · | R 4 x , j ( x ) | e c · | R 4 x , j ( x j , max ) | d x d j = ke c
式中f(x)为相机成像模型下采集的原图像,f4 x(x)表示f(x)对x求四阶导数,ec为未知光源,得到图像的光源信息后,利用Von Kries对角模型将图像颜色校正到标准白光下。
CN201610537416.6A 2016-07-10 2016-07-10 一种图像色彩恒常性的光照估计方法 Active CN106204481B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610537416.6A CN106204481B (zh) 2016-07-10 2016-07-10 一种图像色彩恒常性的光照估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610537416.6A CN106204481B (zh) 2016-07-10 2016-07-10 一种图像色彩恒常性的光照估计方法

Publications (2)

Publication Number Publication Date
CN106204481A true CN106204481A (zh) 2016-12-07
CN106204481B CN106204481B (zh) 2019-06-25

Family

ID=57474273

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610537416.6A Active CN106204481B (zh) 2016-07-10 2016-07-10 一种图像色彩恒常性的光照估计方法

Country Status (1)

Country Link
CN (1) CN106204481B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090021809A1 (en) * 2007-07-20 2009-01-22 Rajeev Ramanath System and Method for Maintaining Hue Constancy
CN101489022A (zh) * 2008-01-16 2009-07-22 延世大学工业学术合作社 色彩恢复方法和系统
CN101706964A (zh) * 2009-08-27 2010-05-12 北京交通大学 一种基于图像导数结构的颜色恒常性计算方法及系统
CN102073995A (zh) * 2010-12-30 2011-05-25 上海交通大学 基于纹理金字塔与正则化局部回归的色彩恒常方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090021809A1 (en) * 2007-07-20 2009-01-22 Rajeev Ramanath System and Method for Maintaining Hue Constancy
CN101489022A (zh) * 2008-01-16 2009-07-22 延世大学工业学术合作社 色彩恢复方法和系统
CN101706964A (zh) * 2009-08-27 2010-05-12 北京交通大学 一种基于图像导数结构的颜色恒常性计算方法及系统
CN102073995A (zh) * 2010-12-30 2011-05-25 上海交通大学 基于纹理金字塔与正则化局部回归的色彩恒常方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ARJAN GIJSENIJ ET AL.: "Color Constancy for Multiple Light Sources", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
JOOST VAN DE WEIJER ET AL.: "Edge-Based Color Constancy", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
杜馨瑜 等: "基于图像导数框架和非负稀疏编码的颜色恒常计算方法", 《电子学报》 *
董俊鹏 等: "基于改进图像导数框架的颜色恒常性计算方法", 《计算机工程与应用》 *
董俊鹏: "基于光照分析的颜色恒常性算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Also Published As

Publication number Publication date
CN106204481B (zh) 2019-06-25

Similar Documents

Publication Publication Date Title
US11810271B2 (en) Domain specific image quality assessment
Berman et al. Air-light estimation using haze-lines
CN101706964B (zh) 一种基于图像导数结构的颜色恒常性计算方法及系统
US7885455B2 (en) Method of combining images of multiple resolutions to produce an enhanced active appearance model
CN101996407B (zh) 一种多相机颜色标定方法
CN103702015B (zh) 一种近红外条件下人脸图像采集系统的曝光控制方法
CN109146897B (zh) 口腔放射图像质量控制方法及装置
CN112884682B (zh) 一种基于匹配与融合的立体图像颜色校正方法及系统
CN107209931A (zh) 颜色校正装置和方法
CN110246580B (zh) 基于神经网络和随机森林的颅侧面影像分析方法和系统
CN108364292B (zh) 一种基于多幅视角图像的光照估计方法
CN102982334B (zh) 基于目标边缘特征与灰度相似性的稀疏视差获取方法
US9185270B2 (en) Ghost artifact detection and removal in HDR image creation using graph based selection of local reference
CN109253722A (zh) 融合语义分割的单目测距系统、方法、设备及存储介质
CN102073995A (zh) 基于纹理金字塔与正则化局部回归的色彩恒常方法
CN103808265A (zh) 油菜叶片及菌核病斑形态同步测量的方法、装置及系统
Eun et al. Oriented tooth localization for periapical dental X-ray images via convolutional neural network
CN104732534B (zh) 一种图像中显著目标的抠取方法及系统
CN105405101B (zh) X射线图像处理方法和装置
CN103637815B (zh) 乳腺自动曝光参考区域的确定方法及系统
CN106204481A (zh) 一种新型图像色彩恒常性的光照估计方法
CN104200469B (zh) 一种视觉智能数控系统的数据融合方法
WO2016058359A1 (zh) 生成三维图像的方法和装置
CN113506230B (zh) 一种基于机器视觉的光伏电站航拍影像匀光处理方法
CN110060212A (zh) 一种基于深度学习的多光谱光度立体表面法向恢复方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201118

Address after: Room 208, 234 Songhu Road, Yangpu District, Shanghai, 200082

Patentee after: SHANGHAI ARTOSYN MICROELECTRONIC Co.,Ltd.

Address before: 200444 Baoshan District Road, Shanghai, No. 99

Patentee before: Shanghai University

CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: Room 501, No.308 Songhu Road, Yangpu District, Shanghai 200082

Patentee after: SHANGHAI ARTOSYN MICROELECTRONIC Co.,Ltd.

Address before: Room 208, 234 Songhu Road, Yangpu District, Shanghai, 200082

Patentee before: SHANGHAI ARTOSYN MICROELECTRONIC Co.,Ltd.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 230088 Building B2, Phase 3, Hefei Innovation Industrial Park, Intersection of Jiangjunling Road and Wanshui Road, High-tech Zone, Hefei City, Anhui Province

Patentee after: Hefei Kuxin Microelectronics Co.,Ltd.

Country or region after: China

Address before: Room 501, No.308 Songhu Road, Yangpu District, Shanghai 200082

Patentee before: SHANGHAI ARTOSYN MICROELECTRONIC Co.,Ltd.

Country or region before: China