CN104103043B - 空间域色域映射方法 - Google Patents
空间域色域映射方法 Download PDFInfo
- Publication number
- CN104103043B CN104103043B CN201410317826.0A CN201410317826A CN104103043B CN 104103043 B CN104103043 B CN 104103043B CN 201410317826 A CN201410317826 A CN 201410317826A CN 104103043 B CN104103043 B CN 104103043B
- Authority
- CN
- China
- Prior art keywords
- image
- gamut
- mapping
- color
- value
- 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
Abstract
本发明提供一种空间域色域映射方法,其特征在于,包括以下步骤:建立源图像和目标设备与CIE LCH颜色空间的映射关系;在CIE LCH颜色空间将源图像映射到目标设备色域,得到初步色域映射图像;将源图像和初步色域映射图像进行分块处理;采用人眼对比度敏感函数分别对分块处理后的源图像和初步色域映射图像进行空间滤波,将源滤波图像的明度值和饱和度值分别减去初步色域映射滤波图像的明度值和饱和度值,得到误差图像;将初步色域映射图像的明度值和饱和度值分别与误差图像的明度值和饱和度值相加,从而得到修正图像;分别建立修正图像和目标设备与CIE LCH颜色空间的映射关系;以及将修正图像映射到目标设备色域,得到最终色域映射图像。
Description
技术领域
本发明涉及图像处理领域,具体涉及一种基于人眼补偿机制的空间域色域映射方法。
背景技术
在跨媒体颜色复制过程中,色域映射一直都是国际上颜色科学与图像复制领域的研究重点。色域是指一个技术系统所能表现的颜色范围,一般分为设备色域和图像色域两种。通常在色彩再现过程中,可能出现色域不兼容的问题。例如,可以在液晶显示器上显示的颜色,可能在印刷机或投影仪上无法表现,这样必然导致颜色信息的丢失。
为了确保在不同色域的目标设备上最佳地再现图像,研究人员提出了一系列的色域映射算法。现有的色域映射算法主要分为两种,第一种是点对点色域映射算法(Point-wise gamut mapping algorithms),第二种是空间域色域映射算法(Spatial gamutmapping algorithms)。
点对点色域映射算法主要分为设备到设备色域映射算法(Device-to-devicegamut mapping algorithms)和图像到设备色域映射算法(Image-to-device gamutmapping algorithms)。然而这些算法都没有考虑临近像素颜色信息的视觉影响,导致映射图像的视觉效果不佳。
空间域色域映射算法主要分为基于视网膜模型(Retinex-based)的色域映射算法、基于视觉感知最小(Image-difference minimization)的色域映射算法和基于频率机制(Frequency-based)的色域映射算法。前两种色域映射算法考虑到了人眼视觉机制,使得复制图像在视觉上忠实于原稿,但是需要大量的计算量,无法满足现实的商业应用。而基于频率机制的色域映射算法将图像的细节部分补偿到映射之后的图像中,方法简单,易于实现,具有很高的商业价值。但是这种算法仅仅是数学意义上的图像细节补偿,没有考虑到实际的人眼视觉系统。
目前,还没有一种考虑人眼视觉系统而计算量较小的色域映射算法。
发明内容
本发明是针对上述问题进行的,目的在于提供一种基于人眼视觉补偿机制、节省计算量和计算时间的空间域色域映射方法。
本发明为解决上述问题,采用了以下的技术方案:
本发明提供一种空间域色域映射方法,其特征在于,包括以下步骤:
分别对源图像和目标设备进行色域边界提取,并分别建立源图像和目标设备与CIE LCH颜色空间的映射关系;
在CIE LCH颜色空间将源图像映射到目标设备色域,得到初步色域映射图像;
将源图像和初步色域映射图像分别进行分块处理;
采用CIE LCH颜色空间下的人眼对比度敏感函数作为空间滤波器,分别对分块处理后的源图像和初步色域映射图像进行空间滤波,得到源滤波图像和初步色域映射滤波图像,将源滤波图像的明度值和饱和度值分别减去初步色域映射滤波图像的明度值和饱和度值,得到误差图像;
将初步色域映射图像的明度值与误差图像的明度值相加,得到修正图像的明度值,将初步色域映射图像的饱和度值与误差图像的饱和度值相加,得到修正图像的饱和度值,根据修正图像的明度值和修正图像的饱和度值对初步色域映射图像进行修正,从而得到修正图像;
分别对修正图像和目标设备进行色域边界提取,建立修正图像和目标设备与CIELCH颜色空间的映射关系;以及
在CIE LCH颜色空间将修正图像映射到目标设备色域,得到最终色域映射图像,
其中,人眼对比度敏感函数为:
CS=k1×[exp(a·f)-exp(b·f)]+k2[1-exp(c·f)],
式中CS是人眼对比度敏感度,f是空间频率,k1、k2、a、b和c是模型系数。
发明的作用与效果
根据本发明所提供的空间域色域映射方法,因为先将源图像映射到目标设备色域得到初步色域映射图像,再利用人眼对比度敏感函数进行空间滤波后的源滤波图像和初步色域映射滤波图像计算误差图像,然后将初步色域映射图像与误差图像相加得到的修正图像映射到目标设备色域,因此最终获得的色域映射图像考虑了人眼视觉系统的特点,充分保留了源图像的细节信息,而且该方法的计算量小,节省计算时间,方便用户使用。
附图说明
图1是本发明的空间域色域方法的流程图;
图2是实施例中多锚点色域映射方法的原理示意图;
图3是明度通道的人眼对比度敏感函数图;以及
图4是饱和度通道的人眼对比度敏感函数图。
具体实施方式
以下结合附图,对本发明提供的空间域色域映射方法作详细阐述。
<实施例>
图1是本发明的空间域色域映射方法的流程图。
如图1所示,本发明所提供的空间域色域映射方法包括以下步骤:
步骤S-01,分别对源图像和目标设备进行色域边界提取,并分别建立源图像和目标设备与CIE LCH颜色空间的映射关系。
步骤S-02,在CIE LCH颜色空间进行图像到目标设备的色域映射,得到初步色域映射图像。在本实施例中,采用多锚点的色域映射算法。
图2是实施例中多锚点色域映射方法的原理示意图。
如图2所示,源图像色域大于目标设备色域,多锚点的色域映射方法将源图像色域按比例压缩到目标设备色域。并且,该方法将目标设备色域分成三个区域,在每个区域采用不同的映射轨迹,分别将相应区域内的源图像色域映射到目标设备色域。
步骤S-03,对源图像和初步色域映射图像分别进行分块处理。在观察图像视场角范围内,视距为w,视角为θ,图像分辨率为RI,则所观察的图像区域宽(高)度为d=2w×tan(θ/2)英寸,其对应的像素行/列数为l=d×RI,根据l的大小对图像进行分块处理。
步骤S-04,采用CIE LCH颜色空间下的人眼对比度敏感函数(CSF)作为空间滤波器,对分块后的源图像和初步色域映射图像进行空间滤波。人眼对比度敏感函数如下:
CS=k1×[exp(a·f)-exp(b·f)]+k2[1-exp(c·f)],
式中CS是人眼对比度敏感度,f是空间频率,k1、k2、a、b和c是模型系数。
图3是明度通道的人眼对比度敏感函数图。
如图3所示,明度通道的CSF对应的模型系数为k1=145.0,k2=-9.0,a=-0.11,b=-1.30,c=-0.08。
图4是饱和度通道的人眼对比度敏感函数图。
如图4所示,在不同的色相角下,饱和度通道的CSF的模型系数各不相同,不同色相角下的模型系数如表1所示。
表1不同色相角下的CSF模型系数
确定视距、视角等观察条件后,对于分块后的源图像和初步色域映射图像的明度值应用统一的CSF空间滤波;对于分块后的源图像和初步色域映射图像的饱和度,首先求出每块图像的色相角均值,对于每块图像,采用与色相角均值相对应的CSF进行滤波,得到源滤波图像和初步色域映射滤波图像。
将将源滤波图像的明度值减去初步色域映射滤波图像的明度值,得到误差图像的明度值;将源滤波图像的饱和度值减去初步色域映射滤波图像的饱和度值,得到误差图像的饱和度值。
步骤S-05,将初步色域映射图像的明度值与误差图像的明度值相加,得到修正图像的明度值;将初步色域映射图像的饱和度值与误差图像的饱和度值相加,得到修正图像的饱和度值,从而得到修正图像。
步骤S-06,对修正图像和目标设备进行色域边界提取,分别建立修正图像和目标设备与CIE LCH颜色空间的映射关系。
步骤S-07,在CIE LCH颜色空间将修正图像映射到目标设备,得到最终色域映射图像,映射方法与步骤S-02相同。
实施例的作用与效果
根据本实施例所提供的空间域色域映射方法,因为先将源图像映射到目标设备色域得到初步色域映射图像,再利用人眼对比度敏感函数进行空间滤波后的源滤波图像和初步色域映射滤波图像计算误差图像,然后将初步色域映射图像与误差图像相加得到的修正图像映射到目标设备色域,因此最终获得的色域映射图像考虑了人眼视觉系统的特点,充分保留了源图像的细节信息和视觉效果,而且该算法的计算量小,节省计算时间,方便用户使用。
以上实施例仅为本发明的优选方案,并不用于限定本发明的保护范围。
Claims (1)
1.一种空间域色域映射方法,其特征在于,包括以下步骤:
分别对源图像和目标设备进行色域边界提取,并分别建立所述源图像和所述目标设备与CIE LCH颜色空间的映射关系;
在CIE LCH颜色空间将所述源图像映射到所述目标设备色域,得到初步色域映射图像;
将所述源图像和所述初步色域映射图像分别进行分块处理;
采用CIE LCH颜色空间下的人眼对比度敏感函数作为空间滤波器,分别对分块处理后的所述源图像和所述初步色域映射图像进行空间滤波,得到源滤波图像和初步色域映射滤波图像,将所述源滤波图像的明度值和饱和度值分别减去所述初步色域映射滤波图像的明度值和饱和度值,得到误差图像;
将所述初步色域映射图像的明度值与所述误差图像的明度值相加,得到修正图像的明度值,将所述初步色域映射图像的饱和度值与所述误差图像的饱和度值相加,得到修正图像的饱和度值,根据所述修正图像的明度值和所述修正图像的饱和度值对所述初步色域映射图像进行修正,从而得到修正图像;
分别对所述修正图像和所述目标设备进行色域边界提取,建立所述修正图像和所述目标设备与CIE LCH颜色空间的映射关系;以及
在CIE LCH颜色空间将所述修正图像映射到目标设备色域,得到最终色域映射图像,
其中,所述人眼对比度敏感函数为:
CS=k1×[exp(a·f)-exp(b·f)]+k2[1-exp(c·f)],
式中CS是人眼对比度敏感度,f是空间频率,k1、k2、a、b和c是模型系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410317826.0A CN104103043B (zh) | 2014-07-04 | 2014-07-04 | 空间域色域映射方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410317826.0A CN104103043B (zh) | 2014-07-04 | 2014-07-04 | 空间域色域映射方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104103043A CN104103043A (zh) | 2014-10-15 |
CN104103043B true CN104103043B (zh) | 2017-02-15 |
Family
ID=51671170
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410317826.0A Expired - Fee Related CN104103043B (zh) | 2014-07-04 | 2014-07-04 | 空间域色域映射方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104103043B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105118026A (zh) * | 2015-07-28 | 2015-12-02 | 小米科技有限责任公司 | 色域模式切换方法及装置 |
CN106408619B (zh) * | 2016-09-13 | 2019-04-26 | 齐鲁工业大学 | 一种基于光谱域来实现跨媒体颜色再现的方法 |
CN107133009B (zh) * | 2017-04-21 | 2020-03-06 | 上海出版印刷高等专科学校 | 一种显示设备色域边界的计算方法 |
CN107738512B (zh) * | 2017-07-19 | 2019-04-02 | 中山火炬职业技术学院 | 基于icc编辑的改进瓷砖喷墨印刷色彩对样方法 |
CN107680049B (zh) * | 2017-09-13 | 2020-09-29 | 齐鲁工业大学 | 一种保持色相恒定的跨媒体色域映射方法 |
CN107492365B (zh) * | 2017-09-27 | 2020-07-03 | 深圳市华星光电半导体显示技术有限公司 | 获得色域映射拟合函数的方法及装置 |
CN107680142B (zh) * | 2017-10-23 | 2019-12-24 | 深圳市华星光电半导体显示技术有限公司 | 改善域外色重叠映射的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1390333A (zh) * | 1999-11-15 | 2003-01-08 | 佳能株式会社 | 嵌入的全域映射算法 |
CN1753452A (zh) * | 2004-09-21 | 2006-03-29 | 三星电子株式会社 | 检测色域边界及使用其来映射色域的设备和方法 |
CN1909594A (zh) * | 2005-08-02 | 2007-02-07 | 佳能株式会社 | 颜色处理方法和设备 |
CN101409066A (zh) * | 2007-10-09 | 2009-04-15 | 三星电子株式会社 | 用于色域外色彩转换的选择性处理的系统和方法 |
CN101523480A (zh) * | 2006-10-12 | 2009-09-02 | 皇家飞利浦电子股份有限公司 | 颜色映射的方法 |
CN102857674A (zh) * | 2011-06-28 | 2013-01-02 | 柯尼卡美能达美国研究所有限公司 | 利用对扩展rgb空间的色调映射处理高动态范围图像的方法 |
-
2014
- 2014-07-04 CN CN201410317826.0A patent/CN104103043B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1390333A (zh) * | 1999-11-15 | 2003-01-08 | 佳能株式会社 | 嵌入的全域映射算法 |
CN1753452A (zh) * | 2004-09-21 | 2006-03-29 | 三星电子株式会社 | 检测色域边界及使用其来映射色域的设备和方法 |
CN1909594A (zh) * | 2005-08-02 | 2007-02-07 | 佳能株式会社 | 颜色处理方法和设备 |
CN101523480A (zh) * | 2006-10-12 | 2009-09-02 | 皇家飞利浦电子股份有限公司 | 颜色映射的方法 |
CN101409066A (zh) * | 2007-10-09 | 2009-04-15 | 三星电子株式会社 | 用于色域外色彩转换的选择性处理的系统和方法 |
CN102857674A (zh) * | 2011-06-28 | 2013-01-02 | 柯尼卡美能达美国研究所有限公司 | 利用对扩展rgb空间的色调映射处理高动态范围图像的方法 |
Non-Patent Citations (3)
Title |
---|
Gamut mapping to preserve spatial luminance variations;Raja Balasubramanian 等;《Journal of Imaging Science and Technology》;20010930;第45卷(第5期);436-443 * |
基于 色空间对比敏感度函数的彩色图像滤波方法;王楠 等;《中国印刷与包装研究》;20140205;第6卷(第1期);25-29 * |
基于高斯拉普拉斯算子的色域映射算法;顾轶凡 等;《包装工程》;20140510;第35卷(第9期);95-98,109 * |
Also Published As
Publication number | Publication date |
---|---|
CN104103043A (zh) | 2014-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104103043B (zh) | 空间域色域映射方法 | |
CN108604133B (zh) | 图像形成的改进 | |
US8462169B2 (en) | Method and system of immersive generation for two-dimension still image and factor dominating method, image content analysis method and scaling parameter prediction method for generating immersive sensation | |
TWI511559B (zh) | 影像處理方法 | |
CN104809694B (zh) | 一种数字图像处理方法和装置 | |
US20120069045A1 (en) | Drawn image sharing apparatus | |
EP2122579A1 (en) | Wide color gamut anaglyphs | |
US7031548B2 (en) | Method and apparatus for filtering noise from a digital image | |
Jang et al. | Adaptive color enhancement based on multi-scaled Retinex using local contrast of the input image | |
US20120113238A1 (en) | Drawn image sharing apparatus, drawn image sharing system, and drawn image sharing method | |
US10620417B2 (en) | Method for generating a reflection-reduced contrast image and corresponding device | |
JP2013013060A (ja) | 拡張rgb空間へのトーンマッピングを用いた高ダイナミックレンジ画像の処理方法 | |
US11328399B2 (en) | Method and apparatus, and storage medium for processing style image | |
US20140098120A1 (en) | Image Color Correction for an Image by a Graphical Processing Unit Looking Up Color Corrections in a Lookup Table | |
WO2022156129A1 (zh) | 图像处理方法、图像处理装置及计算机设备 | |
EP2953340B1 (en) | Image processing device, image processing method, program, and recording medium | |
CN108711142A (zh) | 图像处理方法及图像处理装置 | |
CN103942756B (zh) | 一种深度图后处理滤波的方法 | |
CN111915508A (zh) | 一种面向色觉障碍的图像纹理细节增强方法 | |
CN109636739B (zh) | 图像饱和度增强的细节处理方法及装置 | |
CN110782400A (zh) | 一种自适应的光照均匀实现方法和装置 | |
Kinoshita et al. | Hue-correction scheme considering CIEDE2000 for color-image enhancement including deep-learning-based algorithms | |
CN102510502A (zh) | 用于生成高动态范围图像的方法及系统 | |
CN105208362B (zh) | 基于灰度平衡原理的图像色偏自动校正方法 | |
CN109155845A (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 | ||
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: 20170215 Termination date: 20190704 |
|
CF01 | Termination of patent right due to non-payment of annual fee |