显色指数的计算

显色指数的计算光源显色性定义:是指与参照标准下相比较,一个光源对物体颜色外貌所产生的效果。1965年CIE制定一种评价光源显色性的方法,简称“测验色”法,1974年修订后,正式向国际上推荐使用。此方法是用一个显色指数量值表示光源的显色性。光源的显色指数是待评光源下物体的颜色与参照光源下物体颜色相符程度的度量。为了符合人类长期的照明习惯,CIE规定5000K以下的低色温光源用普郎克辐射体作为参照光源,色温5000K以上的用标准照明体D作为参照光源,设定参照光源的显色指数为100。评价时采用一套14种试验颜色样品,其中1到8用于光源一般显色指数(8个数平均值),各试验色样的数值称之为特殊显色指数。我们平时说的“显色指数”,即是一般显色指数的简称。若某个试验色样在待评光源与参照光源照明下有颜色差Ei那么:特殊显色指数10046iREi；一般显色指数81/8aiRR一、根据待测光源的光功率谱分布,计算待测光源的色度坐标kx，ky，ku,kv及相关色温CT。1、待测光源的色度坐标kx，ky，ku,kv的确定使用光谱仪测出待测光源的光谱功率分布函数()sP，计算光源的三刺激值X，Y，Z：780380()()msXKPxd，780380()()msYKPyd；780380()()msZKPzd;其中:mK为辐射量和光度量之间的比例系数,为常数,等于683lm/W。()x，()y，()z为CIE1931标准色度观察者光谱三刺激值（注：此处的三刺激值可由1931CIE-RGB系统标准色度观察者光谱三刺激值()r，()g，()b来确定。某一波长的光谱刺激()r，()g，()b与光谱色度坐标()r，()g，()b关系如下：rrrgb，ggrgb，bbrgb；某一波长的光谱刺激()r，()g，()b与()x，()y，()z色度坐标关系为：0.49000()0.31000()0.2000()()0.66697()1.13240()1.20063()rgbxrgb，0.17697()0.81240()0.01063()()0.66697()1.13240()1.20063()rgbyrgb，0.00000()0.0100()0.99000()()0.66697()1.13240()1.20063()rgbzrgb，（该三式可由矩阵表示）。在CIE1931标准观察者光谱三刺激值()r，()g，()b中，规定()g与明视觉光谱光效率函数一致，即()g=()V。故使用已知的CIE1931色度图的光谱轨迹色度坐标()x，()y，()z和光谱光效率函数()V就可以球的光谱三刺激值()x，()y，()z()()()()xxVy，()()yV，()()()()zzVy。）在实际测量中,在380nm～780nm区间的光源由此可以得出在x-y色度图中x,y的值XxXYZ,YyXYZ，ZzXYZ以及在CIE1960均匀色度标尺图u-v上色坐标u,v的值4153XuXYZ,6153YvXYZ或42123xuxy，62123yvxy。2、相关色温CT的确定如果一个光源发射光的颜色(即光色,又称色品)与某一温度下的黑体发射光的颜色相同。那么，此时黑体的绝对温度值就叫做该光源的颜色温度（简称色温）。黑体发射光的相对光谱功率分布由普朗克定律给出：2/511(,)(1)cTPTce，其中:T—黑体的绝对温度（K）;-波长（nm）；1c—第一辐射常数，1c=3.7417749x10-16Wm2；l—第二辐射常数，2c=1.4388x10-2m.K。当光源发射光的颜色和黑体不相同时，常用“相关色温”的概念来描述光源的颜色。相关色温的定义是：在某一确定的均匀色度图中，如果一个光源与某一温度下的黑体具有最接近相同的光色，此时黑体的绝对温度值就叫作光源的相关色温。当表示光源光色的色坐标点落在黑体色轨迹上时，说明该光源的光色与某一温度下的黑体的光色相同，计算出的就是光源的色温值。反之，当表示光源光色的色坐标点落在黑体色轨迹以外时，计算出的就是光源的相关色温值。并且,，当表示光源的色坐标点（u,v）在u~v图中偏离黑体色轨迹越远，相关色温的概念越弱。a、计算黑体的色坐标系数u、v值以及选择适当的若干条等相关色温线的斜率在CIE1960UCS均匀色度图中,一种光色对应一个点,有独立的色坐标(u,v)。当黑体的温度从较低的值逐渐升温至K,那么在UCS色度图中,代表黑体光色的色坐标点将会形成一段连续的曲线,称为黑体色轨迹(简称黑体迹)。在均匀色度图中,等相关色温线（下面简称等温线）是一系列垂直于黑体色轨迹（曲线）的直线簇。已知黑体相对光谱分布()P，求得黑体色坐标系数（u,v）等温线的斜率用m表示，它是相关色温值T的函数，且有：1/ml。式中:l为黑体色轨迹与该等温线交点(垂足)处的切线的斜率:''''''''/()3()/2()10()dvdvdTXYXYYZYZldududTXZXZXYXY式中:X、Y、Z是黑体的三刺激值;'X、'Y'Z为黑体三刺激值对于黑体温度T的导数,由下式给出:''/(,)()TXdXdTPTxd，''/(,)()TYdYdTPTyd,''/(,)()TZdZdTPTzd,其中'(,)TPT是(,)TPT对的偏导数。222//'6212(,)(,)//()1cTcTTPTdPTdTccTee。(b)计算光源的(相关)色温如果已知光源色坐标(u,v),便可以计算出光源的(相关)色温。下面是三种(相关)色温的计算方法直接内插法设表示光源的色设表示光源的色品的坐标点(u,v)位于下图中画出的相邻两条等温线T1和T2之间。d1、d2为色品的坐标点到T1、T2的距离,轨迹上两点之间的距离d1+d2近似与(1/T1)-(1/T2)成正比。则光源相关色温Tc近似地可由下式计算得到:112121111cdTTdTT；其中1d，2d由下式得到2()()1iiiiivvmuudm。上式中的mi、ui、vi分别是第i条等温线Ti的斜率与黑体色轨迹交点的色坐标系数。在导出时使用了如下的近似成立的条件:等温线T1和T2之间的黑体色轨迹曲线是一段圆弧，圆心为T1与T2线的交点;夹角12很小；倒色温值是沿弧方向上距离的线性函数。具体计算时，首先要计算出黑体等相关色温线簇与黑体色轨迹交点(垂足)的u、v值。将相关色温线簇的色温值或色温值倒数，及其与黑体色轨迹交点的u、v值做成一表格。根据光源的u、v值查表找到最近邻的两条等温线T1和T2。再由上述公式计算光源的色温值。三角形垂足法首先用与上述方法类似的步骤查表找到与光源色坐标S(u,v)最近邻的两条等温线T1和T2,设两色温线与黑体色轨迹的垂足分别为D1(ub1，vb1)、D2(ub2，vb2)。从光源色温坐标点S(u,v)向D1、D2做垂线，设垂足为D(ub，vb)。设L1=DD1,L2=DD2,根据下面的内插方法可计算出过垂足点D(ub，vb)的等色温线Tc；111221111()cLTTTTL；22111()()Lububvbvb，22222()()Lububvbvb色温逐次逼近法色温逐次逼近法是根据实际情况，凭经验预先设置一起始温度T1(例如做2856K色温实验时可设T1=2000K)。一方面T1的设置不能太接近预期值，需留有一定的余地；另一方面为了节省计算时间，不应偏差太大(如T1=500K)。此外，要预设增减步长STEP(例如STEP=40K)及最小阈值GATE(如GATE=0.0001)。此阈值表示了光源色坐标点(u，v)与待求色温线的接近程度。阈值GATE可根据精度要求适当调节。从T1开始,计算光源色坐标与该色温线的接近程度D1,如不满足所设的最小阈值条件(如D1GATE),则T以步长STEP增减(T2=T1+STEP)。再计算光源色坐标与色温线T2的接近程度……,依此循环进行，直到满足设定的条件。计算期间,步长STEP值应随时进调节,T的增减方式与它和光源色坐标(u，v)的位置有关(图4为此算法的示意图)。色温逐次逼近法不受所查表中等色温线簇间隔大小的影响,计算精度较高。二、由待测光源的光功率谱分布和1-14试验色的光谱辐亮度因数,ki,计算待测光源下1-14号试验色的色度坐标,kix，,kiy；并计算相应的,kiu，,kiv。,,,()()()kikikiXkPxd，,,,()()()kikikiYkPyd，,,,()()()kikikiZkPzd。其中,,100()()kikikPyd。,,,,,4153kikikikikiXuXYZ，,,,,,6153kikikikikiXvXYZ。三、根据Tc选择参照照明体,查表得参照照明体的色度参数ru，rv，rc，rd,以及在参照照明体下试验色的颜色空间坐标,riU,,riV,,riW。在待测光源不高于5000K时,以普朗克辐射体作为基准光源,相对光谱功率分布为2/511(,)(1)cTPTce；当待测光源高于5000K时,以组合昼光做为基准光源,在待测光源相关色温已知时,相对光谱功率分布为01122()()()()SSMSMS，0()S，1()S，2()S计算昼光光谱功率分布用的系数。M1,M2为与光源色坐标相关的量,其量值由下列公式计算给出：11.35151.77035.991140.02410.25620.7341ddddxyMxy，20.030031.442430.07170.02410.25620.7341ddddxyMxy，式中dx，dy为基准光源的CIE1931xy色品坐标值。其值由下列公式给出：23.0002.8700.2750dddyxx，963321010104.60702.96780.099110.244063dcccxTTT40007000cKTK；963321010102.00641.90180.247480.237040dcccxTTT700025000cKTK。ru，rv，由第二步中计算公式计算；rc，rd由下面第四步中公式计算；,riU,,riV,,riW由第五步中公式计算。四、根据转换下列公式,由待测光源色坐标ku，kv求kc，kd以及待测光源下试验色的色坐标,kiu，,kiv求,kic，,kid。并计算试验色的适应色位移',kiu,',kiv。1(410)cuvv，1.7080.4041.481vudv，,,',,,10.8720.404416.5181.481rrririkkkirrririkkcdCdcducdCdcd，',,,5.52016.5181.481kirrririkkvcdCdcd。五、计算待测光源下试验色的,kiU,,kiV,,kiW。13,,',,,',,,25()1713()13()kikikikikirkikikirWYUWuuVWvv六、求,riU,,riV,,riW，,kiU,,kiV,,kiW的色差iE。222,,,,,,()()()irikirikirikiEUUVVWW七、计算特殊显色指数10046iREi；一般显色指数81/8aiRR以上就是用标准方法计算显色指数的主要步骤,其目的就是分别求在待测光