Correlated color temperature

The correlated color temperature (CCT, Tcp) is the temperature of the Planckian radiator whose perceived color most closely resembles that of a given stimulus at the same brightness and under specified viewing conditions

— CIE/IEC 17.4:1987, International Lighting Vocabulary (ISBN 3900734070)
Motivation
Black-body radiators are the reference by which the whiteness of light sources is judged. A black body can be described by its color temperature, whose hues are depicted above. By analogy, nearly Planckian light sources such as certain fluorescent or high-intensity discharge lamps can be judged by their correlated color temperature (CCT), the color temperature of the Planckian radiator that best approximates them. For light source spectra that are not Planckian, color temperature is not a well defined attribute; the concept of correlated color temperature was developed to map such sources as well as possible onto the one-dimensional scale of color temperature, where “as well as possible” is defined in the context 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:

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.
The CIE announced the XYZ color space.
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:

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, Davis, and Judd, with a paper on sensitivity to change in color temperature.

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:


From this, one can find these chromaticities:


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:


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.

The distance from the locus (i.e., degree of departure from a black body) is traditionally indicated in units of {\displaystyle \Delta uv} \Delta uv; positive for points above the locus. This concept of distance has evolved to become Delta 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. The most famous such method is Robertson’s, 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:


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


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


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 sources are nearly white. The CIE recommends that “The concept of correlated color temperature should not be used if the chromaticity of the test source differs more than  from the Planckian radiator.” Beyond a certain value of  , 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) = −449n3 + 3525n2 − 6823.3n + 5520.33,
where n = (x − xe)/(y − ye) is the inverse slope line, and (xe = 0.3320, ye = 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.

A more recent proposal, using exponential terms, considerably extends the applicable range by adding a second epicenter for high color temperatures:

CCT(x,y) = A0 + A1exp(−n/t1) + A2exp(−n/t2) + A3exp(−n/t3),
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 inverse calculation, from color temperature to corresponding chromaticity coordinates, is discussed in Planckian locus.

Source From Wikipedia

Tags: