Correlated color temperature

Correlated color temperature (CCT, T_cp) refers to the "temperature of a Planckian radiator whose perceived color most closely resembles that of a given stimulus at the same brightness and under specified viewing conditions".^[1][2] The SI unit is the Kelvin (K).

Log-log graphs of peak emission wavelength and radiant exitance vs black-body temperature, plotted on the blue line. Red arrows show that 5780 K black bodies have 501 nm peak wavelength and 63.3 MW/m² radiant exitance.

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.^[3][4][5] 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

Judd's (r,g) diagram. The concentric curves indicate the loci of constant purity.

Judd's Maxwell triangle. Planckian locus in gray. Translating from trilinear co-ordinates into Cartesian co-ordinates leads to the next diagram.

Judd's uniform chromaticity space (UCS), with the Planckian locus and the isotherms from 1000 K to 10000 K, perpendicular to the locus. Judd calculated the isotherms in this space before translating them back into the (x,y) chromaticity space, as depicted in the diagram at the top of the article.

Close up of the Planckian locus in the CIE 1960 UCS, with the isotherms in mireds. Note the even spacing of the isotherms when using the reciprocal temperature scale and compare with the similar figure below. The even spacing of the isotherms on the locus implies that the mired scale is a better measure of perceptual color difference than the temperature scale.

The notion of using Planckian radiators as a yardstick against which to judge other light sources is not new.^[6] 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:^[7]

• "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.^[8]
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^[9]

K ΔE = |c₁ − c₂| = max(|r₁ − r₂|, |g₁ − g₂|).

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:^[10]

A difference of one micro-reciprocal-degree (μrd) is fairly representative of the doubtfully perceptible difference under the most favorable conditions of observation.

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,^[7] Davis,^[8] and Judd,^[9] with a paper on sensitivity to change in color temperature.^[11]

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:^[12]

${\displaystyle {\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:^[13]

${\displaystyle 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.^[14] 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

MacAdam's "uniform chromaticity scale" diagram; a simplification of Judd's UCS.

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:^[15]

${\displaystyle 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).^[16] Using other chromaticity spaces, such as u'v', leads to non-standard results that may nevertheless be perceptually meaningful.^[17]

Close up of the CIE 1960 UCS. The isotherms are perpendicular to the Planckian locus, and are drawn to indicate the maximum distance from the locus that the CIE considers the correlated color temperature to be meaningful: Δ_uv = ± 0.05

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.^[18] The most famous such method is Robertson's,^[19] who took advantage of the relatively even spacing of the mired scale (see above) to calculate the CCT T_c using linear interpolation of the isotherm's mired values:^[20]

Computation of the CCT T_c corresponding to the chromaticity coordinate ${\displaystyle \scriptstyle (u_{T},v_{T})}$ in the CIE 1960 UCS.

${\displaystyle {\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 ${\displaystyle T_{i}}$ and ${\displaystyle T_{i+1}}$ are the color temperatures of the look-up isotherms and i is chosen such that ${\displaystyle T_{i}<T_{c}<T_{i+1}}$. (Furthermore, the test chromaticity lies between the only two adjacent lines for which ${\displaystyle d_{i}/d_{i+1}<0}$.)

If the isotherms are tight enough, one can assume ${\displaystyle \theta _{1}/\theta _{2}\approx \sin \theta _{1}/\sin \theta _{2}}$, leading to

${\displaystyle {\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

${\displaystyle d_{i}={\frac {(v_{T}-v_{i})-m_{i}(u_{T}-u_{i})}{\sqrt {1+m_{i}^{2}}}},}$

where ${\displaystyle (u_{i},v_{i})}$ is the chromaticity coordinate of the i-th isotherm on the Planckian locus and m_i is the isotherm's slope. Since it is perpendicular to the locus, it follows that ${\displaystyle m_{i}=-1/l_{i}}$ where l_i is the slope of the locus at ${\displaystyle (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.^[21] 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⁻² from the Planckian radiator."^[22] 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),^[18] McCamy proposed this cubic approximation:^[23]

${\displaystyle CCT(x,y)=449n^{3}+3525n^{2}+6823.3n+5520.33,}$

where n = (x − x_e)/(y_e - y) is the inverse slope line, and (x_e = 0.3320, y_e = 0.1858) 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:^[24]

${\displaystyle 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:

	3–50 kK	50–800 kK
xe	0.3366	0.3356
ye	0.1735	0.1691
A0	−949.86315	36284.48953
A1	6253.80338	0.00228
t1	0.92159	0.07861
A2	28.70599	5.4535×10−36
t2	0.20039	0.01543
A3	0.00004
t3	0.07125

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.^[25] As it does not use one fixed table, it can be applied to any observer color matching function.^[26]

The inverse calculation, from color temperature to corresponding chromaticity coordinates, is discussed in Planckian locus § Approximation.