CN104112081B - 色域边界描述方法 - Google Patents
色域边界描述方法 Download PDFInfo
- Publication number
- CN104112081B CN104112081B CN201410371181.9A CN201410371181A CN104112081B CN 104112081 B CN104112081 B CN 104112081B CN 201410371181 A CN201410371181 A CN 201410371181A CN 104112081 B CN104112081 B CN 104112081B
- Authority
- CN
- China
- Prior art keywords
- point
- value
- gbd
- color
- face
- 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
Landscapes
- Color Image Communication Systems (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
本发明提供一种色域边界描述方法,其特征在于,包括以下步骤:步骤1,利用球面坐标系表示待描述色域;步骤2,将CIELAB标准颜色空间划分为多个分区;步骤3,查找各个分区的最外围颜色样本点GBD点;步骤4,计算待映射颜色点P点的球坐标值;步骤5,计算各GBD点在ψ面上的投影点;步骤6,计算各GBD点在ψ面上的方向角;步骤7,调整各GBD点的方向角和色相角;步骤8,将GBD点总集划分为4个子集;步骤9,在4个子集内分别计算三角形插值的候选GBD点;步骤10,根据最小γ值的GBD点和P点的球坐标值构件颜色点;步骤11,选择3个满足三角形插值的候选GBD点;步骤12,通过三角形插值计算P点对应的色域边界点P*;步骤13,计算其他待映射颜色点对应的色域边界点。
Description
技术领域
本发明涉及计算机图像处理及色彩管理领域,具体涉及一种彩色图像及其显示设备和输出设备的色域边界描述方法。
背景技术
色彩管理是保证色彩信息在跨媒介传输时准确再现的一种技术,其核心是在色彩正确可视化的前提下实现不同颜色空间的匹配,即色域映射。色域映射实施的前提则是色域边界的准确提取与描述。
M.Mahy提出一种基于钮阶堡方程的分区色域描述模型,模型将n色输出设备色域分解成若干个含n-1个基色的分部,对这些分部逐个描述进而完成总体色域的描述,但该方法仅适用于打印输出设备。Raja B.提出一种基于改进凸壳算法的色域描述方法,该方法首先在待描述色域内进行采样,然后使用凸壳算法提取采样点集中处于色域边界的样本点,但由于该算法往往将部分色域内部采样点错判为色域边界点,因此计算出的色域边界与实际色域存在一定的偏差。J.Morovic提出了分区最大边界描述法,将色域边界的描述分为色域边界描述器的计算与线色域边界的计算两步,但该算法在具体实施过程中需要对数据进行两次插值,一是对不存在采样点的色域分区通过插值得到对应分区的边界点数据,二是通过插值得到相应的线色域边界点,此外该算法在求解线色域边界时还需要对插值得到的线色域边界点再进行一次拟合,过多的插值与拟合计算在一定程度上引入了误差,影响了色域边界描述的准确度。
发明内容
本发明是针对上述问题进行的,目的在于提供一种色域边界描述方法,实现在色域映射过程中对彩色图像及其显示和输出设备的色域边界的准确描述。
本发明为实现上述目的,采用了以下的技术方案:
本发明提供一种色域边界描述方法,其特征在于,包括以下步骤:
步骤1,在待描述色域内对颜色样本点进行颜色采样,获得色域颜色样本集,将该色域颜色样本集中的颜色样本值转换到CIELAB标准颜色空间,并通过下列公式计算每个颜色样本点的球面坐标值:
式中,(L,a,b)是颜色样本点对应的CIELAB颜色值,(LE,a,b)是球坐标系原点E对应的CIELAB颜色值,(γ,α,θ)是颜色样本点的球面坐标系值,γ是颜色样本点到球坐标系原点E的距离,α是色相角,在0~360°范围内,θ是色相角为α时的平面内角度,在-90°~90°范围内;
步骤2,采用球坐标表示CIELAB标准颜色空间,通过对球坐标系中α和θ值的均匀划分将CIELAB标准颜色空间划分为多个分区;
步骤3,在步骤1所得的色域颜色样本集中分别查找步骤2所得的各个分区中的最外围颜色样本点,该最外围颜色样本点具有最大的γ值,将各个分区的最外围颜色样本点作为相应的分区的色域边界点GBD点,对于不包含颜色样本点的分区不作处理;
步骤4,将一个待映射颜色点记为P点,计算该P点对应的CIELAB颜色值,记为(Li,ai,bi),通过式1计算P点的球面坐标值,记为(γi,αi,θi);
步骤5,将P点的αi值所确定的色相面记为αi面,将与该αi面相隔180°的色相面记为αi'面,将αi面和αi'面共同组成的球坐标面记为ψ面,采用该ψ面将球坐标系切成两个相等的部分,将步骤3得到的各个GBD点投影到ψ面上,得到对应的投影点GBD'点,将GBD'点的CIELAB颜色值记为(L',a',b'),采用下列等式计算GBD'点在ψ面上的彩度值:
式2
式中,C是GBD'点在ψ面上的彩度值,(a',b')是GBD'点的CIELAB色度值;
步骤6,在ψ面上以球坐标系原点E点为原点,以C为横坐标,以L为纵坐标建立直角坐标系,记为C-L坐标系,并在该C-L坐标系内计算各个GBD点在ψ面上的方向角β,计算公式如下:
β=tan-1((L'-LE)/C) 式3
式中,L'是各个GBD'点的CIELAB亮度值,LE是C-L坐标系原点的CIELAB亮度值,β在0~360°范围内;
步骤7,根据式2和式3计算P点在ψ面上的方向角βi,并以该βi和步骤4计算得到的色相角αi为基准,对各个GBD点的方向角β值及色相角α值进行如下调整:
式4
式中,α、β分别表示各个GBD点的色相角与ψ面上的方向角,αi、βi分别是P点的色相角与ψ面上的方向角,α'、β'分别是调整后的GBD点的色相角与ψ面上的方向角;
步骤8,根据α'、β'值,将所有的GBD点总集划分为四个子集,分别记为子集1、子集2、子集3以及子集4,使子集1所包含的GBD点的α'和β'值均在0~180°之间,子集2所包含的GBD点的α'值在0~180°之间而β'值在180°~360°之间,子集3所包含的GBD点的α'值在180°~360°之间而β'值在0~180°之间,子集4所包含的α'和β'值均在180°~360°之间;
步骤9,在步骤8所得的四个子集中分别计算各个GBD点与P点的球坐标差Δα和Δθ,并选取具有最小(Δα+Δθ)/2值的GBD点作为三角形插值的候选分区色域边界点,共得到四个候选分区色域边界点,Δα和Δθ的计算公式如下:
式5
式中,α、θ是各个GBD点的球坐标值,αi、θi是P点的球坐标值;
步骤10,将步骤9所得的四个候选分区色域边界点中γ最小的记为γmin,以γmin/2和P点的αi和θi值作为球坐标值构建颜色点,记为P'点,则该P'点位于待描述的色域内,且P'点、P点以及球坐标系原点E三点共线;
步骤11,按照如下公式判断任意三个候选分区色域边界点是否满足三角形插值:
式中,L1、a1、b1、L2、a2、b2、L3、a3、b3分别是参与判断的三个候选分区色域边界点的CIELAB颜色值,Li'、ai'、bi'是步骤10得到的P'点的CIELAB颜色值,LE、aE、bE是球坐标系原点E点的CIELAB颜色值,λ、μ、ν是式6计算出的判断系数,当λ≥0,μ≥0,ν≥0,且λ+μ+ν≤1时,P'点位于球坐标系原点E点和所选的三个候选分区色域边界点所构成的四面体内,即E、P'两点连线与三个候选分区色域边界点所构成的空间三角形相交于一点,该交点即为待映射颜色点P点所对应的色域边界点P*,此时所选的三个候选分区色域边界点符合三角形插值的要求;
步骤12,根据步骤11中选取的三个候选分区色域边界点和P'点的CIELAB颜色值,通过三角形插值得到P点所对应的色域边界点P*,插值过程如下:
其中, 式8
式7中,L*、a*、b*为待映射颜色点P点所对应的色域边界点P*的CIELAB颜色值;
步骤13,重复步骤4~12,计算其他所有的待映射颜色点所对应的色域边界点的CIELAB颜色值。
发明的作用与效果
根据本发明所提供的色域边界描述方法,因为在色域映射过程中,采用球坐标系方便、准确地描述颜色样本点和待映射颜色点,并且只需进行一次三角形插值便可以得到待映射颜色点所对应的图像或设备色域边界,因此该方法精度高,计算速度快,并且执行效率高。
附图说明
图1是色域边界描述方法流程图;
图2是实施例中的I23标准图像;
图3是三维空间球坐标系示意图;
图4是ψ面上的C-L直角坐标系示意图;以及
图5是符合插值条件的GBD点与待映射颜色点的空间关系示意图。
具体实施方式
以下结合附图,对本发明所提供的色域边界描述方法作详细阐述。
<实施例>
在本实施例中,待描述色域是HP Designjet Z3200彩色喷墨打印机的设备色域和LCD显示器显示质量评价图形库中的标准图像I23的图像色域。
图1是色域边界描述方法流程图。
如图1所示,色域边界描述方法包括以下步骤:
步骤S-1,本实施例中待描述的HP Designjet Z3200彩色喷墨打印机色域属于输出设备色域,根据颜色分区理论,在该打印机的3个对应分区(CMK,CYK以及MYK)内进行颜色采样,获得设备色域颜色样本集;然后通过Matlab读取打印机ICC文件内置的“AToB1”标签数据,将样本集的网点面积率数值转换为对应的Lab数值,实现样本颜色到标准颜色空间CIELAB的转换。
图2是实施例中的I23标准图像。
如图2所示,待描述的I23标准图的色域属于图像色域,图像大小为512×384,共有196608个颜色像素点,原始颜色模式为RGB模式,在Photoshop图像处理软件中将该图的配置文件设置为“AdobeRGB1998.icc”,并将其颜色转换到CIELAB色空间。
表1所示打印机CMK分区的采样情况,以0.05为步长在C、M、K三个通道上网点面积率为0~1的范围内分别进行等间隔循环采样,共得到9261个采样点颜色,得到各采样点的青红(C)、品红(M)和黑色(K)的网点面积率及相应的Lab颜色值。
表1 CMK分区中采样点网点面积率和对应的Lab值
图3是三维空间球坐标系示意图。
如图3所示,利用球坐标表示CIELAB标准颜色空间,根据如下公式计算打印机色域颜色样本与图像像素颜色对应的球坐标值(γ,α,θ):
式中,(L,a,b)是颜色样本点对应的CIELAB颜色值,(LE,a,b)是球坐标系原点E对应的CIELAB颜色值,(γ,α,θ)是颜色样本点的球面坐标系值,γ是颜色样本点到球坐标系原点E的距离,α是色相角,在0~360°范围内,θ是色相角为α时的平面内角度,在-90°~90°范围内。在本实施例中,E点取值为(50,0,0)。
步骤S-2,分别对α和θ在其取值范围内均匀地划分为16份,从而将CIELAB标准颜色空间分为了256个分区。
步骤S-3,在计算过球坐标值的打印机色域采样点集合与图像像素颜色集合中分别查找步骤2所分的256个分区的色域边界点,记为GBD点,GBD点的α和θ值应该在各个分区的α和θ值范围内,且在该分区内具有最大的γ值。对于不包含采样点的分区则不作处理。
表2 打印机分区色域边界点的lab值及球坐标值
序号 | L值 | a值 | b值 | γ值 | α值 | θ值 |
…… | …… | …… | …… | …… | …… | …… |
33 | 14.39 | 22.21 | 7.88 | 42.70 | 19.53° | -56.51° |
34 | 11.16 | 15.99 | 6.73 | 42.54 | 22.82° | -65.93° |
35 | 21.06 | 13.00 | 14.13 | 34.73 | 47.38° | -56.44° |
…… | …… | …… | …… | …… | …… | …… |
253 | 96.15 | 0.93 | -4.12 | 46.34 | 282.77° | 84.77° |
254 | 93.82 | 4.92 | -6.00 | 44.50 | 309.35° | 79.95° |
255 | 97.79 | 3.79 | -1.79 | 47.97 | 334.77° | 84.99° |
在本实施例中,如表2所示,从打印机色域采样点集合中计算得到了255个GBD点;如表3所示,从图像像素颜色集合中计算得到了235个GBD点。
表3 I23标准图像分区色域边界点的lab值及球坐标值
序号 | L值 | a值 | b值 | γ值 | α值 | θ值 |
…… | …… | …… | …… | …… | …… | …… |
56 | 27.69 | -16.92 | 0.82 | 28.01 | 177.24° | -52.79° |
57 | 26.95 | -15.61 | -1.07 | 27.86 | 183.91° | -55.83° |
58 | 39.55 | -7.22 | -4.51 | 13.48 | 269.29° | -50.86° |
…… | …… | …… | …… | …… | …… | …… |
233 | 98.85 | 1.00 | -1.77 | 48.90 | 299.38° | 87.62° |
234 | 99.79 | 0.59 | -0.32 | 49.80 | 331.79° | 89.23° |
235 | 99.35 | 1.65 | -0.47 | 49.38 | 344.16° | 88.01° |
步骤S-4,在I23标准图中选取一个待映射颜色点,记为P点,计算P点对应的CIELAB颜色值,记为(Li,ai,bi),并通过式1计算P点的球面坐标值,记为(γi,αi,θi)。在本实施例中,选取的P点是I23标准图像中位置坐标为(200,200)的像素点颜色,其对应的CIELAB颜色值为(75.38,-21.02,-6.82),相应的球坐标系值为(33.65,197.99°,48.95°)。
步骤S-5,将P点的αi值(αi=197.99°)所确定的色相面记为αi面,将与该αi面相隔180°的色相面记为αi'面(αi'=17.99°),将αi面和αi'面共同组成的球坐标面记为ψ面,采用该ψ面将球坐标系切成两个相等的部分,将步骤3得到的各个GBD点(如表2、3所示)投影到ψ面上,得到对应的投影点GBD'点,将GBD'点的CIELAB颜色值记为(L',a',b'),采用下列等式计算GBD'点在ψ面上的彩度值:
式2
式中,C是GBD'点在ψ面上的彩度值,(a',b')是GBD'点的CIELAB色度值;
步骤S-6,图4是ψ面上的C-L直角坐标系示意图。
在ψ面上以球坐标系原点E(50,0,0)点为原点,以C为横坐标,以L为纵坐标建立直角坐标系,记为C-L坐标系,如图4所示。并在该C-L坐标系内计算打印机及图像的各个GBD点在ψ面上的方向角β,计算公式如下:
β=tan-1((L'-LE)/C) 式3
式中,L'是各个GBD'点的CIELAB亮度值,LE是C-L坐标系原点的CIELAB亮度值,β在0~360°范围内。
表4 打印机分区色域边界投影后的lab值、彩度值C及方向角β
序号 | L′值 | a′值 | b′值 | C值 | β值 |
…… | …… | …… | …… | …… | …… |
33 | 14.39 | 22.40 | 7.27 | -23.55 | 236.52° |
34 | 11.16 | 16.44 | 5.34 | -17.29 | 246.01° |
35 | 21.06 | 15.91 | 5.17 | -16.73 | 239.97° |
…… | …… | …… | …… | …… | …… |
253 | 96.15 | -0.37 | -0.12 | 0.39 | 89.52° |
254 | 93.82 | 2.69 | 0.87 | -2.83 | 93.69° |
255 | 97.79 | 2.90 | 0.94 | -3.05 | 93.66° |
打印机及图像各个GBD'点的CIELAB颜色值(L',a',b')、彩度值C及GBD在ψ面上的方向角β分别如表4和表5所示。
表5 I23标准图像分区色域边界点投影后的lab值、彩度值C及方向角β
序号 | L′值 | a′值 | b′值 | C值 | β值 |
…… | …… | …… | …… | …… | …… |
56 | 27.69 | -15.06 | -4.89 | 15.84 | 305.37° |
57 | 26.95 | -14.44 | -4.69 | 15.18 | 303.36° |
58 | 39.55 | -7.85 | -2.55 | 8.25 | 308.30° |
…… | …… | …… | …… | …… | …… |
233 | 98.85 | 0.38 | 0.12 | -0.40 | 90.47° |
234 | 99.79 | 0.44 | 0.14 | -0.46 | 90.53° |
235 | 99.35 | 1.36 | 0.44 | -1.43 | 91.66° |
步骤S-7,根据式2和式3计算P点在ψ面上的方向角βi,并以该βi值和色相角αi值为基准,对打印机和图像的各个GBD点的方向角β值及色相角α值进行如下调整:
式4
式中,α、β分别表示各个GBD点的色相角与ψ面上的方向角,αi、βi分别是P点的色相角与ψ面上的方向角,α'、β'分别是调整后的GBD点的色相角与ψ面上的方向角。
在本实施例中,αi=197.99°,βi=48.95°,调整好的打印机及图像的各个GBD点的α和β值如表6所示。
表6 调整后打印机及图像分区色域边界点的α'和β'值
序号 | 打印机α′ | 打印机β′ | 图像α′ | 图像β′ |
…… | …… | …… | …… | …… |
33 | 181.54° | 187.56° | 181.36° | 196.83° |
34 | 184.82° | 197.05° | 195.70° | 198.62° |
35 | 209.39° | 191.01° | 210.41° | 196.60° |
…… | …… | …… | …… | …… |
56 | 338.16° | 262.02° | 339.25° | 256.42° |
57 | 342.34° | 262.07° | 345.92° | 254.41° |
58 | 8.97° | 259.02° | 14.01° | 259.34° |
…… | …… | …… | …… | …… |
步骤S-8,根据如表6所示的调整后的色相角α'和方向角β'值,分别将打印机和图像的GBD点总集分为4个子集,依次记为子集1、子集2、子集3以及子集4。其中子集1所包含的GBD点的α'和β'值均在0~180°之间,子集2所包含的GBD点的α'值在0~180°之间而β'值在180°~360°之间,子集3所包含的GBD点的α'值在180°~360°之间而β'值在0~180°之间,子集4所包含的α'和β'值均在180°~360°之间。本实施例中,打印机子集1~4内的GBD点的数量分别为58、62、70、65,图像的子集1~4内的GBD点的数量分别为48、56、71、60。
步骤S-9,在步骤S-8所得的打印机和图像的GBD点子集中分别计算各个GBD点与P点的球坐标差Δα和Δθ,并选取具有最小(Δα+Δθ)/2值的GBD点作为三角形插值的候选GBD点,共得到四个候选GBD点,Δα和Δθ的计算公式如下:
式5
式中,α、θ是各个GBD点的球坐标值,αi、θi是P点的球坐标值。
本实施例中计算出的打印机和图像候选插值GBD点的CIELAB颜色值及对应的Δα、Δθ和γ值如表7和表8所示。
表7 打印机候选GBD点的CIELAB颜色值及Δα、Δθ和γ值
序号 | 所属子集 | L值 | a值 | b值 | Δα值 | Δθ值 | γ值 |
217 | 1 | 89.21 | -19.70 | -6.79 | 1.02° | 13.06° | 44.40 |
186 | 2 | 79.98 | -31.73 | -27.88 | 23.31° | 13.59° | 51.79 |
201 | 3 | 84.42 | -29.31 | -3.66 | 10.87° | 0.42° | 45.35 |
185 | 4 | 78.68 | -41.25 | -6.86 | 8.54° | 14.51° | 50.71 |
表8 I23标准图像候选GBD点的CIELAB颜色值及Δα、Δθ和γ值
序号 | 所属子集 | L值 | a值 | b值 | Δα值 | Δθ值 | γ值 |
212 | 1 | 92.11 | -8.38 | -3.45 | 4.36° | 28.90° | 43.08 |
183 | 2 | 74.69 | -20.62 | -9.57 | 6.92° | 1.60° | 33.56 |
196 | 3 | 84.45 | -17.80 | -3.14 | 7.80° | 13.36° | 38.90 |
182 | 4 | 79.39 | -26.20 | -7.80 | 1.41° | 1.88° | 40.14 |
步骤S-10,将步骤S-9所得的候选GBD点(如表7、8所示)中γ最小的记为γmin,以γmin/2和P点的αi和θi值作为球坐标值构建颜色点,记为P'点,则该P'点一定位于待描述的色域内,且P'点、P点以及球坐标系原点E三点共线。
本实施例中,打印机候选GBD点的γmin=44.40,构建出的打印机P'点的球坐标值为(22.20,197.99°,48.95°),对应的CIELAB颜色值为(66.74,-13.87,-4.50)。图像候选GBD点的γmin=33.56,构建出的打印机P'点的球坐标值为(16.78,197.99°,48.95°),对应的CIELAB颜色值为(62.66,-10.48,-3.40)。
步骤S-11,分别从如表7和8所示的打印机和图像候选GBD点中选取3个满足三角形插值要求的候选GBD点。
图5是符合插值条件的GBD点与待映射颜色点的空间关系示意图。
如图5所示,当步骤S-10构件的P'点位于球坐标系原点E与所选3个候选GBD点构成的四面体内时,E、P'点连线与所选3个候选GBD点构成的空间三角形相交于一点,该交点即为待映射颜色点P点对应的色域边界点,记为P*。按照如下公式判断任意3个候选GBD点是否满足三角形插值:
式中,L1、a1、b1、L2、a2、b2、L3、a3、b3分别是参与判断的三个候选分区色域边界点的CIELAB颜色值,Li'、ai'、bi'是步骤10得到的P'点的CIELAB颜色值,LE、aE、bE是球坐标系原点E点的CIELAB颜色值,λ、μ、ν是式6计算出的判断系数,当λ≥0,μ≥0,ν≥0,且λ+μ+ν≤1时,P'点位于球坐标系原点E点和所选的3个候选GBD点所构成的四面体内,即E、P'两点连线与3个候选GBD点所构成的空间三角形相交于一点,此时所选的3个候选GBD点符合三角形插值的要求。
本实施例中符合插值条件的3个打印机候选GBD点为表7中子集1、2、3的候选GBD点,符合插值条件的3个图像候选GBD点为表8中子集1、2、4的候选GBD点。
步骤S-12,根据步骤S-11中选取的符合插值条件的候选GBD点和P'点的CIELAB颜色值,通过三角形插值得到P点所对应的色域边界点P*,插值过程如下:
其中, 式8
式7中,L*、a*、b*为待映射颜色点P点所对应的色域边界点P*的CIELAB颜色值。
经计算,实施例所选的P点对应的打印机色域边界点P*的CIELAB颜色值为(84.22,-28.34,-9.20),对应的图像色域边界点P*的CIELAB颜色值为(79.27,-24.24,-7.87)。
步骤S-13,重复步骤S-4~12,计算其他所有的待映射颜色点所对应的打印机和图像色域边界点的CIELAB颜色值。
当然,本实施例仅为本发明的色域边界描述方法的优选案例,并不用于限制本发明的保护范围。
Claims (1)
1.一种色域边界描述方法,其特征在于,包括以下步骤:
步骤1,在待描述色域内对颜色样本点进行颜色采样,获得色域颜色样本集,将该色域颜色样本集中的颜色样本值转换到CIELAB标准颜色空间,并通过下列公式计算每个所述颜色样本点的球面坐标值:
式中,(L,a,b)是所述颜色样本点对应的CIELAB颜色值,(LE,a,b)是球坐标系原点E对应的CIELAB颜色值,(γ,α,θ)是所述颜色样本点的球面坐标系值,γ是所述颜色样本点到球坐标系原点E的距离,α是色相角,在0~360°范围内,θ是色相角为α时的平面内角度,在-90°~90°范围内;
步骤2,采用球坐标表示CIELAB标准颜色空间,通过对球坐标系中α和θ值的均匀划分将CIELAB标准颜色空间划分为多个分区;
步骤3,在所述步骤1所得的所述色域颜色样本集中分别查找所述步骤2所得的各个所述分区中的最外围颜色样本点,该最外围颜色样本点具有最大的γ值,将各个所述分区的所述最外围颜色样本点作为相应的所述分区的色域边界点GBD点,对于不包含所述颜色样本点的所述分区不作处理;
步骤4,将一个待映射颜色点记为P点,计算该P点对应的CIELAB颜色值,记为(Li,ai,bi),通过所述式1计算所述P点的球面坐标值,记为(γi,αi,θi);
步骤5,将所述P点的色相角αi所确定的色相面记为αi面,将与该αi面相隔180°的色相面记为αi'面,将所述αi面和所述αi'面共同组成的球坐标面记为ψ面,采用该ψ面将球坐标系切成两个相等的部分,将所述步骤3得到的各个所述GBD点投影到所述ψ面上,得到对应的投影点GBD'点,将所述GBD'点的CIELAB颜色值记为(L',a',b'),采用下列等式计算所述GBD'点在所述ψ面上的彩度值:
式中,C是所述GBD'点在所述ψ面上的彩度值,(a',b')是所述GBD'点的CIELAB色度值;
步骤6,在所述ψ面上以球坐标系原点E点为原点,以所述C为横坐标,以所述L为纵坐标建立直角坐标系,记为C-L坐标系,并在该C-L坐标系内计算各个所述GBD点在所述ψ面上的方向角β,计算公式如下:β=tan-1((L'-LE)/C) 式3
式中,L'是各个所述GBD'点的CIELAB亮度值,LE是所述C-L坐标系原点的CIELAB亮度值,β在0~360°范围内;
步骤7,根据所述式2和所述式3计算所述P点在所述ψ面上的方向角βi,并以该βi和所述步骤4计算得到的所述色相角αi为基准,对各个所述GBD点的所述方向角β值及色相角α值进行如下调整:
式中,α、β分别表示各个所述GBD点的色相角与所述ψ面上的方向角,αi、βi分别是所述P点的色相角与所述ψ面上的方向角,α'、β'分别是调整后的所述GBD点的色相角与所述ψ面上的方向角;
步骤8,根据所述α'、β'值,将所有的所述GBD点总集划分为四个子集,分别记为子集1、子集2、子集3以及子集4,使所述子集1所包含的GBD点的α'和β'值均在0~180°之间,所述子集2所包含的GBD点的α'值在0~180°之间而β'值在180°~360°之间,所述子集3所包含的GBD点的α'值在180°~360°之间而β'值在0~180°之间,所述子集4所包含的α'和β'值均在180°~360°之间;
步骤9,在所述步骤8所得的四个所述子集中分别计算各个所述GBD点与所述P点的球坐标差Δα和Δθ,并选取具有最小(Δα+Δθ)/2值的GBD点作为三角形插值的候选分区色域边界点,共得到四个候选分区色域边界点,所述Δα和Δθ的计算公式如下:
式中,α是各个所述GBD点的色相角、θ是色相角为α时的平面内角度,αi是所述P点的色相角、θi是所述P点的平面内角度;
步骤10,将所述步骤9所得的四个所述候选分区色域边界点中γ最小的记为γmin,以和所述P点的αi和θi值作为球坐标值构建颜色点,记为P'点,则该P'点位于所述待描述的色域内,且所述P'点、所述P点以及球坐标系原点E三点共线;
步骤11,按照如下公式判断任意三个所述候选分区色域边界点是否满足三角形插值:
式中,L1、a1、b1、L2、a2、b2、L3、a3、b3分别是参与判断的三个所述候选分区色域边界点的CIELAB颜色值,Li'、ai'、bi'是所述步骤10得到的所述P'点的CIELAB颜色值,LE、aE、bE是球坐标系原点E点的CIELAB颜色值,λ、μ、ν是所述式6计算出的判断系数,当λ≥0,μ≥0,ν≥0,且λ+μ+ν≤1时,所述P'点位于球坐标系原点E点和所选的三个所述候选分区色域边界点所构成的四面体内,即E、P'两点连线与三个所述候选分区色域边界点所构成的空间三角形相交于一点,该交点即为所述待映射颜色点P点所对应的色域边界点P*,此时所选的三个所述候选分区色域边界点符合三角形插值的要求;
步骤12,根据所述步骤11中选取的三个所述候选分区色域边界点和所述P'点的CIELAB颜色值,通过三角形插值得到所述P点所对应的色域边界点P*,插值过程如下:
其中,
所述式7中,L*、a*、b*为所述待映射颜色点P点所对应的色域边界点P*的CIELAB颜色值;
步骤13,重复所述步骤4~12,计算其他所有的待映射颜色点所对应的色域边界点的CIELAB颜色值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410371181.9A CN104112081B (zh) | 2014-07-31 | 2014-07-31 | 色域边界描述方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410371181.9A CN104112081B (zh) | 2014-07-31 | 2014-07-31 | 色域边界描述方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104112081A CN104112081A (zh) | 2014-10-22 |
CN104112081B true CN104112081B (zh) | 2017-02-08 |
Family
ID=51708868
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410371181.9A Expired - Fee Related CN104112081B (zh) | 2014-07-31 | 2014-07-31 | 色域边界描述方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104112081B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825020A (zh) * | 2016-03-23 | 2016-08-03 | 天津师范大学 | 三维可感知色域计算方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105654455B (zh) * | 2014-11-11 | 2018-06-15 | 曲阜师范大学 | 一种“图像到设备”的色域映射算法 |
CN104702815A (zh) * | 2015-03-26 | 2015-06-10 | 智川 | 一种彩色输出设备色域可视化方法 |
EP3110124A1 (en) * | 2015-06-25 | 2016-12-28 | Thomson Licensing | Hue changing color gamut mapping |
US9948831B2 (en) * | 2015-08-31 | 2018-04-17 | Kyocera Document Solutions Inc. | In/out-of-gamut determination system for acquiring cross-section of color gamut of image forming apparatus taken on specific plane |
CN106408619B (zh) * | 2016-09-13 | 2019-04-26 | 齐鲁工业大学 | 一种基于光谱域来实现跨媒体颜色再现的方法 |
CN106600659B (zh) * | 2016-12-05 | 2019-06-04 | 浙江工商大学 | 基于随机采样点的多样化颜色表绘制方法 |
CN107680141B (zh) * | 2017-09-21 | 2018-11-23 | 深圳市华星光电半导体显示技术有限公司 | 一种色域映射方法及其色域映射装置 |
CN107888893A (zh) * | 2017-11-07 | 2018-04-06 | 深圳市华星光电半导体显示技术有限公司 | 一种色域映射方法及色域映射装置 |
CN109345457B (zh) * | 2018-10-24 | 2022-09-27 | 河南工程学院 | 一种用于分区最大化色域边界描述器的空区插值方法 |
CN112907721B (zh) * | 2021-03-18 | 2022-06-07 | 武汉纺织大学 | 一种色纺纱线原色体系构建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1753452A (zh) * | 2004-09-21 | 2006-03-29 | 三星电子株式会社 | 检测色域边界及使用其来映射色域的设备和方法 |
-
2014
- 2014-07-31 CN CN201410371181.9A patent/CN104112081B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1753452A (zh) * | 2004-09-21 | 2006-03-29 | 三星电子株式会社 | 检测色域边界及使用其来映射色域的设备和方法 |
Non-Patent Citations (4)
Title |
---|
"Calculating Medium and Image Gamut Boundaries for Gamut Mapping";Morovic J 等;《COLOR research and application》;20001231;第25卷(第6期);全文 * |
"Colorimetric calibration in electronic imaging devices using a look-up-table model and interpolations";Po-Chieh Hung;《Journal of Electronic Imaging》;19930131;第2卷(第1期);全文 * |
"基于区域分割法的图像色域边界描述";冀利利等;《包装工程》;20140531;第35卷(第9期);全文 * |
"彩色复制中的色域映射方法";王淑真等;《北京印刷学院学报》;20060630;第14卷(第3期);全文 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825020A (zh) * | 2016-03-23 | 2016-08-03 | 天津师范大学 | 三维可感知色域计算方法 |
CN105825020B (zh) * | 2016-03-23 | 2019-04-16 | 天津师范大学 | 三维可感知色域计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104112081A (zh) | 2014-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104112081B (zh) | 色域边界描述方法 | |
US9661187B1 (en) | Color gamut mapping method based on color gamut of source image | |
JP4990507B2 (ja) | 映像の主要色成分調整装置、映像の主要色成分調整方法及び映像の主要色成分調整プログラムを記録するコンピュータで読み取り可能な記録媒体 | |
US20160035066A1 (en) | Method and apparatus for image color enhancement | |
US7663788B2 (en) | Image correcting apparatus and method, and image correction program | |
US9616659B2 (en) | Image processing apparatus that determines ejection execution direction of print head | |
EP3327412B1 (en) | Inspection apparatus, inspection method, and non-transitory computer readable medium for storing inspection program of colorimetric value | |
US20050276473A1 (en) | Apparatus and method of detecting color gamut in color device and calculating color space inverse transform function | |
CN102075667B (zh) | 一种基于查表法的反向色彩空间转换方法 | |
US8373721B2 (en) | Method of realism assessment of an image composite | |
JP2014146909A (ja) | 画像処理装置、画像処理方法およびプログラム | |
CN110888612B (zh) | 打印目标色域的评估方法、装置、设备及介质 | |
TWI532384B (zh) | 色彩調整裝置與色彩調整方法 | |
US8593682B2 (en) | Color conversion system and color conversion processing method | |
EP1483897A1 (en) | Choice of chromophores in two color imaging systems | |
JP2011228792A (ja) | 画像処理装置 | |
US9716811B2 (en) | Method and apparatus for obtaining cross-section of color gamut, and recording medium | |
JPH10200777A (ja) | 画像処理装置、画像処理方法および画像処理プログラムを記録した媒体 | |
CN105719244A (zh) | 图像增强处理方法及装置 | |
CN113287162B (zh) | 用于处理图像的方法及其存储介质、设备 | |
CN103079078A (zh) | 彩色图像的颜色数据管理方法 | |
KR100401132B1 (ko) | 모니터 알지비 대 프린터 알지비 매핑방법 | |
JP5020926B2 (ja) | 画像処理装置 | |
US9948831B2 (en) | In/out-of-gamut determination system for acquiring cross-section of color gamut of image forming apparatus taken on specific plane | |
US20070201072A1 (en) | Print data forming apparatus |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170208 Termination date: 20190731 |