CN103931169B - 一种通用的色域映射及色彩管理方法 - Google Patents
一种通用的色域映射及色彩管理方法 Download PDFInfo
- Publication number
- CN103931169B CN103931169B CN201180042671.9A CN201180042671A CN103931169B CN 103931169 B CN103931169 B CN 103931169B CN 201180042671 A CN201180042671 A CN 201180042671A CN 103931169 B CN103931169 B CN 103931169B
- Authority
- CN
- China
- Prior art keywords
- color
- equation
- grey
- rgb
- liu shi
- 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
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/46—Colour picture communication systems
- H04N1/56—Processing of colour picture signals
- H04N1/60—Colour correction or control
- H04N1/6058—Reduction of colour to a range of reproducible colours, e.g. to ink- reproducible colour gamut
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
- G09G5/00—Control arrangements or circuits for visual indicators common to cathode-ray tube indicators and other visual indicators
- G09G5/02—Control arrangements or circuits for visual indicators common to cathode-ray tube indicators and other visual indicators characterised by the way in which colour is displayed
- G09G5/026—Control of mixing and/or overlay of colours in general
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
- G09G2320/00—Control of display operating conditions
- G09G2320/06—Adjustment of display parameters
- G09G2320/0666—Adjustment of display parameters for control of colour parameters, e.g. colour temperature
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
- G09G2320/00—Control of display operating conditions
- G09G2320/06—Adjustment of display parameters
- G09G2320/0673—Adjustment of display parameters for control of gamma adjustment, e.g. selecting another gamma curve
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
- G09G2340/00—Aspects of display data processing
- G09G2340/06—Colour space transformation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Image Processing (AREA)
Abstract
通用的色域映射及色彩管理方法,包括采用结构相同的标定色靶,采用统一的原理和方法进行颜色空间的坐标变换和设备之间的色域映射;设置通道基色量参数并使通道基色量、基准基色量和驱动参数之间形成可逆的幂函数关系;创建生成纯净灰色梯尺的方法,并以此生成灰核参数、灰色平衡幂函数;对图像的灰阶进行伽玛校正;创建了刘氏Dlxlyl特性文件连接颜色空间;创建能够同时完成色域映射、伽玛校正并能够用快速解析算法实施的颜色空间映射方法;创建把刘氏色域映射方程得到的基色量数值精确转换到驱动输入数值的方法。该方法排除了红移干扰、合成和分解颜色时三基色能够保持自身独立性和通道独立性、优先复制灰色成分、按照亮度独立原则传送亮度和色度信息,具有良好的通用性。
Description
一、技术领域
本发明是跨媒体传送图像信息和按感知意图精确重现彩色图像的一种实用的、全新的技术,主要应用对象是图像输入设备(例如数字照相机、家用数码摄像机、扫描仪、电视摄像机、PHOTOCD)、图像显示设备(例如CRT、PDP、LCD、LED显示器)和图像输出设备(例如彩色打印机、多色胶印机、远程图像传送、图像移动通讯、网络图像交换等)和用它们构成的彩色管理***、计算机图像***、多媒体电视***、图像发送和接受***等的设计制造领域,它为前述各类实用***的软件生产和硬件制造提供了一种通用的完全创新的技术支持途径。
二、背景技术
在已有的色彩管理***中,例行地的选择CIE LAB(CIECAM02)作为PCS(特性文件连接空间),然而,CIE LAB和CIECAM02等表示颜色的***仍然存在不可忽视的误差;而CMM(色彩管理模块)在将CIE XYZ转换成RGB或CMYK颜色空间的数据时,由于设备和媒体对X值、Y值和Z值的衰减并不一致,使得CIE XYZ与RGB或CMYK之间成为复杂的非线性关系,但是到目前为止解决这种非线性关系的方法是不统一、不精确的,其症结在于没有解决三基色配色时自身属性的独立性、通道的独立性、灰色成分的独立性,更没有考虑到基色的红移特性等技术问题,影响所及,人们不得不采用查找表方法解决实际遇到的困难,查找表方法不仅繁琐,而且转换结果缺乏唯一性;在彩色图像的再现意图方面,由于缺乏通用的色域映射方法,实现跨媒体色域映射是困难的,直到现在也没有真正达到“所见即所得”的目标。在这样的技术背景之下,本发明另辟蹊径,建立了一种通用的跨媒体的色域映射方法,它使得整个色彩管理***的主要环节能够采用统一的原理、统一的方法和统一的数学描述模型,在具体运用这些原理和方法解决实际问题时,把提高色域映***确性的措施放在建立特性化文件阶段,使得预测色彩的准确性得到保证,在实时工作时,则采用新的色域映射方法和不高于二次方的解析计算方法,这使得色域映射的精确性和生产效率同时得到兼顾;这种方法是透明的,过程是规律的,结果是唯一的,增加了企业参与的机会和充分发挥相关设备潜能的可能性。
对检索有用的背景技术:
(1)J.S.Arney,P.G.Engeldrum,and H.Zeng,An expanded Murray-Davis modelof Tone reproduction in halftone imaging,JOURNAL OF IMAGING SCIENCE ANDTECHNOLOGY,Volume 40,Number 8,November/December 1995.
(2)J.S.Arney and Charles D.Armey.Modeling the Yule-Nielsen HalftoneEffect,Volume 39,Number 3,May/June 1996.
(3)大田 登 著.色彩工学,刘中本译,西安交通大学出版社1997,p49-51.
(4)Roy S.Berns,Ricardo J.Motta and Mark E.Gorzynski,CRT ColorimetryPart I Theory and Practice,Wiley&Sons Ltd,Volume 18,Number 5,October 1983p299-314.
(5)Tian minbo and Yefeng,Principle and Technology of TFT LCD,SciencePress China 2010年p109.
(6)羊晶亮,电视原理,中国人民邮电出版社1988年6月,p100,p137,p241.
(7)童林夙 屠彦Daniel den Engelsen编著,数字电视显示技术,电子工业出版社2007年9月p220.
(8)廖宁放 石俊生 吴文敏 编著,数字图文图像颜色管理***概论,北京理工大学出版社2009年5月p360-364.
(9)Yang Fujia,原子物理学,中国高等教育出版社,2010年pg.29-108.
(10)Neugebauer H.E.J,Die Theoretischen Grundlagen desMehrfarbendruckes,Z.wiss.photogr36,p73-89(1937).
(11)Mahy,M.(1993).Calculation of color gamut based on the Neugebauermodel.Color Res.Appl.,22,365-74.
(12)James R.Huntsman,A SIMMPL E ALTERNATIVE TO THE NEUGEBAUEREQATIONS FOR COMPUTING DOT AREA AND DOT GAIN COLORIMETRICALLY,1991-5(2)TAGAProceedings,p173-195.
(13)Karl J.Heuberger Zhou Mo Jing,Color Transformation from RGB toCMY,UGRA-Bericht 105/1Juni 1991,Switzerland.
(14)Liu Shichang,Optische Phaenomene beim Rasterdruck und die genaueFarbtransformation von XYZ nach CMY,Die FARBE,BAND 46,2002/2003.HEFT 1-2,p3-29.
(15)Schlaepfer,K.and Widmer.E.(1993)which color gamut can be achievedin multicolor printing and in television?Proc.TAGA,1993,41-49.
(16)Kurt Schlaepfer,Classification of color gamuts of printingprocesses,UGRA-Bericht 73/12March 2000,Switzerland.
(17)E.Jensen,Grey Component Removal,prepress Bul.,1985.9/10.
(18)Jan Morovic and M.Ronnier Luo,Developing algorithms for universalcolor gamut mapping,COLOR IMAGING Vision and Technology,JOHN WILEY&SONS LTDp267-280.
(19)Liu Shicheng,Travis Liu,Liu Yun,Image Transmission and DisplayMethod Complying with Chromaticity and Visual Fidelity Principle,PCT/CN2012/000327.
(20)Adobe RGB(1998)Color Image Encoding,May 2005.
三、发明内容
本发明之所以能够把一个具有创新特点的、方法统一的色域映射***建立起来,就在于解决了如下关键性的问题:
(1)在用三基色复制色彩的整个流程中,始终保持基色属性的不变性是最为根本的问题,然而已有的Marry-Davis公式、Yule-Nielson公式、GOG色度预测模型等准确度都是不够的;
(2)在颜色合成空间内三基色在满足自身色相属性不变性的同时,每种基色在每个通道内还必须保持空间的独立性,但是已有的颜色转换方法和数学模型,例如Neugebauer方程和常用的色空间矩阵转换方程等都不能真正满足这样的要求;
(3)彩色图像复制应当优先确保灰色阶调复制是正确的,这涉及到灰色平衡、灰色平衡曲线设置、视觉适应和对媒体衰减效应的伽玛校正问题,但是已有的技术存在灰色和彩色复制互相牵制的缺陷,不能做到灰色阶调的优先独立复制,在色域的某个区域会产生明显的色彩复制误差;
(4)作为PCS特性文件连接空间应当是精确的,但是已有的特性文件连接空间并不是完全均匀的,本发明用新创的Dlxlyl特性文件连接空间解决了色域映射和颜色传送的准确性的问题;
(5)由于已有的技术不能解决转换模型的复杂度和算法效率之间的矛盾,因而不得不采用查找表方法解决问题,但是查找表方法难以做到图像映射的唯一性和色域映射的连贯性;随着时间的推移,当设备老化、特性偏移的时候,用户难以主动采取纠正措施;
(6)面对跨媒体色彩输入设备、输出设备的多样性,已有的色域映射方法缺乏统一的处理原则和方法,这使色域映射技术本身在原理方面就存在产生随机误差和***误差的因素。
本发明在逐一解决上述关键难点的基础上创立的一种全新的色彩映射技术,最终使被映射的颜色保持原有的色相、原有的色度坐标和按照感知意图映射的亮度。由于涉及到较多的非传统的新东西,为了降低理解的难度,下面首先逐项列出映射***中各类设备共用的基本方法和实现这些方法的数学模型,而后把这些基础性发明串联起来构成一个完整的色彩映射***。
声明:本说明书中出现的数学模型采用统一的命名和标示符号,为了避免不必要的重复解释,只在该数学符号首次出现的时候注明含义,敬请读者注意。
1.一种为输入、显示和输出设备所通用的、驱动数值也对应相等的色靶结构
用途和目的:本步骤主要解决***各类设备的统一标定问题。为标定输入设备、显示设备和输出设备生成一种具有共性的互相联通的色靶结构,在不同设备之间用色相上具有对应关系、驱动数值相同的颜色样本进行标定,其目的是让彩色管理***中的各类设备有一个统一的输入标准,获得连续性和继承性良好的映射效果(15)(16(18)。
色靶的结构和生成步骤:
给定基础输入数据:在为输入设备、显示设备和打印输出设备建立特性化文件时,用来实际测量颜色样本的三刺激值,色靶上的颜色样本包括:
(1)单色的三基色梯尺:从最小值0到最大值20之间把驱动数值设定为21级,即让梯级数i=21,对于显示器来说,就是让三基色[r,g,b]的21个驱动数据dri、dgi、dbi都是:0.00,12.75,25.50,38.25,51.00,63.75,76.50,89.25,102.00,114.75,127.50,140.25,153.00,165.75,178.50,191.25,204.00,216.75,229.50,242.25,255.00,在显示器上逐个显示驱动数值生成的颜色并进行实测;对于数码照相机和电视摄像机之类的标定色靶来说,按照与此同样的数据制作成相纸色靶;对于打印机和扫描仪色靶来说,就是让三基色[c,m,y]的21个驱动数据dci、dmi、dyi是:0%=0.00/255,5%=12.75/255,10%=25.50/255,15%=38.25/255,20%=51.00/255,25%=63.75/255,30%=76.50/255,35%=89.25/255,40%=102.00/255,45%=114.75/255,50%=127.50/255,55%=140.25/255,60%=153.00/255,65%=165.75/255,70%=178.50/255,75%=191.25/255,80%=204.00/255,85%=216.75/255,90%=229.50/255,95%=242.25/255,100%=255/255;
对于打印机色靶来说,还应加上一条21级的单色黑墨梯尺,其驱动变量用dki表示;
(2)二次颜色样本:对于显示器来说,它们是用(dr255+dg255)、(dr255+db255)、(dg255+db255)三对驱动数值显示的三个二次色;对于打印机和扫描仪来说,它们是用(dc100%+dm100%)、(dm100%+dy100%)、(dc100%+dy100%)合成的三个二次色;对CMYK四色打印来说还需要增加(dc100%+dk100%)、(dm100%+dk100%)、(dy100%+dk100%)合成的三个二次色;
(3)三次颜色样本:对于显示器来说,它们是用具有相同(dri+dgi+dbi)驱动数值、从0到255顺序显示的基于加色法的灰色样本序列;对于打印机或扫描仪来说,它们是用具有相同驱动数值(dci+dmi+dyi)、从0到100%顺序显示的灰色样本序列,虽然打印机和扫描仪色靶的结构设计是一样的,但制作扫描色靶和打印机色靶的材质和样本颜色的三刺激值并不对应相等,在利用色靶建立特性文件时,还必须分别实测样本色的三刺激值,扫描仪的标定色靶可以用相纸模拟制作;对CMYK四色打印来说还需要增加(dc100%+dm100%+dk100%)、(dm100%+dy100%+dk100%)、(dc100%+dy100%+dk100%)合成的三个三次色及一个用(dc100%+dm100%+dy100%+dk100%)生成的叠印色;请注意:字符d是表示梯尺驱动变量的参数,因为习惯上用CMYK表示打印设备的输入数值,用RGB表示显示和扫描类设备的驱动数值,所以我们在此约定dci、dmi、dyi、dki、dri、dgi、dbi分别和C、M、Y、K、R、G、B等效,前者表示变量,后者表示具体数据。
对于扫描仪来说,不仅需要用分光光度计实测样本色的CIEXYZ数值,而且还需要借助软件获取每个样本色的足够像素的RGB平均值,目的是为计算扫描色的CIEXYZ数值做好准备。
2.为输入设备、显示设备和输出设备所通用的一种保持基色通道独立性的方法
目的:在本发明的色域映射方法中这是一个关键步骤,这种方法能够保证参与配色的每种基色成分保持恒定的色相。保持基色的通道独立性是精确合成颜色和实施色域映射的基础,已有的技术,例如Marry-Davis公式、Yule-Nielson公式和Gain-Offset-Gamma公式都不能确保基色以恒定的色相参加配色,例如在计算网点面积时用Marry-Davis公式和Yule-Nielson公式算出的网点面积数值实际是在几何扩大和光学扩大之间游走的模糊量,本方法的目标在于改变这种状况(1)(2)(3)(6)。
方法:当驱动数值在0-255变化时,让生成的各个基色量都和单位基色量保持一致的色相,如果用数学方法描述这种方法,那么本发明是利用创建的刘氏基色嵌位方程及其衍生的刘氏基色量公式对被媒体干扰的基色的色相、彩度和亮度进行嵌位,从而使基色的色相具有独立性。
保持基色通道独立性的三个步骤:首先以减色基色嵌位方程及其基准基色量公式为例进行解说。
第一步,在基色色靶上用分光光度计实测打印样本色的三刺激值,设实测值为XYZ,同时也测得白点的三刺激值Xw、Yw、Zw和基色实地的三刺激值Xs、Ys、Zs,在此我们规定基色实地的三刺激值Xs、Ys、Zs是该基色的单位基色量,然后用如下所示的刘氏基色嵌位方程表示该样本色的嵌位基色量at、嵌位亮度Yt和色貌保持系数λ等三个变量之间的关系(19):
这个格式的嵌位方程对打印机色靶、扫描仪色靶、常白型显示器的色靶都适用,在不同的场合可以分别称为打印机基色嵌位方程、扫描仪基色嵌位方程和常白型显示器基色嵌位方程。
第二步,求解上面的刘氏基色嵌位方程,可以得到该样本色的嵌位亮度Yt,Yt的解析表达式如下所示,根据下式可以计算出样本色的嵌位亮度Yt的数值:
第三步,将上面算得的嵌位亮度Yt的数值代入刘氏基准基色量公式,算出样本色的基准基色量a,刘氏基准基色量公式如右所示:
由上面的公式可以看出,设单位基色量的数值为1,那么基色量(在印刷工业中俗称网点面积)是分布在0~1的数。在不同场合应用上面所述的刘氏基色嵌位方程和刘氏基准基色量公式时,应辨别清楚公式中参数的具体含义,例如,对于打印机和显示器来说,基色量a代表的分别是青c、品红m、黄y基色的基色量,对于显示器来说代表的分别是红r、绿g、蓝b基色的基色量,即:
c、m、y就是在印刷工业中所称的网点面积;如果是对常白型LCD、LED显示器进行标定计算,基准基色量a代表的则是红r、绿g、蓝b等基准基色量,电视机显示器所用的常黑型基色嵌位方程和基色量公式的格式参见我们的国际申请号为PCT/2011/000327的专利说明书(19)。
在为扫描仪建立特性文件时,不仅需要在扫描色靶上实际测量出三基色梯尺的三刺激值[Xi,Yi,Zi],而且还需要用软件的数据采集模块读取梯尺上的样本颜色的扫描三刺激值[Ri,Gi,Bi],因为[Xi,Yi,Zi]和[Ri,Gi,Bi]是在不同的颜色空间内对同一样本色的描述,所以只需把刘氏嵌位方程稍加改动,即把其中的字符X、Y、Z改写成R、B、G就可以了,改写后的方程具有如下所示的格式:
显然这是一种扫描仪特有的情况,为了便于区分,本发明把它称为RGB扫描嵌位方程,而把前者称为XYZ扫描嵌位方程,在对扫描样本色进行RGB-CMY色空间转换时,要使用这个方程。通过RGB扫描嵌位方程可以计算出样本色的嵌位亮度Gt和基准基色量a,计算公式如下:
本项发明所带来的利益:由基色嵌位方程算出的嵌位亮度Yt排除了‘红移’现象所带来的干扰,这种干扰对于计算基色量的正确性是有害的,但是用刘氏基色量公式确定的基色量a排除了这样的有害干扰,使得基准基色量a所表示的颜色和单位基色量保持一致的色相,同时该基色量所具有的三刺激值不再等于实测三刺激值X、Y、Z,而被修改成为X、Yt、Z。在下面的内容中,基准基色量将作为配色或者三基色配色方程的独立参数使用,它为三基色配色方法走向规范化、***化、精确化提供了一种准确方法。
3.确保基色在配色空间的三个通道内具有独立性的方法及刘氏配色方程
方法:在已有的技术中,用Neugebauer方程实施三基色配色是的最典型方法(10),但是这些方法在三通道中,基色参数的独立性很差,在通道间发生互相串扰,计算误差甚至超过10。从1937年至今,已提出过许多改良方案,但都达不到期望的精确度(11)(12)(13)(14)。时下流行的以插值逼近算法为基础的查找表方法,只能算作模拟方法,不能获得确定性的、准确的色域映射结果。本发明为每个基色在每个配色通道设置‘通道基色量’参数,而‘通道基色量’进一步可以精确地表示为基准基色量参数的函数,与此同时,基色量参数和生成该基色的驱动参数之间又存在精确的、可逆的转换关系,以此为基础构建的配色方程可以保持三基色在配色空间内参加配色的独立性。
刘氏配色方程的种类:
1)基于减色法复制的刘氏减色配色方程:在本发明中,刘氏减色配色方程有两个用途:一是为刘氏四色配色方程的基色cmy建立灰平衡方程式,二是作为扫描仪的颜色预测方程使用,为了区别起见,在前一种情况下使用刘氏减色配色方程时可称为三色打印配色方程,在后一种情况下可称为扫描颜色预测方程,刘氏减色配色方程的格式如下:
在上面的方程中,左端的XYZ表示待匹配颜色的三刺激值,三刺激值XwYwZw、XcYcZc、XmYmZm、XyYyZy、XrYrZr、XgYgZg、XbYbZb、XsYsZs分别表示在色靶上实测8个实地色(白、青、品红、黄、红、绿、蓝、三基色合成黑点)的三刺激值;应当注意,如果涉及的是扫描仪问题,那么就要在扫描仪色靶上实测者8个实地色的三刺激值,如果涉及的是印刷机问题,那么就要在印刷机色靶上实测这8个实地色的三刺激值,请读者随时注意,下面不再重述同类问题;变量yx、yy、yz、mx、my、mz、cx、cy、cz分别表示黄、品红、青的通道基色量。通道基色量参数cx、cy、cz是基准基色量c的函数,mx、my、mz是基准基色量m的函数,yx、yy、yz是基准基色量y的函数,函数格式如下:
请注意:基色量参数c、m和y又进一步是驱动参数dd、dm、dy的函数,即存在如下的函数关系:
对此三式求反解得:
2)刘氏四色配色方程:这种配色方程只是上述刘氏减色配色方程的扩展形式,对于打印机和印刷机来说,cmyk四色复制是标准复制方法,当k=0或者当k为已知数时,四色配色方程退化成普通的刘氏减色预测方程,又返回到标准的三基色复制工艺。刘氏四色配色方程的格式如下:
变量yx、yy、yz、mx、my、mz、cx、cy、cz的含义和函数格式与刘氏减色配色方程完全相同,即有:
基色量参数c、m和y又进一步是驱动参数dd、dm、dy的函数,即存在如下的函数关系:
对此三式求反解得:
方程中相关的色靶实测值也和三色打印配色方程中的一致,kdd表示黑墨的灰成分取代量,假定黑墨的基准基色量用k表示,那么本发明让kdd表示四色复制中的灰色成份取代量,后面将有详细的叙述。扩展式中虽然包含c、m、y、k四个参数,可是由于kdd是在解方程之前预先给定的数值,所以四色配色方程实际上仍然是只有三个变量c、m、y的静定的三基色方程,在用它作标定计算时,可用迭代方法解这个方程。
3)基于加色法复制的刘氏常白型显示器配色方程(计算机用):为了给常白型CRT、PDP、LCD、LED显示器建立特性文件,也就是说,为了获取常白型显示器灰色平衡幂函数式,就需要使用如下的基于加色法复制的常白型刘氏加色配色方程:
在上面的方程中:X、Y、Z表示待匹配颜色的三刺激值,XwYwZw、XkYkZk分别表示显示器白点和黑点的实测三刺激值;XrYrZr、XgYgZg、XbYbZb分别表示驱动数值dr、dg、db取最大值时红、绿、蓝基色的实测三刺激值;
Xc、Yc、Zc对应于G、B取最大值时,由(G+B)同时驱动所得到的青色的三刺激值;
Xm、Ym、Zm对应于R、B取最大值时,由(R+B)同时驱动所得到的品红色的三刺激值;
Xy、Yy、Zy对应于R、G取最大值时,由(R+G)同时驱动所得到的黄色的刺激值;
方程右端的变量参数rx、ry、rz、gx、gy、gz、bx、by、bz被命名为通道基色量,拿它们去匹配方程左端的X、Y、Z刺激值,从这个角度看,通道基色量具有‘通道独立的特性’。然而,通道基色量并不是简单变量,而分别是基色量r、g和b的函数,函数格式是:
请注意:基色量参数r、g和b又进一步是驱动参数dr、dg、db的函数,即存在如下的函数关系:
对此三式求反解得:
4)基于加色法复制的刘氏常黑型显示器配色方程:为了给常黑型电视显示器建立特性文件,也就是说,为了获取常黑型电视显示器灰色平衡幂函数式,就需要使用如下的基于加色法复制的常黑型刘氏加色配色方程:其格式如下所示,
在我们的国际申请号为PCT/2011/000327的专利说明书中,提供过一个与此格式相同的方程,但在当时被称为‘灰色标定方程’,常黑型显示器配色方程就是那里的‘灰色标定方程’。
5)刘氏RGB扫描分色方程:在为扫描仪建立特性文件时,要在RGB颜色空间为扫描仪的灰色平衡函数式作标定计算,需要使用如下所示的RGB扫描分色方程,RGB扫描分色方程是利用CCD获取的RGB三刺激值对被扫描梯尺的青、品红、黄三种基色量进行计算,为了和用刘氏减色配色方程预测扫描梯尺的XYZ问题相区别,在RGB扫描分色方程中使用c′、m′、y′表示青、品红、黄三种基色量参数,并用yx′、yy′、yz′、mx′、my′、mz′、cx′、cy′、cz′表示通道基色量。实际上,基于减色法复制的XYZ刘氏配色方程和RGB扫描分色方程就像是用千克和磅描述同一个物重那样,是在RGB、XYZ色空间用不同的坐标体系描述色靶上的青、品红和黄三种基色,为了便于区别这两种在不同色空间算出的基色量,本发明把在XYZ色空间算出的基色量c、m、y叫做扫描基色量,把在RGB色空间算出的基色量c′、m′、y′叫做孪生扫描基色量。刘氏RGB扫描分色方程的格式如下:
RGB扫描分色方程中的通道基色量是基准基色量的函数,其函数格式和XYZ扫描配色方程相同:
刘氏配色方程在本色域映射方法中的重要地位:为达到灰色平衡的目标,在为扫描仪、打印机、显示器建立特性化文件时,必须利用刘氏配色方程对它们的灰色平衡函数式进行特性化标定,以便获得灰色平衡函数式中的系数数值,以便把得到的数值放在特性文件之中;因为刘氏配色方程需要用迭代方法求解,并且得到的解是精确解,所以在建立特性化文件的阶段用刘氏配色方程对灰色平衡多项式或者幂函数式进行标定计算是非常合适的,这种方法也为后面使用刘氏色域映射方程来快速完成色域映射计算提供了途径。
4.特性化刘氏配色方程的方法
特性化的目的:观察刘氏配色方程中的通道基色量函数可知,通道基色量是基准基色量的幂函数,多项式中包含的常系数是需要通过特性化获得的数值。
实施步骤:对于上列四种刘氏配色方程来说,特性化的步骤是一致的:
1)为了获取三基色梯尺的三刺激值,分别在色靶梯尺的样本色上测量XYZ和RGB三刺激值,我们约定用下标o,p,qw,qk区分扫描仪、打印机、常白型计算机显示器、常黑型电视显示器的相关技术数据,那么可以得到15个三刺激值数组:
对偶的两种扫描仪数据:[Xoci,Yoci,Zoci],[Xomi,Yomi,Zomi],[Xoyi,Yoyi,Zoyi]和[Roci、Goci、Boci],[Romi、Gomi、Bomi],[Royi、Goyi、Boyi];
打印机类:[Xpci,Ypci,Zpci],[Xpmi,Ypmi,Zpmi],[Xpyi,Ypyi,Zpyi];
常白型计算机显示器类:[Xqwri,Yqwri,Zqwri],[Xqwgi,Yqwgi,Zqwgi],[Xqwbi,Yqwbi,Zqwbi];
常黑型电视显示器类:[Xqkri,Yqkri,Zqkri],[Xqkgi,Yqkgi,Zqkgi],[Xqkbi,Yqkbi,Zqkbi];
2)根据以上15组三刺激值数组,用基色嵌位亮度模型分别为扫描仪、打印机、常白型计算机显示器、常黑型电视显示器计算嵌位亮度数值[Ytoci,Ytomi,Ytoyi]和[Gtoci,Gtomi,Gtoyi],[Ytpci,Ytpmi,Ytpyi],[Ytqwri,Ytqwgi,Ytqwbi],[Ytokri,Ytqkgi,Ytqkbi];
3)将算得的嵌位亮度数值代入基准基色量公式计算基准基色量:
4)用下列模型为刘氏减色配色方程、刘氏四色配色方程和显示器配色方程计算通道基色量cxi、cyi,czi,mxi,myi,mzi,yxi,yyi,yzi:
上式中,cxi、cyi、czi、mxi、myi、mzi、yxi、yyi、yzi分别表示青、品红、黄油墨在X、Y、Z通道的通道基色量;
5)用下列模型为RGB扫描分色方程计算通道基色量[cxi′,cyi′,czi′],[mxi′,myi′,mzi′],[yxi′,yyi′,yzi′]:
上式中,[cxi′,cyi′,czi′]、[mxi′,myi′,mzi′]和[yxi′,yyi′,yzi′]分别表示青、品红、黄油墨在R、G、B通道的通道基色量;
6)用下列模型为常白型显示配色方程计算通道基色量[rxi,ryi,rzi],[gxi,gyi,gzi],[bxi,byi,bzi]:
上式中,函数[rxi,ryi,rzi]、[gxi,gyi,gzi],[bxi,byi,bzi]分别表示红、绿、蓝基色在X、Y、Z通道的通道基色量;
7)用下列模型为常黑型显示配色方程计算通道基色量[rxi,ryi,rzi],[gxi,gyi,gzi],[bxi,byi,bzi]:
上式中,函数[rxi,ryi,rzi]、[gxi,gyi,gzi]、[bxi,byi,bzi]分别表示红、绿、蓝基色在X、Y、Z通道的通道基色量。
8)用曲线拟合方法为基于减色法复制的、用于打印机或者扫描仪的刘氏减色配色方程构造通道基色量函数:将青、品红、黄基色的基准基色量数组[ci,mi,yi]分别与对应的通道网点面积率数组[cxi,cyi,czi]、[mxi,myi,mzi]、[yxi,yyi,yzi]进行曲线拟合,即得如下所示的通道基色量的函数表达式:
9)用曲线拟合方法为常白型显示配色方程构造通道基色量函数:将红、绿、蓝基色的基准基色量数组[ri,gi,bb]分别与对应的通道基色量数组[rxi,ryi,rzi]、[gxj,gyj,gzj]、[bxj,byj,bzj]进行曲线拟合,即得通道基色量的函数表达式,对于红基色来说,可得如下的通道基色量函数:
10)用曲线拟合方法为数字照相机或者电视摄像机配色方程构造通道基色量函数:(为常黑型显示配色方程构造通道基色量函数(详见我们的PCT/2011/000327专利说明书。):分别将红、绿、蓝基色的基准基色量数组[rj,gj,bj]分别与对应的通道基色量数组[rxj,ryj,rzj][gxj,gyj,gzj][bxj,byj,bzj]进行曲线拟合,即得如下所示的通道基色量的函数:
通过以上数据拟合步骤,获得了通道多项式中的所有系数值,特性化任务即告完成。
5.一种为扫描仪、打印机、常白型显示器和常黑型显示器生成纯净灰色梯尺的方法
目的:灰色阶调复制质量是决定彩色图像复制质量的首要质量指标,为扫描仪、打印机和显示器各生成一条纯净灰色梯尺的理想三刺激值并以此作为图像灰色成分复制的基础,是本发明的一个重要举措,所谓纯净灰色梯尺是指该梯尺的灰色三刺激值已经被排除了红移成分。
实施步骤:
1)为扫描仪、打印机和显示器生成纯净灰色梯尺:在前述为扫描仪、打印机制作的色靶上以及在显示器上逐一显示的梯尺样本色上,分别实测复色梯尺(即用等驱动数值量生成的近似灰色)的亮度值,分别得初始亮度数组[Yoai]、[Go′ai]、[Ypai]、[Yqwai]、[Yqkai],它们不是纯净的亮度数值;
2)将初始亮度数组[Go′ai]、[Yoai]、[Ypai]、[Yqwai]、[Yqkai]转换成初始密度数组[Doai]、[Do′ai]、[Dpai][Dqwai]、[Dqkai],即让:[Doai]=lg(Yow/Yoai),[Do′ai]=lg(Gow/Go′ai),[Dpai]=lg(Ypw/Ypai),[Dqwai]=lg(Ywq/Yqwai),[Dqkai]=lg(Ywq/Yqkai);
3)对初始密度数组[Doai]、[Do′ai]、[Dpai]、[Dqwi]、[Dqki]进行归一化处理得归一化的初始密度数组,用[Dobi]、[Do′bi]、[Dpbi]、[Dqwbi]、[Dqkbi]表示归一化的初始密度数组,则有:[Dobi]=[Doai]/Doamax,[Do′bi]=[Do′ai]/Do′amax,[Dpbi]=[Dpai]/Dpamax,[Dqwbi]=[Dqwai]/Dqwamax,[Dqkbi]=[Dqkai]/Dqkamax,式中的Doamax、Do′amax、Dpamax、Dqwamax、Dqkamax分别是数组Doai]、[Do′ai]、[Dpai]、[Dqwai]、[Dqkai]中的最大值;4)
4)将基色梯尺的驱动数值进行归一化处理,得归一化的驱动数组[di],字符d是表示梯尺驱动数值的参数,d既代表驱动参数CMYK也代表驱动参数RGB的输入值;
5)以归一化的驱动数组[di]为自变量数组,分别以归一化初始密度数组[Do′bi]、[Dobi]、[Dpbi]、[Dqwbi]、[Dqkbi]为因变量进行幂函数拟合,得灰色梯尺的归一化初始密度模型:
Dob=d^γo,Do′b=d^γo′,Dpb=d^γp,Dqwb=d^γqw,Dqkb=d^γqk;
6)分别将Dob、Do′b、Dpb、Dqb、Dqb进行去归一化处理,并用字符Do、Do′、Dp、Dqw、Dqk表示去归一化初始密度数组,可得:[Doi]=[doi]^γo×Doamax,[Do′]=[doi]^γo′×Do′amax,[Dpi]=[doi]^γp×Dpamax,[Dqwi]=[doi]^γqw×Dqwamax,[Dqki]=[doi]^γqk×Dqkamax;
请注意:虽然Do、Do′、Dp、Dqw、Dqk称为去归一化初始密度数组,但它们的数值和初始密度数组[Doai]、[Do′ai]、[Dpai]、[Dqwi]、[Dqki]并不相等,因为实施前面的数据拟合步骤的目的就是为了改造初始密度数组,起到正本清源的作用,Do和Do′、Dp、Dqw、Dqk才是我们所需要的扫描仪、打印机、常白型显示器和常黑型显示器的纯净灰色阶调密度,色靶梯尺的灰色阶调已经设定为i级,所以扫描仪、打印机、常白型显示器和常黑型显示器的纯净灰色密度数组分别是[Doi]和[Do′i],[Dpi],[Dqwi],[Dqki];
7)分别为扫描仪、打印机、常白型显示器、常黑型显示器生成理想灰色亮度数组[Yoi]和[Goi],[Ypi],[Yqwi],[Yqki],即让:[Yoi]=Ywo/(10^Doi)和[Goi]=Gwo/(10^Do′i),[Ypi]=Ywp/(10^Dpi),[Yqwi]=Yww/(10^Dqwi),[Yqki]=Ywk/(10^Dqki);
8)为扫描仪、打印机和显示器计算参考白点的色度坐标值:
把实测扫描色靶基底白色的三刺激值Xow、Yow、Zow作为扫描仪在XYZ色空间内白点的三刺激值、把实测扫描色靶基底白色的三刺激值Row、Gow、Bow作为扫描仪在RGB色空间内白点的三刺激值、把实测打印机色靶基底白色的三刺激值Xpw、Ypw、Zpw作为打印机在XYZ色空间内的白点三刺激值、把D65光源的色度坐标值作为显示器参考白点的三刺激值,可得三刺激值RowGowBow,设下标o、p、q分别表示扫描仪、打印机和显示器等四种设备,参考白点的色度坐标可用如下模型计算出来:
xow=Xow/(Xow+Yow+Zow),yow=Ycw/(Xow+Yow+Zow)
row=Row/(Row+Gow+Bow),gow=Gow/(Row+Gow+Bow)
xpw=Xpw/(Xpw+Ypw+Zpw),ypw=Ypw/(Xpw+Ypw+Zpw)
xqw=Xqk=0.3127,yqw=yqk=0.3290
9)利用7)中给出的纯净灰色亮度数组[Yoi]、[Goi]、[Ypi]、[Yqwi]、[Yqki]和色度坐标值xow、yow、row、gow、xpw、ypw、xqw、yqw、xqk、yqk为扫描仪、打印机、常白型显示器、常黑型显示器分别构造一条纯净灰色梯尺,每条纯净灰色梯尺的三刺激值如下:
扫描仪Yoi=Yoi,
Roi=(row/gow)·Goi,Goi=Goi,Boi=(1-row-gow)·Goi;
打印机Ypi=Ypi,
常白型显示器Yqwi=Yqwi,
常黑型显示器Yqki=Yqki,
6.一种特性化纯净灰色梯尺组分基色量的方法
目的:用刘氏配色方程对刘氏纯净灰色梯尺特性化是为了把纯净灰色梯尺的组分基色量算出来并表示成纯净灰色密度参数的函数。灰色平衡是决定彩色图像复制质量的重要指标,在灰色平衡的状态下得到的基色量函数式可以推算出任一颜色的‘灰核’组分。可以这样理解‘灰核’的含义,例如利用三基色c、m、y匹配了一个颜色,三个基色量并不相等,其中具有最小基准基色量的基色就是形成该颜色灰色成分的‘灰核’,灰核是三个基色量中的最小值,但是它和另外两个基色相结合,形成一个复合颜色的灰色核心。本发明借助‘灰核’的概念达到了把一个颜色快速、精确地分解成为三种基色成分的目的,同时也达到了优先复制灰色成分的目的。
方法步骤:
1)在扫描仪的RGB色空间为扫描仪建立灰色平衡幂函数式c′(Drgb)、g′(Drgb)、b′(Drgb),方法是:用RGB扫描分色方程对扫描仪的灰色梯尺数组[Roi,Goi,Boi]进行配色计算,得扫描仪灰色梯尺的基准基色量数组[ci′,mi′,yi′];再利用灰色梯尺亮度数组[Goi]计算灰色梯尺的密度数组序列[Drgb],方法是:将亮度数组[Goi]代入密度公式Drgb=lg(Gow/Goi)即可算得密度数组[Drgbi];最后以[Drgb]为自变量、分别以孪生基色量[c′i,m′i,y′i]为因变量进行曲线拟合,最终得扫描仪的灰色平衡幂函数式如下:
在后面的扫描仪分色嵌位方程中,孪生基色量cdd′、mdd′、ydd′将作为灰核使用。
2)为扫描仪建立从基色量c′,m′,y′向基色量c,m,y转换的c′m′y′-cmy幂函数式,方法是:用刘氏减色配色方程对数组[Xoi,Yoi,Zoi]进行配色计算,得到扫描仪在CMY色空间组成灰色梯尺的基准基色量数组[ci,mi,yi];然后分别以[ci′],[mi′],[yi′]为自变量、分别以[ci],[mi],[yi]为因变量进行曲线拟合,得到c′m′y′-cmy的幂函数转换式:
3)为计算机使用的常白型显示器和电视机使用的常黑型显示器建立灰色平衡基色量幂函数式:计算机用的常白型显示器与电视机使用的常黑型显示器所使用的配色方程仅存在微小的差别,前后两种情况分别加下标w或者下标k加以区分。
分别用常白型显示器的配色方程和常黑型显示器的配色方程对灰色梯尺数组[Xqwi,Yqwi,Zqwi],[Xqki,Yqki,Zqki]进行配色计算,得灰色梯尺的基色量数组[rqwi,gqwi,bqwi],[rqki,gqki,bqki],然后分别以纯净灰色密度数组[Dqwi],[Dqki]为自变量,分别以[rqwi,gqwi,bqwi],[rqki,gqki,bqki]为因变量进行曲线拟合,得常白型显示器和常黑型显示器灰色平衡基色量函数式如下:
求上式的反函数得:
rqw、gqw、bqw、rqk、gqk、bqk是根据纯净中性灰色算得的组分基色量,为了表明这一特殊性质,统称rqw、gqw、bqw和rqk、gqk、bqk为灰平衡基色量,其中有一个将被用作构成彩色的灰色核心,也就是说,在一个理想灰色的三个组分基色量中,只有用作灰色核心的那个组分基色量才被称为‘灰核’。
4)为三色打印机建立灰色平衡幂函数式:
用刘氏三色配色方程对灰色梯尺数组[Xpi,Ypi,Zpi]进行配色计算,得灰色梯尺的基色量数组[cpi,mpi,ypi],然后分别以纯净灰色密度数组[Dpi]为自变量,分别以[cpi,mpi,ypi]为因变量进行曲线拟合,得三色打印机灰色平衡基色量幂函数式如下:
求上式的反函数得:
上面的灰色平衡幂函数称为三色打印机灰色平衡幂函数式,cp,mp,yp是根据纯净中性灰色[Xpi,Ypi,Zpi]算得的组分基色量,用刘氏减色配色方程算出的cp,mp,yp并不作为‘灰核’使用,而是在刘氏四色映射方程中,为计算cmy的驱动输入值服务的。在本发明中,四色打印机的纯净中性灰色也是[Xpi,Ypi,Zpi],但是[Xpi,Ypi,Zpi]的组分基色量还包括灰色取代参数kp。
5)为四色打印机(或四色胶印机)建立灰色平衡多项式:
(1)决定灰成分取代量:首先在色靶上测量黑墨印刷梯尺,记录样本色的三刺激值;接着用嵌位亮度公式算出黑基色的嵌位亮度数组[Ytki];然后根据黑基色的嵌位亮度值Ytki用如下公式计算黑基色的基准基色量数组ki:即让ki=(Ywp-Ytki)/(Ywp-Ysk);
(2)以归一化驱动数组[dki]为自变量数组,以[ki]为因变量数组进行曲线拟合,得函数:k=dk^γk;
(3)最后,灰色成分取代量用kdd表示,让kdd=Q(dk^γk)n,Q是一个预先设定的比例常数,用来控制最大黑版量,n是根据黑版阶调长短决定的指数,只要给定n的数值,就可方便的调节黑色阶调的长度;请注意,灰色成分取代量kp是经过净化的黑色成分,称为纯净的灰色成分取代量;
(4)求解与[ki]相配套的纯净灰色的基色量数组[ci]、[mi]、[yi]:
将纯净灰色梯尺的i组三刺激值[Xpi、Ypi、Zpi]逐组放在四色配色方程的左端和同步将[kpi]逐一代入刘氏四色配色方程,那么四色配色方程蜕变成只有未知量c、m、y的静定的刘氏配色方程,与此同时,逐一用迭代法解方程,那么得到基准基色量数组[ci]、[mi]、[yi];
(5)以纯净灰色密度数组[Dpi]为自变量、分别以[cpi]、[mpi]、[ypi]、[kpi]为因变量进行曲线拟合,可得
四色打印的灰色平衡多项式,基色量cp、mp、yp及灰色成份取代量kp如下所示:
cp=a0+a1Dp+a2Dp 2+a3Dp 3+…,mp=b0+b1Dp+b2Dp 2+b3Dp 3+…
yp=c0+c1Dp+c2Dp 2+c3Dp 3+…,kp=d0+d1Dp+d2Dp 2+d3Dp 3+…
用上面函数求得的cp、mp、yp是根据视觉中性灰色算得的组分基色量,是在灰平衡条件下算得的基色量数据,为了表明这一特殊性质,统称cp、mp、yp为打印色空间的灰色平衡基色量,在打印色空间内,cp、mp、yp中必有一个灰色平衡基色量和kp一起形成匹配视觉彩色的灰色核心,这个灰色平衡基色量简称‘灰核’。
7.一种对图像的灰色阶调进行伽玛校正的方法
目的:由于再现图像的媒体和设备总是导致显示和打印图像的色调向偏暗的一端移动,为了校正这种弊病,必须对图像的色调进行总体性的伽玛校正,本发明采取了与现有技术完全不同的伽玛校正方法,它使得伽玛校正和色域映射得以同步进行(15)(16)(17)(18)(20)。
步骤:
1)把纯净灰色密度参数表示成为驱动参数的函数:第5项第6)段中给出的打印机、常白型显示器和常黑型显示器的纯净灰色密度数组分别是[Dpi]、[Dqwi]、[Dqki],将[Dpi]、[Dqwi]、[Dqki]进行归一化处理,处理之后仍然用[Dpi]、[Dqwi]、[Dqki]表示,再把它表示成归一化驱动参数[di]的函数:分别以[di]为自变量、分别以纯净灰色密度数值[Dpi]、[Dqwi]、[Dqki]为因变量进行曲线拟合,得纯净灰色密度的幂函数表达式如下:
2)把伽玛校正密度表示成纯净灰色密度的函数:即把校正后的密度Dp′、Dqw′、Dqk′分别表示成原始密度Dp、Dqw、Dqk的函数,即让
以上三式就是计算伽玛校正密度的刘氏伽玛校正公式。
8.在XYZ空间为扫描仪建立Dlxlyl特性文件连接颜色空间的方法:
目的:数码照相机和扫描仪等颜色输入设备,总是要把获取的颜色信息向显示器、打印机之类的输出设备传送。已有的方法是选择CIE LAB(CIECAM02)作为PCS特性文件连接空间,这种借用均匀颜色空间作为连接纽带的方法是正确的,问题是此类匀色空间还没有达到理想的均匀,作为PCS特性文件连接空间使用,导致的颜色转换误差还比较明显,而且步骤繁琐;再一个问题涉及颜色数值的传送方法,跨媒体远距离传送图像已非常普遍,移动通讯、数字电视和地空图像通讯都打破了一时一地进行色彩管理的格局。但是,传统的构成亮度信号和色差信号的方法在非线性的情况下并不满足恒亮度原则,明显影响电视图像细节的重现质量。针对这样的情况,本发明提出一种全新的、兼有多种功能的Dlxlyl特性文件连接空间,下面介绍在XYZ空间内为扫描仪建立Dlxlyl特性文件连接颜色空间的方法(15)(16)(17)(18)(20)。
方法:Dlxlyl颜色空间是利用刘氏颜色分割方程作为工具建立起来的。该方程的基本原理是:虽然任何一个颜色都是基于三基色合成原理生成的,但是改变一下观察问题的角度,也可以认为:任何一个颜色都是由两部分组成的,一是比例为p的灰色成分,二是比例为(1-p)的两个基色生成的间色成分,即具有[1个视感灰色成分+2个基色成分]的结构特点,前后两种思考方法是等效的,乍一看,后面的思考方式似乎使问题变得更加复杂,但是这种‘复杂’只是表面上的,它不仅能带来许多好处,而且可以使本来影响计算机工作效率的高阶算法变成为算法效率高的解析算法,能够达到事半功倍的效果。刘氏颜色分割方程可以划分成为扫描仪建立的刘氏扫描色分割方程和为数码相机、数码摄像机和电视摄像机建立的刘氏拍摄色分割方程两种,因为扫描仪获取的颜色要向打印设备色域传送,所以要用一个基于减色法复制的刘氏颜色分割方程。因为数码相机获取的颜色一般要向常白型显示器传送,所以要用一个基于常白型配色方程的刘氏颜色分割方程,电视摄像机获取的颜色要向常黑型显示器传送,所以要用一个基于常黑型配色方程的刘氏颜色分割方程。每种刘氏颜色分割方程都包含三个子类型。
1)在XYZ空间内为扫描仪建立Dlxlyl特性文件连接颜色空间的方法
第一步,为扫描仪构建刘氏颜色分割方程:
扫描仪用刘氏颜色分割方程的格式:三个子类型的方程式如下所示:
刘氏颜色分割方程是由三个基本方程组合起来工作的一组方程,三者之间的差别之处在于:
一个是[1个青基色+1个品红基色+1个视感灰色]的分割方程,
一个是[1个青基色+1个黄基色+1个视感灰色]的分割方程,
一个是[1个品红基色+1个黄基色+1个视感灰色]的分割方程,
为简便起见,可把这三种情况分别简记为CMK、CYK和MYK型分割方程,这三个方程貌似复杂,实际上是简单的二次方程式,可以用简单的解析方法进行快速计算。在上面的方程中:X、Y、Z是由扫描仪获取的三刺激值,Xw、Yw、Zw表示在扫描仪色靶上实测的白色实地的三刺激值,Xc、Yc、Zc表示青色实地的实测三刺激值,Xm、Ym、Zm表示品红色实地的实测三刺激值,Xy、Yy、Zy表示黄色实地的实测三刺激值,Xsk、Ysk、Zsk表示由青、品红和黄三基色实地合成的、被称为‘黑点’的三刺激值,并且和Xw、Yw、Zw的色度坐标是一致的。通过程序的判断,可以从三个分割方程中选出一个来执行颜色分割任务并算得XYZ中的灰色量p。
在刘氏分割方程中,让方程右端方括号中的数值等于在扫描仪色靶上实测的白色实地的三刺激值Xw、Yw、Zw代替,那么,刘氏分割方程蜕变成如下的刘氏分割方程的灰色梯尺格式:
Xv=Xw(1-p)+pXsk
Yv=Yw(1-p)+pYsk
Zv=Zw(1-p)+pZsk
如果将[pi]逐一代到刘氏分割方程灰色梯尺格式的右端,那么生成的三刺激值数组[Rvi、Gvi、Bvi]等效于纯净灰色梯尺的三刺激值[Roi、Goi、Boi];
第二步,利用上述知识构建Dlxlyl特性文件连接空间:
(1)根据刘氏分割方程左端的XYZ值计算输入颜色XYZ的色度坐标值:Dlxlyl特性文件连接空间所需要的xl、yl就是根据这组XYZ值、按照下式计算出来的:
(2)应用刘氏颜色分割方程对XYZ进行分割,即通过解刘氏颜色分割方程算出灰色量p的数值;
(3)将算得的p值代入刘氏分割方程灰色梯尺格式中的亮度式,计算颜色XlYlZl中灰阶的亮度值:
Yl=Yv=Yw(1-p)+pYsk;
(4)把灰阶的亮度值Yl换算成灰色密度值Dl,即:Dl=Do=lg(Ywo/Yl);
经过上面4个步骤,就把在XYZ色空间内的任一颜色XYZ或者扫描XYZ颜色变换成了特性文件连接空间的数值Doxoyo。
2)在RGB颜色空间为扫描仪建立Dlxlyl特性文件连接颜色空间的方法:
第一步,为扫描仪构建刘氏颜色分割方程:
在RGB色空间里扫描仪用刘氏颜色分割方程的格式:三个子类型的方程式如下所示:
在上面的方程中,让方程右端方括号中的数值等于在扫描仪色靶上实测的白色实地的三刺激值Rw、Gw、Bw,那么,刘氏分割方程演变成如下的刘氏分割方程灰色梯尺格式:
Rv=Rw·(1-p)+p·Rsk
Gv=Gw·(1-p)+p·Gsk
Bv=Bw·(1-p)+p·Bsk
如果将[pi]逐一代到刘氏分割方程灰色梯尺格式的左端,那么生成的三刺激值数组[Rvi、Gvi、Bvi]等效于纯净灰色梯尺的三刺激值[Roi、Goi、Boi];
第二步,计算灰色成分Rv、Gv、Bv的灰色密度值:
(1)应用刘氏颜色分割方程对RGB进行分割,即通过解刘氏颜色分割方程算出灰色量p的数值;
(2)将算得的p值代入刘氏分割方程灰色梯尺格式中的亮度式,计算颜色RvGvBv中灰阶的亮度值:
Gl=Gv=Gw(1-p)+pGsk;
(3)把灰阶的亮度值Gl换算成灰色密度值Drgb,即:Drgb=Dl=lg(Gwo/Gv);Dlxlyl特性文件连接颜色空间带来的利益:(1)从XYZ到Dlxlyl的转换是精确的,不会给后续的映射变换带来误差;(2)Dlxlyl特性文件连接颜色空间是以刘氏颜色分割方程为工具建立起来的,其实质是把一个用XYZ表示的颜色分解成为以Dl表示的亮度信号和以xlyl表示的色度信号,像电视、卫星通讯这样的场合,由于刘氏颜色分割方程p的特殊作用,可以实现恒亮度、恒色度传送彩色图像信号、可以在节省带宽的情况下以高压缩比无损传送彩色图像信号;(3)在色彩管理***中,这种Dlxlyl特性文件连接颜色空间是实现通用颜色映射方法的一条通途。
9.在RGB色空间快速计算基准基色量c′m′y′的方法及刘氏三色嵌位方程
目的:已有的把扫描的RGB数据转换到成CIEXYZ的方法是利用回归分析方法建立转换多项式,不仅转换误差较大,而且回归分析方法是一种数学逼近方法,不涉及色域映射这样的视觉心理学问题,并不遵循色域映射的规律。本发明从色度学入手,首先在RGB色空间把RGB数据转换成的基准色量c′m′y′,然后把基准基色量c′m′y′转换成为XYZ色空间的扫描基色量cmy,最后把得到的数据cmy代入扫描仪颜色预测方程,就可以把扫描得到的RGB数据准确地转换成CIEXYZ数据。在这几个步骤中,把RGB数据转换成的孪生基色量c′m′y′的方法是借助于刘氏三色嵌位方程完成的(6)(9)(14)。
刘氏三色嵌位方程的格式:刘氏三色嵌位方程继承了刘氏分割方程的类别属性,也分为CMK、CYK、MYK三种子类型,至于选用哪一个刘氏三色嵌位方程对颜色RGB进行转换计算,是由RGB中的最小值决定的,下面列出这三种子类型的刘氏三色嵌位方程:
上式中:λ是刘氏色貌保持参数,cdd′、mdd′、ydd′是来自于扫描仪灰色平衡幂函数式的灰核参数,如果选用CMK型刘氏三色嵌位方程对颜色RGB进行转换计算,那么只需计算灰核ydd的数值,然后把ydd动态地代入CMK型刘氏三色嵌位方程,使该方程变成只有未知量λ、c′、m′的二次方程,所以可以用解析算法快速求解。同样重要的是:由于刘氏色貌保持参数λ和灰核ydd的联合作用,CMK型三色嵌位方程能够使基色量c′、m′、y′=ydd获得基准基色量的特性。如果选用CYK和MYK型刘氏三色嵌位方程,那么方程就分别与灰核mdd和cdd配套使用。
10.把扫描颜色从RGB色空间转换到CIEXYZ色空间的方法
目的:扫描仪是通过CCD感像器件获取色靶图像反射的RGB数据,可以认为扫描仪的CCD感像器件对感受光量具有线性相应,然而被扫描的色靶产生的XYZ数据却是非线性的,借助于下面叙述的方法可以使扫描的RGB数据依据色度学原理完成RGB-CIEXYZ色空间的转换(13)(14)(15)(16)(17)(18)(19)(20)。
方法和步骤:
1)对扫描仪获取的R,G,B三刺激值进行判断,根据R,G,B中的最小数值选择对颜色RGB进行分割的方程格式。选择原则是:如果R是最小值,那么选用GBK型分割方程;如果G是最小值,那么选用RBK型分割方程;如果B是最小值,那么选用RGK型分割方程。为免赘述,下面仅以选用了GBK型分割方程的流程进行解说,另外两种情况依此类推;
2)将扫描仪获取的RGB三刺激值放在GBK型颜色分割方程左端,解方程后得到黑色量p的数值;
3)将黑色量p代入刘氏分割方程的灰色梯尺格式,可以算得灰色灰色亮度刺激值Gl,即有:
Gl=Gw(1-p)+p Gsk;
4)将Gl换算成密度值Drgb:即Drgb=lg(Gwo/Gl),式中的Gwo是在扫描色靶白色面上测得的亮度数值;
5)将Drgb代入扫描仪灰色平衡幂函数式,算得灰核cdd′的数值;
6)将算出的灰核数值cdd′代入GBK型刘氏三色嵌位方程,算得颜色RGB的组分基色量c′、m′、y′;
7)将算出的基色量c′、m′、y′代入扫描仪的c′m′y′-cmy转换式,算得基准基色量c、m、y的数值;
8)将算出的基准基色量c、m、y代进扫描仪配色方程,即可算得目标值XYZ。
11.为数码照相机(或者数码电视摄像机)建立Dlxlyl特性文件连接颜色空间的方法
目的:数码照相机和扫描仪等图像输入设备可以认为是由CCD感像器件的RGB-CIEXYZ坐标变换模块和处理色靶非线性反射的模块组成的,这样处理可以使问题变得准确、简单和符合映射需要(13)(14)(15)(16)(17)(18)(19)(20)。
方法步骤:
数码照相机和数码电视摄像机所用的CCD元件虽有区别,但是可以用一个格式相同的刘氏颜色分割方程生成Dlxlyl特性文件连接颜色空间,生成的步骤如下:
1)在综合考虑当前多媒体设备的普遍需要之后,可以为配色的RGB三基色制定一组普遍认可的色度坐标数值,让白点的色度坐标都等于基准白D65的色度坐标,建立实施RGB-XYZ坐标变换的矩阵方程式和亮度方程式如下所示,因为这种方法是对同一颜色在不同颜色度量***中所进行的几何线性变换,所以这样的变换结果是正确的:
将CCD感像元件获取的任意一组RGB数值代入上面的RGB-XYZ矩阵方程式就可以完成RGB-XYZ色空间坐标数据转换;
2)为标准显示器构造刘氏分割方程的灰色梯尺格式如下所示:
Xgreyi=Xsk·(1-Y)+Y·Xsw
Ygreyi=Ysk·(1-Y)+Y·Ysw
Zgreyi=Zsk·(1-Y)+Y·Zsw;
上式中:XgreyYgreyZgrey表示标准显示器纯净灰色的三刺激值,XwYwZw和XkYkZk分别是标准显示器白点和黑点的三刺激值。
3)把灰阶的亮度值Ygrey换算成灰色密度值Dl,即:Dl=Dn=lg(Yw/Ygrey);
4)根据RGB-XYZ矩阵方程式左端的XYZ值计算颜色XYZ的色度坐标值:
Dlxlyl特性文件连接空间所需要的xl、yl就是根据这组XYZ值、按照下式计算出来的:
xl=xn=X/(X+Y+Z),yl=yn=Y/(X+Y+Z);
经过上面4个步骤,就把扫描XYZ颜色变换成了特性文件连接空间的数值Dnxnyn。
12.一种将特性文件连接空间的参数Dl映射成为目标设备色域的灰色阶调参数方法
目的:建立把特性文件连接空间的参数Dl和目标色域的灰色阶调密度实现无缝连接的通用方法。
方法:在前面的工作中,已经得到数码相机(或数码摄像机和电视摄像机)、扫描仪、显示器和打印机(或印刷机)在设备色域内的灰色阶调密度数组[Dni]、[Doi]、[Dpi]、[Dqwi]、[Dqki],利用这些数组之间的映射关系,就可以得到以输入设备灰色阶调密度数组为自变量、以输出设备灰色阶调密度数组为因变量的密度映射函数,下面叙述这些密度映射函数(9)(14)(15)(16)(17)(18)(19)(20)。
1)将数码相机的灰色阶调密度数组[Dni]映射到常白型显示器(计算机显示器):以数码相机的灰色阶调密度数组[Dni]为自变量数组,以常白型显示器的设备灰色阶调密度数组[Dqwi]为因变量数组进行幂函数拟合,得数码相机-计算机显示器灰密度映射函数式:
2)将电视摄像机的灰色阶调密度数组[Dni]映射到常黑型显示器(电视机显示器)色域中去:以电视摄像机的灰色阶调密度数组[Dni]为自变量数组,以常黑型显示设备灰色阶调密度数组[Dqki]为因变量数组进行幂函数拟合,得电视摄像机-电视机显示器灰密度映射函数式:
3)将扫描仪的灰色阶调密度数组[Doi]映射到打印机(印刷机)设备色域中去:以扫描仪的的灰色阶调密度数组[Doi]为自变量数组,以打印设备灰色阶调密度数组[Dpi]为因变量数组进行幂函数拟合,得扫描仪-打印机(印刷机)灰密度映射函数式:
4)将打印机(印刷机)灰色阶调密度数组[Dpi]映射到显示设备色域中去:以打印机(印刷机)的灰色阶调密度数组[Dpi]为自变量数组,以常白型显示器灰色阶调密度数组[Dqwi]为因变量数组进行幂函数拟合,得打印机(印刷机)-显示器灰密度映射函数式:
13.一种将源设备获取的颜色映射到目标设备色域的方法
目的:建立一种把特性文件连接空间的颜色Dlxlyl无损传递到目标色域的方法(20)。
1)将扫描仪获取的颜色映射到打印设备色域的方法
方法:扫描仪和打印设备之间的色域映射是通过传递特性文件连接空间的Do、xo、yo的数值到打印设备色域和借助于刘氏扫描仪-打印机映射方程实现的。
步骤:
第一步,根据扫描得到的RGB数值,选择刘氏扫描仪-打印机映射方程的子类型:刘氏扫描仪-打印机映射方程具有CMK、CYK、MYK三个子类型,并且和刘氏分割方程的三个子类型一一对应的,下面依次列出这三个子类型的方程式:
为简明起见,下面仅以第一组MYK型方程为例解说。从Dlxlyl特性文件连接空间传来的扫描仪色度坐标数据xo和yo放在方程的左端,参数Yp是未知的亮度参数,Yp实际上具有嵌位亮度Yt=Yp的特性,称为刘氏嵌位亮度,方程左端表示经过扫描-打印映射方程处理之后,扫描输出色Doxoyo将被打印成一个新的颜色,如果这个新颜色是Xp、Yp、Zp,那么XpYpZp实际就是方程左端代表的数值,CMYK就是显示颜色[Xp,Yp,Zp]所需要的驱动数值。即有:Xp=(xo/yo)Yp、Yp=Yp、Zp=[(1-xo-yo)/yo]Yp。
第二步,把扫描仪色空间的灰色阶调密度Do映射成为打印设备色域的灰色阶调密度Dp:
(1)让(2)计算伽玛校正密度:让
(3)用伽玛校正密度Dp′取代下列打印机的灰色平衡多项式中的Dp,生成如下所示的刘氏四色打印反伽玛灰色平衡多项式,可用它计算伽玛校正后的灰色平衡组分色量cp、mp、yp和灰色取代参数kp,并且把cp、mp、yp、kp作为可以选用的灰核使用:
cp=cdd=a0+a1Dp′+a2Dp′2+a3Dp′3+…,mp=mdd=b0+b1Dp′+b2Dp′2+b3Dp′3+…
yp=ydd=c0+c1Dp′+c2Dp′2+c3Dp′3+…,kp=kdd=d0+d1Dp′+d2Dp′2+d3Dp′3+…
根据所选子类型的不同,只需计算kdd,cdd,或者kdd,mdd,或者kdd,ydd;
(4)根据程序选定的刘氏扫描仪-打印机映射方程子类型,将kdd和cdd(或者mdd、或者ydd)数值代入已选定的四色映射方程的右端;
由此可见,方程右边出现的两个新符号cdd、kdd,不是变量,而是由cmyk灰色平衡多项式动态传来的两个数值,对于MYK型方程来说,cdd=c,来源于打印机灰色平衡方程式。但灰色平衡多项式中的自变量最终变成了Dp,由此可见,来自Doxoyo特性文件连接空间的灰色量参数Do是通过Dp、cdd和刘氏扫描-打印映射方程发生联系的。kdd是被转换颜色的灰成分取代量,方程的变量是m、y和Yp,在Yp、m和y的协同作用下,匹配出方程左端的颜色XpYpZp。和刘氏基色嵌位方程的原理一样,变量Yp还起到色貌保持系数的作用,其功能是确保得到的新颜色XpYpZp和扫描仪输出的颜色XYZ保持视觉一致的色貌,通过MYK子类型刘氏扫描-打印映射方程计算出来的三个基色量c=cdd、m、y和cdd是经过伽玛校正的三个基准基色量。
第三步,根据刘氏扫描仪-打印机映射方程解得的三个基色量c=cdd、m、y计算生成颜色Xp、Yp、Zp所需要的驱动数值CMYK:
在方程中,基准基色量cmyk进一步是驱动输入数值dcdmdydk的函数,由此可得把打印机基准基色量cmyk数值转换成驱动输入值CMYK的方法是:
2)将数码相机获取的颜色XYZ映射到常白型显示器色域的方法
方法:数码相机和显示设备之间的色域映射是通过把特性文件连接空间的参数Dn、xn和yn传送到显示设备并通过数码相机-显示器映射方程完成。
步骤:
第一步,把数码相机色空间的灰色阶调密度Dn映射成为显示设备色域的灰色阶调密度Dq,即:(1)让(2)计算伽玛校正密度:让
(3)将伽玛校正密度Dqw′代入常白型显示器的灰平衡组分基色量幂函数式,计算作为灰核利用的灰色平衡组分基色量和灰色取代参数kp:
第二步,把输入设备特性文件连接空间的色度坐标xn和yn映射到常白型显示器并通过数码相机-显示器映射方程完成映射过程:
刘氏数码相机-显示器映射方程也具有rddgb、rgddb、rgbdd三个子类型,并且和刘氏分割方程的三个子类型是一一对应的,下面依次列出这三个子类型的方程式:
根据程序选定的rddgb、rgddb、rgbdd子类型,将第一步得到的rdd、gdd、bdd代入方程的右端,将xnyn和刘氏嵌位亮度Yqw放在方程的左端,即可解得经过伽玛校正和色域映射处理的基准基色量r、g、b和亮度值Yqw。
上式中,驱动输入数值drdgdb进一步是基准基色量rgb的函数,由此可得把常白型显示器基准基色量rgb数值转换成驱动输入值RGB的方法是:
通过刘氏映射方程三个不同子类型计算出来的三个基色量都是经过伽玛校正的基色量。
设Xqw=(xn/yn)Yl、Yqw=Yl、Zqw=[(1-xn-yn)/yn]Yl,则最终显示的颜色就是Xqw、Yqw、Zqw。
RGB就是显示颜色Xqw、Yqw、Zqw所需要的驱动数值。
3)将电视摄像机获取的颜色XYZ映射到常黑型电视显示器色域的方法
方法:电视摄像机和电视显示设备之间的色域映射是通过把特性文件连接空间的参数Dn、xn和yn传送到显示设备并通过电视摄像机-显示器映射方程完成。
步骤:
第一步,把电视摄像机色空间的灰色阶调密度Dn映射成为显示设备色域的灰色阶调密度Dqk,即:
(1)让(2)计算伽玛校正密度:让
(3)将伽玛校正密度Dqk′代入常黑型显示器的灰平衡组分基色量幂函数式,计算作为灰核利用的灰色平衡组分基色量和灰色取代参数kdd:
第二步,把输入设备特性文件连接空间的色度坐标xn和yn映射到电视显示器并用刘氏映射方程完成映射过程,刘氏四色色域映射方程也具有rddgb、rgddb、rgbdd三个子类型,并且和刘氏分割方程的三个子类型是一一对应的,下面依次列出这三个子类型的方程式:
根据程序选定的rddgb、rgddb、rgbdd子类型,将第一步得到的rdd、gdd、bdd代入方程的右端,将参数xn和yn的值和氏嵌位亮度Yqk放在方程的左端,即可解得经过伽玛校正和色域映射处理的基准基色量r、g、b和亮度值Yqk。
上式中,驱动输入数值drdgdb进一步是基准基色量rgb的函数,由此可得把电视显示器基准基色量rgb数值转换成驱动输入值RGB的方法是
通过刘氏综合映射方程三个不同子类型计算出来的三个基色量都是经过伽玛校正的基色量。
设Xqk=(xn/yn)Yn、Yqk=Yn、Zqk=[(1-xn-yn)/yn]Yn,那么,最终显示的颜色就是Xqk、Yqk、Zqk。RGB就是显示颜色Xqk、Yqk、Zqk所需要的驱动数值。
4)将打印设备获取的颜色XYZ映射到常白型显示器色域的方法
方法:从打印设备到显示设备之间实施色域映射是通过特性文件连接空间的密度参数Dp和目标色域的灰色密度参数Dpw的映射及色度坐标参数xp和yp的传递并通过打印机-显示器映射方程完成。
步骤:
第一步,把打印设备色空间的灰色阶调密度Dp映射成为显示设备色域的灰色阶调密度Dq,即:
(1)让
(2)计算伽玛校正密度:让
将伽玛校正密度Do′代入打印机的灰平衡组分基色量幂函数式,计算作为灰核利用的灰色平衡组分基色量和灰色取代参数kp:
根据所选打印机-常白型显示器映射方程型式,只需选择计算其中的rdd,或者gdd,或者bdd;
(3)将算得的rdd,或gdd,或bdd数值代入已选定的打印机-常白型显示器映射方程的右端;
第二步,把打印设备特性文件连接空间的色度坐标xp和yp映射到显示设备色域,方法是利用刘氏打印机-常白型显示器色域映射方程作为工具完成映射过程,的三个子类型rddgb、rgddb、rgbdd的方程式如下:
上式中,驱动输入数值drdgdb进一步是基准基色量rgb的函数,由此可得把常白型显示器基准基色量数值rgb转换成驱动输入值RGB的方法是:
通过刘氏综合映射方程三个不同子类型计算出来的三个基色量都是经过伽玛校正的基色量。设:Xqw=(xn/yn)Yt、Yqw=Yt、Zqw=[(1-xn-yn)/yn]Yt,那么,最终显示的颜色就是Xqw、Yqw、Zqw。RGB就是显示颜色Xqw、Yqw、Zqw所需要的驱动数值。
四、附图说明
Fig.1从扫描仪到打印机设备颜色空间和经由打印机设备颜色空间再到显示器设备颜色空间实施色域映射的流程原理图:从扫描仪颜色预测→将预测色转换到刘氏Dlxlyl特性文件连接空间→到打印机设备颜色空间和经由打印机设备颜色空间→显示器设备颜色空间实施色域映射的流程原理图;
Fig.2把数码相机获取的颜色RGB转换为XYZ数值和根据XYZ数值计算计算机显示器驱动数值RGB的流程原理图;
Fig.3把数码摄像机和电视摄像机获取的颜色RGB转换为XYZ数值和根据XYZ数值计算电视显示器驱动数值RGB的流程原理图;
五、具体实施方式
1.下面参照Fig.1,以映射流程“从扫描仪根据扫描数值RGB预测颜色XYZ→把预测色XYZ转换到刘氏Dlxlyl特性文件连接空间→把颜色XYZ分色成为驱动输入数值CMYK→在常白型显示器上生成预打样彩色图像”说明本发明在色域映射中的通用性、创造性、和工作效率:
1)用本发明提供的方法将R,G,B转换成三刺激值X,Y,Z;
2)对扫描仪CCD元件获取的RGB三刺激值进行判断,根据R,G,B中的最小数值选择对颜色RGB进行分割的方程格式。选择原则是:
如果B是最小值,则选用MYK型刘氏扫描分割方程;
如果R是最小值,则选用CYK型刘氏扫描分割方程;
如果G是最小值,则选用CMK型刘氏扫描分割方程。
为免赘述,下面仅以选用了MYK型分割方程的流程进行解说,另外两种情况依此类推;
3)将目标值XYZ放在MYK型颜色分割方程的左端,解方程后得到黑色量p的数值;
4)黑色量p代入刘氏扫描分割方程的灰色梯尺格式,可以算得灰色灰色亮度刺激值Yl,即有:
Yl=Ywo(1-p)+p Ysk
5)将Yl换算成密度值Do:即让Do=Dl=lg(Ywo/Yl),式中的Ywo是在扫描色靶白色面上测得的亮度数值;
6)把算得的XYZ三刺激值换算成色度坐标值xo和yo,即让:
通过上述步骤,RGB色空间-CIEXYZ色空间-Dlxlyl特性文件连接色空间的转换即告结束。接着的步骤就是向打印机模块传送Do、xo、yo数据了;
7)把扫描仪特性文件连接空间的Do值映射到打印机空间,即让
8)根据Dp计算伽玛校正密度:让
9)将伽玛校正密度Dp′代入打印机的灰色平衡多项式,计算作为灰核利用的组分基色量和灰色取代参数kdd和cdd,即让:
cp=cdd=a0+a1Dp′+a2Dp′2+a3Dp′3+…,kp=kdd=d0+d1Dp′+d2Dp′2+d3Dp′3+…注意:根据所选四色映射方程型式,只需计算kdd,cdd,或者kdd,mdd,或者kdd,ydd;
10)将算得的kdd,cdd代入已选定的MYK型四色映射方程的右端,算出基色量
k=kdd,c=cdd,m和y;
11)把基色量k,c,m和y进一步转换成驱动数值CMYK;
12)计算扫描仪在映射到打印机色域的颜色XpYpZp:
Xp=(xo/yo)Yp,Yp=Yp,Zp=[(1-xo-yo)/yo]Yp;
13)将打印机(印刷机)灰色阶调密度数组[Dp]映射到显示设备色域中去,即让:
14)计算密度值Dq的伽玛校正值Dq′,即让
15)将显示器灰色平衡多项式中的变量Dq用Dq′数据代入,算出灰核cdd的数值;
16)将算得的三刺激值XpYpZp换算成色度坐标值xp和yp,即:
17)将cdd,xp,yp同时代入显示器的RGK型加色映射方程,解此方程可以算得基色量r,g,b=bdd,以及被显示颜色的三刺激值XqYqZq;
18)为了在显示器上显示目标颜色XqYqZq,进一步的工作是把算得的基色量rgb转换到它们的驱动色空间去,即让:
用RGB值驱动显示器,就会显示一个待印刷的预打样图像。
2.下面参照Fig.2,以映射流程“数码照相机根据CCD感知的RGB数值预测颜色XYZ→把预测色XYZ转换到刘氏Dlxlyl特性文件连接空间→把颜色XYZ分色成为驱动输入数值RGB→在显示器上生成彩色图像”说明本发明在色域映射中的通用性、创造性和实时转换效率:
1)将R,G,B代入RGB-XYZ色坐标变换式把RGB转换成为XYZ;
2)将R,G,B代入亮度方程式计算出亮度数值Y;
3)将亮度Y代入刘氏分割方程的灰色梯尺格式算出灰阶亮度Ygrey;
4)把亮度数值换算成灰色密度数值Dl=Dn=lg(Yw/Yn);
5)将Dn映射到显示器色域的密度Dqw=Dn^γqw;对Dqw进行伽玛正,得Dqw′;
6)由Dqw′计算灰核rdd并将灰核rdd代入GBK型显示器映射方程;
7)把颜色的XYZ三刺激值换算成色度坐标值xn和yn,即让:
8)将xn和yn代入GBK型显示器映射方程;
9)解GBK型显示器映射方程:得基色量rgb和预示的显示颜色XqwYqwZqw;
10)从基色量空间的rgb转换到驱动数值RGB,用驱动数值RGB驱动显示器,在显示器上就显示出数码相机所拍摄的彩色图像。
3.下面参照Fig.3,以映射流程“电视摄像机或数码摄像机(DV)根据CCD感知数值RGB预测颜色XYZ→把预测色XYZ转换到刘氏Dlxlyl特性文件连接空间→把颜色XYZ分色成为电视机驱动输入数值RGB→在显示器上生成电视图像”说明本发明在色域映射中的通用性、创造性和高转换效率:
1)将R,G,B代入RGB-XYZ色坐标变换式,把RGB转换成为XYZ;
2)将R,G,B代入亮度方程式计算出亮度数值Y;
3)将亮度Y代入刘氏分割方程的灰色梯尺格式算出灰阶亮度Ygrey;
4)把亮度数值换算成灰色密度数值Dl=Dn=lg(Yw/Yn);
5)将Dn映射到显示器色域的密度Dqk=Dn^γqk;对Dqk进行伽玛正,得Dqk′;
6)由Dqk′计算灰核rdd并将灰核rdd代入GBK型显示器映射方程;
7)把颜色的XYZ三刺激值换算成色度坐标值xn和yn,即让:
8)将xn和yn代入GBK型电视显示器映射方程;
9)解GBK型电视显示器映射方程:得基色量rgb和预示的显示颜色XqkYqkZqk;
10)从基色量空间的rgb转换到驱动数值RGB,用驱动数值RGB驱动电视显示器,在显示器上就显示出电视摄像机所拍摄的彩色图像。
比较Fig.2和Fig.3可知,除了黑、白场颠倒之外,两者并无本质上的不同。
Claims (10)
1.一种通用的色域映射及色彩管理方法,其特征在于:
第一,提出了一种《为输入、显示和输出设备所通用的、驱动数值也对应相等的色靶结构》;该色靶上分布的颜色样本包括:(1)单色的三基色梯尺:对于显示器来说,就是给定三基色rgb的21个驱动数据dri、dgi、dbi,对于数码照相机和电视摄像机之类的标定色靶来说,按照与此同样的数据制作成相纸色靶,对于打印机和扫描仪的标定色靶来说,就是给定三基色cmy和单色黑的21个驱动数据dci、dmi、dyi、dki;利用上述驱动数据生成样本颜色并据此样本实测颜色的三刺激值;(2)色靶上的样本颜色还包括二次颜色样本和三次颜色样本的驱动数据;实测的样本颜色三刺激值和已知的驱动变量[di]将是实现第二款所述刘氏基色嵌位方程、第三款所述刘氏配色方程、第四款所述特性化刘氏配色方程的方法所依赖的基本数据;
第二,提出了《为输入设备、显示设备和输出设备所通用的保持基色通道独立性的方法》;
第一步,根据减色法原理在XYZ和RGB颜色空间分别建立刘氏基色嵌位方程及其衍生的参数模型:刘氏基色嵌位方程及其衍生的参数模型是对前述色靶结构的数学描述,在XYZ和RGB颜色空间对前述色靶结构进行数学描述的刘氏基色嵌位方程如下所示:
左侧的刘氏基色嵌位方程是在XYZ颜色空间为扫描仪、打印机、常白型显示器建立的基色嵌位方程,右侧的刘氏基色嵌位方程是在RGB颜色空间为扫描仪建立的基色嵌位方程,变量at表示用等能三基色[ri,gi,bi]合成的单位白色量的三刺激值[Xs,Ys,Zs]和[Rs,Gs,Bs],即:白点的三刺激值分别是[Xw,Yw,Zw]和[Rw,Gw,Bw];
第二步,保持基色通道独立性的三个步骤:
分别求解XYZ和RGB颜色空间的刘氏基色嵌位方程,可得到计算嵌位亮度Yt和Gt的函数式如下所示:
第三步,把嵌位基准基色量at变成为如下所示的基准基色量函数a:
三次灰色样本序列分别是由等能驱动数据(dri+dgi+dbi)、(dci+dmi+dyi)分别用加色法和减色法生成的中性颜色XsYsZs、RsGsBs;根据电视机的色感匹配原理:目标色的亮度等于各原色光亮度相加的总和,于是衍生出如下所示的基准基色量公式:
适用于打印机类设备:
适用于常白型显示器:
根据第三步所述基准基色量a,即可衍生出RGB色空间的基准基色量公式如下所示:
适用于扫描输入设备:
第三,建立了《确保基色在配色空间的三个通道内具有独立性的方法及刘氏配色方程》:RGB和XYZ色空间的配色方程是根据红、绿、蓝、黄、品红、青、白、黑等8种组分色元的常量三刺激值为入射三刺激值RGB和XYZ匹配组分基色量[r,g,b]和[c,m,y]的,为了使RGB和XYZ的组分三基色量[r,g,b]和[c,m,y]在三通道内保持通道独立性和空间独立性,不同种类的刘氏配色方程都设置9个通道基色量参数并把通道基色量参数设置成基准基色量[r,g,b]和[c,m,y]的幂函数,基准基色量又被设置成为驱动变量[dr,dg,db]和[dc,dm,dy]的幂函数,从而使入射光的三刺激值RGB和XYZ变成12维的空间矢量,利用幂函数指数的伽玛校正功能够把三基色[r,g,b]和[c,m,y]的独立特性同时传递给9个通道基色量和3个驱动变量;
第四,提出了特性化刘氏配色方程的方法:
第一步,为了获取三基色梯尺的三刺激值,根据权利要求1第一款所述《一种为输入、显示和输出设备所通用的、驱动数值也对应相等的色靶结构》,分别在单色三基色梯尺上实测样本色的三刺激值XYZ和RGB,约定用下标o,p,qw,qk区分扫描仪、打印机、常白型计算机显示器、常黑型电视机显示器的三刺激值数据,那么可以得到如下15组三刺激值数据:
对偶的两种扫描仪数据:在RGB色空间实测的三刺激值数据:[Roci,Goci,Boci],[Romi,Gomi,Bomi],[Royi,Goyi,Boyi];在XYZ色空间实测的三刺激值数据:[Xoci,Yoci,Zoci],[Xomi,Yomi,Zomi],[Xoyi,Yoyi,Zoyi];打印机类、常白型计算机显示器类、常黑型电视机显示器类在XYZ色空间实测的三刺激值数据分别如下所示:[Xpci,Ypci,Zpci],[Xpmi,Ypmi,Zpmi],[Xpyi,Ypyi,Zpyi];[Xqwri,Yqwri,Zqwri],[Xqwgi,Yqwgi,Zqwgi],[Xqwbi,Yqwbi,Zqwbi];[Xqkri,Yqkri,Zqkri],[Xqkgi,Yqkgi,Zqkgi],[Xqkbi,Yqkbi,Zqkbi];
在本步骤中:扫描仪在RGB和XYZ色空间具有三组对偶的实测三刺激值,打印机、常白型和常黑型显示器在XYZ色空间各有三组实测的三刺激值,四类设备共计拥有15组实测的三刺激值数据;
第二步,根据以上15组三刺激值,应用权利要求1第二款所述《为输入设备、显示设备和输出设备所通用的一种保持基色通道独立性的方法》中给出的基准基色量公式分别为扫描仪、打印机、常白型计算机显示器、常黑型电视机显示器计算嵌位亮度数据:
[Gtoci,Gtomi,Gtoyi]、[Ytoci,Ytomi,Ytoyi]、[Ytpci,Ytpmi,Ytpyi]、[Ytqwri,Ytqwgi,Ytqwbi]、[Ytqkri,Ytqkgi,Ytqkbi];
第三步,将算得的嵌位亮度数值代入下列基准基色量公式,把基准基色量计算出来:
第四步,用下列模型为刘氏减色配色方程、刘氏四色配色方程以及刘氏扫描仪颜色预测方程计算通道基色量:
第五步,用下列模型为刘氏RGB扫描分色方程计算通道基色量[cxi′,cyi′,czi′],[mxi′,myi′,mzi′],[yxi′,yyi′,yzi′]:
第六步,用下列模型为刘氏常白型显示器配色方程计算通道基色量[rxi,ryi,rzi],[gxi,gyi,gzi],[bxi,byi,bzi]:
第七步,用下列模型为刘氏常黑型显示器配色方程计算通道基色量[rxi,ryi,rzi],[gxi,gyi,gzi],[bxi,byi,bzi]:
第八步,用曲线拟合方法为基于减色法复制的、用于打印机、扫描仪的刘氏减色配色方程构造通道基色量函数:将基准基色量数组[ci,mi,yi]分别与通道基色量数组[cxi,cyi,czi]、[mxi,myi,mzi]、[yxi,yyi,yzi]进行曲线拟合,即得通道基色量的幂函数表达式:
第九步,用曲线拟合方法为刘氏常白型显示配色方程构造通道基色量函数:将三基色的基准基色量数组[rj]、[gj]、[bj]分别与对应的通道基色量数组[rxj,ryj,rzj]、[gxj,gyj,gzj]、[bxj,byj,rzj]进行曲线拟合,即得通道基色量的幂函数表达式:
第十步,用曲线拟合方法为数字照相机或者电视摄像机配色方程构造通道基色量函数:
第五,提出了《一种为扫描仪、打印机、常白型显示器和常黑型显示器生成纯净灰色梯尺的方法》:它是在第一款所述为扫描仪、打印机以及显示器上逐一显示的近似灰色梯尺样本色上,实测五组灰色梯尺亮度数组:[Go′ai]、[Yoai]、[Ypai]、[Yqwai]、[Yqkai],把此5个亮度数组作为初始亮度数组,然后通过九个转换步骤为四类色彩复制设备再现的三刺激值排除了红移误差,藉此生成纯净灰色密度数组[Do′i]、[Doi]、[Dpi]、[Dqwi]、[Dqki],为应用第六款所述的方法把RGB颜色空间的纯净灰色阶调密度Do′和XYZ颜色空间的纯净灰色阶调密度Do、Dp、Dqw、Dqk分解成第六款所述‘组分灰核’创造了条件;本发明首先定义了命名为“灰核”的新概念:即:当用三基色[c,m,y]匹配一个中性灰色时,三基色量并不相等,其中具有最小基准基色量的基色就是该灰色成分的灰核;
第六,提出了《一种特性化纯净灰色梯尺组分基色量的方法》:该方法首先用刘氏配色方程把五个纯净灰色梯尺三刺激值数组[Roi,Goi,Boi]、[Xoi,Yoi,Zoi]、[Xpi,Ypi,Zpi]、[Xqwi,Yqwi,Zqwi]、[Xqki,Yqki,Zqki]分解成为组分‘灰核’[c′i,m′i,y′i]=[cddi′,mddi′,yddi′]、[ci,mi,yi]、[cpi,mpi,ypi]、[rqwi,gqwi,bqwi]、[rqki,gqki,bqki],分别以纯净灰色密度数组[Do′i]=[Drgbi]、[Doi]、[Dpi]、[Dqwi]、[Dqki]为自变量,以组分‘灰核’数组为因变量建立组分‘灰核’与纯净灰色密度之间的灰色平衡幂函数式,利用灰色平衡幂函数式可以实现密度与组分‘灰核’的双向转换;鉴于多媒体介质的光吸收作用导致灰色阶调偏暗,必须采用第七款所述的伽玛校正方法提高灰色密度的亮度;
第七,提出了《一种对图像的灰色阶调进行伽玛校正的方法》;该方法首先把打印机、常白型显示器、常黑型显示器的纯净灰色密度[Dp,Dqw,Dqk]用曲线拟合方法变成以驱动变量[di]为自变量的幂函数,然后利用驱动变量[di]对3类图像输出设备中的复制密度[Dp,Dqw,Dqk]实施十进制平方反比伽玛校正,藉此把XYZ色空间的纯净灰色密度[Dp,Dqw,Dqk]的亮度提高到RGB色空间的灰色密度[Dp′,Dqw′,Dqk′];接着用第八款所述的方法在RGB和XYZ色空间之间双向传递[Dl,xl,yl]特性文件数据;
第八,提出了《在XYZ空间为扫描仪建立Dlxlyl特性文件连接颜色空间的方法》;该方法是利用刘氏颜色分割方程及其灰色梯尺格式作为工具创建Dlxlyl特性文件的,该特性文件利用密度参数Dl和色度坐标[xl,yl]数据在RGB和XYZ色空间之间双向传递亮度和色度坐标信息,扫描仪用刘氏颜色分割方程在XYZ颜色空间的三个子格式如下所示:
扫描仪用刘氏颜色分割方程在XYZ颜色空间的灰色梯尺格式如下所示:
Xv=Xw(1-p)+pXsk
Yv=Yw(1-p)+pYsk
Zv=Zw(1-p)+pZsk
扫描仪用刘氏颜色分割方程在RGB色空间的三个子格式如下所示:
扫描仪用刘氏颜色分割方程在RGB色空间的灰色梯尺格式如下所示:
Rv=Rw·(1-p)+p·Rsk,
Gv=Gw·(1-p)+p·Gsk,
Bv=Bw·(1-p)+p·Bsk;
从第九款所述的方法开始,将要把前述的色域映射方法扩展到在RGB和XYZ之间双向交换三刺激值数据;
第九,提出了《在RGB色空间快速计算基准基色量c′m′y′的方法及刘氏三色嵌位方程》:第一步,建立刘氏三色嵌位方程:本发明给出的刘氏三色嵌位方程是‘刘氏扫描仪颜色预测方程’和‘刘氏RGB扫描分色方程’的交集,具有如下所示的CMK、CYK、MYK三种联立的子格式:
方程中的已知量包括:基色红、绿、蓝、黄、品红、青的实地三刺激值RrGrBr,RgGgBg,RbGbBb,RcGcBc,RmGmBm,RyGyBy、白点的三刺激值RwGwBw、由三基色c′m′y′合成的黑色实地三刺激值RskGskBsk以及根据第六款所述密度Do’=Drgb的组分灰核[cdd′,mdd′,ydd′];要计算未知量是RGB色空间的三基色[c′,m′,y′]和波长参数λ,至于选用哪一个刘氏三色嵌位方程对颜色RGB进行转换计算,是由RGB中的最小值决定的,本发明借助刘氏三色嵌位方程把RGB数据转换成基准基色量数据c′,m′,y′;
第二步,把RGB色空间的基色量c′m′y′转换成为XYZ色空间的扫描基色量cmy,再把得到的cmy数据代入扫描仪颜色预测方程的9个通道基色量幂函数式,进而把9个通道基色量函数代入该方程算得经过灰色标定的CIEXYZ数据;位于转换链上的每个环节都继承刘氏三色嵌位方程的三条运动路经;用第十款所述的方法能够为RGB→XYZ转换方法给出三条确定的光路;
第十,提出了《把扫描颜色从RGB色空间转换到CIEXYZ色空间的方法》:该方法通过刘氏颜色分割方程把刘氏三色嵌位方程的三种子格式[GBK,RBK,RGK]与RGB→XYZ转换流程联系起来,建立了一个[GBK,RBK,RGK]→[p,Gl,Drgb]→[cdd′,mdd′,ydd′]→[c′,m′,y′]→[c,m,y]→[X,Y,Z]转换链;上述从RGB色空间到XYZ色空间转换的矩阵方程式如下所示:
在上述RGB→XYZ转换矩阵方程式等号右侧的方阵是第四款第五步所述RGB扫描分色方程的九个通道基色量函数,它们和白色量p是等效的,与RGB→XYZ色空间转换矩阵方程式对偶的XYZ→RGB色空间转换矩阵方程式如下所示:
在上述色空间转换的矩阵方程式等号右端的方阵是第四款第四步所述刘氏扫描仪颜色预测方程的9个通道基色量函数,它们和密度Dl是等效的,嵌位亮度[Gtc,Gtm,Gty]和[Ytc,Ytm,Yty]在刘氏三色嵌位方程的交集是灰核数值[cdd′,mdd′,ydd′],从而为第十一款的Dlxlyl特性文件连接颜色空间提供了依据;
第十一,提出了《为数码照相机或者数码电视摄像机建立Dlxlyl特性文件连接颜色空间的方法》;该方法首先用三刺激值XYZ中的亮度值Y作为刘氏颜色分割方程灰色梯尺格式的自变量计算纯净灰色的三刺激值XgreyYgreyZgrey,然后对灰色亮度值Ygrey实施对数变换,算得与Ygrey对应的灰色密度Dl=Dn,根据已知输入三刺激值XYZ计算其色度坐标[xl,yl]=[xn,yn],然后向接收端传送特性文件数据Dlxlyy=Doxoyo,接收端将需要应用第十二款所述的伽玛校正方法对接收到的密度数据Dl进行伽玛校正,通过伽玛校正把密度数据Dl的亮度提升为亮度较高的RGB颜色空间的密度数据;
第十二,提出了《一种将特性文件连接空间的参数Dl映射成为目标设备色域的灰色阶调参数的方法》;在第五款中,已经得到数码相机、电视摄像机、扫描仪、打印机与印刷机、显示器在设备色域内的纯净灰色阶调密度数组[Do′i]、[Doi]、[Dpi]、[Dqwi]、[Dqki],接收端用接收到的密度数组作为自变量、用目标设备密度数组作为因变量建立幂函数关系,然后再用第十三款所述的方法将源设备获取的密度数据映射为目标设备的全息彩色数据;
第十三,提出了《一种将源设备获取的颜色映射到目标设备色域的方法》:
1)将扫描仪获取的颜色映射到打印设备色域的方法:
方法:扫描仪和打印设备之间的色域映射是通过传递特性文件连接空间的[Do,xo,yo]数值到打印设备色域和借助于刘氏扫描仪-打印机映射方程实现的;
第一步,根据扫描得到的RGB数值,选择刘氏扫描仪-打印机映射方程的子类型:刘氏扫描仪-打印机映射方程具有CMK、CYK、MYK三个子类型,下面依次列出这三个子类型的方程式:
第二步,把扫描仪色空间的灰色阶调密度Do映射成为打印设备色域的灰色阶调密度Dp,方法是:(1)让(2)计算伽玛校正密度:让(3)用伽玛校正密度Dp′取代打印机灰色平衡多项式中的Dp,生成如下所示的刘氏四色打印反伽玛灰色平衡多项式,应用此多项式可计算伽玛校正后的灰色平衡组分基色量cp、mp、yp和灰色取代参数kp,并且把cp、mp、yp、kp作为可以选用的灰核参数使用:
cp=cdd=a0+a1Dp′+a2Dp′2+a3Dp′3+…,mp=mdd=b0+b1Dp′+b2Dp′2+b3Dp′3+…
yp=ydd=c0+c1Dp′+c2Dp′2+c3Dp′3+…,kp=kdd=d0+d1Dp′+d2Dp′2+d3Dp′3+…
(4)根据所选子类型的不同,只需把kdd,cdd,或者kdd,mdd,或者kdd,ydd计算出来,提供给刘氏扫描仪-打印机映射方程,该方程就退化成三元二次方程;本步骤算出的密度参数Dp和色度坐标[xo,yo]是和特性化文件Dlxlyl等效的参数;黑色参数[kdd,kp]所代表的电子流携带着量子[ε-e,ε+e]把CMK、CYK、MYK三个子类型以通道共享的方式贯通;
第三步,根据刘氏扫描仪-打印机映射方程解得的三组基色量[c=cdd,m,y]、[m=mdd,c,y]、[y=ydd,m,c]计算匹配颜色[Xp,Yp,Zp]所需要的驱动数值CMYK,在方程中,基准基色量cmyk进一步是驱动输入数值dcdmdydk的函数,由此可得把打印机基准基色量cmyk数值转换成驱动输入值CMYK的方法:
2)将数码相机获取的颜色XYZ映射到常白型显示器色域的方法:
方法:在数码相机和显示设备之间的色域映射是通过把特性文件连接空间的参数[Dn,xn,yn]传送到显示设备并通过刘氏数码相机-常白型显示器映射方程完成的;
第一步,把数码相机色空间的灰色阶调密度Dn映射成为显示设备色域的灰色阶调密度Dqw,即:(1)让(2)计算伽玛校正密度:即让(3)将伽玛校正密度Dqw′代入常白型显示器的灰色平衡组分基色量幂函数式,计算作为灰核利用的灰色平衡组分基色量rdd、gdd、bdd和灰色取代参数kp:
第二步,把输入设备特性文件连接空间的色度坐标xn和yn映射到常白型显示器并通过刘氏数码相机-常白型显示器映射方程完成映射过程:刘氏数码相机-常白型显示器映射方程也具有rddgb、rgddb、rgbdd三个子类型,并且和刘氏颜色分割方程的三个子类型是一一对应的,下面依次列出这三个子类型的方程式:
根据程序选定的rddgb、rgddb、rgbdd子类型,将第一步得到的rdd、gdd、bdd代入方程等号的右端,将[xn,yn]和嵌位亮度Yqw放在方程的左端,即可解得经过伽玛校正或者色域映射处理的基准基色量r、g、b和亮度值Yqw;上式中,驱动输入数值[dr,dg,db]进一步是基准基色量[r,g,b]的函数,由此得到把常白型显示器基准基色量[r,g,b]数值转换成驱动输入值RGB的方法是:
利用刘氏数码相机-常白型显示器映射方程三个不同子类型计算出来的三个基色量都是经过伽玛校正后的基色量;鉴于Xqw=(xn/yn)Yqw、Yqw=Yqw、Zqw=[(1-xn-yn)/yn]Yqw,则最终显示的颜色就是Xqw、Yqw、Zqw,RGB就是显示颜色Xqw、Yqw、Zqw所需要的驱动数值;
3)将电视摄像机获取的颜色XYZ映射到常黑型电视显示器色域的方法:电视摄像机和电视显示设备之间的色域映射是通过把特性文件连接空间的参数Dn、xn和yn传送到显示设备并通过刘氏电视摄像机-常黑型显示器映射方程完成的;其步骤是:
第一步,把电视摄像机色空间的灰色阶调密度Dn映射成为常黑型电视显示器色域的灰色阶调密度Dqk,即:(1)让(2)计算伽玛校正密度:让(3)将伽玛校正密度Dqw′代入常黑型显示器的灰平衡组分基色量幂函数式,计算作为灰核利用的灰色平衡组分基色量rdd、gdd、bdd和灰色取代参数kdd:
第二步,把输入设备特性文件连接空间的色度坐标xn和yn映射到常黑型电视显示器并用刘氏电视摄像机-常黑型显示器映射方程完成映射过程,该映射方程具有rddgb、rgddb、rgbdd三个子类型,并且和刘氏颜色分割方程的三个子类型是一一对应的,下面依次列出这三个子类型的方程式:
根据程序选定的rddgb、rgddb、rgbdd子类型,将第一步得到的rdd、gdd、bdd代入方程等号的右端,将参数xn和yn的值和嵌位亮度参数Yqk放在方程的左端,即可解得经过伽玛校正或者色域映射处理的基准基色量r、g、b和亮度值Yqk;上式中,驱动输入数值drdgdb进一步是基准基色量[r,g,b]的函数,由此可得把电视显示器基准基色量[r,g,b]数值转换成驱动输入值RGB的方法是:
观察上式可知:鉴于:Xqk=(xn/yn)Yqk、Yqk=Yqk、Zqk=[(1-xn-yn)/yn]Yqk,利用刘氏电视摄像机-常黑型显示器映射方程三个不同子类型计算出来的三个基色量都是经过伽玛校正后的基色量;
4)将打印设备获取的颜色XYZ映射到常白型显示器色域的方法:
从打印设备到显示设备的映射是通过把特性文件连接空间的参数Dp、xp和yp传送到显示设备并通过刘氏打印机-常白型显示器映射方程完成的;实施步骤如下:
第一步,把打印设备色空间的灰色阶调密度Dp映射成为显示设备色域的灰色阶调密度Dqw,即:(1)让(2)计算伽玛校正密度:让:将伽玛校正密度Dqw′代入打印机的灰色平衡组分基色量幂函数式,计算作为灰核利用的灰色平衡组分基色量rdd、gdd、bdd和灰色取代参数kp:
(3)将算得的rdd,或gdd,或bdd数值代入已选定的刘氏打印机-常白型显示器映射方程等号的右端;
第二步,把打印设备特性文件连接空间的色度坐标xp和yp映射到显示设备色域,方法是利用刘氏打印机-常白型显示器色域映射方程作为工具完成映射过程,刘氏打印机-常白型显示器色域映射方程的三个子类型rddgb、rgddb、rgbdd的方程式如下:
上式中,驱动输入数值drdgdb进一步是基准基色量rgb的函数,由此可得把常白型显示器基准基色量数值rgb转换成驱动输入值RGB的方法是:
通过刘氏打印机-常白型显示器色域映射方程三个不同子类型计算出来的三个基色量都是经过伽玛校正的基色量;设:Xqw=(xn/yn)Yt、Yqw=Yt、Zqw=[(1-xn-yn)/yn]Yt,那么,最终显示的颜色就是Xqw、Yqw、Zqw,RGB就是显示颜色Xqw、Yqw、Zqw所需要的驱动数值。
2.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下方法实施第一款所述《一种为输入、显示和输出设备所通用的、驱动数值也对应相等的色靶结构》:
第一步:生成单色的三基色梯尺:从最小值0到最大值20把驱动数值[di]设为i=21级,对于显示器来说,让三基色[ri,gi,bi]的驱动数据[dri,dgi,dbi]分别等于:0.00,12.75,25.50,38.25,51.00,63.75,76.50,89.25,102.00,114.75,127.50,140.25,153.00,165.75,178.50,191.25,204.00,216.75,229.50,242.25,255.00,利用给定的驱动数据在显示器上逐个显示样本颜色、实测和记录样本色的三刺激值数据;标定数码照相机和电视摄像机的色靶是根据与此同样的数据制作成相纸色靶;对于打印机和扫描仪色靶来说,就是让三基色[ci,mi,yi]的21个驱动数据[dci,dmi,dyi]分别等于:0%=0.00/255=0,5%=12.75/255=0.05,10%=25.50/255=0.1,15%=38.25/255=0.15,20%=51.00/255=0.20,25%=63.75/255=0.25,30%=76.50/255=0.30,35%=89.25/255=0.35,40%=102.00/255=0.40,45%=114.75/255=0.45,50%=127.50/255=0.50,55%=140.25/255=0.55,60%=153.00/255=0.60,65%=165.75/255=0.65,70%=178.50/255=0.70,75%=191.25/255=0.75,80%=204.00/255=0.80,85%=216.75/255=0.85,90%=229.50/255=0.90,95%=242.25/255=0.95,100%=255/255=1;
第二步:生成二次颜色样本:对于电视显示器上显示的幅射光来说,是根据加色法原理生成的:它们是用等能驱动数值(dr255+dg255)、(dr255+db255)、(dg255+db255)合成的三个二次色:对于打印机和扫描仪来说,它们是用(dc100%+dm100%)、(dm100%+dy100%)、(dc100%+dy100%)合成的三个二次色;对于CMYK四色打印来说还需要增加(dc100%+dk100%)、(dm100%+dk100%)、(dy100%+dk100%)等三个二次色,藉此生成、并且实测上述二次颜色样本的三刺激值,以备在后述刘氏配色方程中作为已知常量数据使用;
第三步:生成三次颜色样本:对于显示器来说,它们是用具有相同(dri+dgi+dbi)驱动数据、从0到255顺序显示的基于加色法原理的灰色样本序列;对于打印机或扫描仪来说,它们是用具有相同驱动数据(dci+dmi+dyi)、从0到100%顺序显示的灰色样本序列,常白型显示器是用等能驱动数据(dri+dgi+dbi)相加所生成的灰色样本序列:当dri=dgi=dbi=dmax时,由等能的(drmax+dgmax+dbmax)生成的颜色是由单位刺激值[Rr,Gr,Br]、[Rg,Gg,Bg]、[Rb,Gb,Bb]叠加合成的中性灰色RsGsBs=RwGwBw;对于打印机或扫描仪来说,它们是用等能三基色数据(dci+dmi+dyi)叠加、从0到dmax顺序显示的灰色样本序列,被合成后的三次色亮度降低,当驱动数据dci=dmi=dyi=100%=1时:(dcmax+dmmax+dymax)是由等能基色单位三刺激值[Xc,Yc,Zc]、[Xm,Ym,Zm]、[Xy,Yy,Zy]相互叠加合成中性灰色XsYsZs,三刺激值RsGsBs和XsYsZs将提供给权利要求1第二款所述的刘氏基色嵌位方程使用;在观测颜色时,本方法采用标准D65日光照明,根据电视机原理:虽然与7种驱动数据[dri,dgi,dbi,dci,dmi,dyi,dki]相对应的7种基本色[R,G,B,C,M,Y,K]的色度不会受非线性的影响,可是亮度是不等的。
3.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的方法实施权利要求1第三款所述的《确保基色在配色空间的三个通道内具有独立性的方法及刘氏配色方程》:
第一,基于减色法复制的刘氏减色配色方程:又称三色打印方程,其格式如下所示:
刘氏减色配色方程具有两个用途:一是作为三色打印配色方程使用,二是作为扫描仪的颜色预测方程使用;该方程左端的三刺激值[X,Y,Z]是待匹配颜色的已知三刺激值,三刺激值XwYwZw、XcYcZc、XmYmZm、XyYyZy、XrYrZr、XgYgZg、XbYbZb、XsYsZs分别表示在色靶上实测的8个实地色的三刺激值;三刺激值[Xw,Yw,Zw]是背景白色的三刺激值,三刺激值[Xs,Ys,Zs]是等量三基色[c,m,y]合成的黑色三刺激值,白色和黑色三刺激值分别和刘氏基色嵌位方程中的三刺激值[Xw,Yw,Zw]和[Xs,Ys,Zs]相同,刘氏减色配色方程设置9个未知变量yx、yy、yz、mx、my、mz、cx、cy、cz,它们分别是X、Y、Z三个通道公用的通道基色量,通道基色量分别是以基准基色量[c,m,y]为自变量的幂函数,幂函数格式如下所示:
上式中的基准基色量[c,m,y]又进一步是以驱动参数[dc,dm,dy]为自变量的幂函数,即:
对此三式求反解得:
刘氏减色配色方程是以刘氏基色嵌位方程及其衍生公式为基础、在原始Neugebauer方程中加入9个通道基色量参数[cx,cy,cz,mx,my,mz,yx,yy,yz]构成的方程,通道基色量参数分别是以3个基准基色量[c,m,y]为自变量的幂函数,9个幂函数指数分别是[γxc,γyc,γzc,γxm,γym,γzm,γxy,γyy,γzy];方程求解的目标值是基准基色量[c,m,y],基准基色量[c,m,y]分别是以驱动参数[dc,dm,dy]为自变量的幂函数,幂函数指数分别是[γc,γm,γy],该方程本身是一个高次方程,需要根据方程左端的实测灰色三刺激值XYZ和实测的8组基色的实测三刺激值XwYwZw,XcYcZc,XmYmZm,XyYyZy,XrYrZr,XgYgZg,XbYbZb,XsYsZs应用权利要求1第四款所述《特性化刘氏配色方程的方法》决定9个通道基色量参数[cx,cy,cz,mx,my,mz,yx,yx,yx]、9个幂函数指数[γxc,γyc,γzc,γxm,γym,γzm,γxy,γyy,γzy]、3个基准基色量[c,m,y]以及3个驱动参数[dc,dm,dy]的数值;
第二,刘氏四色配色方程:和前述三色打印方程一样,刘氏四色配色方程也设置9个通道基色量参数[cx,cy,cz,mx,my,mz,yx,yy,yz]、12个幂函数指数[γxc,γyc,γzc,γxm,γym,γzm,γxy,γyy,γzy,γc,γm,γy]以及3个数字驱动参数[dc,dm,dy],其幂函数格式如下所示:
对左式求反解可得:
刘氏四色配色方程继承三色打印方程中的八种实测三刺激值[Xw,Yw,Zw]、[Xr,Yr,Zr]、[Xg,Yg,Zg]、[Xb,Yb,Zb]、[Xc,Yc,Zc]、[Xm,Ym,Zm]、[Xy,Yy,Zy]、[Xs,Ys,Zs];不同的是:在刘氏四色配色方程中增加了黑色变量k和kdd,由于黑墨变量的加入,而使方程增加了8组实测三刺激值:[Xk,Yk,Zk]、[Xck,Yck,Zck]、[Xmk,Ymk,Zmk]、[Xyk,Yyk,Zyk]、[Xrk,Yrk,Zrk]、[Xgk,Ygk,Zgk]、[Xbk,Ybk,Zbk]、[Xsk,Ysk,Zsk],共计16组实测三刺激值,等号右端的黑色变量kdd是在解方程之前预先用基准基色量公式给定的数值,因此刘氏四色配色方程仍然是只有三个变量c、m、y的三元二次方程,可用解析方法求解,刘氏四色配色方程的完整格式如下所示:
第三,基于加色法复制的刘氏常白型显示器配色方程:其格式如下所示:
和前述刘氏三色打印配色方程相比较,刘氏常白型显示器配色方程也继承了三色打印方程中的八种实测三刺激值数据[Xw,Yw,Zw]、[Xr,Yr,Zr]、[Xg,Yg,Zg]、[Xb,Yb,Zb]、[Xc,Yc,Zc]、[Xm,Ym,Zm]、[Xy,Yy,Zy]、[Xs,Ys,Zs],同样设置9个通道基色量参数[rxr,ryr,rzr,gxg,gyg,gzg,bxb,byb,bzb]、9个幂函数指数[γxr,γyr,γzr,γxg,γyg,γzg,γxb,γyb,γzb,γc,γm,γy]以及3个数字驱动参数[drd,dgd,dbd],不同之处在于:刘氏常白型显示器配色方程的通道基色量和驱动变量分别是以基准基色量[r,g,b]为自变量的幂函数,本方程是以加色法原理为三刺激值XYZ匹配基准基色量[r,g,b]的配色方程,而不再是基于减色法原理为三刺激值XYZ匹配基准基色量[c,m,y]的配色方程,其格式如下所示:
上述不同之处表明:打印机是在三色打印配色方程的[c,m,y]减色空间工作的,刘氏常白型显示器配色方程式是在[r,g,b]加色空间工作的;
第四,基于加色法复制的刘氏常黑型显示器配色方程:其格式如下所示:
刘氏常黑型显示器和刘氏常白型显示器的区别在于:前者是在黑色背景幅射三刺激值[Xk,Yk,Zk]的基础上观测到的三刺激值XYZ,其本质是以普朗克黑体幅射为背景观测得到的三刺激值XYZ;后者是在白色背景幅射三刺激值[Xw,Yw,Zw]的基础上观测到的三刺激值XYZ;
第五,刘氏RGB扫描分色方程:其格式如下所示:
刘氏RGB扫描分色方程中的通道基色量是以基准基色量[c′,m′,y′]为自变量的幂函数,其函数格式和三色打印方程相同:
刘氏RGB扫描分色方程左端的三刺激值[R,G,B]是待匹配颜色的已知三刺激值,三刺激值RwGwBw、RcGcBc、RmGmBm、RyGyBy、RrGrBr、RgGgBg、RbGbBb、RsGsBs分别表示在反射式色靶上、即打印机所用反射式色靶上获得的8个实地色的三刺激值;刘氏RGB扫描分色方程是以RGB色空间的刘氏基色嵌位方程及其衍生公式为基础建立的配色方程,方程设置9个通道基色量参数[cx′,cy′,cz′,mx′,my′,mz′,yx′,yy′,yz′]、9个幂函数指数[γxc′,γyc′,γzc′,γxm′,γym′,γzm′,γxy′,γyy′,γzy′];三刺激值[Rw,Gw,Bw]是背景白色的三刺激值,三刺激值[Rs,Gs,Bs]是等量三基色[c′,m′,y′]合成的黑色三刺激值,并且和RGB色空间刘氏基色嵌位方程中的三刺激值[Rw,Gw,Bw]和[Rs,Gs,Bs]是一致的。
4.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施权利要求1第五款所述的《一种为扫描仪、打印机、常白型显示器和常黑型显示器生成纯净灰色梯尺的方法》:
第一步,为扫描仪、打印机和显示器生成纯净灰色梯尺:根据权利要求1第一款所述的《为输入、显示和输出设备所通用的、驱动数值也对应相等的色靶结构》,在为扫描仪、打印机制作的色靶上以及在显示器上逐一显示的梯尺样本色上,分别实测复色梯尺样本的亮度值,得到5个初始亮度数组[Go′ai]、[Yoai]、[Ypai]、[Yqwai]、[Yqkai];
第二步,将初始亮度数组[Go′ai]、[Yoai]、[Ypai]、[Yqwai]、[Yqkai]转换成初始密度数组[Do′ai]、[Doai]、[Dpai]、[Dqwai]、[Dqkai],即让:[Do′ai]=lg(Gow/Go′ai),[Doai]=lg(Yow/Yoai),[Dpai]=lg(Ypw/Ypai),[Dqwai] =lg(Ywq/Yqwai),[Dqkai]=lg(Ywq/Yqkai);
第三步,将初始密度数组[Do′ai]、[Doai]、[Dpai]、[Dqwi]、[Dqki]分别除以各数组中的最大值Do′amax、Doamax、Dpamax、Dqwamax、Dqkamax,从而得到归一化的初始密度数组:[Do′bi]=[Do′ai]/Do′amax,[Dobi]=[Doai]/Doamax,[Dpbi]=[Dpai]/Dpamax,[Dqwbi]=[Dqwai]/Dqwamax,[Dqkbi]=[Dqkai]/Dqkamax;
第四步,将基色梯尺的驱动数值进行归一化处理,得归一化的驱动数组[di];根据权利要求1第一款第一步所述:字符[di]是基色[p,r,g,b]的公用驱动数据,变量[di]使白色量[pi]的三基色[r,g,b]、[c,m,y]在[ri,gi,bi]和[ci,mi,yi]颜色空间内分别获得准确的位置;
第五步,以归一化驱动数组[di]为自变量数组,分别以归一化初始密度数组[Do′bi]、[Dobi]、[Dpbi]、[Dqwbi]、[Dqkbi]为因变量进行幂函数拟合,得灰色梯尺的归一化初始密度模型:
Dob=d^γo,Do′b=d^γo′,Dpb=d^γp,Dqwb=d^γqw,Dqkb=d^γqk;
第六步,为了得到扫描仪、打印机、常白型显示器和常黑型显示器的真正纯净灰色密度数组[Do′i]、[Doi]、[Dpi]、[Dqwi]、[Dqki],应当对上面的归一化初始密度分别进行去归一化处理,为此可分别用字符Do′、Do、Dp、Dqw、Dqk表示去归一化后的初始密度数组,那么可得:
[Doi′]=[doi]^γo′×Do′amax,[Doi]=[di]^γo×Doamax,[Dpi]=[di]^γp×Dpamax,[Dqwi]=[di]^γqw×Dqwamax,[Dqki]=[di]^γqk×Dqkamax;
第七步,为扫描仪、打印机、常白型显示器、常黑型显示器分别生成理想灰色亮度数组[Goi]、[Yoi]、[Ypi]、[Yqwi]、[Yqki]:即让:[Goi]=Gwo/(10^Do′i),[Yoi]=Ywo/(10^Doi),[Ypi]=Ywp/(10^Dpi),[Yqwi]=Yww/(10^Dqwi),[Yqki]=Ywk/(10^Dqki);
第八步,为扫描仪、打印机、常白型显示器、常黑型显示器计算参考白点的色度坐标值:用下标w、k分别表示扫描仪在XYZ、RGB颜色空间的色度坐标,设下标o、p、qw、qk分别表示扫描仪、打印机、常白型显示器、常黑型显示器等四种设备的色度坐标,计算色度坐标的公式如下所示:
row=Row/(Row+Gow+Bow),gow=Gow/(Row+Gow+Bow),
xow=Xow/(Xow+Yow+Zow),yow=Yow/(Xow+Yow+Zow),
xpw=Xpw/(Xpw+Ypw+Zpw),ypw=Ypw/(Xpw+Ypw+Zpw),
xqw=xqk=0.3127,yqw=yqk=0.3290;
第九步,利用第七步给出的纯净灰色亮度数组[Yoi]、[Goi]、[Ypi]、[Yqwi]、[Yqki]和第八步给出的色度坐标值[xow,yow]、[row,gow]、[xpw,ypw]、[xqw,yqw]、[xqk,yqk]分别为扫描仪、打印机、常白型显示器、常黑型显示器构造一条纯净灰色梯尺,每条纯净灰色梯尺的三刺激值如下所示:
扫描仪在RGB色空间的纯净灰色梯尺数组:Roi=(row/gow)·Goi,Goi=Goi,Boi=(1-row-gow)·Goi;
扫描仪在XYZ色空间的纯净灰色梯尺数组:
打印机在XYZ色空间的的纯净灰色梯尺数组:
常白型显示器在XYZ色空间的纯净灰色梯尺数组:
常黑型显示器的灰色梯尺数组:
5.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施第六款所述的《一种特性化纯净灰色梯尺组分基色量的方法》:
为了对纯净灰色梯尺上的组分基色量实施特性化处理,本发明定义了命名为“灰核”的新概念:即:当用三基色[c,m,y]匹配了一个中性灰色时,三个基色量并不相等,其中具有最小基准基色量的基色就是形成该颜色灰色成分的灰核,虽然‘灰核’是三个基色量中的最小值,但是它和另外两个基色相结合,形成为一个复合颜色的灰色核心,借助‘灰核’概念能够使一个未知颜色XYZ达到理想的‘灰色平衡’,这是准确、快速分离中性和彩色成分的关键步骤;
第一步,在RGB色空间为扫描仪建立灰色平衡幂函数式c′(Drgb)、m′(Drgb)、y′(Drgb):方法是:用RGB扫描分色方程对扫描仪的灰色梯尺数组[RoiGoiBoi]进行配色计算,得扫描仪灰色梯尺的基准基色量数组[ci′,mi′,yi′];再将灰色梯尺亮度数组[Goi]代入公式Drgb=lg(Gow/Goi)算得密度数组[Drgbi];然后以[Drgbi]为自变量、分别以孪生基色量[ci′,mi′,yi′]为因变量进行曲线拟合,最终得扫描仪的灰色平衡幂函数式如下:
在后面叙述的扫描仪分色嵌位方程中,将把孪生的三基色量cdd′、mdd′、ydd′作为灰核使用,鉴于di=λmax时,等能三基色[ri,gi,bi]合成为白色量p,密度Drgb和白色量p是反函数,因此三基色量[cdd′,mdd′,ydd′]和密度Drgb是以相同频率(1/γmax)振动的共振体[Drgb,cdd′,mdd′,ydd′];
第二步,为扫描仪建立从孪生基色量[c′,m′,y′]向减色法三基色量[c,m,y]转换的c′m′y′-cmy幂函数式:方法是:用刘氏减色配色方程对灰色三刺激值数组[Xoi,Yoi,Zoi]进行配色计算,即可算得再现灰色三刺激值[Xoi,Yoi,Zoi]的基准基色量数组[ci,mi,yi];然后分别以[ci′],[mi′],[yi′]为自变量、分别以[ci],[mi],[yi]为因变量进行曲线拟合,即得到c′m′y′-cmy的幂函数转换式:
刘氏RGB扫描分色方程和刘氏减色配色方程同属减色复制方程,其白色背景三刺激值[Rw,Gw,Bw]和[Xw,Yw,Zw]给三基色量[c′,m′,y′]和[c,m,y]带来的红移误差已被刘氏基色嵌位方程排除,都具有基准基色量的属性,因此[c′,m′,y′]和[c,m,y]具有相同的振动频率和色相;
第三步,为计算机的常白型显示器和电视机的常黑型显示器建立灰色平衡基色量幂函数式:方法是:分别用常白型显示器的配色方程和常黑型显示器的配色方程对灰色梯尺数组[Xqwi,Yqwi,Zqwi]、[Xqki,Yqki,Zqki]进行配色计算,可得到灰色梯尺的三基色量数组[rqwi,gqwi,bqwi]、[rqki,gqki,bqki],然后分别以纯净灰色密度数组[Dqwi]、[Dqki]为自变量,分别以三基色量数组[rqwi,gqwi,bqwi],[rqki,gqki,bqki]为因变量进行曲线拟合,即得常白型显示器和常黑型显示器的灰色平衡基色量幂函数式如下:
求上式的反函数可得:
第四步,为三色打印机建立灰色平衡幂函数式:
方法:用刘氏三色配色方程对灰色梯尺数组[Xpi,Ypi,Zpi]进行配色计算,得到三基色量数组[cpi,mpi,ypi],然后分别以纯净灰色密度数组[Dpi]为自变量,分别以[cpi,mpi,ypi]为因变量进行曲线拟合,得三色打印机的灰色平衡基色量幂函数式如下所示:
第五步,为四色打印机或四色胶印机建立灰色平衡多项式:方法是:(1)决定灰成分取代量:首先在色靶上测量黑墨印刷梯尺,记录样本色的三刺激值;接着用嵌位亮度公式算出黑基色的嵌位亮度数组[Ytki];然后根据黑基色的嵌位亮度值Ytki用如下公式计算黑基色的基准基色量数组ki:即让ki=(Ywp-Ytki)/(Ywp-Ysk);(2)以归一化驱动数组[dki]为自变量数组,以[ki]为因变量数组进行曲线拟合,得函数:k=dk^γk;(3)用kdd表示灰色成分取代量,让kdd=Q(dk^γk)n,Q是一个预先设定的比例常数,用来控制最大黑版量,n是根据黑版阶调长短决定的指数,只要给定n的数值,就可方便的调节黑色阶调的长度;灰色成分取代量kp是经过净化的黑色成分,称为纯净的灰色成分取代量;(4)求解与[ki]相配套的纯净灰色的基色量数组[ci]、[mi]、[yi]:将纯净灰色梯尺的三刺激值[Xpi,Ypi,Zpi]逐组放在四色配色方程的左端和同步将[kpi]逐一代入刘氏四色配色方程,那么四色配色方程蜕变成只有未知量[c,m,y]的静定的刘氏配色方程,与此同时,逐一用迭代法解方程,那么得到基准基色量数组[ci,mi,yi];(5)以纯净灰色密度数组[Dpi]为自变量、分别以[cpi]、[mpi]、[ypi]、[kpi]为因变量进行曲线拟合,可得四色打印机的灰色平衡多项式,三基色量[cp,mp,yp]及其灰色成份取代量kp如下所示:
cp=a0+a1Dp+a2Dp 2+a3Dp 3+…,mp=b0+b1Dp+b2Dp 2+b3Dp 3+…
yp=co+c1Dp+c2Dp 2+c3Dp 3+…,kp=d0+d1Dp+d2Dp 2+d3Dp 3+…。
6.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施第七款所述的《一种对图像的阶调进行伽玛校正的方法》:
第一步,把纯净灰色密度参数表示成为驱动参数的函数:根据权利要求1第五款所述‘一种为扫描仪、打印机、常白型显示器和常黑型显示器生成纯净灰色梯尺的方法’给出的打印机、常白型显示器和常黑型显示器的纯净灰色密度数组分别是[Dpi]、[Dqwi]、[Dqki],将[Dpi]、[Dqwi]、[Dqki]进行去归一化处理,处理之后仍然用[Dpi]、[Dqwi]、[Dqki]表示,再把它表示成归一化驱动参数[di]的函数:分别以[di]为自变量、分别以纯净灰色密度数值[Dpi]、[Dqwi]、[Dqki]为因变量进行曲线拟合,得纯净灰色密度的幂函数表达式如下所示:
第二步,把伽玛校正密度表示成纯净灰色密度的函数:即让伽玛校正密度Dp′、Dqw′、Dqk′分别等于原始密度Dp、Dqw、Dqk为自变量的函数:
7.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施第八款所述的《一种在XYZ色空间为扫描仪建立Dlxlyl特性文件连接颜色空间的方法》:
1)在XYZ空间内为扫描仪建立Dlxlyl特性文件连接颜色空间的方法
第一步,为扫描仪构建刘氏颜色分割方程:扫描仪用刘氏颜色分割方程的三个子格式CMK、CYK和MYK如下所示:
用色靶上的白色成分[Xw,Yw,Zw]取代比例为(1-p)的彩色成分,即得如下所示的灰色梯尺格式:
Xv=Xw(1-p)+pXsk
Yv=Yw(1-p)+pYsk;
Zv=Zw(1-p)+pZsk
将根据刘氏颜色分割方程解得的[pi]逐一代到刘氏颜色分割方程灰色梯尺格式的右端,那么生成的三刺激值数组[Xvi,Yvi,Zvi]等效于扫描仪纯净灰色梯尺的三刺激值[Xoi,Yoi,Zoi];
第二步,利用上述知识构建Dlxlyl特性文件连接空间:
(1)根据刘氏颜色分割方程左端输入的XYZ值计算灰色三刺激值XYZ的色度坐标值:
(2)应用刘氏颜色分割方程对XYZ进行分割,即通过解刘氏颜色分割方程计算灰色量p的数值;
(3)将算得的p值代入刘氏颜色分割方程灰色梯尺格式中的亮度式,计算刘氏颜色分割方程左端三刺激值XYZ中灰阶的亮度值:Yl=Yv=Yw(1-p)+pYsk;
(4)把灰阶的亮度值Yl换算成灰色密度值Dl,即:Dl=Do=Dxyz=lg(Ywo/Yo);
经过上面4个步骤,就把在XYZ色空间内的任一颜色XYZ或者扫描XYZ颜色变换成了特性文件连接空间的数值Doxoyo;
2)在RGB颜色空间为扫描仪建立Dlxlyl特性文件连接颜色空间的方法:
第一步,为扫描仪构建刘氏颜色分割方程:与前述CMK、CYK和MYK三个子类型相对应,RGB色空间扫描仪用刘氏颜色分割方程的三种格式如下所示:
用白点三刺激值[Rw,Gw,Bw]取代方程右端方括号中的彩色成分,上列扫描仪刘氏颜色分割方程即被演变成如下所示的刘氏颜色分割方程灰色梯尺格式:
Rv=Rw·(1-p)+p·Rsk,
Gv=Gw·(1-p)+p·Gsk,
Bv=Bw·(1-p)+p·Bsk;
如果把根据刘氏颜色分割方程解得的[pi]逐一代到刘氏颜色分割方程灰色梯尺格式的右端,那么生成的三刺激值数组[Rvi,Gvi,Bvi]等效于纯净灰色梯尺的三刺激值[Roi,Goi,Boi];
第二步,计算灰色成分Rv、Gv、Bv的灰色密度值:
方法是:(1)应用刘氏颜色分割方程对RGB进行分割,即:通过解刘氏颜色分割方程算出灰色量p的数值;(2)在RGB和XYZ色空间的白色量p同为灰色色相,具有相等的色度坐标值[xl,yl];(3)将算得的p值代入前述刘氏颜色分割方程灰色梯尺格式中的亮度式,计算颜色RvGvBv中灰阶的亮度值:Gv=Gl=Gw(1-p)+pGsk;(4)把灰阶的亮度值Gl换算成灰色密度值Drgb,即:Drgb=Do′=Dl=lg(Gwo/Gv);
经过上面4个步骤,就把XYZ和RGB色空间内的任一颜色变换成了特性文件连接空间的数值Dlxlyl。
8.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施第十款所述的《把扫描颜色从RGB色空间转换到CIE XYZ色空间的方法》:
第一步,对扫描仪获取的RGB三刺激值进行判断,即根据三刺激值[R,G,B]中的最小数值选择对颜色RGB进行分割的子格式;选用原则是:如果R是最小值,那么选用GBK型分割方程;如果G是最小值,那么选用RBK型刘氏颜色分割方程;如果B是最小值,那么选用RGK型刘氏颜色分割方程;
第二步,根据第一步选出的子格式,将扫描仪获取的RGB三刺激值放在所选刘氏颜色分割方程左端,解方程后得到黑色量p的数值;
第三步,将黑色量p代入刘氏颜色分割方程的灰色梯尺格式,可以算得灰色亮度刺激值Gl,即:Gl=Gw(1-p)+pGsk;
第四步,将亮度Gl换算成密度值Drgb:即让Drgb=Dl=lg(Gwo/Gl);
第五步,将Drgb代入权利要求1第六款所述的《一种特性化纯净灰色梯尺组分基色量的方法》第一步给出的扫描仪灰色平衡幂函数式:
算得灰核[cdd′,mdd′,ydd′]的数值;
第六步,将算出的灰核数值[cdd′,mdd′,ydd′]代入刘氏三色嵌位方程可得颜色RGB的组分基色量[c′,m′,y′];
第七步,将算出的基色量c′、m′、y′代入扫描仪的c′m′y′-cmy转换式算得基准基色量[c,m,y]的数值;
第八步,将算出的基准基色量c、m、y代进刘氏减色配色方程,即可算得目标值XYZ。
9.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施第十一款所述《为数码照相机或者数码电视摄像机建立Dlxlyl特性文件连接颜色空间的方法》:
第一步,在综合考虑当前多媒体设备的普遍需要之后,可以为待匹配三刺激值RGB制定一组普遍认可的色度坐标值,让白点的色度坐标等于D65白色照明的色度坐标,建立如下所示的从RGB色空间到XYZ色空间转换的矩阵方程式:
利用上述从RGB色空间到XYZ色空间转换的矩阵方程式,就可以把扫描输出的三刺激值RGB转换成为CIE XYZ三刺激值;所谓“普遍认可的色度坐标数值”是指符合CIE 1931 XYZ国际坐标制的[Y,x,y]的色度坐标值,该色度坐标值和本发明所述特性文件连接空间[Dl,xl,yl]构成确定性的色域映射关系,从下述步骤中就可以得知[Y,x,y]的色度坐标值与[Dl,xl,yl]特性文件的区别;
第二步,为标准显示器构造刘氏颜色分割方程灰色梯尺格式如下所示:
Xgreyi=Xsk·(1-Y)+Y·Xsw
Ygreyi=Ysk·(1-Y)+Y·Ysw
Zgreyi=Zsk·(1-Y)+Y·Zsw
上式中:等号左侧的灰色三刺激值XgreyYgreyZgrey是根据标准显示器白点的三刺激值XswYswZsw和黑点的三刺激值XskYskZsk计算出来的纯净灰色的三刺激值;
第三步,把纯净灰色的三刺激值XgreyYgreyZgrey中的灰色亮度值Ygrey换算成灰色密度值Dl,即:将灰色亮度值Ygrey代入灰色密度函数Dl:Dl=Dn=lg(Yw/Ygrey),即可算出特性文件连接空间Dlxlyl中的密度值Dl;
第四步,根据第一步得到的CIE XYZ三刺激值计算颜色XYZ的色度坐标值:Dlxlyl特性文件连接空间所需要的色度坐标值[xl,yl]就是根据三刺激值XYZ按照下式计算出来的:
xl=xn=X/(X+Y+Z),yl=yn=Y/(X+Y+Z);
经过上面4个步骤,就把扫描仪得到的三刺激值[R,G,B]颜色变换成了特性文件连接空间的数值Dnxnyn。
10.根据权利要求1所述的一种通用的色域映射及色彩管理方法,其特征是:按照如下的步骤实施第十二款所述《一种将特性文件连接空间的参数Dl映射成为目标设备色域的灰色阶调参数方法》:
第一步,将数码相机的灰色阶调密度数组[Dni]映射到常白型显示器色域中去:以数码相机的灰色阶调密度数组[Dni]为自变量数组,以常白型显示器的灰色阶调密度数组[Dqwi]为因变量数组进行幂函数拟合,即得常白型显示器上显示的密度Dqw:
第二步,将电视摄像机的灰色阶调密度数组[Dni]映射到常黑型电视机显示器色域中去:以电视摄像机的灰色阶调密度数组[Dni]为自变量数组,以常黑型显示设备灰色阶调密度数组[Dqki]为因变量数组进行幂函数拟合,即得常黑型电视机显示器上实际显示的密度值
第三步,将扫描仪的灰色阶调密度数组[Doi]映射到打印机设备色域中去:以扫描仪的灰色阶调密度数组[Doi]为自变量数组,以打印设备灰色阶调密度数组[Dpi]为因变量数组进行幂函数拟合,即得扫描仪-打印机(印刷机)灰密度映射函数式:
第四步,将打印机(印刷机)灰色阶调密度数组[Dpi]映射到显示设备色域中去:以打印机-印刷机的灰色阶调密度数组[Dpi]为自变量数组,以常白型显示器灰色阶调密度数组[Dqwi]为因变量数组进行幂函数拟合,即得打印机(印刷机)-常白型显示器灰密度映射函数式:
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2011/001729 WO2013056387A1 (zh) | 2011-10-17 | 2011-10-17 | 一种通用的色域映射及色彩管理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103931169A CN103931169A (zh) | 2014-07-16 |
CN103931169B true CN103931169B (zh) | 2018-02-02 |
Family
ID=48140277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180042671.9A Active CN103931169B (zh) | 2011-10-17 | 2011-10-17 | 一种通用的色域映射及色彩管理方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20140152687A1 (zh) |
CN (1) | CN103931169B (zh) |
WO (1) | WO2013056387A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110490328A (zh) * | 2019-08-06 | 2019-11-22 | 刘世昌 | 以统一场论和全息电视为基础实现量子计算机的方法 |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160053977A1 (en) | 2008-09-24 | 2016-02-25 | B/E Aerospace, Inc. | Flexible led lighting element |
EP2767144B1 (en) * | 2011-10-12 | 2017-01-11 | B/E Aerospace, Inc. | Methods, apparatus and articles of manufacture to calibrate lighting units |
US8953876B2 (en) * | 2012-08-22 | 2015-02-10 | Facebook, Inc. | Creation of a color profile of an image |
US9697796B2 (en) * | 2015-03-26 | 2017-07-04 | Intel Corporation | Adaptive linear luma domain video pipeline architecture |
CN104766585B (zh) * | 2015-04-29 | 2017-08-25 | 深圳市华星光电技术有限公司 | 设定液晶面板成像时的像素的灰阶值的方法 |
CN104965970B (zh) * | 2015-04-30 | 2017-09-12 | 华文包装股份有限公司 | 用于包装印刷多介质数字打样的专色确定方法 |
CN106937018B (zh) * | 2017-01-12 | 2019-03-05 | 浙江大学 | 基于rbf神经网络用于纺织品喷墨印染的色彩映射方法 |
AU2018230927B2 (en) * | 2017-03-06 | 2020-09-24 | E Ink Corporation | Method for rendering color images |
CN107093031A (zh) * | 2017-05-10 | 2017-08-25 | 广东溢达纺织有限公司 | 颜色数据互联网管理方法与*** |
CN107316602B (zh) | 2017-08-28 | 2019-11-19 | 京东方科技集团股份有限公司 | N基色显示屏的显示控制方法、装置及显示器件 |
CN108205812B (zh) * | 2017-11-22 | 2021-12-17 | 广东工业大学 | 一种颜料色彩混合比例匹配的方法 |
CN109410876B (zh) * | 2018-12-17 | 2021-01-26 | 惠科股份有限公司 | 三色数据到四色数据的转换方法及装置 |
CN110490981B (zh) * | 2019-08-14 | 2020-05-12 | 愉悦家纺有限公司 | 八基色hsb颜色空间的网格化模型及其离散色谱构建方法 |
CN112885300B (zh) * | 2019-11-29 | 2024-04-05 | 美国像素公司 | 使用多个非线性模型的面板校准 |
CN113129222A (zh) * | 2020-01-13 | 2021-07-16 | 华为技术有限公司 | 颜色阴影校正方法、终端设备及计算机可读存储介质 |
TWI745871B (zh) * | 2020-02-24 | 2021-11-11 | 佳世達科技股份有限公司 | 顯示裝置及顏色調整方法 |
JP7428090B2 (ja) | 2020-06-25 | 2024-02-06 | セイコーエプソン株式会社 | 校正装置、校正方法、校正プログラム、及び分光カメラ |
US11615763B2 (en) | 2020-11-05 | 2023-03-28 | Samsung Electronics Co., Ltd. | Color gamut compression and extension |
CN114792363B (zh) * | 2022-04-19 | 2023-07-11 | 江南大学 | 三原色纤维构建的全色域网格化混色模型构建方法及彩色纺纱方法 |
CN114905851B (zh) * | 2022-07-15 | 2022-11-11 | 浙江工业大学 | 一种印刷色密度特征曲线采集方法 |
CN115174881B (zh) * | 2022-07-15 | 2024-02-13 | 深圳市火乐科技发展有限公司 | 色域映射方法、装置、投影设备及存储介质 |
CN117253459B (zh) * | 2023-11-20 | 2024-01-09 | 汕头市超声仪器研究所股份有限公司 | 一种对液晶显示器进行Gamma调校的方法 |
CN117395825B (zh) * | 2023-11-24 | 2024-06-14 | 北京同步风云科技有限公司 | 一种基于实时led亮度和色差的校正软件控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5991456A (en) * | 1996-05-29 | 1999-11-23 | Science And Technology Corporation | Method of improving a digital image |
CN101137070A (zh) * | 2006-08-28 | 2008-03-05 | 华为技术有限公司 | 视频通讯中视频输入设备伽玛特性校正方法和装置 |
CN101421758A (zh) * | 2006-04-18 | 2009-04-29 | 国立大学法人东京工业大学 | 图像变换装置和图像变换程序 |
CN101635858A (zh) * | 2008-07-23 | 2010-01-27 | 胜华科技股份有限公司 | 色彩校正方法与应用其的集成型芯片 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243414A (en) * | 1991-07-29 | 1993-09-07 | Tektronix, Inc. | Color processing system |
US7728845B2 (en) * | 1996-02-26 | 2010-06-01 | Rah Color Technologies Llc | Color calibration of color image rendering devices |
JP2000278546A (ja) * | 1999-01-22 | 2000-10-06 | Sony Corp | 画像処理装置及び画像処理方法、色域変換テーブル作成装置及び色域変換テーブル作成方法、画像処理プログラムを記録した記録媒体、並びに色域変換テーブル作成プログラムを記録した記録媒体 |
KR20060120643A (ko) * | 2003-09-11 | 2006-11-27 | 마쯔시다덴기산교 가부시키가이샤 | 화상처리 장치, 화상처리 방법 및 화상처리 프로그램 |
US6980219B2 (en) * | 2003-10-21 | 2005-12-27 | Clairvoyante, Inc | Hue angle calculation system and methods |
KR20050112319A (ko) * | 2004-05-25 | 2005-11-30 | 삼성전자주식회사 | 원형 색역 사상에 따른 색변환 장치 및 그 방법 |
JP4958638B2 (ja) * | 2007-05-28 | 2012-06-20 | キヤノン株式会社 | 色処理装置および方法 |
US7859721B2 (en) * | 2008-01-15 | 2010-12-28 | Adobe Systems Incorporated | Smoothly changing gray levels in color transforms |
JP4957676B2 (ja) * | 2008-03-19 | 2012-06-20 | セイコーエプソン株式会社 | 画像データ解析装置、画像データ解析方法、およびプログラム |
CN101621700B (zh) * | 2008-10-21 | 2011-04-13 | 青岛海信电器股份有限公司 | 一种多媒体设备的色域匹配方法及电视机 |
US9049410B2 (en) * | 2009-12-23 | 2015-06-02 | Samsung Display Co., Ltd. | Color correction to compensate for displays' luminance and chrominance transfer characteristics |
-
2011
- 2011-10-17 US US13/818,303 patent/US20140152687A1/en not_active Abandoned
- 2011-10-17 WO PCT/CN2011/001729 patent/WO2013056387A1/zh active Application Filing
- 2011-10-17 CN CN201180042671.9A patent/CN103931169B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5991456A (en) * | 1996-05-29 | 1999-11-23 | Science And Technology Corporation | Method of improving a digital image |
CN101421758A (zh) * | 2006-04-18 | 2009-04-29 | 国立大学法人东京工业大学 | 图像变换装置和图像变换程序 |
CN101137070A (zh) * | 2006-08-28 | 2008-03-05 | 华为技术有限公司 | 视频通讯中视频输入设备伽玛特性校正方法和装置 |
CN101635858A (zh) * | 2008-07-23 | 2010-01-27 | 胜华科技股份有限公司 | 色彩校正方法与应用其的集成型芯片 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110490328A (zh) * | 2019-08-06 | 2019-11-22 | 刘世昌 | 以统一场论和全息电视为基础实现量子计算机的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103931169A (zh) | 2014-07-16 |
US20140152687A1 (en) | 2014-06-05 |
WO2013056387A1 (zh) | 2013-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103931169B (zh) | 一种通用的色域映射及色彩管理方法 | |
JP7075687B2 (ja) | 色変換システム、色変換装置及び色変換方法 | |
Hardeberg | Acquisition and reproduction of color images: colorimetric and multispectral approaches | |
KR100255408B1 (ko) | 색상재현방법 | |
CN1323373C (zh) | 投影机的投影平面色彩校正方法、*** | |
JP6687387B2 (ja) | 拡張色域 | |
US20040056965A1 (en) | Method for color correction of digital images | |
CN102870418B (zh) | 按色度和视觉保真原则传送及显示图像的方法 | |
CN101449589A (zh) | 彩色图像显示设备和色彩变换设备 | |
JP5188082B2 (ja) | 画像出力機器の色変換定義の作成方法、作成装置及びそのプログラム | |
WO2005057941A1 (ja) | プロジェクタの色補正方法 | |
CN104956653A (zh) | 颜色分离 | |
WO2006020632A2 (en) | Digital color fidelity | |
CN104010176B (zh) | 一种利用人眼实现从显示器到打印机颜色校正方法 | |
Stone | Representing colors as three numbers [color graphics] | |
JP2010268294A (ja) | カラーマネジメントシステム及びプログラム | |
JPH10164381A (ja) | 色変換方法 | |
Deshpande | N-colour separation methods for accurate reproduction of spot colours | |
Gofaizen et al. | Adaptation to observation conditions in television systems by means of signal correction | |
JP3427615B2 (ja) | 色再現方法 | |
US6909523B2 (en) | Information color system and printing method in the information color system | |
JPH10164352A (ja) | 画像信号処理方法 | |
Yanzhe | A novel algorithm based on improved BP neural network and its application in color management | |
Vladimirova | COLOR MANAGEMENT TECHNIQUES AND COLOR RENDERING IN THE VISUAL COMMUNICATION DESIGN | |
Kim et al. | Applying CIECAM97s in a color management system |
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 |