Correlated color temperature

Motivation

Black-body radiators are the reference by which the whiteness of light sources is judged. A black body is characterized by its temperature and emits light of a specific hue, which is referred to as color temperature. In practice, light sources that approximate Planckian radiators, such as certain fluorescent or high-intensity discharge lamps, are assessed based on their CCT, which is the temperature of a Planckian radiator whose color most closely resembles that of the light source. For light sources that do not follow the Planckian distribution, aligning them with a black body is not straightforward; thus, the concept of CCT is extended to represent these sources as accurately as possible on a one-dimensional color temperature scale, where "as accurately as possible" is determined within the framework of an objective color space.

Background

The notion of using Planckian radiators as a yardstick against which to judge other light sources is not new. In 1923, writing about "grading of illuminants with reference to quality of color ... the temperature of the source as an index of the quality of color", Priest essentially described CCT as we understand it today, going so far as to use the term "apparent color temperature", and astutely recognized three cases:

• "Those for which the spectral distribution of energy is identical with that given by the Planckian formula."
• "Those for which the spectral distribution of energy is not identical with that given by the Planckian formula, but still is of such a form that the quality of the color evoked is the same as would be evoked by the energy from a Planckian radiator at the given color temperature."
• "Those for which the spectral distribution of energy is such that the color can be matched only approximately by a stimulus of the Planckian form of spectral distribution."

Several important developments occurred in 1931. In chronological order:

1. Raymond Davis published a paper on "correlated color temperature" (his term). Referring to the Planckian locus on the r-g diagram, he defined the CCT as the average of the "primary component temperatures" (RGB CCTs), using trilinear coordinates.
2. The CIE announced the XYZ color space.
3. Deane B. Judd published a paper on the nature of "least perceptible differences" with respect to chromatic stimuli. By empirical means he determined that the difference in sensation, which he termed ΔE for a "discriminatory step between colors ... Empfindung" (German for sensation) was proportional to the distance of the colors on the chromaticity diagram. Referring to the (r,g) chromaticity diagram depicted aside, he hypothesized that

:::K ΔE = |c1 − c2| = max(|r1 − r2|, |g1 − g2|).

These developments paved the way for the development of new chromaticity spaces that are more suited to estimating correlated color temperatures and chromaticity differences. Bridging the concepts of color difference and color temperature, Priest made the observation that the eye is sensitive to constant differences in "reciprocal" temperature:

Priest proposed to use "the scale of temperature as a scale for arranging the chromaticities of the several illuminants in a serial order". Over the next few years, Judd published three more significant papers:

The first verified the findings of Priest,

The second proposed a new chromaticity space, guided by a principle that has become the holy grail of color spaces: perceptual uniformity (chromaticity distance should be commensurate with perceptual difference). By means of a projective transformation, Judd found a more "uniform chromaticity space" (UCS) in which to find the CCT. Judd determined the "nearest color temperature" by simply finding the point on the Planckian locus nearest to the chromaticity of the stimulus on Maxwell's color triangle, depicted aside. The transformation matrix he used to convert X,Y,Z tristimulus values to R,G,B coordinates was:

:\begin{bmatrix} R \ G \ B \end{bmatrix} = \begin{bmatrix} 3.1956 & 2.4478 & -0.1434 \ -2.5455 & 7.0492 & 0.9963 \ 0.0000 & 0.0000 & 1.0000 \end{bmatrix} \begin{bmatrix} X \ Y \ Z \end{bmatrix}.

From this, one can find these chromaticities:

:u=\frac{0.4661x+0.1593y}{y-0.15735x+0.2424}, \quad v=\frac{0.6581y}{y-0.15735x+0.2424}.

The third depicted the locus of the isothermal chromaticities on the CIE 1931 x,y chromaticity diagram. Since the isothermal points formed normals on his UCS diagram, transformation back into the xy plane revealed them still to be lines, but no longer perpendicular to the locus.

Calculation

Judd's idea of determining the nearest point to the Planckian locus on a uniform chromaticity space is current. In 1937, MacAdam suggested a "modified uniform chromaticity scale diagram", based on certain simplifying geometrical considerations:

:u = \frac{4x}{-2x+12y+3}, \quad v = \frac{6y}{-2x+12y+3}.

This (u,v) chromaticity space became the CIE 1960 color space, which is still used to calculate the CCT (even though MacAdam did not devise it with this purpose in mind). Using other chromaticity spaces, such as ''u'v''', leads to non-standard results that may nevertheless be perceptually meaningful.{{cite journal

The distance from the locus (i.e., degree of departure from a black body) is traditionally indicated in units of Δuv; positive for points above the locus. This concept of distance has evolved to become CIELAB ΔE*, which continues to be used today.

Robertson's method

Before the advent of powerful personal computers, it was common to estimate the correlated color temperature by way of interpolation from look-up tables and charts. who took advantage of the relatively even spacing of the mired scale (see above) to calculate the CCT Tc using linear interpolation of the isotherm's mired values:

:\frac{1}{T_c} = \frac{1}{T_i} + \frac{\theta_1}{\theta_1+\theta_2} \left( \frac{1}{T_{i+1}} - \frac{1}{T_i} \right),

where T_i and T_{i+1} are the color temperatures of the look-up isotherms and i is chosen such that T_i . (Furthermore, the test chromaticity lies between the only two adjacent lines for which d_i/d_{i+1} .)

If the isotherms are tight enough, one can assume \theta_1/\theta_2 \approx \sin \theta_1/\sin \theta_2, leading to

:\frac{1}{T_c} = \frac{1}{T_i} + \frac{d_i}{d_i-d_{i+1}} \left( \frac{1}{T_{i+1}} - \frac{1}{T_i} \right).

The distance of the test point to the i-th isotherm is given by

:d_i = \frac{ (v_T-v_i)-m_i (u_T-u_i) }{\sqrt {1+m_i^2}},

where (u_i, v_i) is the chromaticity coordinate of the i-th isotherm on the Planckian locus and mi is the isotherm's slope. Since it is perpendicular to the locus, it follows that m_i = -1/l_i where li is the slope of the locus at (u_i, v_i).

Precautions

Although the CCT can be calculated for any chromaticity coordinate, the result is meaningful only if the light source somewhat approximates a Planckian radiator. The CIE recommends that "The concept of correlated color temperature should not be used if the chromaticity of the test source differs more than Δuv = 5×10−2 from the Planckian radiator." Beyond a certain value of Δuv, a chromaticity co-ordinate may be equidistant to two points on the locus, causing ambiguity in the CCT.

Approximation

If a narrow range of color temperatures is considered—those encapsulating daylight being the most practical case—one can approximate the Planckian locus in order to calculate the CCT in terms of chromaticity coordinates. Following Kelly's observation that the isotherms intersect in the purple region near (x = 0.325, y = 0.154), McCamy proposed this cubic approximation:

:CCT(x, y) = 449 n^3 + 3525 n^2 + 6823.3 n + 5520.33,

where is the inverse slope line, and is the "epicenter"; quite close to the intersection point mentioned by Kelly. The maximum absolute error for color temperatures ranging from 2856 K (illuminant A) to 6504 K (D65) is under 2 K.

Hernández-André's 1999 proposal, using exponential terms, considerably extends the applicable range by adding a second epicenter for high color temperatures:

: CCT(x, y) = A_0 + A_1 \exp(-n/t_1) + A_2 \exp(-n/t_2) + A_3 \exp(-n/t_3)

where n is as before and the other constants are defined below:

The author suggests that one use the low-temperature equation to determine whether the higher-temperature parameters are needed.

Ohno (2013) proposes an accurate combined method based on a lookup table, a "parabolic" search, and a "triangular" search. The paper stresses the importance of also returning the Δuv value for evalulation of light sources. As it does not use one fixed table, it can be applied to any observer color matching function.

The inverse calculation, from color temperature to corresponding chromaticity coordinates, is discussed in .

References

References

1. [http://www.cie.co.at/publ/abst/17-4-89.html CIE/IEC 17.4:1987 International Lighting Vocabulary] {{Webarchive. link. (2010-02-27 ({{ISBN). 3900734070)
2. (December 2001). "The concept of correlated colour temperature revisited". Color Research & Application.
3. (19 March 2007). "The Color White".
4. (30 May 2024). "What is CCT?".
5. (12 July 2021). "Effects of illuminance and correlated color temperature of indoor light on emotion perception".
6. Hyde, Edward P.. (June 1911). "A New Determination of the Selective Radiation from Tantalum (abstract)". The American Physical Society.
7. Priest, Irwin G.. (1923). "The colorimetry and photometry of daylight ·and incandescent illuminants by the method of rotatory dispersion". [[JOSA]].
8. Davis, Raymond. (1931). "A Correlated Color Temperature for Illuminants". Bureau of Standards Journal of Research.
9. Judd, Deane B.. (1931). "Chromaticity sensibility to stimulus differences". [[JOSA]].
10. Priest, Irwin G.. (February 1933). "A proposed scale for use in specifying the chromaticity of incandescent illuminants and various phases of daylight". [[JOSA]].
11. Judd, Deane B.. (January 1933). "Sensibility to Color-Temperature Change as a Function of Temperature". [[JOSA]].
12. Judd, Deane B.. (January 1935). "A Maxwell Triangle Yielding Uniform Chromaticity Scales". [[JOSA]].
13. OSA Committee on Colorimetry. (November 1944). "Quantitative data and methods for colorimetry". [[JOSA]].
14. Judd, Deane B.. (November 1936). "Estimation of Chromaticity Differences and Nearest Color Temperatures on the Standard 1931 I.C.I. Colorimetric Coordinate System". [[JOSA]].
15. MacAdam, David L.. (August 1937). "Projective transformations of I.C.I. color specifications". [[JOSA]].
16. [http://www.delta.dk/C1256ED600446B80/sysOakFil/i102/$File/I102%20Correlated%20Colour%20Temperature.pdf The CIE definition of correlated color temperature (removed)] {{webarchive. link. (2009-02-05)
17. Robertson, Alan R.. (November 1968). "Computation of Correlated Color Temperature and Distribution Temperature". [[JOSA]].
18. [http://www.brucelindbloom.com/index.html?Eqn_XYZ_to_T.html ANSI C implementation] {{webarchive. link. (2008-04-22 , Bruce Lindbloom)
19. Walter, Wolfgang. (February 1992). "Determination of correlated color temperature based on a color-appearance model". Color Research & Application.
20. Schanda, János. (2007). "Colorimetry: Understanding the CIE System". [[Wiley Interscience]].
21. Kelly, Kenneth L.. (August 1963). "Lines of Constant Correlated Color Temperature Based on MacAdam's (u,v) Uniform Chromaticity Transformation of the CIE Diagram". JOSA.
22. McCamy, Calvin S.. (April 1992). "Correlated color temperature as an explicit function of chromaticity coordinates". Color Research & Application.
23. (September 20, 1999). "Calculating Correlated Color Temperatures Across the Entire Gamut of Daylight and Skylight Chromaticities". Applied Optics.
24. (2 January 2014). "Practical Use and Calculation of CCT and Duv". LEUKOS.
25. "colour.temperature.uv_to_CCT_Ohno2013 — Colour 0.4.4 documentation".