def test_luminous_efficacy(self):
"""
Tests :func:`colour.colorimetry.photometry.luminous_efficacy`
definition.
"""
self.assertAlmostEqual(
luminous_efficacy(
ILLUMINANTS_RELATIVE_SPDS.get('F2').clone().normalise()),
336.83937176,
places=7)
self.assertAlmostEqual(
luminous_efficacy(LIGHT_SOURCES_RELATIVE_SPDS.get(
'Neodimium Incandescent')),
136.21708032,
places=7)
self.assertAlmostEqual(
luminous_efficacy(LIGHT_SOURCES_RELATIVE_SPDS.get(
'F32T8/TL841 (Triphosphor)')),
348.88267549,
places=7)
spd = zeros_spd()
spd[555] = 1
self.assertAlmostEqual(
luminous_efficacy(spd),
683.00000000,
places=7)
def RGB_to_spectral_Smits1999(RGB):
"""
Recovers the spectral power distribution of given *RGB* colourspace matrix
using Smits (1999) method.
Parameters
----------
RGB : array_like, (3,)
*RGB* colourspace matrix.
Returns
-------
SpectralPowerDistribution
Recovered spectral power distribution.
Examples
--------
>>> RGB = np.array([0.02144962, 0.13154603, 0.09287601])
>>> RGB_to_spectral_Smits1999(RGB) # doctest: +ELLIPSIS
<...SpectralPowerDistribution object at 0x...>
"""
white_spd = SMITS_1999_SPDS.get('white').clone()
cyan_spd = SMITS_1999_SPDS.get('cyan').clone()
magenta_spd = SMITS_1999_SPDS.get('magenta').clone()
yellow_spd = SMITS_1999_SPDS.get('yellow').clone()
red_spd = SMITS_1999_SPDS.get('red').clone()
green_spd = SMITS_1999_SPDS.get('green').clone()
blue_spd = SMITS_1999_SPDS.get('blue').clone()
R, G, B = np.ravel(RGB)
spd = zeros_spd(SMITS_1999_SPDS.get('white').shape)
if R <= G and R <= B:
spd += white_spd * R
if G <= B:
spd += cyan_spd * (G - R)
spd += blue_spd * (B - G)
else:
spd += cyan_spd * (B - R)
spd += green_spd * (G - B)
elif G <= R and G <= B:
spd += white_spd * G
if R <= B:
spd += magenta_spd * (R - G)
spd += blue_spd * (B - R)
else:
spd += magenta_spd * (B - G)
spd += red_spd * (R - B)
else:
spd += white_spd * B
if R <= G:
spd += yellow_spd * (R - B)
spd += green_spd * (G - R)
else:
spd += yellow_spd * (G - B)
spd += red_spd * (R - G)
return spd
def RGB_to_spectral_Smits1999(RGB):
"""
Recovers the spectral power distribution of given *RGB* colourspace array
using Smits (1999) method.
Parameters
----------
RGB : array_like, (3,)
*RGB* colourspace array.
Returns
-------
SpectralPowerDistribution
Recovered spectral power distribution.
Examples
--------
>>> RGB = np.array([0.02144962, 0.13154603, 0.09287601])
>>> RGB_to_spectral_Smits1999(RGB) # doctest: +ELLIPSIS
<...SpectralPowerDistribution object at 0x...>
"""
white_spd = SMITS_1999_SPDS.get('white').clone()
cyan_spd = SMITS_1999_SPDS.get('cyan').clone()
magenta_spd = SMITS_1999_SPDS.get('magenta').clone()
yellow_spd = SMITS_1999_SPDS.get('yellow').clone()
red_spd = SMITS_1999_SPDS.get('red').clone()
green_spd = SMITS_1999_SPDS.get('green').clone()
blue_spd = SMITS_1999_SPDS.get('blue').clone()
R, G, B = np.ravel(RGB)
spd = zeros_spd(SMITS_1999_SPDS.get('white').shape)
if R <= G and R <= B:
spd += white_spd * R
if G <= B:
spd += cyan_spd * (G - R)
spd += blue_spd * (B - G)
else:
spd += cyan_spd * (B - R)
spd += green_spd * (G - B)
elif G <= R and G <= B:
spd += white_spd * G
if R <= B:
spd += magenta_spd * (R - G)
spd += blue_spd * (B - R)
else:
spd += magenta_spd * (B - G)
spd += red_spd * (R - B)
else:
spd += white_spd * B
if R <= G:
spd += yellow_spd * (R - B)
spd += green_spd * (G - R)
else:
spd += yellow_spd * (G - B)
spd += red_spd * (R - G)
return spd
def RGB_to_spectral_Smits1999(RGB):
"""
Recovers the spectral power distribution of given *RGB* colourspace array
using *Smits (1999)* method.
Parameters
----------
RGB : array_like, (3,)
*RGB* colourspace array.
Returns
-------
SpectralPowerDistribution
Recovered spectral power distribution.
Examples
--------
>>> RGB = np.array([0.02144962, 0.13154603, 0.09287601])
>>> print(RGB_to_spectral_Smits1999(RGB)) # doctest: +ELLIPSIS
SpectralPowerDistribution('0 Constant', (380.0, 720.0, 37.7777777...))
"""
white_spd = SMITS_1999_SPDS['white'].clone()
cyan_spd = SMITS_1999_SPDS['cyan'].clone()
magenta_spd = SMITS_1999_SPDS['magenta'].clone()
yellow_spd = SMITS_1999_SPDS['yellow'].clone()
red_spd = SMITS_1999_SPDS['red'].clone()
green_spd = SMITS_1999_SPDS['green'].clone()
blue_spd = SMITS_1999_SPDS['blue'].clone()
R, G, B = np.ravel(RGB)
spd = zeros_spd(SMITS_1999_SPDS['white'].shape)
if R <= G and R <= B:
spd += white_spd * R
if G <= B:
spd += cyan_spd * (G - R)
spd += blue_spd * (B - G)
else:
spd += cyan_spd * (B - R)
spd += green_spd * (G - B)
elif G <= R and G <= B:
spd += white_spd * G
if R <= B:
spd += magenta_spd * (R - G)
spd += blue_spd * (B - R)
else:
spd += magenta_spd * (B - G)
spd += red_spd * (R - B)
else:
spd += white_spd * B
if R <= G:
spd += yellow_spd * (R - B)
spd += green_spd * (G - R)
else:
spd += yellow_spd * (G - B)
spd += red_spd * (R - G)
return spd
def test_luminous_efficacy(self):
"""
Tests :func:`colour.colorimetry.photometry.luminous_efficacy`
definition.
"""
self.assertAlmostEqual(luminous_efficacy(
ILLUMINANTS_RELATIVE_SPDS.get('F2').clone().normalise()),
336.83937176,
places=7)
self.assertAlmostEqual(luminous_efficacy(
LIGHT_SOURCES_RELATIVE_SPDS.get('Neodimium Incandescent')),
136.21708032,
places=7)
self.assertAlmostEqual(luminous_efficacy(
LIGHT_SOURCES_RELATIVE_SPDS.get('F32T8/TL841 (Triphosphor)')),
348.88267549,
places=7)
spd = zeros_spd()
spd[555] = 1
self.assertAlmostEqual(luminous_efficacy(spd), 683.00000000, places=7)