Practical Use and Calculation of CCT and Duv 2013
Practical Use and Calculation of CCT and Duv 2013
To cite this article: Yoshi Ohno (2014) Practical Use and Calculation of CCT and Duv, LEUKOS: The Journal of the Illuminating
Engineering Society of North America, 10:1, 47-55, DOI: 10.1080/15502724.2014.839020
Taylor & Francis makes every effort to ensure the accuracy of all the information (the “Content”) contained
in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no
representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the
Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and
are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and
should be independently verified with primary sources of information. Taylor and Francis shall not be liable for
any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever
or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of
the Content.
This article may be used for research, teaching, and private study purposes. Any substantial or systematic
reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any
form to anyone is expressly forbidden. Terms & Conditions of access and use can be found at http://
www.tandfonline.com/page/terms-and-conditions
LEUKOS, 10:47–55, 2014
ISSN: 1550-2724 print / 1550-2716 online
DOI: 10.1080/15502724.2014.839020
1. INTRODUCTION
The chromaticity is one of the critical parameters for light sources for general light-
ing and is normally specified with chromaticity coordinates CIE (x, y) or (u , v ).
However, these two numbers do not provide the color information intuitively. For
practical purposes, correlated color temperature (CCT) is commonly used to pro-
vide the chromaticity information of general illumination sources. CCT, however,
provides only one dimension of the chromaticity, and there is another dimension,
which is the position of chromaticity with respect to Planckian locus. For this pur-
pose, “duv” or similar terms as the distance from the Planckian locus have been
used in some parts of the industry, but these had not been officially defined in any
standards. Recently, Duv was defined in an American National Standards Institute
standard [ANSI 2008].
For color quality information of lighting products (for example, in data sheets),
only CCT and the color rendering index [CIE 1995] are often used, and Duv is
not shown though it is also critical. The chromaticity specification for fluorescent
lamps [ANSI 2001] specifying (x, y) coordinates has long been used, and that for
solid-state lighting products, specifying the ranges of Duv, is also available [ANSI
This article not subject to US copyright
law. 2008]. However, not all products in the market meet these requirements. It is a
Received 6 April 2013; revised 22 August
problem that some products do not meet the requirements in these standards because
2013; accepted 23 August 2013. inappropriate Duv may cause rejection by customers. It has also been anecdotally
Address correspondence to Yoshi reported that chromaticity points slightly below the Planckian locus are preferred.
Ohno, National Institute of Standards Duv is important for the users’ acceptance and preference of light sources. Related
and Technology, 100 Bureau Dr. MS
8442, Gaithersburg, MD 20899-8442. research on white point for light sources has been reported [Rea and Freyssinier-Nova
E-mail: [email protected] 2013].
47
A Duv value provides information on the distance
and direction of a color shift from the Planckian locus c1 λ−5 c
2
−1
Le,λ (λ, T ) = exp −1 (1)
(yellowish/greenish or pinkish), whereas the “duv” or other π λT
term often showed only the distance and no informa-
where L is the spectral radiance, T is the temperature (K),
tion on the direction of the shift. With the use of CCT
and λ is the wavelength in the medium. Note that the
and Duv, the two numbers can provide full information
refractive index of the medium, n, which is normally used
on white light chromaticity of light in an intuitive man-
in the Planck’s equation for radiometric applications, is not
ner, and these two numbers can replace (x, y) or (u ,
used for the calculation of CCT. c 1 , c 2 are the radiation
v ) in test reports and specification data. In this arti-
constants; the value of c 1 (3.741 7749 × 10−16 W·m2 ) is
cle, the use of CCT and Duv is proposed to express
not relevant to the value of CCT, and only c 2 is relevant.
chromaticity of white light sources for lighting, whereas
The value of c 2 used in the calculation of CCT should fol-
color differences should use distances in (u , v ) coordi-
low that adopted in the current International Temperature
nates. Because Duv is not widely known, its definition
Scale (ITS-90) [Preston-Thomas 1990], namely, c 2 =
together with that of CCT is first introduced in sec-
1.4388 × 10−2 m·K. See Appendix E of CIE 15:2004
Downloaded by [University of California, San Francisco] at 14:01 28 October 2014
48 Y. Ohno
Fig. 1 The ANSI specification [ANSI 2008] plotted on the (u, v ) diagram (left) and on a CCT-Duv chart (right) (color figure available
online).
This is an attempt to make such a chart for the whole (2) to calculate Duv. If the value of CCT is not available,
Downloaded by [University of California, San Francisco] at 14:01 28 October 2014
CCT range. In this case, a log scale is used for CCT to CCT needs to be calculated first from (x, y) or (u , v ) to
make the color space reasonably uniform, as shown in the use this method.
sizes of the seven-step MacAdam ellipses [DOE 2008] also
plotted. In this example, the CCT scale is reversed so
that the direction of the CCT scale is the same as that 4.2. Simple Approximation Formula
on the chromaticity diagrams. Charts with the CCT scale
This simple method can be used if the value of CCT is
not reversed can also be used. Such a chart can be used,
not available, and Duv needs to be calculated from (x,
for example, for quality control, by plotting (T cp , Duv ) of
y) or (u , v ) only and without using the full calculation
products measured.
procedures for CCT as described in section 5. First, the
A further consideration is that reporting Duv values
chromaticity coordinate of the light source (x, y) or
with many zeros after the decimal point is not convenient.
(u , v ) needs to be converted to the 1960
For more practical use of this quantity, a further idea would
(u, v)2 :
be a new quantity (for example, called “tint”) defined as
1000 times Duv, so that, for example, Duv −0.003 would
u = 4x/(−2x + 12y + 3) u = u
be “tint” −3.
or (3)
Section 4.2. However, if Duv needs to be calculated sepa- LFP = (u − 0.292)2 + (v − 0.240)2
rately from calculation of the CCT, the following methods
can be used. u − 0.292
a = arccos
LFP
If the value of CCT as well as (x, y) or (u , v ) of the light Duv = LFP − LBB
source is known, Duv can be calculated simply as follows. (4)
First, obtain the spectral distribution of the Planckian
radiation at temperature T of the CCT using (1). Then, where k6 = −0.00616793, k5 = 0.0893944,
from the Planckian spectrum, calculate the chromaticity k4 = −0.5179722, k3 = 1.5317403, k2 = −2.4243787,
coordinate (u0 , v0 ) using any program available. Then, use k1 = 1.925865, and k0 = −0.471106.
Fig. 2 Errors in Duv using the approximation polynomial (color figure available online).
Ti (K) ui vi di (example) i
Fig. 5 Errors in CCT [K] using the triangular solution with 1% CT-step table. Errors are in absolute value (color figure available online).
52 Y. Ohno
Fig. 6 Errors of the triangular solution at Duv = 0, −0.002, and +0.002 using a 1% CT-step Planckian (u, v ) table (a) before and (b) after
the correction is made (color figure available online).
d (T ) = a T 2 + b T + c, (9)
where
Fig. 7 Example of parabolic fit for d m-1 , d m , d m-1 (color figure
X = (Tm+1 − Tm ) (Tm−1 − Tm+1 ) (Tm − Tm−1 ) available online).
dm+1 (Tm − Tm−1 ) · Tm−1 · Tm ] · X −1 . SIGN (z) = 1 for z ≥ 0 and SIGN (x) = −1 for z < 0.
The CCT of the input light is determined as the color tem- vTx is the v coordinate of the Planck’s radiation at temper-
perature T = Tx where the parabolic function reaches the ature Tx and may be calculated from the Planck’s equation
minimum and is obtained as or using the formula for vTx in (8) if the triangular solution
is also calculated.
b The accuracy of this method has been evaluated for the
Tx = − ∵ d (T ) = 2aTx + b = 0 (10) example in Table 1 and the results are shown in Fig. 8.
2a
It is clearly shown that this method works very well for a
The same correction factor for the nonlinearity effect is wide range of Duv, but errors are very large around Duv
applied as = 0 (right on the Planckian locus). The reason for this
is explained in Fig. 9, which shows the case for CCT =
Tx, cor = Tx × 0.99991. 4000 K and Duv = 0. The parabolic minimum (solution)
Practical Use and Calculation of CCT and Duv 53
Fig. 8
Downloaded by [University of California, San Francisco] at 14:01 28 October 2014
Errors in CCT [K] using the parabolic solution with the 1% CT-step Planckian (u, v ) table. Errors are in absolute values (color
figure available online).
Fig. 9 The case of parabolic fit for 4000 K and Duv = 0. The dashed line shows plots of actual distances. In this case the distances do
not follow a parabolic curve (color figure available online).
does not reach distance 0, though it should. In this case, within 0.5 K in the 2000 to 20,000 K range and −0.03
the actual distances follow V-shaped linear lines (dashed < Duv < 0.03 and the maximum error is 0.8 K (which
line) rather than a parabolic curve, causing the error in occurs at conditions other than shown in Fig. 10).
the obtained CCT. This occurs only at or very near the
Planckian locus.
5.5. Further Improvements
The error analyses presented above are for the cases using
5.4. Combined Solution
the 1% CT-step Planckian (u, v) table shown in Table 1,
The parabolic solution seems to work accurately except on consisting of 303 rows of data. The accuracy obtained in
or near the Planckian locus. To solve this problem, the tri- this condition is considered to be sufficient for most prac-
angular solution and the parabolic solution are combined tical applications. However, when necessary, the accuracy
so that the results of the triangular solution are used for can be easily improved by using a larger table. For exam-
|Duv | < 0.002 and the parabolic solution is used for other ple, by using a 0.25% CT-step table (1210 rows), the errors
regions. The threshold was determined where the errors will be reduced to one quarter of the results shown in
of both methods are at similar level. Figure 10 shows the Fig. 10, less than 0.2 K in the same ranges of CCT and
results of this combined solution. The errors are mostly Duv. By using a table with much smaller CCT intervals,
54 Y. Ohno
Fig. 10 Errors in CCT [K] using the combined solution (1% CT-step table). Errors are in absolute value (color figure available online).
Downloaded by [University of California, San Francisco] at 14:01 28 October 2014