CN105509901B - 发光体相对色温的高精度计量方法 - Google Patents
发光体相对色温的高精度计量方法 Download PDFInfo
- Publication number
- CN105509901B CN105509901B CN201511008796.6A CN201511008796A CN105509901B CN 105509901 B CN105509901 B CN 105509901B CN 201511008796 A CN201511008796 A CN 201511008796A CN 105509901 B CN105509901 B CN 105509901B
- Authority
- CN
- China
- Prior art keywords
- colour temperature
- calculate
- use formula
- formula
- illuminator
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 230000002596 correlated effect Effects 0.000 claims description 8
- 230000000875 corresponding effect Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 230000005457 Black-body radiation Effects 0.000 claims description 2
- 238000009529 body temperature measurement Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/60—Radiation pyrometry, e.g. infrared or optical thermometry using determination of colour temperature
Abstract
本发明公开了一种发光体相对色温的高精度计量方法,包括以下步骤:I)采用分光谱仪或光/颜色传感器测量出光源或发光体的CIE1960UCS色度坐标:uw,vw;II)采用迭代方法并以下式计量相对色温:其中j=0,1,..,直到Tj+1和Tj之间的绝对值小于10‑3为止,此时,Tj+1为测量的发光体的相对色温。本发明的优点是计量精度高,适用相对色温范围广。
Description
技术领域
本发明属于发光体如光源的相对色温测量装置的一部分,如分光亮度计或分光辐射度计等,具体是发光体相对色温的高精度计量方法。
背景技术
色温是描述光源或发光体的一个重要指标,在颜色测量中有重要意义。光源的色温定义为与光源辐射颜色相同的黑体所对应的色温。但实际上完全与光源的颜色相匹配的黑体可能不存在,这就延伸出相对色温的概念。目前,根据国际照明委员会(CIE)的定义,在CIE1960 UCS色度坐标u,v下,最接近于给定光源色度坐标uw,vw的黑体辐射坐标点u(T),v(T)所对应的色温T为光源的相对色温。
因此测量光源的相对色温可转化为求距离的平方函数f(T)关于T的最小值问题。这里函数f(T)由下式定义:
(1)
因此测量光源的相对色温变为求解一个非线性优化问题,相当复杂。
现有技术中人们测量相对色温的基本步骤为:
I)测量光源的色度坐标uw,vw(或xw,yw);
II) 通过色度坐标预测出相对色温。
其中步骤I)光源的色度坐标uw,vw的获取办法:
1)光源的色度坐标uw,vw可以通过分光光谱仪测量并计算出;
2)通过任何光传感器检测出CIE三刺激值的办法给出光源的色度坐标uw,vw;
3)理论上给定光的相对光谱功率分布函数,其色度坐标uw,vw可以计算出。
步骤II)中,当光源的色度坐标已知,便可通过算法,如Robertson方法,McCamy方法,或Qiu方法等预测出相对色温,其计算出的光源或发光体的相对色温(记为CCT)的显示表达式,可统一写成,这里的表达式q视方法不同而不同, 如Robertson方法,McCamy方法或Qiu方法都属这类;或者,还可以采用迭代方法预估出相关色温,而迭代方法中需要的函数值以及相关的一阶导数和二阶导数都采用近似的手段获得,Krystek和Gardner方法属于这类。但这些方法都只能给出近似的相对色温,精度不佳,因此有必要对此进行改进。
发明内容
本发明的目的是为了克服现有技术存在的缺点和不足,而提供一种计量精度高的发光体相对色温的高精度计量方法。
为实现上述目的,本发明的技术方案包括以下步骤:
I) 采用分光光谱仪或光/颜色传感器测量出光源或发光体的色度坐标:uw,vw;以最接近于坐标u w ,v w 的黑体辐射坐标u(T),v(T)所对应的色温T为光源的相对色温,被测光源的相对色温T可转化为求以上两坐标点之间距离的平方函数f(T)的最小值问题,即:
(1)
II) 采用迭代方法并以下式计量相对色温:
(2)
直到Tj+1和Tj 之间的绝对值小于10-3 为止,此时,Tj+1为测量的发光体的相对色温;
(3)
(4)
式(1)-(4)中的u(T),v(T)由式(5)-(9)计算出:
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
进一步设置是式(2)的迭代方法是:
首先预置参数:设定允许精度Tol=0.001; 最大允许迭代次数ITMAX=10; 迭代次数计数器IT=1; 然后用 Qiu方法根据已经测得的(uw,vw)预测CCT的初始估值T;
该Qiu方法为:令P=( uw,vw),P1=(0.2861884,0.246725)和P2=(0.328151,0.1333451);d1, d2分别是P和P1及P和P2的距离,是P和P1连线与u轴的负方向所成的夹角;是P和P2连线与u轴的负方向所成的夹角,预测估值T由式(18)-(22)计算出:
(18)
(19)
(20)
(21)
(22)
步骤1:用式(9)计算X(T),Y(T),Z(T)
步骤2:用式(5)计算u(T),v(T);
步骤3:用式(10),(11)计算X’(T),Y’(T),Z’(T);
步骤4:用式(13),(14)计算X’’(T),Y’’(T),Z’’(T);
步骤5:用式(12)计算u’(T),v’(T);
步骤6:用式(15)计算P(T),Q(T);
步骤7:用式(16)计算P’(T),Q’(T);
步骤8:用式(12)计算u’’(T),v’’(T);
步骤8:用式(3)计算f’(T);
步骤9:用式(14)计算f’’(T);
步骤10:DT=f’(T)/f’’(T);
步骤12:若DT的绝对值不超过Tol,则接受T1相关色温,并停止计算,否则更新
T=T1;
IT=IT+1;
步骤13:若IT不超过ITMAX转到步骤1,并依次计算,否则接受当前的T为相关色温,并停止计算。
本发明的优点是解决了两个技术问题:(1)本方法预测的相对色温与理论值保持一致,解决了现有相对色温测量中普遍存在的误差问题;(2)本方法的适用相对色温范围为500K~30000K,解决了现有色温测量中适用范围的限制问题。
下面结合说明书附图和具体实施方式对本发明做进一步介绍。
附图说明
图1 本发明的原理图。
具体实施方式
下面通过实施例对本发明进行具体的描述,只用于对本发明进行进一步说明,不能理解为对本发明保护范围的限定,该领域的技术工程师可根据上述发明的内容对本发明作出一些非本质的改进和调整。
实施例1:
用分光光谱仪实测在360nm到830nm的日光光谱,1nm间隔,计算获得色度坐标(uw,vw)=(0.1982, 0.3104);(xw,yw)=(0.3109, 0.3246); 理论相对色温值(四舍五入)为6642K;
基于本方法预测的相关色温为:6642 K; 若用其他方法: Qiu, Javier et al,McCamy,Robertson,Gardener,和Krestek可分别获得如下结果:6640,6640,6642,6639,6646,6644。可见除了本方法和McCamy方法外,其他方法均不能给出准确结果;
例2:
采用CIE标准钨丝灯光源(CIE A光源),通过其360nm到830nm,1nm间隔的光谱,计算获得色度坐标(uw,vw)=(0.2560, 0.3495);(xw,yw)=(0.4476, 0.4074); 理论相对色温值(四舍五入)为2856 K;
基于本方法预测的相关色温为:2856 K; 若用其他方法: Qiu, Javier et al,McCamy, Robertson, Gardener, 和Krestek可分别获得如下结果:2854, 2790, 2857,2855, 2858, 2855。可见除了本方法外,其他方法均不能给出准确结果;
例3:
用分光光谱仪实测在360nm到830nm的办公室灯光光谱,1nm间隔,计算获得色度坐标(uw,vw)=(0.1977, 0.3225);(xw,yw)=( 0.3267, 0.3552); 理论相对色温值(四舍五入)为5740K;
基于本方法预测的相关色温为:5740 K; 若用其他方法: Qiu, Javier et al,McCamy, Robertson, Gardener, 和Krestek可分别获得如下结果:5738, 5724, 5739,5738, 5734, 5743。可见除了本方法外,其他方法均不能给出准确结果。
本发明与已有方法的适用范围和计算精度比较如下表所示:
方法 | 适用范围 | 例1 | 例2 | 例3 |
Qiu | 2500 K~10000 K | 6640 | 2854 | 5738 |
Javier et al | 3000K~ 8×105 K(分两段) | 6640 | 2790 | 5724 |
McCamy | 2000K ~12500K | 6642 | 2857 | 5739 |
Robertson | 取决于查找表范围 | 6639 | 2855 | 5738 |
Gardener | 1000K~10000 K | 6646 | 2858 | 5734 |
Krestek | 1000K~15000K | 6644 | 2855 | 5743 |
本方法 | 500K~30000K | 6642 | 2856 | 5740 |
理论CCT | 所有范围 | 6642 | 2856 | 5740 |
Claims (2)
1.一种发光体相对色温的高精度计量方法,其特征在于包括以下步骤:
采用分光谱仪或光/颜色传感器测量出发光体的CIE1960UCS色度坐标:uw,vw;以最接近于坐标uw,vw的黑体辐射坐标u(T),v(T)所对应的色温T为光源的相对色温,被测光源的相对色温T可转化为求以上两坐标点之间距离的平方函数f(T)的最小值问题,即:
(1)
采用迭代方法并以下式计量相对色温:
(2)
直到Tj+1和Tj 之间的绝对值小于10-3 为止,此时,Tj+1为测量的发光体的相对色温;
(3)
(4)
式(1)-(4)中的u(T),v(T)由式(5)-(9)计算出:
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)。
2.根据权利要求1所述的一种发光体相对色温的高精度计量方法,其特征在于:
式(2)的迭代方法是:
预置参数:设定允许精度Tol=0.001; 最大允许迭代次数ITMAX=10; 迭代次数计数器IT=1;然后用 Qiu方法根据已经测得的(uw,vw)预测CCT的估值T;该Qiu方法预测CCT的估值T的过程为:令P=( uw,vw), P1=(0.2861884, 0.246725) 和 P2 =(0.328151, 0.1333451);
(18)
(19)
(20)
(21)
(22)
步骤1:用式(9)计算X(T),Y(T),Z(T)
步骤2:用式(5)计算u(T),v(T);
步骤3:用式(10),(11)计算X’(T),Y’(T),Z’(T);
步骤4:用式(13),(14)计算X’’(T),Y’’(T),Z’’(T);
步骤5:用式(12)计算u’(T),v’(T);
步骤6:用式(15)计算P(T),Q(T);
步骤7:用式(16)计算P’(T),Q’(T);
步骤8:用式(12)计算u’’(T),v’’(T);
步骤8:用式(3)计算f’(T);
步骤9:用式(14)计算f’’(T);
步骤10:DT=f’(T)/f’’(T);
步骤12:若DT的绝对值不超过Tol,则接受T1相关色温,并停止计算,否则更新
T=T1;
IT=IT+1;
步骤13:若IT不超过ITMAX转到步骤1,并依次计算,否则接受当前的T为相关色温,并停止计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511008796.6A CN105509901B (zh) | 2015-12-29 | 2015-12-29 | 发光体相对色温的高精度计量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511008796.6A CN105509901B (zh) | 2015-12-29 | 2015-12-29 | 发光体相对色温的高精度计量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105509901A CN105509901A (zh) | 2016-04-20 |
CN105509901B true CN105509901B (zh) | 2018-11-02 |
Family
ID=55718048
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511008796.6A Active CN105509901B (zh) | 2015-12-29 | 2015-12-29 | 发光体相对色温的高精度计量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105509901B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108917932B (zh) * | 2018-07-16 | 2020-09-11 | 宁波公牛光电科技有限公司 | 一种合成光谱的评价方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102903343A (zh) * | 2012-09-26 | 2013-01-30 | 苏州佳世达电通有限公司 | 快速调整色温方法 |
CN103048050A (zh) * | 2012-07-30 | 2013-04-17 | 中国科学院自动化研究所 | 一种测量目标真实温度的辐射测温方法和仪器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009056665A1 (de) * | 2009-12-02 | 2011-06-09 | Osram Opto Semiconductors Gmbh | Verfahren zur Klasseneinteilung strahlungsemittierender, optoelektronischer Halbleiterbauelemente |
-
2015
- 2015-12-29 CN CN201511008796.6A patent/CN105509901B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103048050A (zh) * | 2012-07-30 | 2013-04-17 | 中国科学院自动化研究所 | 一种测量目标真实温度的辐射测温方法和仪器 |
CN102903343A (zh) * | 2012-09-26 | 2013-01-30 | 苏州佳世达电通有限公司 | 快速调整色温方法 |
Non-Patent Citations (2)
Title |
---|
LED光源相关色温计算新方法;林岳等;《计量学报》;20100930;第31卷(第5期);第436-439页 * |
光源相关色温计算方法的讨论;代彩红等;《计量学报》;20000731;第21卷(第3期);第183-188页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105509901A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102067883B1 (ko) | 분광 복사계를 교정하는 방법 | |
US8144316B2 (en) | Instrument and method for measuring total luminous flux of luminous elements | |
EP3136065B1 (en) | Method for measuring actual temperature of flame using all information about radiation spectrum and measurement system | |
CN102959370B (zh) | 积分球光度计及其测量方法 | |
CN104062010B (zh) | 一种优化定标算法的分光光源颜色照度测量仪器 | |
CN105509901B (zh) | 发光体相对色温的高精度计量方法 | |
CN106908150B (zh) | 嵌入式系统相关色温快速计算方法 | |
CN105675148A (zh) | 一种标准光源对色观察箱的光源校准检测方法 | |
Yuan et al. | LED-based spectrally tunable light source with optimized f itting | |
CN106908149A (zh) | 一种机器人物体颜色识别系统及方法 | |
Krüger et al. | Spectral mismatch correction factor estimation for white LED spectra based on the photometer’s f1′ value | |
CN116256069A (zh) | 一种双波长温度测量系统定标方法 | |
Winter et al. | Convenient integrating sphere scanner for accurate luminous flux measurements | |
Pavlova et al. | Optimization of the methods for measuring color characteristics of light-emitting diodes in laboratory conditions | |
KR101557769B1 (ko) | 복사조도 측정장치 및 그 측정방법 | |
JP2010169427A (ja) | 測色方法及び測色装置 | |
Pousset et al. | Uncertainty evaluation for measurement of LED colour by Monte Carlo simulations | |
CN105092088A (zh) | 一种宽带荧光光谱测温方法 | |
CN112014764A (zh) | Led灯具加速老化测试中光通量维持率修正的方法及装置 | |
CN204206520U (zh) | 一种色温可调黑体光源 | |
CN107748010B (zh) | 一种探测器光谱测量方法 | |
Achmadi et al. | The effect of glass neutral density filter on illuminance measurement error | |
CN109341856B (zh) | 光谱不对称单色led的光谱分布函数拟合方法 | |
Yuan et al. | Integral colorimeter based on compound LED illumination | |
CN114486183B (zh) | 一种灯具配光测试方法以及灯具配光测量系统 |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20160420 Assignee: Intelligent lock Research Institute of Wenzhou University Assignor: Wenzhou University Contract record no.: X2020330000086 Denomination of invention: High precision measurement method of relative color temperature of luminescent body Granted publication date: 20181102 License type: Common License Record date: 20201030 |