CN106908150B - 嵌入式系统相关色温快速计算方法 - Google Patents

嵌入式系统相关色温快速计算方法 Download PDF

Info

Publication number
CN106908150B
CN106908150B CN201710072609.3A CN201710072609A CN106908150B CN 106908150 B CN106908150 B CN 106908150B CN 201710072609 A CN201710072609 A CN 201710072609A CN 106908150 B CN106908150 B CN 106908150B
Authority
CN
China
Prior art keywords
color temperature
values
calculating
light source
value
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
Application number
CN201710072609.3A
Other languages
English (en)
Other versions
CN106908150A (zh
Inventor
叶新怡
章致好
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HANGZHOU XINYE OPTOELECTRONIC ENGINEERING CO LTD
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201710072609.3A priority Critical patent/CN106908150B/zh
Publication of CN106908150A publication Critical patent/CN106908150A/zh
Application granted granted Critical
Publication of CN106908150B publication Critical patent/CN106908150B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J3/52—Measurement of colour; Colour measuring devices, e.g. colorimeters using colour charts
    • 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
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J2003/467—Colour computing

Abstract

本发明提供了一种嵌入式系统相关色温快速计算方法。所述方法中的u、v值表，是由1000K至4999K色温的u、v值表和5000K至25000K色温的u、v值表，再添加900K色温和26000K色温的u、v值合并完成，其中900K色温和26000K色温的u、v值用于色温区域划分；该表中的u、v值直接存储于便携式测量仪器中微处理器的闪存中；在针对被测光源的相关色温计算中，可以直接获取表中数据，加快运算速度。

Description

嵌入式系统相关色温快速计算方法
技术领域
本发明涉及嵌入式系统相关色温快速计算方法。
背景技术
光源在温度T时所呈现的颜色与在某一温度Tc时的黑体颜色相同时，则称Tc为该光源的色温度。如果光源在温度T时的相对光谱功率分布所决定的色度坐标，不在色品图的黑体温度轨迹上，而在此轨迹的附近。这时，用色度坐标与其最靠近的黑体温度表示该光源的温度，称之为该光源的相关色温。
光源(相关)色温的计算方法
由光源的三刺激值X、Y、Z，计算光源在CIE1931-XYZ色品图中的x、y：
计算CIE1960UCS均匀色度坐标系中的色坐标u、v值：
在CIE1960UCS均匀色度图中，每种光源对应一个独立的色坐标(u，v)。当黑体的温度从较低的值逐渐升温至∞K，在UCS色度图中，代表黑体光色的色坐标点将会形成一段连续的曲线(如图1所示)，称为黑体色轨迹。在均匀色度图中，等相关色温线是一系列垂直于黑体色轨迹(曲线)的直线簇。
相关色温计算方法较多，主要相关论文如下：
《光源相关色温计算方法的讨论》比较了内插法、三角形垂足法、逐次逼近法计算相关色温的精度；其中逐次逼近法是精度最高的相关色温计算方法；
《光源相关色温算法的比较研究》比较了三角垂足插值法、黑体轨迹的Chebyshev法、模拟黑体轨迹弧线法、McCamy近似公式法计算相关色温的精度；
《光源相关色温计算方法的研究》利用遗传模拟退火算法得到的计算色温和相关色温的经验公式为，是一种近似的计算方法；
《关于光源颜色温度标定方法的讨论》主要关注颜色温度标准灯及其标定和测试方法；
《白光LED光学仿真中的相关色温计算》采用逐点逼近法计算相关色温，重点关注TracePro光学软件对白光LED的模拟仿真；
《LED光源相关色温计算新方法》通过等点间隔和等几何距离间隔，相关色温计算速度较快且精度高；
《二分法优化计算LED光源相关色温》采用多重二分的计算方法快速计算相关色温。
在光源相关色温测量中，积分式测量方法是应用光学滤光片组将探测器的光谱响应曲线匹配成CIE推荐的XYZ光谱三刺激值曲线，探测器可以对被测光源的光谱组成进行快速积分测量。积分式测量仪器特点是系统简易，性价比高，测量速度快，可靠性高。作为相关色温的便携式测量仪器采用积分式测量方法，获得被测光源的XYZ，进而计算相关色温。便携式测量仪器采用嵌入式系统设计，采用单片机、MCU等而不是PC作为核心处理器，因此计算能力受限，对于较为复杂的计算过程，则需要较长的时间，不适用便携式测量仪器的快速测量功能。因此研究适应于嵌入式系统的快速且高精度的相关色温计算方法，具有重要意义和应用价值。
发明内容
本发明所要解决的技术问题是提供一种能够快速获得相关色温的嵌入式系统相关色温快速计算方法。为此，本发明采用以下技术方案：
1、嵌入式系统相关色温快速计算方法，其特征在于它包括以下步骤：
1.1建立黑体和标准照明体D在不同色温时的u、v值表
以黑体和标准照明体D为标准光源，计算黑体和标准照明体D在不同色温时的u、v值，并保存为色温的u、v值表。计算过程中将色温分为五种情况，分别为小于1000K，大于等于1000K小于5000K，大于等于5000K小于7000K，大于等于7000K小于等于25000K，大于25000K；由计算机编程计算后获得从1000K到25000K色温的u、v值表，间隔为1K；
1.1.1黑体的u、v值计算包括以下步骤
(1)、由普朗克定律计算得到的黑体发射光的相对光谱功率：
其中：T为黑体的绝对温度(K)；λ为波长(nm)；c1为第一辐射常数，c1＝3.741844×10-12W·cm2；c2为第二辐射常数，c2＝1.438833cm·K；
(2)、由黑体相对光谱功率分布P(λ)，根据下式计算黑体的三刺激值X、Y、Z：
其中Δλ＝5nm；
(3)、由(3)式计算CIE1931-XYZ色品图中的色品坐标x、y，并由(4)式计算CIE1960UCS均匀色度坐标系中的色品坐标u、v值；
(4)、对于按上述过程计算黑体1000K到4999K色温的三刺激值X、Y、Z，x、y、u、v，间隔为1K，建立黑体从1000K到4999K色温的u、v值表，间隔为1K；
1.1.2标准照明体D的u、v值计算包括以下步骤
(1)、针对5000K到25000K色温，采用标准照明体D的色品坐标，当色温在5000K到6999K之间：
(2)、当色温在7000K到25000K之间：
(3)、针对5000K到25000K色温，采用(7)式计算u、v值：
(4)、建立标准照明体D从5000K到25000K色温的u、v值表，间隔为1K；
1.1.3存储u、v值表
嵌入式系统相关色温快速计算方法中的u、v值表，是由1.1.1建立的1000K至4999K色温的u、v值表和由1.1.2建立的5000K至25000K色温的u、v值表，再添加900K色温和26000K色温的u、v值合并完成；其中900K色温和26000K色温的u、v值用于色温区域划分；该表中的u、v值直接存储于便携式测量仪器中微处理器的闪存(Flash Memory)中；在针对被测光源的相关色温计算中，可以直接获取表中数据，加快运算速度；
1.2计算被测光源的相关色温
1.2.1计算被测光源的色品坐标
对于一个被测光源的色品坐标x0、y0，由(7)式计算u0、v0值；
1.2.2色温分区
(1)、分别计算被测光源的u0、v0和色温u、v值表中900K及1000K色温的u、v值距离平方E900K和E1000K：
ET＝(u0-uT)*(v0-vT) (8)
其中uT、vT分别表示某色温的u、v值；ET表示被测光源的u、v值和u、v值表中某色温的u、v值的距离平方，如果E900K<E1000K，则被测光源的相关色温设定为0；
(2)、由(8)式分别计算被测光源的u0、v0和色温的u、v值表中25000K及26000K色温的u、v值距离平方E25000K及E26000K；如果E26000K<E25000K，则被测光源的相关色温设定为0；
(3)、从色温u、v值表中获得5000K及5001K色温的u、v值：u5000K＝＝0.209144，v5000K＝0.325443，u5001K＝0.209133，v5001K＝0.325433；用式(8)计算被测光源的u0、v0和色温u、v值表中的5000K及5001K的u、v值距离平方E5000K及E5001K，如果E5000K<E5001K，设定为被测光源的相关色温处在低色温区域，否则设定为高色温区域；
1.2.3相关色温快速计算
采用McCamy近似公式计算被测光源的相关色温T的近似值：
根据下式计算该相关色温值的麦勒德值(mired)：
mired＝1000000/T (10)
根据所设定的高、低色温区域，低色温区域的相关色温最大值为4999K，高色温区域的相关色温最小值为5000K；在所述麦勒德值正负1个麦勒德的范围内，获取色温表中所有u、v值，并分别计算距离平方，求出最短距离，则该色温值为被测光源的相关色温。
本发明具有以下优点：
1、针对光源的相关色温测量，划定常用光源的相关色温区间，对于超出计算范围的以0K表示。在计算范围内，分别以黑体和标准照明D作为参考标准，并以5000K为参考标准的划分阈值，提高了相关色温计算精度；
2、由预处理提供的u、v表直接存储在Flash中，间隔为1K，无需便携式仪器重新计算，提高了计算速度并保证了计算精度；
3、采用距离平方，而不是距离计算，在保证相同效果的前提下，降低计算量；
4、综合比较了各种近似算法，采用近似精度较高且较为简洁的McCamy近似公式估算被测光源的相关色温，大致确定相关色温范围；而后采用等色温间隔的麦勒德值划定精确计算范围，最终采用全查表计算法获得最终结果。其计算速度极快且精度为1K，对便携式测量仪器的微处理器内存要求低，达到便携式测量仪器的速度及精度要求。
附图说明
图1为CIE1960UCS色度图。
具体实施方式
嵌入式系统相关色温快速计算方法如下：
1、针对被测光源的相关色温一般在1000K与25000K之间，将分为五种情形，分别为T<1000K，1000K≤T<5000K，5000K≤T<7000K，7000K≤T≤25000K，T>25000K；
2、计算黑体轨迹，由普朗克定律计算得到的黑体发射光的相对光谱功率：
其中：T为黑体的绝对温度(K)；λ为波长(nm)；c1为第一辐射常数，c1＝3.741844×10-12W·cm2；c2为第二辐射常数，c2＝1.438833cm·K；
3、由黑体相对光谱功率分布P(λ)，由下式计算黑体的三刺激值X、Y、Z：
其中Δλ＝5nm；
4、由(3)式计算CIE1931-XYZ色品图中的x、y，并由(4)式计算CIE1960UCS均匀色度坐标系中的色坐标u、v值；
5、按上述过程计算黑体1000K到4999K色温的三刺激值X、Y、Z，x、y、u、v，间隔为1K。建立黑体从1000K到4999K色温的uv值表，间隔为1K。
6、针对5000K到25000K色温，采用标准照明体D的色品坐标，当色温在5000K到6999K之间：
7、当色温在7000K到25000K之间：
8、计算uv值：
9、计算并建立标准照明体D从5000K到25000K色温的uv值表，间隔为1K；
10、按黑体相对光谱功率计算900K色温的uv值，并计算标准照明体D26000K色温的uv值，保存在色温uv表中；
以上为相关色温的预处理程序，由PC计算后获得uv值保存为色温的u、v值表，该表中的uv值直接存储于便携式测量仪器中微处理器的闪存中；在针对被测光源的相关色温计算中，可以直接获取表中数据，加快运算速度。表1、2、3显示了该u、v值表的一部分。
11、对于一个被测光源的色品坐标x0、y0，由(7)式计算u0、v0值；
12、分别计算u0、v0和表中900K及1000K色温的距离平方：
ET＝(u0-uT)*(v0-vT) (8)
如果E900K<E1000K，则被测光源的相关色温设定为0；
13、由(8)式分别计算u0、v0和表中25000K及26000K色温的距离平方。如果E26000K<E25000K，则被测光源的相关色温设定为0；
14、从uv表中获得5000K及5001K色温的uv值：u5000K＝＝0.209144，v5000K＝0.325443，u5001K＝0.209133，v5001K＝0.325433。计算E5000K及E5001K。如果E5000K<E5001K，设定为低色温区域，否则设定为高色温区域。
15、采用McCamy近似公式计算被测光源的相关色温近似值：
16、根据下式计算该相关色温值的麦勒德值：
mired＝1000000/T (10)
17、根据14中所设定的高低色温区域，低色温区域的最大值为4999K，高色温区域的最小值为5000K。在16中所计算的麦勒德值正负1个麦勒德的范围内，获取uv表中所有uv值，并分别计算距离平方，求出最短距离，则该色温值为被测光源的相关色温。
表1为u、v值表中色温900K和1000K-1500K的那部分表。
表2为u、v值表中色温4999K-5500K的那部分表。
表3为u、v值表中色温24499K-25000K和26000K的那部分表。
表1
表2
表3

Claims (1)

1.嵌入式系统相关色温快速计算方法，其特征在于它包括以下步骤：
1.1建立黑体和标准照明体D在不同色温时的u、v值表
以黑体和标准照明体D为标准光源，计算黑体和标准照明体D在不同色温时的u、v值，并保存为色温的u、v值表；计算过程中将色温分为五种情况，分别为小于1000K，大于等于1000K小于5000K，大于等于5000K小于7000K，大于等于7000K小于等于25000K，大于25000K；由计算机编程计算后获得从1000K到25000K色温的u、v值表，间隔为1K；
1.1.1黑体的u、v值计算包括以下步骤
(1)、由普朗克定律计算得到的黑体发射光的相对光谱功率：
<mrow> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>,</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>c</mi> <mn>1</mn> </msub> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>5</mn> </mrow> </msup> <msup> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mrow> <msub> <mi>c</mi> <mn>2</mn> </msub> <mo>/</mo> <mi>&amp;lambda;</mi> <mi>T</mi> </mrow> </msup> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中：T为黑体的绝对温度(K)；λ为波长(nm)；c1为第一辐射常数，c1＝3.741844×10- 12W·cm2；c2为第二辐射常数，c2＝1.438833cm·K；
(2)、由黑体相对光谱功率分布P(λ)，根据下式计算黑体的三刺激值X、Y、Z：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>X</mi> <mo>=</mo> <mo>&amp;Integral;</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;lambda;</mi> <mo>=</mo> <mo>&amp;Sigma;</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>&amp;Delta;</mi> <mi>&amp;lambda;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>Y</mi> <mo>=</mo> <mo>&amp;Integral;</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mover> <mi>y</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;lambda;</mi> <mo>=</mo> <mo>&amp;Sigma;</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mover> <mi>y</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>&amp;Delta;</mi> <mi>&amp;lambda;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>Z</mi> <mo>=</mo> <mo>&amp;Integral;</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mover> <mi>z</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;lambda;</mi> <mo>=</mo> <mo>&amp;Sigma;</mo> <mi>P</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mover> <mi>z</mi> <mo>&amp;OverBar;</mo> </mover> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>&amp;Delta;</mi> <mi>&amp;lambda;</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中Δλ＝5nm；
(3)、由(3)式计算CIE1931-XYZ色品图中的色品坐标x、y，并由(4)式计算CIE1960UCS均匀色度坐标系中的色品坐标u、v值；
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mi>X</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>+</mo> <mi>Y</mi> <mo>+</mo> <mi>Z</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mi>Y</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>+</mo> <mi>Y</mi> <mo>+</mo> <mi>Z</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>z</mi> <mo>=</mo> <mi>Z</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>+</mo> <mi>Y</mi> <mo>+</mo> <mi>Z</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>u</mi> <mo>=</mo> <mn>6</mn> <mi>X</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>+</mo> <mn>15</mn> <mi>Y</mi> <mo>+</mo> <mn>3</mn> <mi>Z</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>v</mi> <mo>=</mo> <mn>6</mn> <mi>Y</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>+</mo> <mn>15</mn> <mi>Y</mi> <mo>+</mo> <mn>3</mn> <mi>Z</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
(4)、对于按上述过程计算黑体1000K到4999K色温的三刺激值X、Y、Z，x、y、u、v，间隔为1K，建立黑体从1000K到4999K色温的u、v值表，间隔为1K；
1.1.2标准照明体D的u、v值计算包括以下步骤
(1)、针对5000K到25000K色温，采用标准照明体D的色品坐标，当色温在5000K到6999K之间：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mo>-</mo> <mn>4.6070</mn> <mfrac> <msup> <mn>10</mn> <mn>9</mn> </msup> <msup> <mi>T</mi> <mn>3</mn> </msup> </mfrac> <mo>+</mo> <mn>2.9678</mn> <mfrac> <msup> <mn>10</mn> <mn>6</mn> </msup> <msup> <mi>T</mi> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mn>0.09911</mn> <mfrac> <msup> <mn>10</mn> <mn>3</mn> </msup> <mi>T</mi> </mfrac> <mo>+</mo> <mn>0.244063</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mo>-</mo> <mn>3.000</mn> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <mn>2.870</mn> <mi>x</mi> <mo>-</mo> <mn>0.275</mn> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
(2)、当色温在7000K到25000K之间：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>x</mi> <mo>=</mo> <mo>-</mo> <mn>2.0064</mn> <mfrac> <msup> <mn>10</mn> <mn>9</mn> </msup> <msup> <mi>T</mi> <mn>3</mn> </msup> </mfrac> <mo>+</mo> <mn>1.9018</mn> <mfrac> <msup> <mn>10</mn> <mn>6</mn> </msup> <msup> <mi>T</mi> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mn>0.24748</mn> <mfrac> <msup> <mn>10</mn> <mn>3</mn> </msup> <mi>T</mi> </mfrac> <mo>+</mo> <mn>0.237040</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>y</mi> <mo>=</mo> <mo>-</mo> <mn>3.000</mn> <msup> <mi>x</mi> <mn>2</mn> </msup> <mo>+</mo> <mn>2.870</mn> <mi>x</mi> <mo>-</mo> <mn>0.275</mn> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
(3)、针对5000K到25000K色温，采用(7)式计算u、v值：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>u</mi> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mi>x</mi> </mrow> <mrow> <mo>-</mo> <mn>2</mn> <mi>x</mi> <mo>+</mo> <mn>12</mn> <mi>y</mi> <mo>+</mo> <mn>3</mn> </mrow> </mfrac> </mrow> </mtd> <mtd> <mrow> <mi>v</mi> <mo>=</mo> <mfrac> <mrow> <mn>6</mn> <mi>y</mi> </mrow> <mrow> <mo>-</mo> <mn>2</mn> <mi>x</mi> <mo>+</mo> <mn>12</mn> <mi>y</mi> <mo>+</mo> <mn>3</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
(4)、建立标准照明体D从5000K到25000K色温的u、v值表，间隔为1K；
1.1.3存储u、v值表
嵌入式系统相关色温快速计算方法中的u、v值表，是由1.1.1建立的1000K至4999K色温的u、v值表和由1.1.2建立的5000K至25000K色温的u、v值表，再添加900K色温和26000K色温的u、v值合并完成；其中900K色温和26000K色温的u、v值用于色温区域划分；该表中的u、v值直接存储于便携式测量仪器中微处理器的闪存(Flash Memory)中；在针对被测光源的相关色温计算中，可以直接获取表中数据，加快运算速度；
1.2计算被测光源的相关色温
1.2.1计算被测光源的色品坐标
对于一个被测光源的色品坐标x0、y0，由(7)式计算u0、v0值；
1.2.2色温分区
(1)、分别计算被测光源的u0、v0和色温u、v值表中900K及1000K色温的u、v值距离平方E900K和E1000K：
ET＝(u0-uT)*(v0-vT) (8)
其中uT、vT分别表示某色温的u、v值；ET表示被测光源的u、v值和u、v值表中某色温的u、v值的距离平方，如果E900K<E1000K，则被测光源的相关色温设定为0；
(2)、由(8)式分别计算被测光源的u0、v0和色温的u、v值表中25000K及26000K色温的u、v值距离平方E25000K及E26000K；如果E26000K<E25000K，则被测光源的相关色温设定为0；
(3)、从色温u、v值表中获得5000K及5001K色温的u、v值：u5000K＝0.209144，v5000K＝0.325443，u5001K＝0.209133，v5001K＝0.325433；用式(8)计算被测光源的u0、v0和色温u、v值表中的5000K及5001K的u、v值距离平方E5000K及E5001K，如果E5000K<E5001K，设定为被测光源的相关色温处在低色温区域，否则设定为高色温区域；
1.2.3相关色温快速计算
采用McCamy近似公式计算被测光源的相关色温T的近似值：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>T</mi> <mo>=</mo> <mo>-</mo> <mn>437</mn> <msup> <mi>n</mi> <mn>3</mn> </msup> <mo>+</mo> <mn>3601</mn> <msup> <mi>n</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>6861</mn> <mi>n</mi> <mo>+</mo> <mn>5514.31</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>n</mi> <mo>=</mo> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>0.3320</mn> <mo>)</mo> </mrow> <mo>/</mo> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <mn>0.1858</mn> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
根据下式计算该相关色温值的麦勒德值(mired)：
mired＝1000000/T (10)
根据所设定的高、低色温区域，低色温区域的相关色温最大值为4999K，高色温区域的相关色温最小值为5000K；在所述麦勒德值正负1个麦勒德的范围内，获取色温表中所有u、v值，并分别计算距离平方，求出最短距离，则该色温值为被测光源的相关色温。
CN201710072609.3A 2017-02-10 2017-02-10 嵌入式系统相关色温快速计算方法 Active CN106908150B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710072609.3A CN106908150B (zh) 2017-02-10 2017-02-10 嵌入式系统相关色温快速计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710072609.3A CN106908150B (zh) 2017-02-10 2017-02-10 嵌入式系统相关色温快速计算方法

Publications (2)

Publication Number Publication Date
CN106908150A CN106908150A (zh) 2017-06-30
CN106908150B true CN106908150B (zh) 2018-03-27

Family

ID=59209198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710072609.3A Active CN106908150B (zh) 2017-02-10 2017-02-10 嵌入式系统相关色温快速计算方法

Country Status (1)

Country Link
CN (1) CN106908150B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107368461B (zh) * 2017-07-11 2020-12-22 电子科技大学中山学院 光源相关色温和光源Duv的计算方法及系统
CN110060308B (zh) * 2019-03-28 2021-02-02 杭州电子科技大学 一种基于光源颜色分布限制的颜色恒常性方法
CN113324664A (zh) * 2021-04-20 2021-08-31 安徽曼德克环境科技有限公司 一种色温原理检测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105606230A (zh) * 2016-01-26 2016-05-25 浙江大学 一种精确获取光源相关色温的方法及系统
CN105934020A (zh) * 2016-04-27 2016-09-07 浙江大学 一种多色led匹配光谱和照度的方法
CN105973572A (zh) * 2016-04-27 2016-09-28 浙江大学 一种多色led实现光源最佳显色性的光谱优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105606230A (zh) * 2016-01-26 2016-05-25 浙江大学 一种精确获取光源相关色温的方法及系统
CN105934020A (zh) * 2016-04-27 2016-09-07 浙江大学 一种多色led匹配光谱和照度的方法
CN105973572A (zh) * 2016-04-27 2016-09-28 浙江大学 一种多色led实现光源最佳显色性的光谱优化方法

Also Published As

Publication number Publication date
CN106908150A (zh) 2017-06-30

Similar Documents

Publication Publication Date Title
CN106908150B (zh) 嵌入式系统相关色温快速计算方法
JP2014225463A5 (zh)
JP2014225464A5 (zh)
US9565734B1 (en) System and method for rapidly generating color models for LED-based lamps
US9144140B1 (en) System and method for controlling a plurality of light fixture outputs
ES2886112T3 (es) Procedimiento de determinación de la reflectancia de un objeto y dispositivo asociado
JP2013201132A5 (zh)
JP2013179060A5 (zh)
CN114746728A (zh) 环境光源分类
WO2021103954A1 (zh) 一种光源信息的确定方法、装置、电子设备和存储介质
CN112683338B (zh) 多参数同步测量方法、装置及系统
Wego Accuracy simulation of an LED based spectrophotometer
TW201421183A (zh) 最佳視覺照明裝置及方法
RU2017120364A (ru) Показатель для определения качества цвета
TWI476383B (zh) 以黑體輻射光譜相似性分析光源品質之方法
JP5016428B2 (ja) Ｌｅｄ照明装置
Burgos-Fernandez et al. Spectrally tunable light source based on light-emitting diodes for custom lighting solutions
CN113465742B (zh) 一种基于照度优化的白光光源光照颜色分辨能力量化方法及系统
JP5297821B2 (ja) 測色方法及び測色装置
CN115597726A (zh) 色温检测的方法、装置和电子设备
CN105509901B (zh) 发光体相对色温的高精度计量方法
Pavlova et al. Optimization of the methods for measuring color characteristics of light-emitting diodes in laboratory conditions
CN107368461B (zh) 光源相关色温和光源Duv的计算方法及系统
Durmus Optimizing a Three-channel sensor spectral sensitivity using a genetic algorithm
US20210037622A1 (en) Control design for perceptually uniform color tuning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240325

Address after: 310013 room 1801, Kangxin business building, No. 583, Wensan Road, Xihu District, Hangzhou City, Zhejiang Province

Patentee after: Hangzhou Xinye Optoelectronic Engineering Co.,Ltd.

Country or region after: China

Address before: 310030 room 1002, unit 1, building 26, Zijin Wenyuan, Xihu District, Hangzhou, Zhejiang

Patentee before: Ye Xinyi

Country or region before: China

Patentee before: Zhang Zhihao

TR01 Transfer of patent right