The derivation of the true colour of natural waters is based on the calculation of the Tristimulus values that are the three primaries (X, Y, Z) that specify a colour stimulus of the human eye [6 ,14 ]. Suppose the radiation spectrum that comes from the water is given by I that is a function of wavelength (λ), then the tristimulus values are given by:
The CIE 1931 standard colourimetric two degree Colour Matching Functions (CMFs) are presented by x (red), y (green) and z (blue). These serve as weighting functions for the determination of the tristimulus values. The intensity I can be replaced by I = E × R, the product of the illumination E times the remote sensing reflectance of water (R) [10]. For notation purposes we introduce the symbol T that represent the three tristimulus values (X, Y, Z) and that represents the three CMFs:
To further simplify the calculations the illumination E is taken as a constant, independent of wavelength, and the remote sensing reflectance is assumed to be corrected for the surface effects (Fresnel reflectance, foam, capillary waves). We refer to the standard books on water remote sensing by Mobley [11 ] and Kirk [10 ]. Thus RRS (λ) describes the intrinsic colour of the water, independent of air-water interface effects or illumination effects. Finally, because the integrals cannot be solved analytically, T can be written as the summation:
Note that the summation is taken between 400 and 710 nm. This will be discussed in more detail below. Also E is taken as unity and y is the product of the remote sensing reflectance times the CMF weighting functions:
Because ocean colour satellites do not provide full-spectral coverage, the y-spectrum must be first reconstructed by linear interpolation, based on the remote sensing reflection measured at the spectral bands (b). The contribution to T of a small interval of the spectrum between wavelengths L1 and L2 can be approximated by the trapezium rule (Figure 1 ):
To calculate y (Equation (4)), R(λ) = RRS(λ) must be retrieved from the values at the satellite bands b1 and b2. This can be done by linear interpolation at wavelength L1 and L2:
If Equations (4)–(6) are combined we find:
with:
Rewriting Equation (7) to an expression that is linear in the satellite bands Rb1 and Rb2 we find:
This implies that if we have the measured R values at b1 and b2, we can estimate ΔT between those bands (Equations (3), (8) and (9)) as a linear combination of those two, because for every wavelength interval we can calculate A, B and C and know at wavelengths L1 and L2 from [6 ]. Once the tristimulus values T (X, Y, Z) have been calculated, the three values are normalized and the colour is expressed in the coordinates:
The white point has the coordinates x = y = 1/3. In the (x, y) chromaticity plane, the coordinates are transformed to polar coordinates with respect to the white point and the hue angle is derived. The hue angle (α) lies between the vector to a point with coordinates (x − xw, y − yw) and the positive x-axis (at y − yw = 0), giving higher angles in an anti-clockwise direction (seeFigure 2 ).
All calculations in this paper were made with the ATAN2 function (four-quadrant inverse tangent) and the derived angles (in radians) are multiplied by 180/π to get the angles in degrees. In the development of the water hue angle concept, Wernand et al. [12 (link)] used (αM) for the hue angle derived with the FUME algorithm for MERIS, while Novoa et al. [15 (link)] introduced (αW) for the hue angle of water. In this manuscript we will refrain from indices and use (α) as the hue angle that represents the “true” or “intrinsic” colour of a natural water, which can be approximated by satellite remote sensing reflectance measurements.
The CIE 1931 standard colourimetric two degree Colour Matching Functions (CMFs) are presented by x (red), y (green) and z (blue). These serve as weighting functions for the determination of the tristimulus values. The intensity I can be replaced by I = E × R, the product of the illumination E times the remote sensing reflectance of water (R) [10]. For notation purposes we introduce the symbol T that represent the three tristimulus values (X, Y, Z) and that represents the three CMFs:
To further simplify the calculations the illumination E is taken as a constant, independent of wavelength, and the remote sensing reflectance is assumed to be corrected for the surface effects (Fresnel reflectance, foam, capillary waves). We refer to the standard books on water remote sensing by Mobley [11 ] and Kirk [10 ]. Thus RRS (λ) describes the intrinsic colour of the water, independent of air-water interface effects or illumination effects. Finally, because the integrals cannot be solved analytically, T can be written as the summation:
Note that the summation is taken between 400 and 710 nm. This will be discussed in more detail below. Also E is taken as unity and y is the product of the remote sensing reflectance times the CMF weighting functions:
Because ocean colour satellites do not provide full-spectral coverage, the y-spectrum must be first reconstructed by linear interpolation, based on the remote sensing reflection measured at the spectral bands (b). The contribution to T of a small interval of the spectrum between wavelengths L1 and L2 can be approximated by the trapezium rule (
To calculate y (Equation (4)), R(λ) = RRS(λ) must be retrieved from the values at the satellite bands b1 and b2. This can be done by linear interpolation at wavelength L1 and L2:
If Equations (4)–(6) are combined we find:
with:
Rewriting Equation (7) to an expression that is linear in the satellite bands Rb1 and Rb2 we find:
This implies that if we have the measured R values at b1 and b2, we can estimate ΔT between those bands (Equations (3), (8) and (9)) as a linear combination of those two, because for every wavelength interval we can calculate A, B and C and know at wavelengths L1 and L2 from [6 ]. Once the tristimulus values T (X, Y, Z) have been calculated, the three values are normalized and the colour is expressed in the coordinates:
The white point has the coordinates x = y = 1/3. In the (x, y) chromaticity plane, the coordinates are transformed to polar coordinates with respect to the white point and the hue angle is derived. The hue angle (α) lies between the vector to a point with coordinates (x − xw, y − yw) and the positive x-axis (at y − yw = 0), giving higher angles in an anti-clockwise direction (see
All calculations in this paper were made with the ATAN2 function (four-quadrant inverse tangent) and the derived angles (in radians) are multiplied by 180/π to get the angles in degrees. In the development of the water hue angle concept, Wernand et al. [12 (link)] used (αM) for the hue angle derived with the FUME algorithm for MERIS, while Novoa et al. [15 (link)] introduced (αW) for the hue angle of water. In this manuscript we will refrain from indices and use (α) as the hue angle that represents the “true” or “intrinsic” colour of a natural water, which can be approximated by satellite remote sensing reflectance measurements.
Full text: Click here