GB_T 43789.31-2024 电子纸显示器件 第3-1部分：光学性能测试方法.docx

ICS31.120CCS1.53GB中华人民共和国国家标准GB/T43789.312024/IEC62679-3-1：2014电子纸显示器件第3-1部分：光学性能测试方法Electronicpaperdisplays-PartS-IQpticalmeasuringmethods(IEC62679-3-1：2014,IDT)2024-07-01实施2024-03-15发布国家市场监督管理总局国家标准化管理委员会前言引言1范困2规范性引用文件32端珞语A标准刃I试冬件43标准照明条件44武I©十4546Wjyi式1义胴IK5士去ct¼5JQ视角范困5Il9射录A（密科性）日光色域的计U方法A1S66A.3计算QE1.AB立体色域的表面细分法I11111123911111212151718182226272831343636363843本文件按照GB/T1.12020标准化工作导则第1部分:标准化文件的结构和起草规则的规定起草.本文件是GB/T43789电子纸显示器件的第3-1部分.GB/T43789已经发布了以下部分：一第1-1部分:术谙；一第3T部分：光学性能测试方法；一第3-2部分：光电性能测试方法.本文件等同采用IEC62679-3-1:20144电子纸显示器件第31部分:光学性能测试方法).请注意本文件的某些内容可能涉及专利.本文件的发布机构不承担识别专利的责任.本文件由中华人民共和国工业和信息化部提出.本文件由全国电子显示器件标准化技术委员会(SAC/TC547)归.本文件起草单位：广州奥贯电子科技股份有限公司.中国电子技术标准化研究院、东南大学.上海唯视锐光电技术有限公司、浙江智爱科技有限公司.浙江智慧照明技术有限公司、京东方科技集团股份有限公司.Ta华星光电技术有限公司.颜色空间科技有限公司、厦门市计量检定测试院.广州赛西标准检源研究院有限公司、合阳维信诺科技有限公司、义乌涌越光电科技有限公司、四川赛狄佶息技术股份公司.本文件壬靛草人:王直杜.曲典禄、李晓茹、赵英、王飞S1.李晓华、王蔚生、李俊凯、牟同升、张志刚、黄卫东、旗听、阮育娇、黄艺宾、吴杜雄'勉痂、高裕弟、肖红.为了满足我国电子纸显示器件产品的生产制造、检验和出口货易需求，更好地促进相关领域的交流和技术合作,有必要制定电子纸显示器杵系列标准.G8/T43789电子纸显示据件拟由以下部分构成：一第1-1部分:术语.目的在于确定电子纸显示器件的基本术语.一第1-2部分:总规范.目的在于确定电子纸显示翳件检脆的总规范.一第2部分,基本额定值和特性.目的在于确定评估电子纸显示器件的光学性能.电光特性及环境测试的额定值和特性.一第3-1部分:光学性能测试方法.目的在于确定电子纸显示器件光学性能的标准测试条件和测试方法.一第3-2部分：光电性烧测试方法.目的在于确定电子纸显示提件光电性畿的标准测试条件和测试方法.一第33部分：带集成照明单元显示器件的光学测试方法.目的在于确定带集成照典单元的电子纸显示器件的光学测试方法.一第4-2部分:环境试版方法.目的在于确定电子纸显示器件的环境试脸方法.一第Sr部分：空间频率下电子堆显示的可读性.目的在于斓定电子纸空间频率下的可读性评估标准.电子纸显示器件第3-1部分：光学性能测试方法1必图本文件磁述了电子纸显示器件(EPD)光学性能的标注测试条件和测试方法.本文件适用于反射式单色或彩色的无源字段或有源矩阵型电子纸显示器件应用.对于包含集成式光源组件(I1.U)的电子纸显示器件，在测试时，集成式光源组件需关闭.本文件不适用于貂过三种颜色的彩色体系的产品.2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款.其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本(包括所有的修改单)适用于本文件.IEC62679-1-1电子堆显示器件第1-1S8：iS(ElectronicpaperdisplaysPartl-1.Termi-nology)注：GB/T43789.112024胆子母属示Zf件«11部分：术语UEC62679l1.2014.1DT)IEC60050(所有部分)国际电工词汇(InternationaIeIeCtrOteChniCalvocabulary)注：GB/T2900(所有部分)电工术的lEC60050(所有部分)IEC61966-2-1多媒体系统和设筋SS色测试和险理第2-1部分:颜色管理»；认RG8色彩空间sRGB(MultimediasystemsandequipmentColourmeasurementandmanagementPart2-1：CoIourmanagement-DefauItRGBcolourspace-sRGB)CIE15比色法(Colorimetry)CIE38材料的辐射和光度特性及只测试(Radiometricandphotometriccharacteristicsofmaterialsandtheirmeasurement)3术语和定义.缩珞语3.1 术语和定义IEC62679-1-1.IEC60050(所有部分)界定的以及下列术语和定义适用于本文杵.3.1.1环境光对比度ambientcontrastratio同时使用半球漫射光源和直射光源腌射显示屏去面枭模拟JS实照明环境时，显示屏的对比度.3.12日光下的显示屏颜色daylightdisplaycolour同时使用半球漫射光源和定向光源以特定几何条件.光谱和光强射于显示屏去面来模拟百实日光照明环境下的显示屏颜色.3.1.3立体色域colourgamutvolume显示颜色的最大可能范围(包括原色.白色W和黑色K的所有可能混合)，是在三维颜色空间如QE1.AB中的体积.3.1.4日光下的立体色域daylightcolourgamutvolume同时使用半球漫射光源和定向光源以特定几何条件、光谙和光强照射于显示屏表面来模拟真实日光照明环境下的立体色域.3.2 缩略语下列缩略语适用于本文件.CCT：相关色温(correlatedcolourtemperature)CIE：国际照明委员会(IntematiOnaIeommi$siononIllumination)C1E1.AB：C1E1976(ab')色空间IClE1976(1.a-b')colourspaceDUT:被测器件(deviceundertest)EPD：电子纸显示器件(electronicpaperdisplay)I1.U：集成光遁组件(如侧光照照的前导光板Mimegratedlightingunit(e.g.anedge-Iitfrontguidepate)JISO：国际标准化组织(IntematiOnalorganiZatiOnfOrStandardiZatiOn)1.ED：发光二极管(Iightemittingdiode)1.MD:光测试仪器(rightmeasuringdevice)RGB:红，绿,三(red,green.blue)SDCM:色容差(standarddeviationOfcolourmatching)SRGB:在IEC61966-2-1定义的标准RGB色空间(StandardRGBcolourSPaCeasdefinedinIEC61966-2-1)4标准测试条件4.1 标准测试环境条件光学和光电测试应在标准环境条件下进行：温度25°C±3cC,相对温度25%8S%,气压86kPa106kPa.当使用不同环境条件时，应在报告中注明.4.2 观察方向坐标系垓观察方向是指观测者看向祓测翳件(DUT)测试点的方向.在测试过程中，光测试设备(1.MD)模拟观测者，在观察方向上对准DUT上的测试点.观索方向由两个角度来表征：倾角6(相对于DUT的去面法线)和旋缝角¢(也称为方位角)，如图1所示.方位用是按逆时针测量方向，它与时怦表面的方向如下:¢=0°是3点钟方向(右)，¢=90°是12点钟方向(上)，¢=180°是9点钟方向(左)，¢=270°是6点钟方向(下).图1由极坐标系统中的倾角和旋转角（方位角）表示的观察方向或测试方向4.3 标准酣明条件4.4 .1EPD测试的一般要求本文件中的EPD是反射型显示翳件.反射型显示器件是一种利用反射光实现信息显示的显示器件.反射型显示器件并不发射任何光线，需有环境光才能看到显示信息，因此反射型显示器件的测试规范中应包括测试期间的照明条件.测试用照明由一个或多个光源组成，每一个光源的光谙分布和光源几何条件都应详细说明.因此,为使测试可再现，显示性能测试应在具体明确规定的照明和检测条件下进行.IlU集成到EPD中来提供补充照明，以补偿不够充足的环境照明.本文件中的测试方法是在I1.U关闭的情况下进行.4.3辎述了用于测试EPD光学性锢指标的标准照明条件.除了标准几何条件外，EPD还可在其他照明和抽测几何条件下进行测试.照明光源需要预热,在整个测试过程中，光源信号应稳定在±5%以内.4.3.2 哨室条件EPD的测试应在受控照明条件下进行.应把不需要的背景照明降到堆低,通常在胞室中照射显示器件.来自昭室的背景照明光谱辐射亮度，即从DUT反射测狎的光谱辐射亮度,应不大于打开照明光源时DUT簿色状态的光滑辐射光度的1/100.如果不能满足该条件，则应扣除肾寰照明在DUT反射的光滑辐射亮度且在报告中注照.另外,如果IMD的灵敏度不足以在这些较低光诺箱射亮度水平下进行测试,则应在掖告中注明IMD的下限.除非另有规定外,标注省累照明条件应为暗室条件.4.3.3 标准环境阳明光谱以下照照条件适用于在环境照明下，对反射式显示器件进行光学及光电测试.环境照明应模拟室内或室外照明条件.通常使用两种几何照明条件的光源姐合来模拟室内环境照明和请空下的室外日光醺明（见参考文献国、）.均匀的半球湿射醺明将用于模拟定向光漉非百射条件下的室内背景照明,例如室内的灯具，或模拟太阳光漫射照明在显示屏.定向光源将模拟室内灯具或阳光声射在显示屏上定向照明的效果.以下障明条件用于模拟在室内和空外观看显示屏的环境.a）室内照明条件均匀的半球漫射照明：使用近似CIE15中定义的QE标准照明体A、CIE标准照明体D6S或ClE标准照明体D50的光源.对于光谱测试，应使用光滑平滑的宽光谱光源（例如近似于ClE标准照明体A）.使用宽光谱光源（例如CIE标准照照体A）测试光滑反射因数，室内亮度和色度可用要求的参考光谱（例如QE标准照明体D65）计算得到.对于室内阅读环境，应使用300IX来计算性地指标（见参考文献3）.为使测试更精确，实际测试半球漫反射系数时可在较高的照度进行，然后根据实际光环境的照度值对结果按比例缩小.定向照明：应使用与半球漫射照明相同光谱的光源.对于室内阅读环境且显示屏处于垂直方向，在计算室内光度和色度时，应使用定向光源入射到显示屏去面的照度值200Ix.为使测试更精确，实际测试反射系数时可能需要更高的照度等级，然后根据实际光环境的照度值对结果按比例缩小,定向光通与表面法线的夹角应为4'（即ft-4S1且对向角不超过5,对向角被定义为颁到显示屏测试区域中G的全角度,图，；|算室内照明条件下的环墩光对比度,除以上规定的照度等级外，还可使用其他等级.但是，总照明大约60%是半球漫射照明，而40%是定向照明.b）日光照明条件均匀的半球濯射照明：使用与QE标准照明体D75的B光光诣分布近似的光源（见叁考文献4）.根据实际应用场合，还可使用其他的ClE日光照明体（例如D6S）.对于光谱测试，可使用光潜平滑的宽光谙光源（例如接近于QE标准瞟明体A）冽试光谙反射因数，之后可用ClE标准照明体D7S的光谱计算日光亮度和色度.计算日光亮度和色度时，应使用半球漫射照明（包括镜面反射）里出入射在显示解表面的照度值150001x（见参考文献的、5）,实际半球漫反射系数的测试可在较低的照度下进行，然后根据要求的照度佰对结果按比例放大.定向照明：定向光源应接近CIE日光照明体D50（见参考文献4）.根据实际应用场合,还可使用其他的ClE日光照明体（例如D65）可使用光造平滑的宽带光源（例如接近于QE标准照明体A）测试光诺反射因数，之后可用QE标准照明体DSO的光谙计算日光亮度和色度.在计算日光对比度或色度时，应使用定向光源以倾角&=45。入射到显示屏表面的照度值6S000lx,且IMD应垂直于显示屏表面法线（Cd=。"）（见参考文献（4）、（5）.实际的反射系数测试可在较低的照度等级下进行，然后将结果按比例放大到所需的照度等级.对比度和色度是根据按比例放大的照度等级计算的.定向光源的对向角应约为0.5°.根据测试的光谙反射因数计M日光光度和色度时，应使用ClEIS中列出的CI雨准照明体A、D50,D65和D75的相对光滑分布.其他的QE日光照明体应使用CIE15中定义的适当的特征函数来确定.光源的紫外战区域（380nm）应通过紫外战止波光片迸行阻隔.当使用高的光源照明等级时，宜使用红外做止谑光片以最大程度地减少器件发热.4.3.4 标准照明几何条件4.3.4.1 通则EPD的性能测试应使用三种照明几何条件.4.3.4定义用于实现这些照躬几何条件的标准配置，也可使用其他的照明几何条件.SJ试所使用的照照几何条件的详细信息应在报告中注典.关于如何正造实现这些照明几何条件的更多指南，请查疑SID显示屏测试标次（见参考文献（1）.4.3,42定向照明当光源产生几乎平行的光线入射到DUT上时即为定向照明.光线与光钝的最大偏差取决于光源和测试点的直径，与光轴的最大俏圄角由公式（1）给出:arctan（rm.+r,Jd）<5*（1）式中：rt光源半径；d光源到测试点的距淘；自一测试点半径.光束横截面的照度均匀性应在S%以内.距圄DUT足够远的光源都可当成定向照明（例如太阳和月亮）.当模拟室外定向环境照明（如太阳和月亮）时，光源对向角（从DUT观看）宜£0.5°.当光源尺寸小于光源和样品上测试点之间的距离时，可使用三种不同类型的光源实现定向照明.这些几何条件见图2.平面朗伯源，例如积分球的出口，见图2a）；球形各向同性光源，例如漫反射侦璃球内的白炽灯泡，见图2b）；带透断或反射tft的投影系统，见图2c）.定向照明是通过使用与DUT表面法统形成夹角的小口径（与到测试点的距离相比）光源予以实现.这种定向光源会在DUT上产生光辘.1.MD在光入射平面内以倾斜角仇放翼,并且其测试区域位于光珊的中心.光源和IMD可在一定的倾斜角度范围内进行调节，但IMD应保持在入射平面内（即#180。）.该配置见图3a.图3b）则为其在板坐标系统中的表示.DUT的测试区域为Di在IMD检3!器上所成像的部，.IVfDa)fM(H机图图4环形光照明测试设置示例（其中仇±=35±5e=20°）为确保测试误差在±5%以内，环形灯和IMD的对准精度宜为±0.7当使用环形光模拟室外定向环境照明时，光源（从DUT观看）的对向角2直£0.5。.在这种情况下，直使用光纤环形灯.标引序号说明：2j一环形光源河边；s一环形光法便向；C.一环影光源80离.r.一环形光源半径；O.一环形光源交度，m.一测区域H径：D.一光照射区域H径（光度）；%,一光照测试直径对边.图5环形灯将性的详如示意图与光柏的最大梅阁角由公式（2）给出：arctan-arctanrJ'-"J,)式中：G一环形光源距陶；r,环形光源半径；m,一测量区域m径.因此，环形光直径（D）宜至少是测量区域直径（mJ的6倍.照射区域直径（DJ直至少是测量区域直径（m,）的1.S倍.如果显示屏包含反射表面上面的厚薄膜层，宜注意测试环形光到显示屏最远可见层的距离（U）,而不是显示屏的上表面.如果从显示肝反射出来的光显现不对称的散射图案（例如矩阵散射），则宜使用环形光通进行定向照明测试.如果非对称散射时使用定向光遁，则可能导致结果无法直现（见参考文献（8）.4.3.44半球照明半球照明是均匀入射在显示表面上的全向光.它可用于侵拟室内环境的漫射背景照明或室外太阳光（遮挡百射阳光）的漫射照明在真实的半球照朗情况下，恒定亮度的入射光线在测试区域所有角度内都是均匀分布.这种测试几何条件的两种可能示例见图6,近似理想的半球照明（例如限大到90。侦用的亮度恒定）只能通过积分球来实现，其中测试前口（图6,装置A）或采样（装MtB）的直径比球体的亘径小.如果积分球直径至少是显示屏外对角域的7倍，则可把显示屏安装在球体的中心（图6,装置A）.对于大的显示屏来说，采样球（图6,装窗B）或半球会更合适.无论那种情况，配置均应港况标准比8°到d，:10照明/探测几何条件，其中d是漫射的标准符号.当进行半球照明测试时，应使用以下步骤.a）将显示解窗于积分球/半球的中心或紧靠采样球的采样.b）在测试区域建立所需的照明.色湿和照明光谓可通过位于显示屏测试区域附近的漫反射标准白板（图6,装置A）或临近采样口的采样球球（四6,装置B）的反射光测试得到.测试报告中应注明所使用的光源类型及其相关色温.图6使用装置A和装置B的半球照明测试几何条件示例C）调整1.MD使其与显示屏法线的夹角为8。10。,并通过球内壁中的测试端口对准显示屏的中心.也可通过倾斜积分球内的显示屏以实现所需的IMDI项角.IMD聚焦于显示屏表面，d）测试端口直径应比1.MD透搐的有效孔径大20%到30%.应注意:避免来自光源的百射光或任何表面（非显示屏本身）反射的亮光进入1.MD的透镜，以最大程度减少反射亮度测试的眩光干扰.将1.MD从测试端口往后移到球壁对1.MD不可见的位翼.另外，通常采样口直径要4.5工作标准器4.5.1漫反射标准器漫反射标准白板可用一个漫反射率大于或等于98%的样品，也可用不同灰度的白板,对于已知的探测几何条件和照明光谱及配置，标准白板的照度E按公式(3)计15:R一式中：E漫反射标准白板的照度；Rl”一该测试装置的校准光反射因数；1.l”一漫反射标准板的亮度.当照明装置是均匀的半球照明时，则Rn等于光反射率pg.标准板的光反射率仅在半球照明下校准才有效,如果使用任意角度的定向光源，则不认为其发光反射率为该装置或光谱的光反射率因子.术语“光反射率”和”光反射应分别缩写为“反射率”和"反射因数".4.5.2钺面反射标准器里玻璃(例如BG-100o)或中性高密度吸收沌光片(光学密度24)可用于确定光源亮度1”通过公式(4)计算:I：上(4)1.1式中：1.一光源亮度；1.o,一黑玻埼反射的虚像亮度值；AII一茶俄埔的tft面反射率.当进行确面测试时，探测器聚焦于光源的虚像上.醇破遮可视为是具有4%至5%的低康面反射率的前表面潦.当测试几何条件下，不能亘接测得光源亮度而只能通过设面测试时，则可使用标准黑装璃.在与反射测试相同的数履级情况下,低镜面反射率的黑色玻璃镯用于测试光源亮度.黑色坡璃的镜面反射率受镜面角度、照朗光语以及镜面清洁度的影响.当测试的几何条件改变时应重新校席(见叁考文献(1).4.6标准测试区域位置4.6.1 矩阵显示器可在DUT表面的买个指定位道进行亮度.光谱分布和/或三刺激值的测试.为此，将显示器表面分成25个相同的虚构矩形，见图8.除非另有说照，否则应在每个矩形的中心进行测试.矩形的顶序从显示器的中心开始，开以顺时针的螭旋方式向边缘推进.注意显示器上的测试区域不应亚Si.在水平（三）和垂直(V)方向指定位置的测试场定位偏差应分别在H和V的7%之内.显示器或探测器应沿水平和垂直方向平移，以确保能测试所需的显示位2三,所有测试均垂亘于显示器.任何与上述标准位31不同的情况，应在报告中予以说明.ffl8标准测试位置XS、!)-8-三z4.6.2 段码显示器段码显示器的标准测试位置与上述矩阵显示器所规定的位置相同.但是，对于段码显示器来说,所有测试均应在字段的中心进行，且选中的字段直接近指定矩形的中心.因此，当要求在位SiPi(i=024)进行测试时，宜使用最接近矩形P.中心的字段的几何中心来定位探测器.与上述标准位置有任何不同的情况,应予以详细SS述.测试区域应完全包含在字段之内.5光学测试方法5.1 反射率测试5.1.1 目的此方法的目的是确定EPD在规定照明条件下的反射率指标.5.1.2 测试条件3!试条件如下.a)仪器驱动电遁、驱动信号设备、积分球、采样球和半球、定向光源.对于光谱测试,需要能测试亮度和光谙福亮度的光谙辐亮度计，以及巳知半球光谐漫反射系数和特定几何条件下校准的定向光滑漫反射系数的标次泡反射白板.对于充度测试,需要能海试充度的检测器,以及已知半球漫反射光谱反射因数和特定几何条件下校准的定向光谱反射因数的标准浸反射白板.如果接触式光潜翱射度计与经校准的非接触式光谱辐射度计或光度计相比证实有效时，也可使用接触式光谱辎射度计测试.b)照明条件附录A(资料性)日光色域的计算方法Al目的本方法的目的是描述一种在三维QESB色空间中讨口分散色点的立体色域的过程.A.2计算立体色域的过程A.1描述了计算立体色域的过程.图A,1计算立体色域的分析流程图在必定的日光条件下，至少测试显示屏的红色，嫁色，宦色，青色，品红色，黄色，黑色和白色8个色度.表A.1为使用SRGB基色测试时的示例，其中$RG8基色使用了CIED65白点开将白色亮度(Y)归一化为100%.表A.1SRGBlS色的三色值SS色P0Xqar>bYQ.arvbZQxb红色0.640033041.23921.2641.933绿色0.3000&0035.75871.51711.919蓝色0.1500.06018.0487.21995.053育色0.2250.32953.80678.736106.973BtI色0321015459.28728,48396.986黄色0.419OSOS7699892,78113.853“色0000000.0白色0.31270.329095.046100.000108.906用公式(30)公式(33)将所有颇色点就换为QEIAB籁色空间在C正1.AB颜色空间中设置的SRGB颜色示例如表A.2和图A.2.表A.2QE1.AB,领色空间中表示的SRGBRji色示例毋色(一弊.2KX2HT.nM<-2i.94.4%.Dffft(-4H.I-K.i9om7o6oso4n7/-50!OO段色.b,1.*红色80.10567.22351233绿色-86.18883.18687.737茎色79.1%107,85432303Jl色-48.084-1412891117品红色98.25060.83360.320黄色-21.56194.48897.138烹色000白色00100图A2ClE1.AB色空间中SRGB的立体色域图示将显示色点中包含的所有四面体相加来计算立体色域且与SRGB立体色域相比，该计算的立体色域以体枳百分比形式记录.表A.3为在QE1.AB色空间中计算出的SRGB立体色域的显示示例.表A.3ClE1.AB色空间中SRGB立体色域的示例立体色域总和8.20XKP用对于SRG8的百分比100%A.3计算QElAB立体色域的表面细分法A.3.1目的在QEI931XYZ=色值中，此算法接受指定的任意一组色域尖角情况.最小的颜色集是红色，绿色，蓝色，青色，品红色，黄色，黑色和白色.至少需要八个颜色角条件，才能使XYZ值排列在输入变量P行中，输出值则是计算出的立体色域.A.3.2假设假设CIEXYZ色空间中的色域定义为指定尖角情况的凸包.C】EIAB色空间中的色域将是此凸包,且在QEXYZ空间中通过具有最大光度（视为白点）的尖角情况进行归一化，然后转换为ClE1.AB色空间，其中该色空间不再是完墨的凸面.A.3.3运JI法则运投法则如下.a）获得P中颜色尖角点的凸包.将此包型表面的S0分曲面存储在T中.将总体枳V初始化为0.注1在这种情况下,拐角点是标准的RGBCMYKW88色.b）计算色域中间点的点P的平均值，并存储在Pm中.C）对于T中的每个三角形表面碎片，D令S等于在1.,ab坐标中具有范围大于10的边缘数量.注2使用范国不能被当作长度,因为它们计算时速度更快.2）如果s=o,则计U在表面碎片的顶点和Pm之间定义的体积.将此体积添加到V中.3）如果S=3,则计算CIEXYZ空间中的中点,并将三角形碎片细分为由每个尖角顶点定义的4个子碎片，其中该角顶点具有两个最近的中点和3个中点.对每个三角形子碎片田复计口3次.4）如果S=I或2,则在ClE1.AB中计算只有最大更囹的边缘的CIEXYZ空间中的中点,并沿中点和相对顶点之间的线格三角形碎片细分为2个子碎片.对每个三角形子图块至官计算3次.d）返回现在V中包含的总体积.A.3.4软件示例MATlAB代码提供了立体色域计算的软件使用示例.该代码由两个主要模块执行，下面两个枢GB/T43789.312024/IEC62679-3-1:2014图柢述了这些模块.为了执行Matlab程序，将相应的SRGB数据加载到内存中来执行以下命令：>>P=GetGamutCornersO$RGB)DefaultD65whiteused0000.41240.21260.01930.77000.92780.13850.35760.71520.11920.53810.78741.06970.18050.07220.95050.59290.28480.96990.95051.00001.0891其中数据矩阵对应于以下三色坐标，如下所示：×YkZkXrYRZRXvYvZvXgYgZgXKYkZkXcYcZcXeYeZbXmYmZmXwYWZw然后通过执行以下命令来获得QE1.AB立体色域：>>CIE1.abVo1.subd(P)ans三8.2013x10'ClE1.abVo1.subd.mfutionIvJCIE1.abVoI.sub(P).t11-EachrOfPcontains×ZstimulusvaluesOfgamutcornerpoints.%The3DgamutisdefinedastheconvexhullofthesepointsinXYZspace.%ThesurfaceisrecursivelysubdivideddowntoathresholdscaleinCIE1.AB%andthevolumemadebyeachsurfacetiletoaCentralpointissummedthresh=10；%CIE1.absubdivisionthreshold%GetthehulldefinedbythepointsT=ConvhuIIn(P)；%Getthewhitepoint(takenastheprimarywiththemaximumY)W,i1-max(P(2);W=Pt.,3；%NormalisethegamuttothewhitepointPn=P./(repmat(W.size(P.1).1)；%getthemid-pointPm-mean(Pn)；%addontheCIE1.abpointsPn=Pn×YZ21.ab(Pn;；PmP-.XVZ21.ab(P)：%Calculateandsumthe1.abvolumeofeachSurfacetiletothemid-pointvOiIrsize(TiI).v=v+SubD1.abVol(Pn(T(nl:),:),PmJhresh)；end%sub-functions%XYZ21.abConvertsXYZValuGsarrangedincoIumnsto1.abfunctiontl=XYZ21.ab(t)i=(t>O.008856)；t(P=t(i).A(1/3)；t(-i)三7.787t(-i)÷16116:t=116t(I.2).16,500-(t(：,1)-t(：.2).200w(t(：,2)-t(：.3)；end%RecursivefunctiontodevideuptheSurfacetilethenreturnthevolumefunctionv)>SubD1.abVo(vp.c,th)%GetthemaxGxtentoleachGdg(quckrthanlengthcalculation)mmaxtabv-circshift<vp，1)，：2)；(ounthomanyGdge$haveext«>!$IargerthanthethresholdS=sum(m>th)；if(s-O),%noedgeslarger：returnthevolumev=abs(dt(vp(:.4:6)-repmat(c<1.4：6),3.1)/6)；els,if(s=3),%,aIledgeslarger：dividetileinfourgGtmidpointsip=(vp(：.1：3)÷circshift(vp(：.1：3).1)/2;%CaIcuIateCIE1.abpointsofthemid-pointsip=(ip.XYZ21.ab(ip)：%andCallrecursivelyforeachsub-tileV=SubD1.abVo(vp(1.：)：ip(1：2,：).c,th)：v=vtSubD1.abVol(vp(2,：);ip(2：3,：),c,th)；v.v÷SubD1.abVo(vp(3.：)：ip(1：2：3.：).c,th)；v=v+SubD1.abVo(ip,c,th)：else%oeortwoedgeslarger：SplitthetileontheIargestedge%Shifttheorderso12hastheargestextetm,i=ma×(m)；vp=circshift(vp.2-i)；%calculatethemid-poi11tof1-2andtheCE1.abpointip-(VP11：3)+vp(21：3)2IpIipAY721.ab(ip)：%andCallfecursivelyforthetwosub-tilesv-SubD1.abVo(vp(l3.：);ip.c.th)：v-v»SubD1.abVol(vp(2：3.0;ip,c,th)；endendendGetGamutCorners.mfunction(P三GetGamutCorers(P.wh)%GETPRIMreturnsasetofcolourcornerpointsbasedonastandardgamut%inputstringmusiconIainoneof：%,sRGB*,*Rc709',*EBU'.*NTSC*%optionallyoneof%'D5''D55,D65'.*D75','IHA''川E'ifischar(P)ifnargin<2wh=P:endifstrfind(P.'sRGBl)strfid(P.'Rec709')prim=(0.64.0.33:0.3.0.6:0.15,0.06)；lSifstrfindPJEBU')prim=O.64.0.33:0.29.0.6:0.15.0.06：lseifstrfid(P,NTSCt)prim-0.67.0.33；0.21.0.71;0.14.0,08；elseerror<lnon-validcolourprimaryspecification1)；endP=prim；endIfischar(Wh)ifstrfind(w"D50*)wh=0.3457.0.3585：elseifstrfind(wh,'D55,)wh=0,3324.0.3474；elseifstrfind(wh.1D65t)wh-(0.3127,0.3290；elsfstrfind(wh.'D75')wh=(0.2990.0.3149；lsifstrfid(wh.IIIA,)WMO.44757。40745；elseifstrfindfWhZIIIE*)wh=0.3333.0.3333)；elsewh=0.3127.O3290；dlay('DefauItD65whiteused)：endendwh=(wh.1-sum(wh)*'wh(2)iP=P(P.O；PPrrPn,l»-Pfj1j3)P-KRY(；CHMWP=Io00：P(1.：)um(P(1：2.：)P(2.：)um<P(2J.：)；.P(3,：)；sum(P(«31.：)；sum(P)i；end参考文献(1) InternationaICommitteeforDispayMetrology*InformationDisplayMeasurementsStandard.Ver.1.03,(2012)(2) E.F.Kelley.M.1.indfors.andJ.Penczek.Displaydaylightambientcontr