高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

多基色显示系统基色亮度求解及讨论

姚昞晖 邓林宵 许立新

引用本文:
Citation:

多基色显示系统基色亮度求解及讨论

    通讯作者: 许立新, xulixin@ustc.edu.cn
  • 基金项目:

    国家重点研发计划资助项目 2021YFF0307804

  • 中图分类号: TN27;TN249

Calculation and discussion of the primary luminances in multi-primary display system

    Corresponding author: XU Lixin, xulixin@ustc.edu.cn
  • CLC number: TN27;TN249

  • 摘要: 为了进一步提升显示系统的颜色显示能力, 并解决三基色显示系统理论对多基色显示系统不兼容的问题, 采用参数不等式的方法进行了理论分析, 基于多基色显示系统和三基色显示系统的差异, 求解出多基色显示系统的峰值亮度所对应的解空间, 从理论上推导了多基色显示系统中基色峰值亮度的设置, 并对解空间的等价性和完备性进行了模拟实验验证。结果表明, 四基色、五基色和六基色激光显示系统对应的理论最大色域体积分别是2185100, 2258400和2395800。该方案对多基色显示系统的理论模拟具有重要的指导意义。
  • 图 1  R2G1B1构成的理想三基色显示系统的色域

    Figure 1.  Gamut of the ideal three-primary color display composed of R2, G1, B1

    图 2  a—解空间及对应的k1的取值范围和色域体积  b—当k1=0.34时的色域

    Figure 2.  a—solution space and the range of k1 and the corresponding values of color gamut volume  b—the gamut when k1= 0.34

    图 3  a—解空间及对应的k1k2的取值范围  b—当k1=0.34,k2=0.15时的色域

    Figure 3.  a—solution space and the range of k1 and k2  b—the gamut when k1=0.34 and k2=0.15

    图 4  a—解空间及对应的k1k2k3的取值范围  b—当k1=0.34,k2=0.15,k3=0.01时的色域

    Figure 4.  a—solution space and the range of k1, k2 and k3  b—the gamut when k1=0.34, k2=0.15 and k3=0.01

    表 1  模拟中采用的基色色坐标

    Table 1.  Chromatic coordinates in the simulation

    R1 R2 G1 G2 B1 B2
    coordinates
    3P (0.7080, 0.2920) (0.1700, 0.7970) (0.1310, 0.0460)
    4P (0.7080, 0.2920) (0.1700, 0.7970) (0.0388, 0.8116) (0.1310, 0.0460)
    5P (0.7260, 0.2740) (0.7080, 0.2920) (0.1700, 0.7970) (0.0388, 0.8116) (0.1310, 0.0460)
    6P (0.7260, 0.2740) (0.7080, 0.2920) (0.1700, 0.7970) (0.0388, 0.8116) (0.1310, 0.0460) (0.1611, 0.0138)
    下载: 导出CSV
  • [1] 许祖彦. 大色域显示——新一代显示技术[J]. 物理, 2010, 39(4): 227-231.

    XU Z Y. Large colour gamut display——the new generation of display technique[J]. Physics, 2010, 39(4): 227-231(in Chinese).
    [2] 刘旭. 现代投影显示技术[M]. 杭州: 浙江大学出版社, 2009: 46-57.

    LIU X. Modern projection display technology[M]. Hangzhou: Zhejiang University Press, 2009: 46-57(in Chinese).
    [3] 李亨. 颜色技术原理及其应用[M]. 北京: 科学出版社, 1994: 1-30.

    LI H. Principle of color technology and its application[M]. Beijing: Science Press, 1994: 1-30(in Chinese).
    [4] 郁道银. 工程光学[M]. 第4版. 北京: 机械工业出版社, 2018: 77-106.

    YU D Y. Engineering optics[M]. 4th ed. Beijing: China Machine Press, 2018: 77-106(in Chinese).
    [5] 胡威捷. 现代颜色技术原理及应用[M]. 北京: 北京理工大学出版社, 2007: 31-97.

    HU W J. Modern color science and application[M]. Beijing: Beijing Institute of Technology Press, 2007: 31-97(in Chinese).
    [6]

    FAIRCHILD M D. Color appearance models[M]. Hoboken, USA: John Wiley & Sons, 2013: 53-82.
    [7]

    MAcADAM D L. The theory of the maximum visual efficiency of colored materials[J]. Journal of the Optical Society of America, 1935, 25(8): 249-252. doi: 10.1364/JOSA.25.000249
    [8]

    MAcADAM D L. Maximum visual efficiency of colored materials[J]. Journal of the Optical Society of America, 1935, 25(11): 361-367. doi: 10.1364/JOSA.25.000361
    [9]

    MASAOKA K. Fast and accurate model for optimal color computation[J]. Optics Letters, 2010, 35(12): 2031-2033. doi: 10.1364/OL.35.002031
    [10]

    WANG G, YANG Y H, DONG T H, et al. Macadam's theory in RGB laser display[J]. Chinese Physics, 2016, B28(6): 064209.
    [11] 王聪, 姚昞晖, 马佳飞, 等. 三基色激光显示中白平衡点的选择研究[J]. 激光技术, 2021, 45(4): 411-416.

    WANG C, YAO B H, MA J F, et al. Study on the selection of white equilibrium point in three primary color laser display[J]. Laser Technology, 2021, 45(4): 411-416(in Chinese). 
    [12]

    INTERNATIONAL TECOMMUNICATION UNION-RADIOCOMMUNICATION SECTOR. Parameter values for the HDTV standards for production and international programme exchange[EB/OL]. [2022-11-18]. https://eva.fing.edu.uy/pluginfile.php/207457/mod_resource/content/1/R-REC-BT.709-5.pdf.
    [13]

    INTERNATIONAL TECOMMUNICATION UNION-RADIOCOMMUNICATION SECTOR. Parameter values for ultra-high definition television systems for production and international programme exchange[EB/OL]. [2022-11-18]. https://glenwing.github.io/docs/ITU-R-BT.2020-1.pdf.
    [14] 王贯. 新型显示的立体色域理论[D]. 合肥: 中国科学技术大学, 2018: 47-71.

    WANG G. Stereoscopic color gamut theory of new generation display systems[D]. Hefei: University of Science and Technology of China, 2018: 47-71(in Chinese).
    [15]

    DENG L X, DONG T H, FANG Y W, et al. Speckle reduction in laser projection based on a rotating ball lens[J]. Optics & Laser Technology, 2021, 135: 106686.
    [16]

    SONG H Y, LI H, LIU X. Studies on different primaries for a nearly-ultimate gamut in a laser display[J]. Optics Express, 2018, 26(18): 23436-23448. doi: 10.1364/OE.26.023436
    [17]

    WEN S. Design of relative primary luminances for four-primary displays[J]. Displays, 2005, 26(4/5): 171-176.
    [18]

    ZHU L Q, WANG G, YANG Y H, et al. Six-primary-laser projection display system: Demonstration and stereo color gamut measurement[J]. Optics Express, 2021, 29(26): 43885-43898. doi: 10.1364/OE.446944
    [19]

    YAO B H, ZHU L Q, YANG Y H, et al. General solution to the calculation of peak luminance of primaries in multi-primary display systems[J]. Optics Express, 2021, 30(2): 1036-1055.
    [20]

    MASAOKA K, JIANG F, FAIRCHILD M D, et al. Color gamut of multi-chromatic displays[J]. SID Symposium Digest of Technical Papers, 2019, 50(1): 861-864. doi: 10.1002/sdtp.13058
  • [1] 邓林宵杨雨桦姚昞晖朱立全王贯顾春许立新 . 多基色激光显示系统彩色散斑分布特性的研究. 激光技术, 2023, 47(6): 736-741. doi: 10.7510/jgjs.issn.1001-3806.2023.06.002
    [2] 王聪姚昞晖马佳飞方昱玮顾春许立新王贯 . 三基色激光显示中白平衡点的选择研究. 激光技术, 2021, 45(4): 411-416. doi: 10.7510/jgjs.issn.1001-3806.2021.04.002
    [3] 朱立全姚昞晖邓林宵杨雨桦王贯顾春许立新 . RGGB四基色激光显示系统的研究. 激光技术, 2023, 47(6): 772-777. doi: 10.7510/jgjs.issn.1001-3806.2023.06.007
    [4] 胡菊菊马军山高小燕 . 耦合半导体激光器的多调制延时混沌同步. 激光技术, 2010, 34(4): 440-442,455. doi: 10.3969/j.issn.1001-3806.2010.04.003
    [5] 高胜淼闫珂柱韩培高许春玉王荣新 . 飞秒激光诱导硅材料表面周期结构的研究. 激光技术, 2015, 39(3): 395-398. doi: 10.7510/jgjs.issn.1001-3806.2015.03.025
    [6] 李蕾臧景峰 . 双狭缝扫描法测量激光光束质量. 激光技术, 2015, 39(6): 845-849. doi: 10.7510/jgjs.issn.1001-3806.2015.06.024
    [7] 郑晖林季鹏史斐戴殊韬江雄康治军翁文林文雄 . 倍频过程对激光光束质量及空间分布的影响. 激光技术, 2009, 33(1): 67-70.
    [8] 骆永全张大勇张翠娟罗飞沈志学刘海涛 . 液晶光学器件激光损伤研究. 激光技术, 2010, 34(3): 392-394. doi: 10.3969/j.issn.1001-3806.2010.03.030
    [9] 刘伟孟冬冬范彩连李品蔡野伍彦伟刘会霞 . 激光透射焊接尼龙66光学属性研究. 激光技术, 2016, 40(5): 716-721. doi: 10.7510/jgjs.issn.1001-3806.2016.05.020
    [10] 张靳黄磊王东生殷聪巩马理 . 光学组合半导体激光器输出光束特性研究. 激光技术, 2007, 31(3): 228-231,241.
    [11] 王天明李斌成赵斌兴孙启明 . 高功率激光作用下光学元件非线性热效应研究. 激光技术, 2022, 46(6): 729-735. doi: 10.7510/jgjs.issn.1001-3806.2022.06.003
    [12] 徐成伟姚梅许振领张文攀胡欣王军刘艳芳 . 光学暗室内杂散激光能量分布模型仿真与测量. 激光技术, 2012, 36(1): 141-144. doi: 10.3969/j.issn.1001-3806.2012.01.037
    [13] 王燕陈家璧庄松林 . 激光多普勒效应的实验验证. 激光技术, 2010, 34(2): 170-172. doi: 10.3969/j.issn.1001-3806.2010.02.008
    [14] 刘航李义冯立强 . 红外激光与远紫外激光场驱动H2+辐射谐波. 激光技术, 2018, 42(2): 145-150. doi: 10.7510/jgjs.issn.1001-3806.2018.02.001
    [15] 何宁谭智诚吕杏利蒋红艳 . 基于语音识别控制的激光超声水下遥感研究. 激光技术, 2017, 41(3): 391-395. doi: 10.7510/jgjs.issn.1001-3806.2017.03.017
    [16] 王晓琳贺锋涛贾琼瑶刘佳 . 基于光纤振动的激光散斑控制. 激光技术, 2014, 38(2): 177-180. doi: 10.7510/jgjs.issn.1001-3806.2014.02.007
    [17] 郭商勇胡雄闫召爱程永强郭文杰 . 国外星载激光雷达研究进展. 激光技术, 2016, 40(5): 772-778. doi: 10.7510/jgjs.issn.1001-3806.2016.05.032
    [18] 刘航李义姚震冯立强 . 啁啾激光调控谐波截止能量及强度的研究. 激光技术, 2017, 41(5): 708-711. doi: 10.7510/jgjs.issn.1001-3806.2017.05.018
    [19] 和一凡申雨婷王文霄田友伟 . 激光脉冲初始相位对电子辐射的影响. 激光技术, 2023, 47(1): 103-107. doi: 10.7510/jgjs.issn.1001-3806.2023.01.016
    [20] 江天程湘爱 . 激光辐照3CCD规则分布亮点的研究. 激光技术, 2010, 34(2): 168-169,213. doi: 10.3969/j.issn.1001-3806.2010.02.007
  • 加载中
图(4) / 表(1)
计量
  • 文章访问数:  2145
  • HTML全文浏览量:  1272
  • PDF下载量:  20
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-11-21
  • 录用日期:  2022-12-07
  • 刊出日期:  2024-01-25

多基色显示系统基色亮度求解及讨论

    通讯作者: 许立新, xulixin@ustc.edu.cn
  • 1. 中国科学技术大学 核探测与核电子学国家重点实验室, 合肥 230026, 中国
  • 2. 中国科学技术大学 物理学院 安徽省光电子科学与技术重点实验室, 合肥 230026, 中国
  • 3. 先进激光技术安徽省实验室, 合肥 230026, 中国
基金项目:  国家重点研发计划资助项目 2021YFF0307804

摘要: 为了进一步提升显示系统的颜色显示能力, 并解决三基色显示系统理论对多基色显示系统不兼容的问题, 采用参数不等式的方法进行了理论分析, 基于多基色显示系统和三基色显示系统的差异, 求解出多基色显示系统的峰值亮度所对应的解空间, 从理论上推导了多基色显示系统中基色峰值亮度的设置, 并对解空间的等价性和完备性进行了模拟实验验证。结果表明, 四基色、五基色和六基色激光显示系统对应的理论最大色域体积分别是2185100, 2258400和2395800。该方案对多基色显示系统的理论模拟具有重要的指导意义。

English Abstract

    • 从早期的黑白显示到现在蓬勃发展的激光显示[1],显示一直在人们的生活中扮演着重要的角色,人们对它的要求也越来越高。显示色彩的优劣直接影响显示的效果。对于显示系统,通常使用色域来衡量其显示颜色的能力[2]。色域越大,则显示系统可渲染的颜色越丰富。为了定量表征颜色,国际照明委员会(Commission Internationale de I’Eclairage, CIE)首先提出了CIEXYZ颜色空间,并将其简化为xy色度图[3],其直观的特点让其被广泛应用在显示领域来表征色域。但颜色是一个3维的物理量,2维的xy色品图为了直观而牺牲了亮度属性,这也导致其准确性相对较低[4]。为了更准确地描述人眼感知的颜色,CIE进一步提出了CIEL*a*b*颜色空间,它是目前应用最广泛的均匀颜色空间[5]。立体颜色空间中的色域称为立体色域。在均匀颜色空间中,精确定量描述颜色成为可能[6]。MAcADAM最早于1935年提出了计算色域的方法[7],但由于时代限制,只将色域在xy色品图中进行了表示[8]。MASAOKA改进了他的方法, 并将其应用于立体颜色空间,获得了立体色域[9]。而WANG等人继续改进了MAcADAM的方法,使其适用于当前常见的基色显示系统[10]。WANG等人进一步研究了色温或者说白平衡点色坐标对色域的影响[11]。根据这些研究可以看出,目前的三基色显示系统的理论基础已经较为完备。在此基础上,为了方便不同设备之间色彩还原的准确性,美国国家电视标准委员会(National Television Standards Committee,NTSC)与电影和电视工程师协会等提出了NTSC Rec.709[12]、Rec.2020[13]等色域标准,其对应的色域大小也不尽相同,其中Rec.2020是当今最先进的色域标准。《中国制造2025》中明确提出,2020年激光显示的色域要达到160% NTSC标准所对应的色域,2025年要进一步达到200% NTSC标准所对应的色域,这说明色域有重要意义且具有很高的研究价值。

      随着激光光源的发展,三基色激光显示系统已经达到149% NTSC的平面色域面积[14], 这说明减小基色谱宽是提升色域的重要方法,但基色谱宽过窄会导致激光散斑增大的问题[15],并且三基色显示系统的上限较为明显[16]。为了进一步提升色域,势必要采用多基色显示的方案[17]。随着多基色显示技术的发展,目前已经有相关的研究,成功搭建了六基色显示系统,并使其立体色域的体积超过了180% NTSC[18]。在作者团队前期的研究中,提出了多基色显示系统基色亮度通解的计算方案,利用格拉斯曼定律来获得多基色显示系统对应的亮度解空间,以得到其所有可能的基色搭配,并采用实际的六基色显示系统进行了验证[19],但这些探究尚不完善。本文作者将对之前提出的解空间理论计算的数学形式进行规范化,并对其数学形式上的完备性和等价性进行证明。进一步地,将探究了解空间的几何性质,并采用理想的多基色显示系统进行演示。为简化起见,本文中的讨论仅限于颜色可加性良好的显示系统,为了探究显示系统的基本特性,本文中不考虑不同环境光源对颜色外观带来的颜色变化,即认为显示系统位于一个理想暗室中。

    • 对于显示系统而言,存在一个最基本的约束条件,即所有基色的最高亮度之和为白色,该白色也称为该系统的白平衡点。为了更好地说明多基色显示系统亮度通解的求解,此处对作者团队之前的工作进行了数学形式的优化。根据格拉斯曼定律,一个颜色可以由3种在颜色空间中不共线的颜色混合而得到,利用公式表达,即为:

      $ \left\{\begin{array}{l} X_1+X_2+X_3=X_{\text {mix }} \\ Y_1+Y_2+Y_3=Y_{\text {mix }} \\ Z_1+Z_2+Z_3=Z_{\text {mix }} \end{array}\right. $

      (1)

      $ \left\{\begin{array}{l} \frac{X_1}{x_1}+\frac{X_2}{x_2}+\frac{X_3}{x_3}=\frac{X_{\text {mix }}}{x_{\text {mix }}} \\ \frac{Y_1}{y_1}+\frac{Y_2}{y_2}+\frac{Y_3}{y_3}=\frac{Y_{\text {mix }}}{y_{\text {mix }}} \\ \frac{Z_1}{z_1}+\frac{Z_2}{z_2}+\frac{Z_3}{z_3}=\frac{Z_{\text {mix }}}{z_{\text {mix }}} \end{array}\right. $

      (2)

      式中:XYZ及其下标1, 2, 3分别代表三基色的三刺激值;xyz代表对应的色坐标, 下文同理; 下标mix代表混合后的颜色。如果令混合后的颜色为白色,则最终得到的计算结果就是常见的三基色显示系统。在已知基色和白色的色坐标的前提下,上述6个方程可以简化成含有3个未知数的三方程组,因此对应特定的基色和白色,三基色显示系统只有一组特定的Y刺激值的解。需要说明的是,Y刺激值和亮度存在对应的转化关系,因此获得了Y刺激值即相当于获得了基色的亮度比例。而对于N(N>3)基色显示系统,令j=N-3,公式则变为:

      $ \left\{\begin{array}{l} X_1+X_2+X_3+\sum\limits_j X_{j+3}=X_{w h i t e} \\ Y_1+Y_2+Y_3+\sum\limits_j Y_{j+3}=Y_{\text {white }} \\ Z_1+Z_2+Z_3+\sum\limits_j Z_{j+3}=Z_{w h i t e} \end{array}\right. $

      (3)

      $ \left\{\begin{array}{l} \frac{X_1}{x_1}+\frac{X_2}{x_2}+\frac{X_3}{x_3}+\sum\limits_j \frac{X_{j+3}}{x_{j+3}}=\frac{X_{\text {white }}}{x_{\text {white }}} \\ \frac{Y_1}{y_1}+\frac{Y_2}{y_2}+\frac{Y_3}{y_3}+\sum\limits_j \frac{Y_{j+3}}{y_{j+3}}=\frac{Y_{\text {white }}}{y_{\text {white }}} \\ \frac{Z_1}{z_1}+\frac{Z_2}{z_2}+\frac{Z_3}{z_3}+\sum\limits_j \frac{Z_{j+3}}{z_{j+3}}=\frac{Z_{\text {white }}}{z_{\text {white }}} \end{array}\right. $

      (4)

      式中:下标white代表白色。显然,对于多基色系统,上述方程由于增加了新的未知数,因此变为不定解方程,不存在如同三基色系统的一组定解来获得白平衡,而是有无数种解。这对控制多基色显示系统的白平衡,以及优化系统的基色参数提出了挑战。

      对于一个N(N>3)基色显示系统,假设已知所有基色的色坐标和白色的色坐标。首先任意选取其中3个基色,得到如下刺激值条件限制:

      $ \left\{\begin{array}{l} X_{\mathrm{o}, 1}+X_{\mathrm{o}, 2}+X_{\mathrm{o}, 3}=X_{\text {white }} \\ Y_{\mathrm{o}, 1}+Y_{\mathrm{o}, 2}+Y_{\mathrm{o}, 3}=Y_{\text {white }} \\ Z_{\mathrm{o}, 1}+Z_{\mathrm{o}, 2}+Z_{\mathrm{o}, 3}=Z_{\text {white }} \end{array}\right. $

      (5)

      $ \left\{\begin{array}{l} \frac{X_{\mathrm{o}, 1}}{x_1}+\frac{X_{\mathrm{o}, 2}}{x_2}+\frac{X_{\mathrm{o}, 3}}{x_3}=\frac{X_{\text {white }}}{x_{\text {white }}} \\ \frac{Y_{\mathrm{o}, 1}}{y_1}+\frac{Y_{\mathrm{o}, 2}}{y_2}+\frac{Y_{\mathrm{o}, 3}}{y_3}=\frac{Y_{\text {white }}}{y_{\text {white }}} \\ \frac{Z_{\mathrm{o}, 1}}{z_1}+\frac{Z_{\mathrm{o}, 2}}{z_2}+\frac{Z_{\mathrm{o}, 3}}{z_3}=\frac{Z_{\text {white }}}{z_{\text {white }}} \end{array}\right. $

      (6)

      式中:下标o代表选取的初始三基色混合成白色的情形。计算得到Yo, 1Yo, 2Yo, 3,该结果代表选取的初始三基色混合成白色所需要的Y刺激值。在理论计算中通常令Ywhite=100来简化计算, 再将第j+3 (j≥1)个基色用前3个基色进行表示,有:

      $ \left\{\begin{array}{l} X_{j, 1}+X_{j, 2}+X_{j, 3}=X_{j+3}{ }^{\prime} \\ Y_{j, 1}+Y_{j, 2}+Y_{j, 3}=Y_{j+3}{ }^{\prime} \\ Z_{j, 1}+Z_{j, 2}+Z_{j, 3}=Z_{j+3}{ }^{\prime} \end{array}\right. $

      (7)

      $ \left\{\begin{array}{l} \frac{X_{j, 1}}{x_1}+\frac{X_{j, 2}}{x_2}+\frac{X_{j, 3}}{x_3}=\frac{X_{j+3}^{\prime}}{x_{j+3}} \\ \frac{Y_{j, 1}}{y_1}+\frac{Y_{j, 2}}{y_2}+\frac{Y_{j, 3}}{y_3}=\frac{Y_{j+3}^{\prime}}{y_{j+3}} \\ \frac{Z_{j, 1}}{z_1}+\frac{Z_{j, 2}}{z_2}+\frac{Z_{j, 3}}{z_3}=\frac{Z_{j+3}^{\prime}}{z_{j+3}} \end{array}\right. $

      (8)

      式中:下标j代表初始三基色混合成第j+3个基色的情形;Xj+3′,Yj+3′,Zj+3′代表此时第j+3个基色的三刺激值;xj+3yj+3zj+3代表其对应色坐标,下文同理。假设第j+3个基色的Y刺激值为kjYj+3′,引入kj参数作为变量,并令Yj+3′为常数方便计算,用这6个式子可以计算出Yj, 1Yj, 2Yj, 3。由于多基色系统的定义是每个基色都不能由其它基色混合得到,因此此时Yj, 1Yj, 2Yj, 3必然有负值,但并不影响后续计算, 则此时多基色系统满足:

      $ \begin{gathered} Y_{\mathrm{o}, 1}+Y_{\mathrm{o}, 2}+Y_{\mathrm{o}, 3}+\sum\limits_{j=1}^{N-3} k_j\left[Y_{j+3}^{\prime}-\right. \\ \left.\left(Y_{j, 1}+Y_{j, 2}+Y_{j, 3}\right)\right]=Y_{\text {white }} \end{gathered} $

      (9)

      即:

      $ \begin{gathered} \left(Y_{\mathrm{o}, 1}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 1}\right)+\left(Y_{\mathrm{o}, 2}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 2}\right)+ \\ \left(Y_{\mathrm{o}, 3}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 2}\right)+\sum\limits_{j=1}^{N-3} k_j Y_{j+3}^{\prime}=Y_{\text {white }} \end{gathered} $

      (10)

      上述两个式子的物理意义是: 为了增加新基色的Y刺激值,则必须按照比例降低原有三基色的Y刺激值; 不难看出,式中每一项即对应每一个基色的Y刺激值。对于多基色显示系统,每个基色的Y刺激值都应该大于0,因此有:

      $ \left\{\begin{array}{l} Y_{\mathrm{o}, 1}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 1}>0 \\ Y_{\mathrm{o}, 2}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 2}>0 \\ Y_{\mathrm{o}, 3}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 3}>0 \\ k_j Y_{j+3}^{\prime}>0 \\ k_{j+1} Y_{j+4}^{\prime}>0 \\ \ldots \end{array}\right. $

      (11)

      此不等式组即为N(N>3)基色系统的解空间,根据不等式组的性质,不难得出解空间的维度为N-3, 而此时解空间中的点(k1, k2, …, kj, …)对应的多基色系统的基色亮度解为:

      $ \left\{\begin{array}{l} Y_1=Y_{\mathrm{o}, 1}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 1} \\ Y_2=Y_{\mathrm{o}, 2}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 2} \\ Y_3=Y_{\mathrm{o}, 3}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 3} \\ Y_{j+3}=k_j Y_{j+3}{ }^{\prime} \\ Y_{j+4}=k_{j+1} Y_{j+4}^{\prime} \\ \cdots \end{array}\right. $

      (12)

      式中:Yj+3′都为常数,为了计算方便,通常设置为100。获得了每个基色的Y刺激值后,根据每个基色的色坐标(xi, yi)(其中i=1, 2, …, N)求出所有基色的三刺激值,将显示系统的亮度按照Y刺激值的比例进行设置,即可得到一个多基色显示系统。

      在获得了三刺激值后,在立体颜色空间中,根据参考文献[20]中的计算方法,即得到立体色域,代表此多基色显示系统可以显示颜色的集合,目前最常用的立体颜色空间为CIEL*a*b*颜色空间。对于一组特定的基色和白平衡点参数,比较解空间中不同点对应的立体色域,即可获得某组基色的最优亮度组合;而对比不同的基色和白平衡点参数组合的解空间,即可获得最优的基色组合。

    • 首先,只有获得多基色显示系统的所有可能的参数组合前提下,才能进行最优化。因此有必要证明上述推导中的解空间是完备的,即该解空间包含了该特定基色和白色参数组合下所有的Y刺激值的可能搭配。

      此处证明采用反证法。假设在某基色和白色参数组合下,存在一组Y刺激值(Y1′, Y2′, …, Yj′, …)在解空间中不存在点(k1, k2, …, kj, …)和其对应。将(Y1′, Y2′, …, Yj′, …)代入式(12),由于式(12)的未知数为Yo, 1′,Yo, 2′,Yo, 3′和kj′,未知数数目和方程数一致,且方程互相线性无关,故该方程可必然反解出唯一解(k1′, k2′, …, kj′, …),此点不在解空间内, 故该点必然不满足式(11),即Yj′中存在负值。因为Y刺激值和亮度可按照比例对应,由于负亮度不存在,故不存在这样一组Y刺激值不在式(11)和式(12)的包含范围内。因此解空间是完备的。

      然后需要证明解空间的等价性。从之前的推导中可以看出,解空间的推导基于最先选择的三基色的参数。不同的初始三基色计算出来的解空间的范围显然不同。但是,解空间中的点对应的基色Y刺激值是等价的,即一个初始三基色计算出的解空间中的点(k1, k2, …, kj, …)对应的一组基色Y刺激值,必然可以在另一组初始三基色计算出的解空间找到一个对应点(k1*, k2*, …, kj*, …)。

      对于解空间1,其中的坐标点(k1, k2, …, kj, …)对应的基色Y刺激值如下式所示:

      $ \left\{\begin{array}{l} Y_1=Y_{\mathrm{o}, 1}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 1} \\ Y_2=Y_{\mathrm{o}, 2}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 2} \\ Y_3=Y_{\mathrm{o}, 3}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 3} \\ Y_{j+3}=k_j Y_{j+3}^{\prime} \\ Y_{j+4}=k_{j+1} Y_{j+4}^{\prime} \\ \cdots \end{array}\right. $

      (13)

      其对应的解空间为:

      $ \left\{\begin{array}{l} Y_{\mathrm{o}, 1}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 1}>0 \\ Y_{\mathrm{o}, 2}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 2}>0 \\ Y_{\mathrm{o}, 3}-\sum\limits_{j=1}^{N-3} k_j Y_{j, 3}>0 \\ k_j Y_{j+3}{ }^{\prime}>0 \\ k_{j+1} Y_{j+4}{ }^{\prime}>0 \\ \ldots \end{array}\right. $

      (14)

      对于解空间2,初始三基色更换后,其解空间和基色Y刺激值如下所示:

      $ \left\{\begin{array}{l} Y_{\mathrm{o}, 1}{ }^*-\sum\limits_{j=1}^{N-3} k_j^* Y_{j, 1}{ }^*>0 \\ Y_{\mathrm{o}, 2}{ }^*-\sum\limits_{j=1}^{N-3} k_j^* Y_{j, 2}{ }^*>0 \\ Y_{\mathrm{o}, 3}{ }^*-\sum\limits_{j=1}^{N-3} k_j^* Y_{j, 3}{ }^*>0 \\ k_j{ }^* Y_{j+3}{ }^{\prime \prime}>0 \\ k_{j+1}{ }^* Y_{j+4}{ }^{\prime \prime}>0 \\ \cdots \end{array}\right. $

      (15)

      $ \left\{\begin{array}{l} Y_1{ }^*=Y_{\mathrm{o}, 1}{ }^*-\sum\limits_{j=1}^{N-3} k_j{ }^* Y_{j, 1}{ }^* \\ Y_2{ }^*=Y_{\mathrm{o}, 2}{ }^*-\sum\limits_{j=1}^{N-3} k_j{ }^* Y_{j, 2}{ }^* \\ Y_3{ }^*=Y_{\mathrm{o}, 3}{ }^*-\sum\limits_{j=1}^{N-3} k_j{ }^* Y_{j, 3}{ }^* \\ Y_{j+3}{ }^*=k_j{ }^* Y_{j+3}{ }^{\prime \prime} \\ Y_{j+4}{ }^*=k_{j+1}{ }^* Y_{j+4}{ }^{\prime \prime} \\ \cdots \end{array}\right. $

      (16)

      注意此时的两套符号中,虽然仍为原基色组,但编号顺序发生了变化,这是由于选取的初始三基色不同导致的,也即式(13)和式(14)中的Yi与式(15)及式(16)中的Yi*(其中i=1, 2, …, N)代表不同的基色,相当于交换了顺序,更换了初始三基色,原先用来表征其它基色的颜色变成了被表征颜色, 而Yj+3′和Yj+3对应,Yj+3″和Yj+3*对应,形式的含义和式(11)、式(12)表征的解空间的基本形式一致。根据解空间的完备性,即式(13)和式(15)对应的解空间都是完备的,则解空间1中的基色Y刺激值组,必然满足式(15),因此根据式(16)必然可以反解出解空间2中的唯一点(k1*, k2*, …, kj*, …)。因此不同初始三基色的选择虽然会导致解空间的差异,但是由解空间中的点导出的基色Y刺激值组,是完备并且一致的。由此可以得出结论: 该方案对于多基色系统的基色亮度计算是成功的。

    • 为了进一步说明本文中的算法,现采用模拟样例进行说明。表 1中给出了模拟样例的色坐标,这里从三基色激光显示系统出发,添加新基色扩展为四基色、五基色和六基色系统,观察其解空间的演化。最初的三基色显示系统,选择其色坐标匹配Rec.2020色域标准。白色的色坐标(xwhite, ywhite)均设置为(0.3127, 0.3290),其Y刺激值设置为100。表中3P,4P,5P,6P分别代表三基色、四基色、五基色、六基色激光显示系统,R1R2G1G2B1B2分别代表不同波长的两种红基色、绿基色和蓝基色。

      表 1  模拟中采用的基色色坐标

      Table 1.  Chromatic coordinates in the simulation

      R1 R2 G1 G2 B1 B2
      coordinates
      3P (0.7080, 0.2920) (0.1700, 0.7970) (0.1310, 0.0460)
      4P (0.7080, 0.2920) (0.1700, 0.7970) (0.0388, 0.8116) (0.1310, 0.0460)
      5P (0.7260, 0.2740) (0.7080, 0.2920) (0.1700, 0.7970) (0.0388, 0.8116) (0.1310, 0.0460)
      6P (0.7260, 0.2740) (0.7080, 0.2920) (0.1700, 0.7970) (0.0388, 0.8116) (0.1310, 0.0460) (0.1611, 0.0138)

      对于三基色系统,根据式(5)和式(6)可计算出:

      $ \left\{\begin{array}{l} Y_{\mathrm{o}, R_2}=26.2700 \\ Y_{\mathrm{o}, G_1}=67.7998 \\ Y_{\mathrm{o}, B_1}=5.9302 \end{array}\right. $

      (17)

      根据MASAOKA提出的方案[20]来计算立体色域,在CIEL*a*b*空间中画出其立体色域,如图 1所示。图中, L*为亮度,a*轴正负分别代表红绿色,b*轴正负分别代表黄蓝色, 其立体色域的体积为1854900。

      图  1  R2G1B1构成的理想三基色显示系统的色域

      Figure 1.  Gamut of the ideal three-primary color display composed of R2, G1, B1

      利用式(1)和式(2),将第4个基色G2用初始的三基色R2, G1, B1进行表示,有:

      $ \left\{\begin{array}{l} Y_{1, R_2}=-8.4144 \\ Y_{1, G_1}=107.6332 \\ Y_{1, B_1}=0.7812 \end{array}\right. $

      (18)

      代入式(11),得到四基色显示系统的解空间:

      $ 0<k_1<0.6299 $

      (19)

      图 2a中给出了四基色显示系统的色域体积随k1变化的趋势。利用MASAOKA提出的多基色立体色域计算方法[20]对色域体积进行计算,可以看出,在k1=0.34时,色域体积取最大值2185100,此时色域如图 2b所示。

      图  2  a—解空间及对应的k1的取值范围和色域体积  b—当k1=0.34时的色域

      Figure 2.  a—solution space and the range of k1 and the corresponding values of color gamut volume  b—the gamut when k1= 0.34

      再添加基色R1,同样按照上述步骤,有:

      $ \left\{\begin{array}{l} Y_{2, R_2}=110.1457 \\ Y_{2, G_1}=-10.1692 \\ Y_{2, B_1}=0.0235 \end{array}\right. $

      (20)

      图 3a为五基色显示系统的解空间,在一个平面范围; 图 3bk1=0.34和k2=0.15时的立体色域,色域体积为2258400。

      图  3  a—解空间及对应的k1k2的取值范围  b—当k1=0.34,k2=0.15时的色域

      Figure 3.  a—solution space and the range of k1 and k2  b—the gamut when k1=0.34 and k2=0.15

      对于六基色显示系统,添加基色B2,同样的,有:

      $ \left\{\begin{array}{l} Y_{3, R_2}=119.1512 \\ Y_{3, G_1}=-354.1547 \\ Y_{3, B_1}=335.0035 \end{array}\right. $

      (21)

      其解空间如图 4所示。当k1=0.34、k2=0.15、k3= 0.01时,立体色域如图 4b所示,色域体积为2395800。

      图  4  a—解空间及对应的k1k2k3的取值范围  b—当k1=0.34,k2=0.15,k3=0.01时的色域

      Figure 4.  a—solution space and the range of k1, k2 and k3  b—the gamut when k1=0.34, k2=0.15 and k3=0.01

    • 根据图例中的解空间可以看出,每增加一个基色,解空间就增加一个维度。本文中的解空间获取,是在基色亮度无限制的情况下进行计算的。在实际应用中,考虑到实际光源的亮度和成本限制,必然存在更多的约束条件,比如要求某种基色的亮度低于某个值。在这种情况下,就需要对计算添加更多的约束条件,即增加不等式的数目,但这并不改变解空间的维度,只是改变了解空间的形状。在约束后的解空间中,依旧可以采取同样的方式,来获得需要的基色亮度配比。从上述的计算样例中可以发现,四基色显示系统的解空间,为五基色显示系统的解空间的一条边界,而五基色显示系统的解空间,为六基色显示系统解空间的底面。因为解空间由式(11)构成,其边界应当为其中某些不等式等于0的情形,即对应维数降低的显示系统。因此N基色显示系统的解空间的边界实际上涵盖了(N-1),(N-2)···,一直到三基色,所有以这些参数作为基色的理论显示系统的解空间。本文作者提出的方法成功解决了多基色显示系统中基色亮度选择的问题。结合色域体积的计算,可以对多基色显示系统的光源参数进行适当的优化,以满足实际的显示需求。

    • 基于多基色显示系统和和三基色显示系统的差异,本文作者通过参数不等式的方式,在给定白平衡点和基色色坐标的情况下,理论上推导了多基色显示系统中基色亮度的设置问题。求解出多基色显示系统的亮度所对应的解空间,并对解空间的等价性和完备性进行了证明。同时,本文中采用理想的四基色、五基色和六基色显示系统对整体算法进行了演示,并探究了解空间的相关性质。结合立体色域的计算,通过遍历解空间,可以很容易地得到理论最大色域体积。该算法可以推广至任意种光源的多基色显示系统。结合色域体积的计算,该理论可用于指导基色波长、谱宽和亮度的选择。在未来的工作中,将继续探索基于该方法的颜色管理理论,提升色彩的还原程度和表现力。

参考文献 (20)

目录

    /

    返回文章
    返回