CIECAM02

From Infogalactic: the planetary knowledge core
Jump to: navigation, search
Observing field model. Not drawn to scale.

In colorimetry, CIECAM02 is the color appearance model published in 2002 by the International Commission on Illumination (CIE) Technical Committee 8-01 (Color Appearance Modelling for Color Management Systems) and the successor of CIECAM97s.[1]

The two major parts of the model are its chromatic adaptation transform, CIECAT02, and its equations for calculating mathematical correlates for the six technically defined dimensions of color appearance: brightness (luminance), lightness, colorfulness, chroma, saturation, and hue.

Brightness is the subjective appearance of how bright an object appears given its surroundings and how it is illuminated. Lightness is the subjective appearance of how light a color appears to be. Colorfulness is the degree of difference between a color and grey. Chroma is the colorfulness relative to the brightness of another color that appears white under similar viewing conditions. This allows for the fact that a surface of a given chroma displays increasing colorfulness as the level of illumination increases. Saturation is the colorfulness of a color relative to its own brightness. Hue is the degree to which a stimulus can be described as similar to or different from stimuli that are described as red, green, blue, and yellow, the so-called unique hues. The colors that make up an object’s appearance are best described in terms of lightness and chroma when talking about the colors that make up the object’s surface, and it in terms of brightness, saturation and colorfulness when talking about the light that is emitted by or reflected off the object.

CIECAM02 takes for its input the tristimulus values of the stimulus, the tristimulus values of an adapting white point, adapting background, and surround luminance information, and whether or not observers are discounting the illuminant (color constancy is in effect). The model can be used to predict these appearance attributes or, with forward and reverse implementations for distinct viewing conditions, to compute corresponding colors.

CIECAM02 is used in Windows Vista's Windows Color System.[2]

Viewing conditions

The inner circle is the stimulus, from which the tristimulus values should be measured in CIE XYZ using the 2° standard observer. The intermediate circle is the proximal field, extending out another 2°. The outer circle is the background, reaching out to 10°, from which the relative luminance (Yb) need be measured. If the proximal field is the same color as the background, the background is considered to be adjacent to the stimulus. Beyond the circles which comprise the display field (display area, viewing area) is the surround field (or peripheral area), which can be considered to be the entire room. The totality of the proximal field, background, and surround is called the adapting field (the field of view that supports adaptation—extends to the limit of vision).[3]

When referring to the literature, it is also useful to be aware of the difference between the terms adopted white point (the computational white point) and the adapted white point (the observer white point).[4] The distinction may be important in mixed mode illumination, where psychophysical phenomena come into play. This is a subject of research.

Parameter decision table

CIECAM02 defines three surround(ing)s – average, dim, and dark – with associated parameters defined here for reference in the rest of this article:[5]

Surround
condition
Surround
ratio
F c Nc Application
Average SR > 0.2 1.0 0.69 1.0 Viewing surface colors
Dim 0 < SR < 0.2 0.9 0.59 0.95 Viewing television
Dark SR = 0 0.8 0.525 0.8 Using a projector in a dark room
  • SR = Lsw / Ldw: ratio of the absolute luminance of the reference white (white point) measured in the surround field to the display area. The 0.2 coefficient derives from the "gray world" assumption (~18%–20% reflectivity). It tests whether the surround luminance is darker or brighter than medium gray.
  • F: factor determining degree of adaptation
  • c: impact of surrounding
  • Nc: chromatic induction factor

For intermediate conditions, these values can be linearly interpolated.[5]

The absolute luminance of the adapting field, which is a quantity that will be needed later, should be measured with a photometer. If one is not available, it can be calculated using a reference white:


  L_A = \frac{E_w}{\pi} \frac{Y_b}{Y_w} = \frac{L_W Y_b}{Y_w}

where Yb is the relative luminance of background, the Ew = πLW is the illuminance of the reference white in lux, LW is the absolute luminance of the reference white in cd/m2, and Yw is the relative luminance of the reference white in the adapting field. If unknown, the adapting field can be assumed to have average reflectance ("gray world" assumption): LA = LW / 5.

Note: Care should be taken not to confuse LW, the absolute luminance of the reference white in cd/m2, and Lw the red cone response in the LMS color space.

Chromatic adaptation

Summary

  1. Convert to the "spectrally sharpened" CAT02 LMS space to prepare for adaptation. Spectral sharpening is the transformation of the tristimulus values into new values that would have resulted from a sharper, more concentrated set of spectral sensitivities. It is argued that this aids color constancy, especially in the blue region.(Compare Finlayson et al. 94, Spectral Sharpening:Sensor Transformations for Improved Colour Constancy)
  2. Perform chromatic adaptation using CAT02 (also known as the "modified CMCCAT2000 transform").
  3. Convert to an LMS space closer to the cone fundamentals. It is argued that predicting perceptual attribute correlates is best done in such spaces.[5]
  4. Perform post-adaptation cone response compression.

CAT02

Given a set of tristimulus values in XYZ, the corresponding LMS values can be determined by the MCAT02 transformation matrix (calculated using the CIE 1931 2° standard colorimetric observer).[1] The sample color in the test illuminant is:


  \begin{bmatrix}
    L\\
    M\\
    S
  \end{bmatrix}
  =
  \mathbf{M}_{CAT02}
  \begin{bmatrix}
    X\\
    Y\\
    Z
  \end{bmatrix},\quad
  \mathbf{M}_{CAT02}
  =
  \begin{bmatrix}
     \;\;\,0.7328 & 0.4296 & -0.1624\\
    -0.7036 & 1.6975 &  \;\;\,0.0061\\
     \;\;\,0.0030 & 0.0136 &  \;\;\,0.9834
  \end{bmatrix}

Once in LMS, the white point can be adapted to the desired degree by choosing the parameter D.[3] For the general CAT02, the corresponding color in the reference illuminant is:

\begin{align}
  L_c &= \Big(\frac{Y_w L_{wr}}{Y_{wr} L_w} D + 1-D\Big)L\\
  M_c &=\Big(\frac{Y_w M_{wr}}{Y_{wr} M_w} D + 1-D\Big)M\\
  S_c &= \Big(\frac{Y_w S_{wr}}{Y_{wr} S_w} D + 1-D\Big)S\\
\end{align}
Ciecam02 degree of adaptation.svg

where the Yw / Ywr factor accounts for the two illuminants having the same chromaticity but different reference whites.[6] The subscripts indicate the cone response for white under the test (w) and reference illuminant (wr). The degree of adaptation (discounting) D can be set to zero for no adaptation (stimulus is considered self-luminous) and unity for complete adaptation (color constancy). In practice, it ranges from 0.65 to 1.0, as can be seen from the diagram. Intermediate values can be calculated by:[5]


  D = F \left( 1 - \textstyle{\frac{1}{3.6}} e^{-(L_A + 42) / 92} \right)

where surround F is as defined above and LA is the adapting field luminance in cd/m2.[1]

log-log plot of FL vs. LA (LA ranges from 10−4 to 104, FL ranges from 10−4 to 10). The cube root approximation of FL is 0.1715LA1/3

In CIECAM02, the reference illuminant has equal energy Lwr = Mwr = Swr = 100) and the reference white is the perfect reflecting diffuser (i.e., unity reflectance, and Ywr = 100) hence:

\begin{align}
  L_c &= \Big(\frac{Y_w}{L_w} D + 1-D\Big)L\\
  M_c &=\Big(\frac{Y_w}{M_w} D + 1-D\Big)M\\
  S_c &= \Big(\frac{Y_w}{S_w} D + 1-D\Big)S\\
\end{align}

Furthermore, if the reference white in both illuminants have the Y tristimulus value (Ywr = Yw) then:

\begin{align}
  L_c &= \Big(\frac{L_{wr}}{L_w} D + 1-D\Big)L\\
  M_c &=\Big(\frac{M_{wr}}{M_w} D + 1-D\Big)M\\
  S_c &= \Big(\frac{S_{wr}}{S_w} D + 1-D\Big)S\\
\end{align}

Post-adaptation

After adaptation, the cone responses are converted to the Hunt–Pointer–Estévez space by going to XYZ and back:[5]


  \begin{bmatrix}
    L' \\
    M' \\
    S'
  \end{bmatrix}
  =
  \mathbf{M}_H
  \begin{bmatrix}
    X_c \\
    Y_c \\
    Z_c
  \end{bmatrix}
  =
  \mathbf{M}_H
  \mathbf{M}_{CAT02}^{-1}
  \begin{bmatrix}
    L_c \\
    M_c \\
    S_c
  \end{bmatrix}

  \mathbf{M}_H
  =
  \begin{bmatrix}
     \;\;\,0.38971 & 0.68898 & -0.07868 \\
    -0.22981 & 1.18340 &  \;\;\,0.04641 \\
     \;\;\,0.00000 & 0.00000 &  \;\;\,1.00000
  \end{bmatrix}
log La vs. log L′ for LA = 200 (FL = 1)

Finally, the response is compressed based on the generalized Michaelis–Menten equation (as depicted aside):[5]


  k = \frac{1}{5 L_A + 1}

  F_L = \textstyle{\frac{1}{5}} k^4 \left( 5 L_A \right) + \textstyle{\frac{1}{10}} {(1 - k^4)}^2 {\left( 5 L_A \right)}^{1/3}

FL is the luminance level adaptation factor.

\begin{align}
  L'_a &= \frac{400 {\left(F_L L'/100\right)}^{0.42}}{27.13 + {\left(F_L L'/100\right)}^{0.42}} + 0.1 \\
  M'_a &= \frac{400 {\left(F_L M'/100\right)}^{0.42}}{27.13 + {\left(F_L M'/100\right)}^{0.42}} + 0.1 \\
  S'_a &= \frac{400 {\left(F_L S'/100\right)}^{0.42}}{27.13 + {\left(F_L S'/100\right)}^{0.42}} + 0.1
\end{align}

As previously mentioned, if the luminance level of the background is unknown, it can estimated from the absolute luminance of the white point as LA = LW / 5 using the "medium gray" assumption. (The expression for FL is given in terms of 5LA for convenience.) In photopic conditions, the luminance level adaptation factor (FL) is proportional to the cube root of the luminance of the adapting field (LA). In scotopic conditions, it is proportional to LA (meaning no luminance level adaptation). The photopic threshold is roughly LW = 1 (see FLLA graph above).

Appearance correlates

CIECAM02 defines correlates for yellow-blue, red-green, brightness, and colorfulness. Let us make some preliminary definitions.

\begin{align}
  C_1 &= L^\prime_a - M^\prime_a \\
  C_2 &= M^\prime_a - S^\prime_a \\
  C_3 &= S^\prime_a - L^\prime_a
\end{align}

The correlate for red–green (a) is the magnitude of the departure of C1 from the criterion for unique yellow (C1 = C2 / 11), and the correlate for yellow–blue (b) is based on the mean of the magnitude of the departures of C1 from unique red (C1 = C2) and unique green (C1 = C3).[3]

\begin{align}
  a &= C_1 - \textstyle{\frac{1}{11}}C_2
    &= L^\prime_a - \textstyle{\frac{12}{11}} M^\prime_a + \textstyle{\frac{1}{11}} S^\prime_a \\
  b &= \textstyle{\frac{1}{2}} \left( C_2 - C_1 + C_1 - C_3 \right) / 4.5
    &= \textstyle{\frac{1}{9}} \left( L^\prime_a + M^\prime_a - 2S^\prime_a \right)
\end{align}

The 4.5 factor accounts for the fact that there are fewer cones at shorter wavelengths (the eye is less sensitive to blue). The order of the terms is such that b is positive for yellowish colors (rather than blueish).

The hue angle (h) can be found by converting the rectangular coordinate (a, b) into polar coordinates:


  h = \angle (a, b),\ (0 < h < 360^\circ)

To calculate the eccentricity (et) and hue composition (H), determine which quadrant the hue is in with the aid of the following table. Choose i such that hih′ < hi+1, where h′ = h if h > h1 and h′ = h + 360° otherwise.

Red Yellow Green Blue Red
i 1 2 3 4 5
hi 20.14 90.00 164.25 237.53 380.14
ei 0.8 0.7 1.0 1.2 0.8
Hi 0.0 100.0 200.0 300.0 400.0
\begin{align}
  H &= H_i + \frac{100 (h^\prime - h_i) / e_i}{(h^\prime - h_i) / e_i + (h_{i+1} - h^\prime) / e_{i+1}} \\
  e_t &= \textstyle{\frac{1}{4}} \left[ \cos\left( \textstyle{\frac{\pi}{180}}h + 2\right) + 3.8 \right]
\end{align}

(This is not exactly the same as the eccentricity factor given in the table.)

Calculate the achromatic response A:


  A = (2 L^\prime_a + M^\prime_a + \textstyle{\frac{1}{20}} S^\prime_a - 0.305) N_{bb}

where

\begin{align}
  &N_{bb} = N_{cb} = 0.725 n^{-0.2} \\
  &n = Y_b / Y_w
\end{align}

The correlate of lightness is


  J = 100 \left( A / A_w \right)^{c z}

where c is the impact of surround (see above), and


  z = 1.48 + \sqrt{n}

The correlate of brightness is


  Q = \left(4 / c \right) \sqrt{\textstyle{\frac{1}{100}} J} \left(A_w + 4\right) F_L^{1/4}

Then calculate a temporary quantity t,


  t = \frac{ \textstyle{\frac{50\,000}{13}} N_c N_{cb} e_t \sqrt{a^2+b^2} }
           { L_a^\prime + M_a^\prime + \textstyle{\frac{21}{20}} S_a^\prime }

The correlate of chroma is


  C = t^{0.9} \sqrt {\textstyle{\frac{1}{100}} J} (1.64 - 0.29^n)^{0.73}

The correlate of colorfulness is


  M = C \cdot F_L^{1/4}

The correlate of saturation is


  s = 100 \sqrt {M / Q}

References

  1. 1.0 1.1 1.2 Lua error in package.lua at line 80: module 'strict' not found.
  2. “Windows Color System: The Next Generation Color Management System”. Microsoft white paper. September 13, 2005.
  3. 3.0 3.1 3.2 Lua error in package.lua at line 80: module 'strict' not found.
  4. Lua error in package.lua at line 80: module 'strict' not found.
  5. 5.0 5.1 5.2 5.3 5.4 5.5 Lua error in package.lua at line 80: module 'strict' not found.
  6. Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.

External links