def chromaticity_diagram_colours_CIE1976UCS(
samples=4096,
cmfs='CIE 1931 2 Degree Standard Observer',
antialiasing=True):
"""
Plots the *CIE 1976 UCS Chromaticity Diagram* colours.
Parameters
----------
samples : numeric, optional
Samples count on one axis.
cmfs : unicode, optional
Standard observer colour matching functions used for diagram bounds.
antialiasing : bool, optional
Whether to apply anti-aliasing to the image.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`colour.plotting.render`},
Please refer to the documentation of the previously listed definition.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> chromaticity_diagram_colours_CIE1976UCS() # doctest: +SKIP
"""
cmfs = get_cmfs(cmfs)
illuminant = DEFAULT_PLOTTING_ILLUMINANT
triangulation = Delaunay(Luv_to_uv(XYZ_to_Luv(cmfs.values, illuminant),
illuminant),
qhull_options='Qu QJ')
xx, yy = np.meshgrid(np.linspace(0, 1, samples),
np.linspace(1, 0, samples))
xy = tstack((xx, yy))
mask = (triangulation.find_simplex(xy) < 0).astype(DEFAULT_FLOAT_DTYPE)
if antialiasing:
kernel = np.array([
[0, 1, 0],
[1, 2, 1],
[0, 1, 0],
]).astype(DEFAULT_FLOAT_DTYPE)
kernel /= np.sum(kernel)
mask = convolve(mask, kernel)
mask = 1 - mask[:, :, np.newaxis]
XYZ = xy_to_XYZ(Luv_uv_to_xy(xy))
RGB = normalise_maximum(XYZ_to_sRGB(XYZ, illuminant), axis=-1)
return np.dstack([RGB, mask])
def test_XYZ_to_plotting_colourspace(self):
"""
Tests :func:`colour.plotting.common.XYZ_to_plotting_colourspace`
definition.
"""
XYZ = np.random.random(3)
np.testing.assert_almost_equal(
XYZ_to_sRGB(XYZ), XYZ_to_plotting_colourspace(XYZ), decimal=7)
def single_spd_plot(spd, cmfs='CIE 1931 2 Degree Standard Observer', **kwargs):
"""
Plots given spectral power distribution.
Parameters
----------
spd : SpectralPowerDistribution, optional
Spectral power distribution to plot.
cmfs : unicode
Standard observer colour matching functions used for spectrum creation.
\*\*kwargs : \*\*
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> from colour import SpectralPowerDistribution
>>> data = {400: 0.0641, 420: 0.0645, 440: 0.0562}
>>> spd = SpectralPowerDistribution('Custom', data)
>>> single_spd_plot(spd) # doctest: +SKIP
True
"""
cmfs, name = get_cmfs(cmfs), cmfs
shape = cmfs.shape
spd = spd.clone().interpolate(shape)
wavelengths = shape.range()
colours = []
y1 = []
for wavelength, value in spd:
XYZ = wavelength_to_XYZ(wavelength, cmfs)
colours.append(XYZ_to_sRGB(XYZ))
y1.append(value)
colours = normalise(colours)
settings = {
'title': '"{0}" - {1}'.format(spd.name, cmfs.name),
'x_label': u'Wavelength λ (nm)',
'y_label': 'Spectral Power Distribution',
'x_tighten': True,
'x_ticker': True,
'y_ticker': True
}
settings.update(kwargs)
return colour_parameters_plot([
colour_parameter(x=x[0], y1=x[1], RGB=x[2])
for x in tuple(zip(wavelengths, y1, colours))
], **settings)
def visible_spectrum_plot(cmfs='CIE 1931 2 Degree Standard Observer',
out_of_gamut_clipping=True,
**kwargs):
"""
Plots the visible colours spectrum using given standard observer *CIE XYZ*
colour matching functions.
Parameters
----------
cmfs : unicode, optional
Standard observer colour matching functions used for spectrum creation.
out_of_gamut_clipping : bool, optional
Out of gamut colours will be clipped if *True* otherwise, the colours
will be offset by the absolute minimal colour leading to a rendering on
gray background, less saturated and smoother. [1]_
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> visible_spectrum_plot() # doctest: +SKIP
"""
cmfs = get_cmfs(cmfs)
cmfs = cmfs.clone().align(DEFAULT_SPECTRAL_SHAPE)
wavelengths = cmfs.shape.range()
colours = XYZ_to_sRGB(
wavelength_to_XYZ(wavelengths, cmfs),
ILLUMINANTS['CIE 1931 2 Degree Standard Observer']['E'],
apply_encoding_cctf=False)
if not out_of_gamut_clipping:
colours += np.abs(np.min(colours))
colours = DEFAULT_PLOTTING_ENCODING_CCTF(normalise_maximum(colours))
settings = {
'title': 'The Visible Spectrum - {0}'.format(cmfs.title),
'x_label': 'Wavelength $\\lambda$ (nm)',
'y_label': False,
'x_tighten': True,
'y_ticker': False
}
settings.update(kwargs)
return colour_parameters_plot([
ColourParameter(x=x[0], RGB=x[1])
for x in tuple(zip(wavelengths, colours))
], **settings)
def blackbody_colours_plot(shape=SpectralShape(150, 12500, 50),
cmfs='CIE 1931 2 Degree Standard Observer',
**kwargs):
"""
Plots blackbody colours.
Parameters
----------
shape : SpectralShape, optional
Spectral shape to use as plot boundaries.
cmfs : unicode, optional
Standard observer colour matching functions.
\*\*kwargs : \*\*
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> blackbody_colours_plot() # doctest: +SKIP
True
"""
cmfs, name = get_cmfs(cmfs), cmfs
colours = []
temperatures = []
for temperature in shape:
spd = blackbody_spd(temperature, cmfs.shape)
XYZ = spectral_to_XYZ(spd, cmfs)
RGB = normalise(XYZ_to_sRGB(XYZ / 100))
colours.append(RGB)
temperatures.append(temperature)
settings = {
'title': 'Blackbody Colours',
'x_label': 'Temperature K',
'y_label': '',
'x_tighten': True,
'x_ticker': True,
'y_ticker': False
}
settings.update(kwargs)
return colour_parameters_plot([
colour_parameter(x=x[0], RGB=x[1])
for x in tuple(zip(temperatures, colours))
], **settings)
def visible_spectrum_plot(cmfs='CIE 1931 2 Degree Standard Observer',
**kwargs):
"""
Plots the visible colours spectrum using given standard observer *CIE XYZ*
colour matching functions.
Parameters
----------
cmfs : unicode, optional
Standard observer colour matching functions used for spectrum creation.
\*\*kwargs : \*\*
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> visible_spectrum_plot() # doctest: +SKIP
True
"""
cmfs, name = get_cmfs(cmfs), cmfs
cmfs = cmfs.clone().align(DEFAULT_SPECTRAL_SHAPE)
wavelengths = cmfs.shape.range()
colours = []
for i in wavelengths:
XYZ = wavelength_to_XYZ(i, cmfs)
colours.append(XYZ_to_sRGB(XYZ))
colours = np.array([np.ravel(x) for x in colours])
colours *= 1 / np.max(colours)
colours = np.clip(colours, 0, 1)
settings = {
'title': 'The Visible Spectrum - {0}'.format(name),
'x_label': u'Wavelength λ (nm)',
'x_tighten': True
}
settings.update(kwargs)
return colour_parameters_plot([
colour_parameter(x=x[0], RGB=x[1])
for x in tuple(zip(wavelengths, colours))
], **settings)
def test_XYZ_to_sRGB(self):
"""Test :func:`colour.models.rgb.common.XYZ_to_sRGB` definition."""
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.20654008, 0.12197225, 0.05136952])),
np.array([0.70573936, 0.19248266, 0.22354169]),
decimal=7,
)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.14222010, 0.23042768, 0.10495772])),
np.array([0.25847003, 0.58276102, 0.29718877]),
decimal=7,
)
np.testing.assert_almost_equal(
XYZ_to_sRGB(
np.array([0.07818780, 0.06157201, 0.28099326]),
np.array([0.34570, 0.35850]),
),
np.array([0.09838967, 0.25404426, 0.65130925]),
decimal=7,
)
np.testing.assert_almost_equal(
XYZ_to_sRGB(
np.array([0.00000000, 0.00000000, 0.00000000]),
np.array([0.44757, 0.40745]),
),
np.array([0.00000000, 0.00000000, 0.00000000]),
decimal=7,
)
np.testing.assert_almost_equal(
XYZ_to_sRGB(
np.array([0.20654008, 0.12197225, 0.05136952]),
np.array([0.44757, 0.40745]),
chromatic_adaptation_transform="Bradford",
),
np.array([0.60873814, 0.23259548, 0.43714892]),
decimal=7,
)
np.testing.assert_almost_equal(
XYZ_to_sRGB(
np.array([0.20654008, 0.12197225, 0.05136952]),
apply_cctf_encoding=False,
),
np.array([0.45620520, 0.03081070, 0.04091953]),
decimal=7,
)
def test_XYZ_to_sRGB(self):
"""
Tests :func:`colour.models.rgb.common.XYZ_to_sRGB` definition.
"""
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313])),
np.array([0.17498172, 0.38818743, 0.32159978]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.47097710, 0.34950000, 0.11301649])),
np.array([0.96979222, 0.48891569, 0.30231574]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.00000000, 0.00000000, 0.00000000]),
np.array([0.44757, 0.40745])),
np.array([0.00000000, 0.00000000, 0.00000000]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.11805834, 0.10340000, 0.05150892]),
np.array([0.31270, 0.32900])),
np.array([0.48222001, 0.31654775, 0.22070353]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313]),
chromatic_adaptation_transform='Bradford'),
np.array([0.17498172, 0.38818743, 0.32159978]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313]),
apply_encoding_cctf=False),
np.array([0.02583795, 0.12473949, 0.08439296]),
decimal=7)
def test_XYZ_to_sRGB(self):
"""
Tests :func:`colour.models.common.XYZ_to_sRGB` definition.
"""
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313])),
np.array([0.17501358, 0.38818795, 0.32161955]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.47097710, 0.34950000, 0.11301649])),
np.array([0.96984378, 0.48883420, 0.30229060]),
decimal=7)
np.testing.assert_almost_equal(XYZ_to_sRGB(
np.array([0.00000000, 0.00000000, 0.00000000]),
np.array([0.44757, 0.40745])),
np.array([0., 0., 0.]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.1180583421, 0.10340000, 0.0515089229]),
np.array([0.31271, 0.32902])),
np.array([0.48224885, 0.31651974, 0.22070513]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313]),
chromatic_adaptation_transform='Bradford'),
np.array([0.17501358, 0.38818795, 0.32161955]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313]),
apply_OECF=False),
np.array([0.02584654, 0.12473983, 0.0844036]),
decimal=7)
def test_XYZ_to_sRGB(self):
"""
Tests :func:`colour.models.rgb.common.XYZ_to_sRGB` definition.
"""
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313])),
np.array([0.17498817, 0.38819472, 0.32160312]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.47097710, 0.34950000, 0.11301649])),
np.array([0.96978949, 0.48894663, 0.30232185]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.00000000, 0.00000000, 0.00000000]),
np.array([0.44757, 0.40745])),
np.array([0.00000000, 0.00000000, 0.00000000]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.11805834, 0.10340000, 0.05150892]),
np.array([0.31270, 0.32900])),
np.array([0.48221920, 0.31656152, 0.22070695]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313]),
chromatic_adaptation_transform='Bradford'),
np.array([0.17498817, 0.38819472, 0.32160312]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.07049534, 0.10080000, 0.09558313]),
apply_encoding_cctf=False),
np.array([0.02583969, 0.12474440, 0.08439476]),
decimal=7)
def test_XYZ_to_RGB(self):
"""
Tests :func:`colour.models.common.XYZ_to_sRGB` definition.
"""
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.1180583421, 0.1034, 0.0515089229])),
np.array([0.48224885, 0.31651974, 0.22070513]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.21341854, 0.19387471, 0.16653476])),
np.array([0.59313312, 0.44141487, 0.42141429]),
decimal=7)
np.testing.assert_almost_equal(XYZ_to_sRGB(np.array([0, 0, 0]),
(0.34567, 0.35850)),
np.array([0., 0., 0.]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.1180583421, 0.1034, 0.0515089229]),
illuminant=(0.32168, 0.33767)),
np.array([0.47572655, 0.31766531, 0.23306098]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.1180583421, 0.1034, 0.0515089229]),
chromatic_adaptation_method='Bradford'),
np.array([0.48224885, 0.31651974, 0.22070513]),
decimal=7)
np.testing.assert_almost_equal(
XYZ_to_sRGB(np.array([0.1180583421, 0.1034, 0.0515089229]),
transfer_function=False),
np.array([0.19797725, 0.08168657, 0.03992654]),
decimal=7)
def XYZ_to_sd(
XYZ: ArrayLike,
method: Union[Literal["Jakob 2019", "Mallett 2019", "Meng 2015",
"Otsu 2018", "Smits 1999", ], str, ] = "Meng 2015",
**kwargs: Any,
) -> SpectralDistribution:
"""
Recover the spectral distribution of given *CIE XYZ* tristimulus
values using given method.
Parameters
----------
XYZ
*CIE XYZ* tristimulus values to recover the spectral distribution
from.
method
Computation method.
Other Parameters
----------------
additional_data
{:func:`colour.recovery.XYZ_to_sd_Jakob2019`},
If *True*, ``error`` will be returned alongside ``sd``.
basis_functions
{:func:`colour.recovery.RGB_to_sd_Mallett2019`},
Basis functions for the method. The default is to use the built-in
*sRGB* basis functions, i.e.
:attr:`colour.recovery.MSDS_BASIS_FUNCTIONS_sRGB_MALLETT2019`.
clip
{:func:`colour.recovery.XYZ_to_sd_Otsu2018`},
If *True*, the default, values below zero and above unity in the
recovered spectral distributions will be clipped. This ensures that the
returned reflectance is physical and conserves energy, but will cause
noticeable colour differences in case of very saturated colours.
cmfs
{:func:`colour.recovery.XYZ_to_sd_Meng2015`},
Standard observer colour matching functions.
colourspace
{:func:`colour.recovery.XYZ_to_sd_Jakob2019`},
*RGB* colourspace of the target colour. Note that no chromatic
adaptation is performed between ``illuminant`` and the colourspace
whitepoint.
dataset
{:func:`colour.recovery.XYZ_to_sd_Otsu2018`},
Dataset to use for reconstruction. The default is to use the published
data.
illuminant
{:func:`colour.recovery.XYZ_to_sd_Jakob2019`,
:func:`colour.recovery.XYZ_to_sd_Meng2015`},
Illuminant spectral distribution, default to
*CIE Standard Illuminant D65*.
interval
{:func:`colour.recovery.XYZ_to_sd_Meng2015`},
Wavelength :math:`\\lambda_{i}` range interval in nm. The smaller
``interval`` is, the longer the computations will be.
optimisation_kwargs
{:func:`colour.recovery.XYZ_to_sd_Jakob2019`,
:func:`colour.recovery.XYZ_to_sd_Meng2015`},
Parameters for :func:`scipy.optimize.minimize` and
:func:`colour.recovery.find_coefficients_Jakob2019` definitions.
Returns
-------
:class:`colour.SpectralDistribution`
Recovered spectral distribution.
Notes
-----
+------------+-----------------------+---------------+
| **Domain** | **Scale - Reference** | **Scale - 1** |
+============+=======================+===============+
| ``XYZ`` | [0, 1] | [0, 1] |
+------------+-----------------------+---------------+
- *Smits (1999)* method will internally convert given *CIE XYZ*
tristimulus values to *sRGB* colourspace array assuming equal energy
illuminant *E*.
References
----------
:cite:`Jakob2019`, :cite:`Mallett2019`, :cite:`Meng2015c`,
:cite:`Otsu2018`, :cite:`Smits1999a`
Examples
--------
*Jakob and Hanika (2009)* reflectance recovery:
>>> import numpy as np
>>> from colour import MSDS_CMFS, SDS_ILLUMINANTS, SpectralShape
>>> from colour.colorimetry import sd_to_XYZ_integration
>>> from colour.utilities import numpy_print_options
>>> XYZ = np.array([0.20654008, 0.12197225, 0.05136952])
>>> cmfs = (
... MSDS_CMFS['CIE 1931 2 Degree Standard Observer'].
... copy().align(SpectralShape(360, 780, 10))
... )
>>> illuminant = SDS_ILLUMINANTS['D65'].copy().align(cmfs.shape)
>>> sd = XYZ_to_sd(
... XYZ, method='Jakob 2019', cmfs=cmfs, illuminant=illuminant)
>>> with numpy_print_options(suppress=True):
... sd # doctest: +ELLIPSIS
SpectralDistribution([[ 360. , 0.4893773...],
[ 370. , 0.3258214...],
[ 380. , 0.2147792...],
[ 390. , 0.1482413...],
[ 400. , 0.1086169...],
[ 410. , 0.0841255...],
[ 420. , 0.0683114...],
[ 430. , 0.0577144...],
[ 440. , 0.0504267...],
[ 450. , 0.0453552...],
[ 460. , 0.0418520...],
[ 470. , 0.0395259...],
[ 480. , 0.0381430...],
[ 490. , 0.0375741...],
[ 500. , 0.0377685...],
[ 510. , 0.0387432...],
[ 520. , 0.0405871...],
[ 530. , 0.0434783...],
[ 540. , 0.0477225...],
[ 550. , 0.0538256...],
[ 560. , 0.0626314...],
[ 570. , 0.0755869...],
[ 580. , 0.0952675...],
[ 590. , 0.1264265...],
[ 600. , 0.1779272...],
[ 610. , 0.2649393...],
[ 620. , 0.4039779...],
[ 630. , 0.5832105...],
[ 640. , 0.7445440...],
[ 650. , 0.8499970...],
[ 660. , 0.9094792...],
[ 670. , 0.9425378...],
[ 680. , 0.9616376...],
[ 690. , 0.9732481...],
[ 700. , 0.9806562...],
[ 710. , 0.9855873...],
[ 720. , 0.9889903...],
[ 730. , 0.9914117...],
[ 740. , 0.9931801...],
[ 750. , 0.9945009...],
[ 760. , 0.9955066...],
[ 770. , 0.9962855...],
[ 780. , 0.9968976...]],
interpolator=SpragueInterpolator,
interpolator_kwargs={},
extrapolator=Extrapolator,
extrapolator_kwargs={...})
>>> sd_to_XYZ_integration(sd, cmfs, illuminant) / 100 # doctest: +ELLIPSIS
array([ 0.2066217..., 0.1220128..., 0.0513958...])
*Mallett and Yuksel (2019)* reflectance recovery:
>>> cmfs = (
... MSDS_CMFS['CIE 1931 2 Degree Standard Observer'].
... copy().align(SPECTRAL_SHAPE_sRGB_MALLETT2019)
... )
>>> illuminant = SDS_ILLUMINANTS['D65'].copy().align(cmfs.shape)
>>> sd = XYZ_to_sd(XYZ, method='Mallett 2019')
>>> with numpy_print_options(suppress=True):
... sd # doctest: +ELLIPSIS
SpectralDistribution([[ 380. , 0.1735531...],
[ 385. , 0.1720357...],
[ 390. , 0.1677721...],
[ 395. , 0.1576605...],
[ 400. , 0.1372829...],
[ 405. , 0.1170849...],
[ 410. , 0.0895694...],
[ 415. , 0.0706232...],
[ 420. , 0.0585765...],
[ 425. , 0.0523959...],
[ 430. , 0.0497598...],
[ 435. , 0.0476057...],
[ 440. , 0.0465079...],
[ 445. , 0.0460337...],
[ 450. , 0.0455839...],
[ 455. , 0.0452872...],
[ 460. , 0.0450981...],
[ 465. , 0.0448895...],
[ 470. , 0.0449257...],
[ 475. , 0.0448987...],
[ 480. , 0.0446834...],
[ 485. , 0.0441372...],
[ 490. , 0.0417137...],
[ 495. , 0.0373832...],
[ 500. , 0.0357657...],
[ 505. , 0.0348263...],
[ 510. , 0.0341953...],
[ 515. , 0.0337683...],
[ 520. , 0.0334979...],
[ 525. , 0.0332991...],
[ 530. , 0.0331909...],
[ 535. , 0.0332181...],
[ 540. , 0.0333387...],
[ 545. , 0.0334970...],
[ 550. , 0.0337381...],
[ 555. , 0.0341847...],
[ 560. , 0.0346447...],
[ 565. , 0.0353993...],
[ 570. , 0.0367367...],
[ 575. , 0.0392007...],
[ 580. , 0.0445902...],
[ 585. , 0.0625633...],
[ 590. , 0.2965381...],
[ 595. , 0.4215576...],
[ 600. , 0.4347139...],
[ 605. , 0.4385134...],
[ 610. , 0.4385184...],
[ 615. , 0.4385249...],
[ 620. , 0.4374694...],
[ 625. , 0.4384672...],
[ 630. , 0.4368251...],
[ 635. , 0.4340867...],
[ 640. , 0.4303219...],
[ 645. , 0.4243257...],
[ 650. , 0.4159482...],
[ 655. , 0.4057443...],
[ 660. , 0.3919874...],
[ 665. , 0.3742784...],
[ 670. , 0.3518421...],
[ 675. , 0.3240127...],
[ 680. , 0.2955145...],
[ 685. , 0.2625658...],
[ 690. , 0.2343423...],
[ 695. , 0.2174830...],
[ 700. , 0.2060461...],
[ 705. , 0.1977437...],
[ 710. , 0.1916846...],
[ 715. , 0.1861020...],
[ 720. , 0.1823908...],
[ 725. , 0.1807923...],
[ 730. , 0.1795571...],
[ 735. , 0.1785623...],
[ 740. , 0.1775758...],
[ 745. , 0.1771614...],
[ 750. , 0.1767431...],
[ 755. , 0.1764319...],
[ 760. , 0.1762597...],
[ 765. , 0.1762209...],
[ 770. , 0.1761803...],
[ 775. , 0.1761195...],
[ 780. , 0.1760763...]],
interpolator=SpragueInterpolator,
interpolator_kwargs={},
extrapolator=Extrapolator,
extrapolator_kwargs={...})
>>> sd_to_XYZ_integration(sd, cmfs, illuminant) / 100
... # doctest: +ELLIPSIS
array([ 0.2065436..., 0.1219996..., 0.0513764...])
*Meng (2015)* reflectance recovery:
>>> cmfs = (
... MSDS_CMFS['CIE 1931 2 Degree Standard Observer'].
... copy().align(SpectralShape(360, 780, 10))
... )
>>> illuminant = SDS_ILLUMINANTS['D65'].copy().align(cmfs.shape)
>>> sd = XYZ_to_sd(
... XYZ, method='Meng 2015', cmfs=cmfs, illuminant=illuminant)
>>> with numpy_print_options(suppress=True):
... sd # doctest: +SKIP
SpectralDistribution([[ 360. , 0.0762005...],
[ 370. , 0.0761792...],
[ 380. , 0.0761363...],
[ 390. , 0.0761194...],
[ 400. , 0.0762539...],
[ 410. , 0.0761671...],
[ 420. , 0.0754649...],
[ 430. , 0.0731519...],
[ 440. , 0.0676701...],
[ 450. , 0.0577800...],
[ 460. , 0.0441993...],
[ 470. , 0.0285064...],
[ 480. , 0.0138728...],
[ 490. , 0.0033585...],
[ 500. , 0. ...],
[ 510. , 0. ...],
[ 520. , 0. ...],
[ 530. , 0. ...],
[ 540. , 0.0055767...],
[ 550. , 0.0317581...],
[ 560. , 0.0754491...],
[ 570. , 0.1314115...],
[ 580. , 0.1937649...],
[ 590. , 0.2559311...],
[ 600. , 0.3123173...],
[ 610. , 0.3584966...],
[ 620. , 0.3927335...],
[ 630. , 0.4159458...],
[ 640. , 0.4306660...],
[ 650. , 0.4391040...],
[ 660. , 0.4439497...],
[ 670. , 0.4463618...],
[ 680. , 0.4474625...],
[ 690. , 0.4479868...],
[ 700. , 0.4482116...],
[ 710. , 0.4482800...],
[ 720. , 0.4483472...],
[ 730. , 0.4484251...],
[ 740. , 0.4484633...],
[ 750. , 0.4485071...],
[ 760. , 0.4484969...],
[ 770. , 0.4484853...],
[ 780. , 0.4485134...]],
interpolator=SpragueInterpolator,
interpolator_kwargs={},
extrapolator=Extrapolator,
extrapolator_kwargs={...})
>>> sd_to_XYZ_integration(sd, cmfs, illuminant) / 100 # doctest: +ELLIPSIS
array([ 0.2065400..., 0.1219722..., 0.0513695...])
*Otsu, Yamamoto and Hachisuka (2018)* reflectance recovery:
>>> cmfs = (
... MSDS_CMFS['CIE 1931 2 Degree Standard Observer'].
... copy().align(SPECTRAL_SHAPE_OTSU2018)
... )
>>> illuminant = SDS_ILLUMINANTS['D65'].copy().align(cmfs.shape)
>>> sd = XYZ_to_sd(
... XYZ, method='Otsu 2018', cmfs=cmfs, illuminant=illuminant)
>>> with numpy_print_options(suppress=True):
... sd # doctest: +ELLIPSIS
SpectralDistribution([[ 380. , 0.0601939...],
[ 390. , 0.0568063...],
[ 400. , 0.0517429...],
[ 410. , 0.0495841...],
[ 420. , 0.0502007...],
[ 430. , 0.0506489...],
[ 440. , 0.0510020...],
[ 450. , 0.0493782...],
[ 460. , 0.0468046...],
[ 470. , 0.0437132...],
[ 480. , 0.0416957...],
[ 490. , 0.0403783...],
[ 500. , 0.0405197...],
[ 510. , 0.0406031...],
[ 520. , 0.0416912...],
[ 530. , 0.0430956...],
[ 540. , 0.0444474...],
[ 550. , 0.0459336...],
[ 560. , 0.0507631...],
[ 570. , 0.0628967...],
[ 580. , 0.0844661...],
[ 590. , 0.1334277...],
[ 600. , 0.2262428...],
[ 610. , 0.3599330...],
[ 620. , 0.4885571...],
[ 630. , 0.5752546...],
[ 640. , 0.6193023...],
[ 650. , 0.6450744...],
[ 660. , 0.6610548...],
[ 670. , 0.6688673...],
[ 680. , 0.6795426...],
[ 690. , 0.6887933...],
[ 700. , 0.7003469...],
[ 710. , 0.7084128...],
[ 720. , 0.7154674...],
[ 730. , 0.7234334...]],
interpolator=SpragueInterpolator,
interpolator_kwargs={},
extrapolator=Extrapolator,
extrapolator_kwargs={...})
>>> sd_to_XYZ_integration(sd, cmfs, illuminant) / 100 # doctest: +ELLIPSIS
array([ 0.2065494..., 0.1219712..., 0.0514002...])
*Smits (1999)* reflectance recovery:
>>> cmfs = (
... MSDS_CMFS['CIE 1931 2 Degree Standard Observer'].
... copy().align(SpectralShape(360, 780, 10))
... )
>>> illuminant = SDS_ILLUMINANTS['E'].copy().align(cmfs.shape)
>>> sd = XYZ_to_sd(XYZ, method='Smits 1999')
>>> with numpy_print_options(suppress=True):
... sd # doctest: +ELLIPSIS
SpectralDistribution([[ 380. , 0.0787830...],
[ 417.7778 , 0.0622018...],
[ 455.5556 , 0.0446206...],
[ 493.3333 , 0.0352220...],
[ 531.1111 , 0.0324149...],
[ 568.8889 , 0.0330105...],
[ 606.6667 , 0.3207115...],
[ 644.4444 , 0.3836164...],
[ 682.2222 , 0.3836164...],
[ 720. , 0.3835649...]],
interpolator=LinearInterpolator,
interpolator_kwargs={},
extrapolator=Extrapolator,
extrapolator_kwargs={...})
>>> sd_to_XYZ_integration(sd, cmfs, illuminant) / 100 # doctest: +ELLIPSIS
array([ 0.1894770..., 0.1126470..., 0.0474420...])
"""
a = as_float_array(XYZ)
method = validate_method(method, XYZ_TO_SD_METHODS)
function = XYZ_TO_SD_METHODS[method]
if function is RGB_to_sd_Smits1999:
from colour.recovery.smits1999 import XYZ_to_RGB_Smits1999
a = XYZ_to_RGB_Smits1999(XYZ)
elif function is RGB_to_sd_Mallett2019:
from colour.models import XYZ_to_sRGB
a = XYZ_to_sRGB(XYZ, apply_cctf_encoding=False)
return function(a, **filter_kwargs(function, **kwargs))
def colour_checker_plot(colour_checker='ColorChecker 2005', **kwargs):
"""
Plots given colour checker.
Parameters
----------
colour_checker : unicode, optional
Color checker name.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`boundaries`, :func:`canvas`, :func:`decorate`,
:func:`display`},
Please refer to the documentation of the previously listed definitions.
width : numeric, optional
{:func:`multi_colour_plot`},
Colour polygon width.
height : numeric, optional
{:func:`multi_colour_plot`},
Colour polygon height.
spacing : numeric, optional
{:func:`multi_colour_plot`},
Colour polygons spacing.
across : int, optional
{:func:`multi_colour_plot`},
Colour polygons count per row.
text_display : bool, optional
{:func:`multi_colour_plot`},
Display colour text.
text_size : numeric, optional
{:func:`multi_colour_plot`},
Colour text size.
text_offset : numeric, optional
{:func:`multi_colour_plot`},
Colour text offset.
Returns
-------
Figure
Current figure or None.
Raises
------
KeyError
If the given colour rendition chart is not found in the factory colour
rendition charts.
Examples
--------
>>> colour_checker_plot() # doctest: +SKIP
"""
canvas(**kwargs)
colour_checker, name = COLOURCHECKERS.get(colour_checker), colour_checker
if colour_checker is None:
raise KeyError(('Colour checker "{0}" not found in '
'factory colour checkers: "{1}".').format(
name, sorted(COLOURCHECKERS.keys())))
_name, data, illuminant = colour_checker
colour_parameters = []
for _index, label, xyY in data:
XYZ = xyY_to_XYZ(xyY)
RGB = XYZ_to_sRGB(XYZ, illuminant)
colour_parameters.append(
ColourParameter(label.title(), np.clip(np.ravel(RGB), 0, 1)))
background_colour = '0.1'
width = height = 1.0
spacing = 0.25
across = 6
settings = {
'standalone': False,
'width': width,
'height': height,
'spacing': spacing,
'across': across,
'text_size': 8,
'background_colour': background_colour,
'margins': (-0.125, 0.125, -0.5, 0.125)
}
settings.update(kwargs)
multi_colour_plot(colour_parameters, **settings)
text_x = width * (across / 2) + (across * (spacing / 2)) - spacing / 2
text_y = -(len(colour_parameters) / across + spacing / 2)
pylab.text(
text_x,
text_y,
'{0} - {1} - Colour Rendition Chart'.format(
name, RGB_COLOURSPACES['sRGB'].name),
color='0.95',
clip_on=True,
ha='center')
settings.update({
'title': name,
'facecolor': background_colour,
'edgecolor': None,
'standalone': True
})
boundaries(**settings)
decorate(**settings)
return display(**settings)
def single_spd_plot(spd,
cmfs='CIE 1931 2 Degree Standard Observer',
out_of_gamut_clipping=True,
**kwargs):
"""
Plots given spectral power distribution.
Parameters
----------
spd : SpectralPowerDistribution
Spectral power distribution to plot.
out_of_gamut_clipping : bool, optional
Out of gamut colours will be clipped if *True* otherwise, the colours
will be offset by the absolute minimal colour leading to a rendering on
gray background, less saturated and smoother. [1]_
cmfs : unicode
Standard observer colour matching functions used for spectrum creation.
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> from colour import SpectralPowerDistribution
>>> data = {400: 0.0641, 420: 0.0645, 440: 0.0562}
>>> spd = SpectralPowerDistribution('Custom', data)
>>> single_spd_plot(spd) # doctest: +SKIP
"""
cmfs = get_cmfs(cmfs)
shape = cmfs.shape
spd = spd.clone().interpolate(shape, 'Linear')
wavelengths = spd.wavelengths
values = spd.values
y1 = values
colours = XYZ_to_sRGB(
wavelength_to_XYZ(wavelengths, cmfs),
ILLUMINANTS['CIE 1931 2 Degree Standard Observer']['E'],
apply_encoding_cctf=False)
if not out_of_gamut_clipping:
colours += np.abs(np.min(colours))
colours = DEFAULT_PLOTTING_ENCODING_CCTF(normalise_maximum(colours))
settings = {
'title': '{0} - {1}'.format(spd.title, cmfs.title),
'x_label': 'Wavelength $\\lambda$ (nm)',
'y_label': 'Spectral Power Distribution',
'x_tighten': True
}
settings.update(kwargs)
return colour_parameters_plot([
ColourParameter(x=x[0], y1=x[1], RGB=x[2])
for x in tuple(zip(wavelengths, y1, colours))
], **settings)
def blackbody_spectral_radiance_plot(
temperature=3500,
cmfs='CIE 1931 2 Degree Standard Observer',
blackbody='VY Canis Major',
**kwargs):
"""
Plots given blackbody spectral radiance.
Parameters
----------
temperature : numeric, optional
Blackbody temperature.
cmfs : unicode, optional
Standard observer colour matching functions.
blackbody : unicode, optional
Blackbody name.
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> blackbody_spectral_radiance_plot() # doctest: +SKIP
"""
canvas(**kwargs)
cmfs = get_cmfs(cmfs)
matplotlib.pyplot.subplots_adjust(hspace=0.4)
spd = blackbody_spd(temperature, cmfs.shape)
matplotlib.pyplot.figure(1)
matplotlib.pyplot.subplot(211)
settings = {
'title': '{0} - Spectral Radiance'.format(blackbody),
'y_label': 'W / (sr m$^2$) / m',
'standalone': False
}
settings.update(kwargs)
single_spd_plot(spd, cmfs.name, **settings)
XYZ = spectral_to_XYZ(spd, cmfs)
RGB = normalise_maximum(XYZ_to_sRGB(XYZ / 100))
matplotlib.pyplot.subplot(212)
settings = {
'title': '{0} - Colour'.format(blackbody),
'x_label': '{0}K'.format(temperature),
'y_label': '',
'aspect': None,
'standalone': False
}
single_colour_plot(ColourParameter(name='', RGB=RGB), **settings)
settings = {'standalone': True}
settings.update(kwargs)
boundaries(**settings)
decorate(**settings)
return display(**settings)
def blackbody_colours_plot(shape=SpectralShape(150, 12500, 50),
cmfs='CIE 1931 2 Degree Standard Observer',
**kwargs):
"""
Plots blackbody colours.
Parameters
----------
shape : SpectralShape, optional
Spectral shape to use as plot boundaries.
cmfs : unicode, optional
Standard observer colour matching functions.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`colour.plotting.render`},
Please refer to the documentation of the previously listed definition.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> blackbody_colours_plot() # doctest: +SKIP
"""
axes = canvas(**kwargs).gca()
cmfs = get_cmfs(cmfs)
colours = []
temperatures = []
with suppress_warnings():
for temperature in shape:
spd = blackbody_spd(temperature, cmfs.shape)
XYZ = spectral_to_XYZ(spd, cmfs)
RGB = normalise_maximum(XYZ_to_sRGB(XYZ / 100))
colours.append(RGB)
temperatures.append(temperature)
x_min, x_max = min(temperatures), max(temperatures)
y_min, y_max = 0, 1
axes.bar(x=temperatures,
height=1,
width=shape.interval,
color=colours,
align='edge')
settings = {
'title': 'Blackbody Colours',
'x_label': 'Temperature K',
'y_label': None,
'limits': (x_min, x_max, y_min, y_max),
'x_tighten': True,
'y_tighten': True,
'y_ticker': False
}
settings.update(kwargs)
return render(**settings)
def multi_spd_plot(spds,
cmfs='CIE 1931 2 Degree Standard Observer',
use_spds_colours=False,
normalise_spds_colours=False,
**kwargs):
"""
Plots given spectral power distributions.
Parameters
----------
spds : list, optional
Spectral power distributions to plot.
cmfs : unicode, optional
Standard observer colour matching functions used for spectrum creation.
use_spds_colours : bool, optional
Use spectral power distributions colours.
normalise_spds_colours : bool
Should spectral power distributions colours normalised.
\*\*kwargs : \*\*
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> from colour import SpectralPowerDistribution
>>> data1 = {400: 0.0641, 420: 0.0645, 440: 0.0562}
>>> data2 = {400: 0.134, 420: 0.789, 440: 1.289}
>>> spd1 = SpectralPowerDistribution('Custom1', data1)
>>> spd2 = SpectralPowerDistribution('Custom2', data2)
>>> multi_spd_plot([spd1, spd2]) # doctest: +SKIP
True
"""
cmfs, name = get_cmfs(cmfs), cmfs
if use_spds_colours:
illuminant = ILLUMINANTS_RELATIVE_SPDS.get('D65')
x_limit_min, x_limit_max, y_limit_min, y_limit_max = [], [], [], []
for spd in spds:
wavelengths, values = tuple(zip(*[(key, value) for key, value in spd]))
shape = spd.shape
x_limit_min.append(shape.start)
x_limit_max.append(shape.end)
y_limit_min.append(min(values))
y_limit_max.append(max(values))
matplotlib.pyplot.rc("axes", color_cycle=["r", "g", "b", "y"])
if use_spds_colours:
XYZ = spectral_to_XYZ(spd, cmfs, illuminant) / 100
if normalise_spds_colours:
XYZ = normalise(XYZ, clip=False)
RGB = np.clip(XYZ_to_sRGB(XYZ), 0, 1)
pylab.plot(wavelengths,
values,
color=RGB,
label=spd.name,
linewidth=2)
else:
pylab.plot(wavelengths, values, label=spd.name, linewidth=2)
settings = {
'x_label':
u'Wavelength λ (nm)',
'y_label':
'Spectral Power Distribution',
'x_tighten':
True,
'legend':
True,
'legend_location':
'upper left',
'x_ticker':
True,
'y_ticker':
True,
'limits': [
min(x_limit_min),
max(x_limit_max),
min(y_limit_min),
max(y_limit_max)
]
}
settings.update(kwargs)
bounding_box(**settings)
aspect(**settings)
return display(**settings)
def the_blue_sky_plot(cmfs='CIE 1931 2 Degree Standard Observer', **kwargs):
"""
Plots the blue sky.
Parameters
----------
cmfs : unicode, optional
Standard observer colour matching functions.
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> the_blue_sky_plot() # doctest: +SKIP
True
"""
canvas(**kwargs)
cmfs, name = get_cmfs(cmfs), cmfs
ASTM_G_173_spd = ASTM_G_173_ETR.clone()
rayleigh_spd = rayleigh_scattering_spd()
ASTM_G_173_spd.align(rayleigh_spd.shape)
spd = rayleigh_spd * ASTM_G_173_spd
matplotlib.pyplot.subplots_adjust(hspace=0.4)
matplotlib.pyplot.figure(1)
matplotlib.pyplot.subplot(211)
settings = {
'title': 'The Blue Sky - Synthetic Spectral Power Distribution',
'y_label': u'W / m-2 / nm-1',
'standalone': False
}
settings.update(kwargs)
single_spd_plot(spd, name, **settings)
matplotlib.pyplot.subplot(212)
settings = {
'title':
'The Blue Sky - Colour',
'x_label': ('The sky is blue because molecules in the atmosphere '
'scatter shorter wavelengths more than longer ones.\n'
'The synthetic spectral power distribution is computed as '
'follows: '
'(ASTM G-173 ETR * Standard Air Rayleigh Scattering).'),
'y_label':
'',
'aspect':
None,
'standalone':
False
}
blue_sky_color = XYZ_to_sRGB(spectral_to_XYZ(spd))
single_colour_plot(ColourParameter('', normalise(blue_sky_color)),
**settings)
settings = {'standalone': True}
settings.update(kwargs)
boundaries(**settings)
decorate(**settings)
return display(**settings)
"""
Compute sRGB colour values of CIE 224:2017's 99 test color samples under D65.
"""
from colour.colorimetry import (SpectralShape, SDS_ILLUMINANTS, sd_to_XYZ)
from colour.models import XYZ_to_sRGB
from colour.notation import RGB_to_HEX
from colour.quality.cfi2017 import get_tcs_CFI2017
if __name__ == '__main__':
sds = get_tcs_CFI2017(SpectralShape(380, 780, 5)).to_sds()
D65 = SDS_ILLUMINANTS['D65']
for sd in sds:
XYZ = sd_to_XYZ(sd, illuminant=D65) / 100
RGB = XYZ_to_sRGB(XYZ, cctf_encoding=True)
HEX = RGB_to_HEX(RGB)
print('\'%s\', ' % HEX, end='')
print('')
def the_blue_sky_plot(cmfs='CIE 1931 2 Degree Standard Observer', **kwargs):
"""
Plots the blue sky.
Parameters
----------
cmfs : unicode, optional
Standard observer colour matching functions.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`colour.plotting.render`},
Please refer to the documentation of the previously listed definition.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> the_blue_sky_plot() # doctest: +SKIP
"""
canvas(**kwargs)
cmfs, name = get_cmfs(cmfs), cmfs
ASTM_G_173_spd = ASTM_G_173_ETR.copy()
rayleigh_spd = rayleigh_scattering_spd()
ASTM_G_173_spd.align(rayleigh_spd.shape)
spd = rayleigh_spd * ASTM_G_173_spd
matplotlib.pyplot.subplots_adjust(hspace=0.4)
matplotlib.pyplot.figure(1)
matplotlib.pyplot.subplot(211)
settings = {
'title': 'The Blue Sky - Synthetic Spectral Power Distribution',
'y_label': u'W / m-2 / nm-1',
'standalone': False
}
settings.update(kwargs)
single_spd_plot(spd, name, **settings)
matplotlib.pyplot.subplot(212)
settings = {
'title':
'The Blue Sky - Colour',
'x_label': ('The sky is blue because molecules in the atmosphere '
'scatter shorter wavelengths more than longer ones.\n'
'The synthetic spectral power distribution is computed as '
'follows: '
'(ASTM G-173 ETR * Standard Air Rayleigh Scattering).'),
'y_label':
'',
'aspect':
None,
'standalone':
False
}
blue_sky_color = XYZ_to_sRGB(spectral_to_XYZ(spd))
single_colour_swatch_plot(
ColourSwatch('', normalise_maximum(blue_sky_color)), **settings)
settings = {'standalone': True}
settings.update(kwargs)
return render(**settings)
def generate_documentation_plots(output_directory: str):
"""
Generate documentation plots.
Parameters
----------
output_directory
Output directory.
"""
filter_warnings()
colour_style()
np.random.seed(0)
# *************************************************************************
# "README.rst"
# *************************************************************************
filename = os.path.join(
output_directory, "Examples_Colour_Automatic_Conversion_Graph.png"
)
plot_automatic_colour_conversion_graph(filename)
arguments = {
"tight_layout": True,
"transparent_background": True,
"filename": os.path.join(
output_directory, "Examples_Plotting_Visible_Spectrum.png"
),
}
plt.close(
plot_visible_spectrum(
"CIE 1931 2 Degree Standard Observer", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Illuminant_F1_SD.png"
)
plt.close(plot_single_illuminant_sd("FL1", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Blackbodies.png"
)
blackbody_sds = [
sd_blackbody(i, SpectralShape(0, 10000, 10))
for i in range(1000, 15000, 1000)
]
plt.close(
plot_multi_sds(
blackbody_sds,
y_label="W / (sr m$^2$) / m",
plot_kwargs={"use_sd_colours": True, "normalise_sd_colours": True},
legend_location="upper right",
bounding_box=(0, 1250, 0, 2.5e6),
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Cone_Fundamentals.png"
)
plt.close(
plot_single_cmfs(
"Stockman & Sharpe 2 Degree Cone Fundamentals",
y_label="Sensitivity",
bounding_box=(390, 870, 0, 1.1),
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Luminous_Efficiency.png"
)
plt.close(
plot_multi_sds(
(
sd_mesopic_luminous_efficiency_function(0.2),
SDS_LEFS_PHOTOPIC["CIE 1924 Photopic Standard Observer"],
SDS_LEFS_SCOTOPIC["CIE 1951 Scotopic Standard Observer"],
),
y_label="Luminous Efficiency",
legend_location="upper right",
y_tighten=True,
margins=(0, 0, 0, 0.1),
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_BabelColor_Average.png"
)
plt.close(
plot_multi_sds(
SDS_COLOURCHECKERS["BabelColor Average"].values(),
plot_kwargs={"use_sd_colours": True},
title=("BabelColor Average - " "Spectral Distributions"),
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_ColorChecker_2005.png"
)
plt.close(
plot_single_colour_checker(
"ColorChecker 2005", text_kwargs={"visible": False}, **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Chromaticities_Prediction.png"
)
plt.close(
plot_corresponding_chromaticities_prediction(
2, "Von Kries", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Examples_Plotting_Chromaticities_CIE_1931_Chromaticity_Diagram.png",
)
RGB = np.random.random((32, 32, 3))
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1931(
RGB,
"ITU-R BT.709",
colourspaces=["ACEScg", "S-Gamut"],
show_pointer_gamut=True,
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_CRI.png"
)
plt.close(
plot_single_sd_colour_rendering_index_bars(
SDS_ILLUMINANTS["FL2"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Colour_Rendition_Report.png"
)
plt.close(
plot_single_sd_colour_rendition_report(
SDS_ILLUMINANTS["FL2"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Plot_Visible_Spectrum_Section.png"
)
plt.close(
plot_visible_spectrum_section(
section_colours="RGB", section_opacity=0.15, **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Examples_Plotting_Plot_RGB_Colourspace_Section.png"
)
plt.close(
plot_RGB_colourspace_section(
"sRGB", section_colours="RGB", section_opacity=0.15, **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Examples_Plotting_CCT_CIE_1960_UCS_Chromaticity_Diagram.png",
)
plt.close(
plot_planckian_locus_in_chromaticity_diagram_CIE1960UCS(
["A", "B", "C"], **arguments
)[0]
)
# *************************************************************************
# Documentation
# *************************************************************************
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_CVD_Simulation_Machado2009.png"
)
plt.close(plot_cvd_simulation_Machado2009(RGB, **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Colour_Checker.png"
)
plt.close(plot_single_colour_checker("ColorChecker 2005", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Colour_Checkers.png"
)
plt.close(
plot_multi_colour_checkers(
["ColorChecker 1976", "ColorChecker 2005"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_SD.png"
)
data = {
500: 0.0651,
520: 0.0705,
540: 0.0772,
560: 0.0870,
580: 0.1128,
600: 0.1360,
}
sd = SpectralDistribution(data, name="Custom")
plt.close(plot_single_sd(sd, **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_SDS.png"
)
data_1 = {
500: 0.004900,
510: 0.009300,
520: 0.063270,
530: 0.165500,
540: 0.290400,
550: 0.433450,
560: 0.594500,
}
data_2 = {
500: 0.323000,
510: 0.503000,
520: 0.710000,
530: 0.862000,
540: 0.954000,
550: 0.994950,
560: 0.995000,
}
spd1 = SpectralDistribution(data_1, name="Custom 1")
spd2 = SpectralDistribution(data_2, name="Custom 2")
plt.close(plot_multi_sds([spd1, spd2], **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_CMFS.png"
)
plt.close(
plot_single_cmfs("CIE 1931 2 Degree Standard Observer", **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_CMFS.png"
)
cmfs = (
"CIE 1931 2 Degree Standard Observer",
"CIE 1964 10 Degree Standard Observer",
)
plt.close(plot_multi_cmfs(cmfs, **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Illuminant_SD.png"
)
plt.close(plot_single_illuminant_sd("A", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Illuminant_SDS.png"
)
plt.close(plot_multi_illuminant_sds(["A", "B", "C"], **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Visible_Spectrum.png"
)
plt.close(plot_visible_spectrum(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Lightness_Function.png"
)
plt.close(plot_single_lightness_function("CIE 1976", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Lightness_Functions.png"
)
plt.close(
plot_multi_lightness_functions(
["CIE 1976", "Wyszecki 1963"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Luminance_Function.png"
)
plt.close(plot_single_luminance_function("CIE 1976", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Luminance_Functions.png"
)
plt.close(
plot_multi_luminance_functions(
["CIE 1976", "Newhall 1943"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Blackbody_Spectral_Radiance.png"
)
plt.close(
plot_blackbody_spectral_radiance(
3500, blackbody="VY Canis Major", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Blackbody_Colours.png"
)
plt.close(
plot_blackbody_colours(SpectralShape(150, 12500, 50), **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Colour_Swatch.png"
)
RGB = ColourSwatch((0.45620519, 0.03081071, 0.04091952))
plt.close(plot_single_colour_swatch(RGB, **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Colour_Swatches.png"
)
RGB_1 = ColourSwatch((0.45293517, 0.31732158, 0.26414773))
RGB_2 = ColourSwatch((0.77875824, 0.57726450, 0.50453169))
plt.close(plot_multi_colour_swatches([RGB_1, RGB_2], **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Function.png"
)
plt.close(plot_single_function(lambda x: x ** (1 / 2.2), **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Functions.png"
)
functions = {
"Gamma 2.2": lambda x: x ** (1 / 2.2),
"Gamma 2.4": lambda x: x ** (1 / 2.4),
"Gamma 2.6": lambda x: x ** (1 / 2.6),
}
plt.close(plot_multi_functions(functions, **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Image.png"
)
path = os.path.join(
colour.__path__[0],
"examples",
"plotting",
"resources",
"Ishihara_Colour_Blindness_Test_Plate_3.png",
)
plt.close(plot_image(read_image(str(path)), **arguments)[0])
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Corresponding_Chromaticities_Prediction.png",
)
plt.close(
plot_corresponding_chromaticities_prediction(
1, "Von Kries", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Spectral_Locus.png"
)
plt.close(
plot_spectral_locus(spectral_locus_colours="RGB", **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Chromaticity_Diagram_Colours.png"
)
plt.close(
plot_chromaticity_diagram_colours(diagram_colours="RGB", **arguments)[
0
]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Chromaticity_Diagram.png"
)
plt.close(plot_chromaticity_diagram(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Chromaticity_Diagram_CIE1931.png"
)
plt.close(plot_chromaticity_diagram_CIE1931(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Chromaticity_Diagram_CIE1960UCS.png"
)
plt.close(plot_chromaticity_diagram_CIE1960UCS(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Chromaticity_Diagram_CIE1976UCS.png"
)
plt.close(plot_chromaticity_diagram_CIE1976UCS(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_SDS_In_Chromaticity_Diagram.png"
)
A = SDS_ILLUMINANTS["A"]
D65 = SDS_ILLUMINANTS["D65"]
plt.close(plot_sds_in_chromaticity_diagram([A, D65], **arguments)[0])
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_SDS_In_Chromaticity_Diagram_CIE1931.png",
)
plt.close(
plot_sds_in_chromaticity_diagram_CIE1931([A, D65], **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_SDS_In_Chromaticity_Diagram_CIE1960UCS.png",
)
plt.close(
plot_sds_in_chromaticity_diagram_CIE1960UCS([A, D65], **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_SDS_In_Chromaticity_Diagram_CIE1976UCS.png",
)
plt.close(
plot_sds_in_chromaticity_diagram_CIE1976UCS([A, D65], **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Pointer_Gamut.png"
)
plt.close(plot_pointer_gamut(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Colourspaces_In_Chromaticity_Diagram.png",
)
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram(
["ITU-R BT.709", "ACEScg", "S-Gamut"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Colourspaces_In_Chromaticity_Diagram_CIE1931.png",
)
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram_CIE1931(
["ITU-R BT.709", "ACEScg", "S-Gamut"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Colourspaces_In_"
"Chromaticity_Diagram_CIE1960UCS.png",
)
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram_CIE1960UCS(
["ITU-R BT.709", "ACEScg", "S-Gamut"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Colourspaces_In_"
"Chromaticity_Diagram_CIE1976UCS.png",
)
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram_CIE1976UCS(
["ITU-R BT.709", "ACEScg", "S-Gamut"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Chromaticities_In_" "Chromaticity_Diagram.png",
)
RGB = np.random.random((128, 128, 3))
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram(
RGB, "ITU-R BT.709", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Chromaticities_In_"
"Chromaticity_Diagram_CIE1931.png",
)
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1931(
RGB, "ITU-R BT.709", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Chromaticities_In_"
"Chromaticity_Diagram_CIE1960UCS.png",
)
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1960UCS(
RGB, "ITU-R BT.709", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_RGB_Chromaticities_In_"
"Chromaticity_Diagram_CIE1976UCS.png",
)
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1976UCS(
RGB, "ITU-R BT.709", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Ellipses_MacAdam1942_In_Chromaticity_Diagram.png",
)
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram(**arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Ellipses_MacAdam1942_In_"
"Chromaticity_Diagram_CIE1931.png",
)
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram_CIE1931(**arguments)[
0
]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Ellipses_MacAdam1942_In_"
"Chromaticity_Diagram_CIE1960UCS.png",
)
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram_CIE1960UCS(
**arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Ellipses_MacAdam1942_In_"
"Chromaticity_Diagram_CIE1976UCS.png",
)
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram_CIE1976UCS(
**arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_CCTF.png"
)
plt.close(plot_single_cctf("ITU-R BT.709", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_CCTFs.png"
)
plt.close(plot_multi_cctfs(["ITU-R BT.709", "sRGB"], **arguments)[0])
data = np.array(
[
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.40920000, 0.28120000, 0.30600000]),
np.array(
[
[0.02495100, 0.01908600, 0.02032900],
[0.10944300, 0.06235900, 0.06788100],
[0.27186500, 0.18418700, 0.19565300],
[0.48898900, 0.40749400, 0.44854600],
]
),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.30760000, 0.48280000, 0.42770000]),
np.array(
[
[0.02108000, 0.02989100, 0.02790400],
[0.06194700, 0.11251000, 0.09334400],
[0.15255800, 0.28123300, 0.23234900],
[0.34157700, 0.56681300, 0.47035300],
]
),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.39530000, 0.28120000, 0.18450000]),
np.array(
[
[0.02436400, 0.01908600, 0.01468800],
[0.10331200, 0.06235900, 0.02854600],
[0.26311900, 0.18418700, 0.12109700],
[0.43158700, 0.40749400, 0.39008600],
]
),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.20510000, 0.18420000, 0.57130000]),
np.array(
[
[0.03039800, 0.02989100, 0.06123300],
[0.08870000, 0.08498400, 0.21843500],
[0.18405800, 0.18418700, 0.40111400],
[0.32550100, 0.34047200, 0.50296900],
[0.53826100, 0.56681300, 0.80010400],
]
),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.35770000, 0.28120000, 0.11250000]),
np.array(
[
[0.03678100, 0.02989100, 0.01481100],
[0.17127700, 0.11251000, 0.01229900],
[0.30080900, 0.28123300, 0.21229800],
[0.52976000, 0.40749400, 0.11720000],
]
),
None,
],
]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Constant_Hue_Loci.png"
)
plt.close(plot_constant_hue_loci(data, "IPT", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_Munsell_Value_Function.png"
)
plt.close(plot_single_munsell_value_function("ASTM D1535", **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Multi_Munsell_Value_Functions.png"
)
plt.close(
plot_multi_munsell_value_functions(
["ASTM D1535", "McCamy 1987"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Single_SD_Rayleigh_Scattering.png"
)
plt.close(plot_single_sd_rayleigh_scattering(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_The_Blue_Sky.png"
)
plt.close(plot_the_blue_sky(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Colour_Quality_Bars.png"
)
illuminant = SDS_ILLUMINANTS["FL2"]
light_source = SDS_LIGHT_SOURCES["Kinoton 75P"]
light_source = light_source.copy().align(SpectralShape(360, 830, 1))
cqs_i = colour_quality_scale(illuminant, additional_data=True)
cqs_l = colour_quality_scale(light_source, additional_data=True)
plt.close(plot_colour_quality_bars([cqs_i, cqs_l], **arguments)[0])
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Single_SD_Colour_Rendering_Index_Bars.png",
)
illuminant = SDS_ILLUMINANTS["FL2"]
plt.close(
plot_single_sd_colour_rendering_index_bars(illuminant, **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Multi_SDS_Colour_Rendering_Indexes_Bars.png",
)
light_source = SDS_LIGHT_SOURCES["Kinoton 75P"]
plt.close(
plot_multi_sds_colour_rendering_indexes_bars(
[illuminant, light_source], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Single_SD_Colour_Quality_Scale_Bars.png",
)
illuminant = SDS_ILLUMINANTS["FL2"]
plt.close(
plot_single_sd_colour_quality_scale_bars(illuminant, **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Multi_SDS_Colour_Quality_Scales_Bars.png",
)
light_source = SDS_LIGHT_SOURCES["Kinoton 75P"]
plt.close(
plot_multi_sds_colour_quality_scales_bars(
[illuminant, light_source], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Hull_Section_Colours.png"
)
vertices, faces, _outline = primitive_cube(1, 1, 1, 64, 64, 64)
XYZ_vertices = RGB_to_XYZ(
vertices["position"] + 0.5,
RGB_COLOURSPACE_sRGB.whitepoint,
RGB_COLOURSPACE_sRGB.whitepoint,
RGB_COLOURSPACE_sRGB.matrix_RGB_to_XYZ,
)
hull = trimesh.Trimesh(XYZ_vertices, faces, process=False)
plt.close(
plot_hull_section_colours(hull, section_colours="RGB", **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Hull_Section_Contour.png"
)
plt.close(
plot_hull_section_contour(hull, section_colours="RGB", **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Visible_Spectrum_Section.png"
)
plt.close(
plot_visible_spectrum_section(
section_colours="RGB", section_opacity=0.15, **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_RGB_Colourspace_Section.png"
)
plt.close(
plot_RGB_colourspace_section(
"sRGB", section_colours="RGB", section_opacity=0.15, **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_Planckian_Locus.png"
)
plt.close(
plot_planckian_locus(planckian_locus_colours="RGB", **arguments)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Planckian_Locus_In_Chromaticity_Diagram.png",
)
plt.close(
plot_planckian_locus_in_chromaticity_diagram(
["A", "B", "C"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Planckian_Locus_In_Chromaticity_Diagram_CIE1931.png",
)
plt.close(
plot_planckian_locus_in_chromaticity_diagram_CIE1931(
["A", "B", "C"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Planckian_Locus_In_Chromaticity_Diagram_CIE1960UCS.png",
)
plt.close(
plot_planckian_locus_in_chromaticity_diagram_CIE1960UCS(
["A", "B", "C"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Single_SD_Colour_Rendition_Report_Full.png",
)
plt.close(
plot_single_sd_colour_rendition_report(
SDS_ILLUMINANTS["FL2"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Single_SD_Colour_Rendition_Report_Intermediate.png",
)
plt.close(
plot_single_sd_colour_rendition_report(
SDS_ILLUMINANTS["FL2"], "Intermediate", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory,
"Plotting_Plot_Single_SD_Colour_Rendition_Report_Simple.png",
)
plt.close(
plot_single_sd_colour_rendition_report(
SDS_ILLUMINANTS["FL2"], "Simple", **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_RGB_Colourspaces_Gamuts.png"
)
plt.close(
plot_RGB_colourspaces_gamuts(
["ITU-R BT.709", "ACEScg", "S-Gamut"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_RGB_Colourspaces_Gamuts.png"
)
plt.close(
plot_RGB_colourspaces_gamuts(
["ITU-R BT.709", "ACEScg", "S-Gamut"], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Plotting_Plot_RGB_Scatter.png"
)
plt.close(plot_RGB_scatter(RGB, "ITU-R BT.709", **arguments)[0])
filename = os.path.join(
output_directory, "Plotting_Plot_Colour_Automatic_Conversion_Graph.png"
)
plot_automatic_colour_conversion_graph(filename)
# *************************************************************************
# "tutorial.rst"
# *************************************************************************
arguments["filename"] = os.path.join(
output_directory, "Tutorial_Visible_Spectrum.png"
)
plt.close(plot_visible_spectrum(**arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Tutorial_Sample_SD.png"
)
sample_sd_data = {
380: 0.048,
385: 0.051,
390: 0.055,
395: 0.060,
400: 0.065,
405: 0.068,
410: 0.068,
415: 0.067,
420: 0.064,
425: 0.062,
430: 0.059,
435: 0.057,
440: 0.055,
445: 0.054,
450: 0.053,
455: 0.053,
460: 0.052,
465: 0.052,
470: 0.052,
475: 0.053,
480: 0.054,
485: 0.055,
490: 0.057,
495: 0.059,
500: 0.061,
505: 0.062,
510: 0.065,
515: 0.067,
520: 0.070,
525: 0.072,
530: 0.074,
535: 0.075,
540: 0.076,
545: 0.078,
550: 0.079,
555: 0.082,
560: 0.087,
565: 0.092,
570: 0.100,
575: 0.107,
580: 0.115,
585: 0.122,
590: 0.129,
595: 0.134,
600: 0.138,
605: 0.142,
610: 0.146,
615: 0.150,
620: 0.154,
625: 0.158,
630: 0.163,
635: 0.167,
640: 0.173,
645: 0.180,
650: 0.188,
655: 0.196,
660: 0.204,
665: 0.213,
670: 0.222,
675: 0.231,
680: 0.242,
685: 0.251,
690: 0.261,
695: 0.271,
700: 0.282,
705: 0.294,
710: 0.305,
715: 0.318,
720: 0.334,
725: 0.354,
730: 0.372,
735: 0.392,
740: 0.409,
745: 0.420,
750: 0.436,
755: 0.450,
760: 0.462,
765: 0.465,
770: 0.448,
775: 0.432,
780: 0.421,
}
sd = SpectralDistribution(sample_sd_data, name="Sample")
plt.close(plot_single_sd(sd, **arguments)[0])
arguments["filename"] = os.path.join(
output_directory, "Tutorial_SD_Interpolation.png"
)
sd_copy = sd.copy()
sd_copy.interpolate(SpectralShape(400, 770, 1))
plt.close(
plot_multi_sds(
[sd, sd_copy], bounding_box=[730, 780, 0.25, 0.5], **arguments
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Tutorial_Sample_Swatch.png"
)
sd = SpectralDistribution(sample_sd_data)
cmfs = MSDS_CMFS_STANDARD_OBSERVER["CIE 1931 2 Degree Standard Observer"]
illuminant = SDS_ILLUMINANTS["D65"]
with domain_range_scale("1"):
XYZ = sd_to_XYZ(sd, cmfs, illuminant)
RGB = XYZ_to_sRGB(XYZ)
plt.close(
plot_single_colour_swatch(
ColourSwatch(RGB, "Sample"),
text_kwargs={"size": "x-large"},
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Tutorial_Neutral5.png"
)
patch_name = "neutral 5 (.70 D)"
patch_sd = SDS_COLOURCHECKERS["ColorChecker N Ohta"][patch_name]
with domain_range_scale("1"):
XYZ = sd_to_XYZ(patch_sd, cmfs, illuminant)
RGB = XYZ_to_sRGB(XYZ)
plt.close(
plot_single_colour_swatch(
ColourSwatch(RGB, patch_name.title()),
text_kwargs={"size": "x-large"},
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Tutorial_Colour_Checker.png"
)
plt.close(
plot_single_colour_checker(
colour_checker="ColorChecker 2005",
text_kwargs={"visible": False},
**arguments,
)[0]
)
arguments["filename"] = os.path.join(
output_directory, "Tutorial_CIE_1931_Chromaticity_Diagram.png"
)
xy = XYZ_to_xy(XYZ)
plot_chromaticity_diagram_CIE1931(standalone=False)
x, y = xy
plt.plot(x, y, "o-", color="white")
# Annotating the plot.
plt.annotate(
patch_sd.name.title(),
xy=xy,
xytext=(-50, 30),
textcoords="offset points",
arrowprops=dict(arrowstyle="->", connectionstyle="arc3, rad=-0.2"),
)
plt.close(
render(
standalone=True,
limits=(-0.1, 0.9, -0.1, 0.9),
x_tighten=True,
y_tighten=True,
**arguments,
)[0]
)
# *************************************************************************
# "basics.rst"
# *************************************************************************
arguments["filename"] = os.path.join(
output_directory, "Basics_Logo_Small_001_CIE_XYZ.png"
)
RGB = read_image(os.path.join(output_directory, "Logo_Small_001.png"))[
..., 0:3
]
XYZ = sRGB_to_XYZ(RGB)
plt.close(
plot_image(XYZ, text_kwargs={"text": "sRGB to XYZ"}, **arguments)[0]
)
def generate_documentation_plots(output_directory):
"""
Generates documentation plots.
Parameters
----------
output_directory : unicode
Output directory.
"""
filter_warnings()
colour_style()
np.random.seed(0)
# *************************************************************************
# "README.rst"
# *************************************************************************
filename = os.path.join(output_directory,
'Examples_Colour_Automatic_Conversion_Graph.png')
plot_automatic_colour_conversion_graph(filename)
arguments = {
'tight_layout':
True,
'transparent_background':
True,
'filename':
os.path.join(output_directory,
'Examples_Plotting_Visible_Spectrum.png')
}
plt.close(
plot_visible_spectrum('CIE 1931 2 Degree Standard Observer',
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Examples_Plotting_Illuminant_F1_SD.png')
plt.close(plot_single_illuminant_sd('FL1', **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Examples_Plotting_Blackbodies.png')
blackbody_sds = [
sd_blackbody(i, SpectralShape(0, 10000, 10))
for i in range(1000, 15000, 1000)
]
plt.close(
plot_multi_sds(blackbody_sds,
y_label='W / (sr m$^2$) / m',
use_sds_colours=True,
normalise_sds_colours=True,
legend_location='upper right',
bounding_box=(0, 1250, 0, 2.5e15),
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Examples_Plotting_Cone_Fundamentals.png')
plt.close(
plot_single_cmfs('Stockman & Sharpe 2 Degree Cone Fundamentals',
y_label='Sensitivity',
bounding_box=(390, 870, 0, 1.1),
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Examples_Plotting_Luminous_Efficiency.png')
plt.close(
plot_multi_sds((sd_mesopic_luminous_efficiency_function(0.2),
PHOTOPIC_LEFS['CIE 1924 Photopic Standard Observer'],
SCOTOPIC_LEFS['CIE 1951 Scotopic Standard Observer']),
y_label='Luminous Efficiency',
legend_location='upper right',
y_tighten=True,
margins=(0, 0, 0, .1),
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Examples_Plotting_BabelColor_Average.png')
plt.close(
plot_multi_sds(COLOURCHECKERS_SDS['BabelColor Average'].values(),
use_sds_colours=True,
title=('BabelColor Average - '
'Spectral Distributions'),
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Examples_Plotting_ColorChecker_2005.png')
plt.close(
plot_single_colour_checker('ColorChecker 2005',
text_parameters={'visible': False},
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Examples_Plotting_Chromaticities_Prediction.png')
plt.close(
plot_corresponding_chromaticities_prediction(2, 'Von Kries', 'Bianco',
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Examples_Plotting_CCT_CIE_1960_UCS_Chromaticity_Diagram.png')
plt.close(
plot_planckian_locus_in_chromaticity_diagram_CIE1960UCS(
['A', 'B', 'C'], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Examples_Plotting_Chromaticities_CIE_1931_Chromaticity_Diagram.png')
RGB = np.random.random((32, 32, 3))
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1931(
RGB,
'ITU-R BT.709',
colourspaces=['ACEScg', 'S-Gamut'],
show_pointer_gamut=True,
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Examples_Plotting_CRI.png')
plt.close(
plot_single_sd_colour_rendering_index_bars(ILLUMINANTS_SDS['FL2'],
**arguments)[0])
# *************************************************************************
# Documentation
# *************************************************************************
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_CVD_Simulation_Machado2009.png')
plt.close(plot_cvd_simulation_Machado2009(RGB, **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_Colour_Checker.png')
plt.close(plot_single_colour_checker('ColorChecker 2005', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Multi_Colour_Checkers.png')
plt.close(
plot_multi_colour_checkers(['ColorChecker 1976', 'ColorChecker 2005'],
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Single_SD.png')
data = {
500: 0.0651,
520: 0.0705,
540: 0.0772,
560: 0.0870,
580: 0.1128,
600: 0.1360
}
sd = SpectralDistribution(data, name='Custom')
plt.close(plot_single_sd(sd, **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Multi_SDS.png')
data_1 = {
500: 0.004900,
510: 0.009300,
520: 0.063270,
530: 0.165500,
540: 0.290400,
550: 0.433450,
560: 0.594500
}
data_2 = {
500: 0.323000,
510: 0.503000,
520: 0.710000,
530: 0.862000,
540: 0.954000,
550: 0.994950,
560: 0.995000
}
spd1 = SpectralDistribution(data_1, name='Custom 1')
spd2 = SpectralDistribution(data_2, name='Custom 2')
plt.close(plot_multi_sds([spd1, spd2], **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Single_CMFS.png')
plt.close(
plot_single_cmfs('CIE 1931 2 Degree Standard Observer',
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Multi_CMFS.png')
cmfs = ('CIE 1931 2 Degree Standard Observer',
'CIE 1964 10 Degree Standard Observer')
plt.close(plot_multi_cmfs(cmfs, **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_Illuminant_SD.png')
plt.close(plot_single_illuminant_sd('A', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Multi_Illuminant_SDS.png')
plt.close(plot_multi_illuminant_sds(['A', 'B', 'C'], **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Visible_Spectrum.png')
plt.close(plot_visible_spectrum(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_Lightness_Function.png')
plt.close(plot_single_lightness_function('CIE 1976', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Multi_Lightness_Functions.png')
plt.close(
plot_multi_lightness_functions(['CIE 1976', 'Wyszecki 1963'],
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_Luminance_Function.png')
plt.close(plot_single_luminance_function('CIE 1976', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Multi_Luminance_Functions.png')
plt.close(
plot_multi_luminance_functions(['CIE 1976', 'Newhall 1943'],
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Blackbody_Spectral_Radiance.png')
plt.close(
plot_blackbody_spectral_radiance(3500,
blackbody='VY Canis Major',
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Blackbody_Colours.png')
plt.close(
plot_blackbody_colours(SpectralShape(150, 12500, 50), **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_Colour_Swatch.png')
RGB = ColourSwatch(RGB=(0.45620519, 0.03081071, 0.04091952))
plt.close(plot_single_colour_swatch(RGB, **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Multi_Colour_Swatches.png')
RGB_1 = ColourSwatch(RGB=(0.45293517, 0.31732158, 0.26414773))
RGB_2 = ColourSwatch(RGB=(0.77875824, 0.57726450, 0.50453169))
plt.close(plot_multi_colour_swatches([RGB_1, RGB_2], **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Single_Function.png')
plt.close(plot_single_function(lambda x: x**(1 / 2.2), **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Multi_Functions.png')
functions = {
'Gamma 2.2': lambda x: x**(1 / 2.2),
'Gamma 2.4': lambda x: x**(1 / 2.4),
'Gamma 2.6': lambda x: x**(1 / 2.6),
}
plt.close(plot_multi_functions(functions, **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Image.png')
path = os.path.join(output_directory, 'Logo_Medium_001.png')
plt.close(plot_image(read_image(str(path)), **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Corresponding_Chromaticities_Prediction.png')
plt.close(
plot_corresponding_chromaticities_prediction(1, 'Von Kries', 'CAT02',
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Spectral_Locus.png')
plt.close(
plot_spectral_locus(spectral_locus_colours='RGB', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Chromaticity_Diagram_Colours.png')
plt.close(plot_chromaticity_diagram_colours(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Chromaticity_Diagram.png')
plt.close(plot_chromaticity_diagram(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Chromaticity_Diagram_CIE1931.png')
plt.close(plot_chromaticity_diagram_CIE1931(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Chromaticity_Diagram_CIE1960UCS.png')
plt.close(plot_chromaticity_diagram_CIE1960UCS(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Chromaticity_Diagram_CIE1976UCS.png')
plt.close(plot_chromaticity_diagram_CIE1976UCS(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_SDS_In_Chromaticity_Diagram.png')
A = ILLUMINANTS_SDS['A']
D65 = ILLUMINANTS_SDS['D65']
plt.close(plot_sds_in_chromaticity_diagram([A, D65], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_SDS_In_Chromaticity_Diagram_CIE1931.png')
plt.close(
plot_sds_in_chromaticity_diagram_CIE1931([A, D65], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_SDS_In_Chromaticity_Diagram_CIE1960UCS.png')
plt.close(
plot_sds_in_chromaticity_diagram_CIE1960UCS([A, D65], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_SDS_In_Chromaticity_Diagram_CIE1976UCS.png')
plt.close(
plot_sds_in_chromaticity_diagram_CIE1976UCS([A, D65], **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Pointer_Gamut.png')
plt.close(plot_pointer_gamut(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_RGB_Colourspaces_In_Chromaticity_Diagram.png')
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram(
['ITU-R BT.709', 'ACEScg', 'S-Gamut'], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_RGB_Colourspaces_In_Chromaticity_Diagram_CIE1931.png')
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram_CIE1931(
['ITU-R BT.709', 'ACEScg', 'S-Gamut'], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Colourspaces_In_'
'Chromaticity_Diagram_CIE1960UCS.png')
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram_CIE1960UCS(
['ITU-R BT.709', 'ACEScg', 'S-Gamut'], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Colourspaces_In_'
'Chromaticity_Diagram_CIE1976UCS.png')
plt.close(
plot_RGB_colourspaces_in_chromaticity_diagram_CIE1976UCS(
['ITU-R BT.709', 'ACEScg', 'S-Gamut'], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Chromaticities_In_'
'Chromaticity_Diagram.png')
RGB = np.random.random((128, 128, 3))
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram(
RGB, 'ITU-R BT.709', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Chromaticities_In_'
'Chromaticity_Diagram_CIE1931.png')
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1931(
RGB, 'ITU-R BT.709', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Chromaticities_In_'
'Chromaticity_Diagram_CIE1960UCS.png')
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1960UCS(
RGB, 'ITU-R BT.709', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Chromaticities_In_'
'Chromaticity_Diagram_CIE1976UCS.png')
plt.close(
plot_RGB_chromaticities_in_chromaticity_diagram_CIE1976UCS(
RGB, 'ITU-R BT.709', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Ellipses_MacAdam1942_In_Chromaticity_Diagram.png')
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Ellipses_MacAdam1942_In_'
'Chromaticity_Diagram_CIE1931.png')
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram_CIE1931(
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Ellipses_MacAdam1942_In_'
'Chromaticity_Diagram_CIE1960UCS.png')
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram_CIE1960UCS(
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Ellipses_MacAdam1942_In_'
'Chromaticity_Diagram_CIE1976UCS.png')
plt.close(
plot_ellipses_MacAdam1942_in_chromaticity_diagram_CIE1976UCS(
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Single_CCTF.png')
plt.close(plot_single_cctf('ITU-R BT.709', **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Multi_CCTFs.png')
plt.close(plot_multi_cctfs(['ITU-R BT.709', 'sRGB'], **arguments)[0])
data = np.array([
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.40920000, 0.28120000, 0.30600000]),
np.array([
[0.02495100, 0.01908600, 0.02032900],
[0.10944300, 0.06235900, 0.06788100],
[0.27186500, 0.18418700, 0.19565300],
[0.48898900, 0.40749400, 0.44854600],
]),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.30760000, 0.48280000, 0.42770000]),
np.array([
[0.02108000, 0.02989100, 0.02790400],
[0.06194700, 0.11251000, 0.09334400],
[0.15255800, 0.28123300, 0.23234900],
[0.34157700, 0.56681300, 0.47035300],
]),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.39530000, 0.28120000, 0.18450000]),
np.array([
[0.02436400, 0.01908600, 0.01468800],
[0.10331200, 0.06235900, 0.02854600],
[0.26311900, 0.18418700, 0.12109700],
[0.43158700, 0.40749400, 0.39008600],
]),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.20510000, 0.18420000, 0.57130000]),
np.array([
[0.03039800, 0.02989100, 0.06123300],
[0.08870000, 0.08498400, 0.21843500],
[0.18405800, 0.18418700, 0.40111400],
[0.32550100, 0.34047200, 0.50296900],
[0.53826100, 0.56681300, 0.80010400],
]),
None,
],
[
None,
np.array([0.95010000, 1.00000000, 1.08810000]),
np.array([0.35770000, 0.28120000, 0.11250000]),
np.array([
[0.03678100, 0.02989100, 0.01481100],
[0.17127700, 0.11251000, 0.01229900],
[0.30080900, 0.28123300, 0.21229800],
[0.52976000, 0.40749400, 0.11720000],
]),
None,
],
])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Constant_Hue_Loci.png')
plt.close(plot_constant_hue_loci(data, 'IPT', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_Munsell_Value_Function.png')
plt.close(plot_single_munsell_value_function('ASTM D1535', **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Multi_Munsell_Value_Functions.png')
plt.close(
plot_multi_munsell_value_functions(['ASTM D1535', 'McCamy 1987'],
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Single_SD_Rayleigh_Scattering.png')
plt.close(plot_single_sd_rayleigh_scattering(**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_The_Blue_Sky.png')
plt.close(plot_the_blue_sky(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Colour_Quality_Bars.png')
illuminant = ILLUMINANTS_SDS['FL2']
light_source = LIGHT_SOURCES_SDS['Kinoton 75P']
light_source = light_source.copy().align(SpectralShape(360, 830, 1))
cqs_i = colour_quality_scale(illuminant, additional_data=True)
cqs_l = colour_quality_scale(light_source, additional_data=True)
plt.close(plot_colour_quality_bars([cqs_i, cqs_l], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Single_SD_Colour_Rendering_Index_Bars.png')
illuminant = ILLUMINANTS_SDS['FL2']
plt.close(
plot_single_sd_colour_rendering_index_bars(illuminant, **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Multi_SDS_Colour_Rendering_Indexes_Bars.png')
light_source = LIGHT_SOURCES_SDS['Kinoton 75P']
plt.close(
plot_multi_sds_colour_rendering_indexes_bars(
[illuminant, light_source], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Single_SD_Colour_Quality_Scale_Bars.png')
illuminant = ILLUMINANTS_SDS['FL2']
plt.close(
plot_single_sd_colour_quality_scale_bars(illuminant, **arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Multi_SDS_Colour_Quality_Scales_Bars.png')
light_source = LIGHT_SOURCES_SDS['Kinoton 75P']
plt.close(
plot_multi_sds_colour_quality_scales_bars([illuminant, light_source],
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_Planckian_Locus.png')
plt.close(plot_planckian_locus(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Planckian_Locus_CIE1931.png')
plt.close(plot_planckian_locus_CIE1931(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_Planckian_Locus_CIE1960UCS.png')
plt.close(plot_planckian_locus_CIE1960UCS(**arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Planckian_Locus_In_Chromaticity_Diagram.png')
plt.close(
plot_planckian_locus_in_chromaticity_diagram(['A', 'B', 'C'],
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Planckian_Locus_In_Chromaticity_Diagram_CIE1931.png')
plt.close(
plot_planckian_locus_in_chromaticity_diagram_CIE1931(['A', 'B', 'C'],
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory,
'Plotting_Plot_Planckian_Locus_In_Chromaticity_Diagram_CIE1960UCS.png')
plt.close(
plot_planckian_locus_in_chromaticity_diagram_CIE1960UCS(
['A', 'B', 'C'], **arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Colourspaces_Gamuts.png')
plt.close(
plot_RGB_colourspaces_gamuts(['ITU-R BT.709', 'ACEScg', 'S-Gamut'],
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Plotting_Plot_RGB_Colourspaces_Gamuts.png')
plt.close(
plot_RGB_colourspaces_gamuts(['ITU-R BT.709', 'ACEScg', 'S-Gamut'],
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Plotting_Plot_RGB_Scatter.png')
plt.close(plot_RGB_scatter(RGB, 'ITU-R BT.709', **arguments)[0])
filename = os.path.join(
output_directory,
'Plotting_Plot_Colour_Automatic_Conversion_Graph.png')
plot_automatic_colour_conversion_graph(filename)
# *************************************************************************
# "tutorial.rst"
# *************************************************************************
arguments['filename'] = os.path.join(output_directory,
'Tutorial_Visible_Spectrum.png')
plt.close(plot_visible_spectrum(**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Tutorial_Sample_SD.png')
sample_sd_data = {
380: 0.048,
385: 0.051,
390: 0.055,
395: 0.060,
400: 0.065,
405: 0.068,
410: 0.068,
415: 0.067,
420: 0.064,
425: 0.062,
430: 0.059,
435: 0.057,
440: 0.055,
445: 0.054,
450: 0.053,
455: 0.053,
460: 0.052,
465: 0.052,
470: 0.052,
475: 0.053,
480: 0.054,
485: 0.055,
490: 0.057,
495: 0.059,
500: 0.061,
505: 0.062,
510: 0.065,
515: 0.067,
520: 0.070,
525: 0.072,
530: 0.074,
535: 0.075,
540: 0.076,
545: 0.078,
550: 0.079,
555: 0.082,
560: 0.087,
565: 0.092,
570: 0.100,
575: 0.107,
580: 0.115,
585: 0.122,
590: 0.129,
595: 0.134,
600: 0.138,
605: 0.142,
610: 0.146,
615: 0.150,
620: 0.154,
625: 0.158,
630: 0.163,
635: 0.167,
640: 0.173,
645: 0.180,
650: 0.188,
655: 0.196,
660: 0.204,
665: 0.213,
670: 0.222,
675: 0.231,
680: 0.242,
685: 0.251,
690: 0.261,
695: 0.271,
700: 0.282,
705: 0.294,
710: 0.305,
715: 0.318,
720: 0.334,
725: 0.354,
730: 0.372,
735: 0.392,
740: 0.409,
745: 0.420,
750: 0.436,
755: 0.450,
760: 0.462,
765: 0.465,
770: 0.448,
775: 0.432,
780: 0.421
}
sd = SpectralDistribution(sample_sd_data, name='Sample')
plt.close(plot_single_sd(sd, **arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Tutorial_SD_Interpolation.png')
sd_copy = sd.copy()
sd_copy.interpolate(SpectralShape(400, 770, 1))
plt.close(
plot_multi_sds([sd, sd_copy],
bounding_box=[730, 780, 0.25, 0.5],
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Tutorial_Sample_Swatch.png')
sd = SpectralDistribution(sample_sd_data)
cmfs = STANDARD_OBSERVERS_CMFS['CIE 1931 2 Degree Standard Observer']
illuminant = ILLUMINANTS_SDS['D65']
with domain_range_scale('1'):
XYZ = sd_to_XYZ(sd, cmfs, illuminant)
RGB = XYZ_to_sRGB(XYZ)
plt.close(
plot_single_colour_swatch(ColourSwatch('Sample', RGB),
text_parameters={'size': 'x-large'},
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Tutorial_Neutral5.png')
patch_name = 'neutral 5 (.70 D)'
patch_sd = COLOURCHECKERS_SDS['ColorChecker N Ohta'][patch_name]
with domain_range_scale('1'):
XYZ = sd_to_XYZ(patch_sd, cmfs, illuminant)
RGB = XYZ_to_sRGB(XYZ)
plt.close(
plot_single_colour_swatch(ColourSwatch(patch_name.title(), RGB),
text_parameters={'size': 'x-large'},
**arguments)[0])
arguments['filename'] = os.path.join(output_directory,
'Tutorial_Colour_Checker.png')
plt.close(
plot_single_colour_checker(colour_checker='ColorChecker 2005',
text_parameters={'visible': False},
**arguments)[0])
arguments['filename'] = os.path.join(
output_directory, 'Tutorial_CIE_1931_Chromaticity_Diagram.png')
xy = XYZ_to_xy(XYZ)
plot_chromaticity_diagram_CIE1931(standalone=False)
x, y = xy
plt.plot(x, y, 'o-', color='white')
# Annotating the plot.
plt.annotate(patch_sd.name.title(),
xy=xy,
xytext=(-50, 30),
textcoords='offset points',
arrowprops=dict(arrowstyle='->',
connectionstyle='arc3, rad=-0.2'))
plt.close(
render(standalone=True,
limits=(-0.1, 0.9, -0.1, 0.9),
x_tighten=True,
y_tighten=True,
**arguments)[0])
# *************************************************************************
# "basics.rst"
# *************************************************************************
arguments['filename'] = os.path.join(output_directory,
'Basics_Logo_Small_001_CIE_XYZ.png')
RGB = read_image(os.path.join(output_directory, 'Logo_Small_001.png'))[...,
0:3]
XYZ = sRGB_to_XYZ(RGB)
plt.close(
plot_image(XYZ, text_parameters={'text': 'sRGB to XYZ'},
**arguments)[0])
def colour_quality_bars_plot(specifications,
labels=True,
hatching=None,
hatching_repeat=1,
**kwargs):
"""
Plots the colour quality data of given illuminants or light sources colour
quality specifications.
Parameters
----------
specifications : array_like
Array of illuminants or light sources colour quality specifications.
labels : bool, optional
Add labels above bars.
hatching : bool or None, optional
Use hatching for the bars.
hatching_repeat : int, optional
Hatching pattern repeat.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`boundaries`, :func:`canvas`, :func:`decorate`,
:func:`display`},
Please refer to the documentation of the previously listed definitions.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> from colour import (
... ILLUMINANTS_RELATIVE_SPDS,
... LIGHT_SOURCES_RELATIVE_SPDS,
... SpectralShape)
>>> illuminant = ILLUMINANTS_RELATIVE_SPDS['F2']
>>> light_source = LIGHT_SOURCES_RELATIVE_SPDS['Kinoton 75P']
>>> light_source = light_source.clone().align(SpectralShape(360, 830, 1))
>>> cqs_i = colour_quality_scale(illuminant, additional_data=True)
>>> cqs_l = colour_quality_scale(light_source, additional_data=True)
>>> colour_quality_bars_plot([cqs_i, cqs_l]) # doctest: +SKIP
"""
settings = {'figure_size': (DEFAULT_FIGURE_WIDTH, DEFAULT_FIGURE_WIDTH)}
settings.update(kwargs)
canvas(**settings)
bar_width = 0.5
y_ticks_interval = 10
count_s, count_Q_as = len(specifications), 0
patterns = cycle(DEFAULT_HATCH_PATTERNS)
if hatching is None:
hatching = False if count_s == 1 else True
for i, specification in enumerate(specifications):
Q_a, Q_as, colorimetry_data = (specification.Q_a, specification.Q_as,
specification.colorimetry_data)
count_Q_as = len(Q_as)
colours = (
[[1] * 3] +
[np.clip(XYZ_to_sRGB(x.XYZ), 0, 1) for x in colorimetry_data[0]])
x = (i + np.arange(0, (count_Q_as + 1) * (count_s + 1), (count_s + 1),
dtype=np.float_)) * bar_width
y = [s[1].Q_a for s in sorted(Q_as.items(), key=lambda s: s[0])]
y = np.array([Q_a] + list(y))
if np.sign(np.min(y)) < 0:
warning(
('"{0}" spectral distribution has negative "Q_a" value(s), '
'using absolute value(s) '
'for plotting purpose!'.format(specification.name)))
y = np.abs(y)
bars = pylab.bar(x,
y,
color=colours,
width=bar_width,
hatch=(next(patterns) *
hatching_repeat if hatching else None),
label=specification.name)
if labels:
label_rectangles(
bars,
rotation='horizontal' if count_s == 1 else 'vertical',
offset=(0 if count_s == 1 else 3 / 100 * count_s + 65 / 1000,
0.025),
text_size=-5 / 7 * count_s + 12.5)
pylab.axhline(y=100, color='black', linestyle='--')
pylab.xticks(
(np.arange(0, (count_Q_as + 1) * (count_s + 1), (count_s + 1),
dtype=np.float_) * bar_width + (count_s * bar_width / 2)),
['Qa'] +
['Q{0}'.format(index + 1) for index in range(0, count_Q_as + 1, 1)])
pylab.yticks(range(0, 100 + y_ticks_interval, y_ticks_interval))
settings.update({
'title':
'Colour Quality',
'legend':
hatching,
'x_tighten':
True,
'y_tighten':
True,
'limits': (-bar_width, ((count_Q_as + 1) * (count_s + 1)) / 2, 0, 120),
'aspect':
1 / (120 / (bar_width + len(Q_as) + bar_width * 2))
})
settings.update(kwargs)
boundaries(**settings)
decorate(**settings)
return display(**settings)
def multi_spd_plot(spds,
cmfs='CIE 1931 2 Degree Standard Observer',
use_spds_colours=False,
normalise_spds_colours=False,
**kwargs):
"""
Plots given spectral power distributions.
Parameters
----------
spds : list
Spectral power distributions to plot.
cmfs : unicode, optional
Standard observer colour matching functions used for spectrum creation.
use_spds_colours : bool, optional
Whether to use spectral power distributions colours.
normalise_spds_colours : bool
Whether to normalise spectral power distributions colours.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`colour.plotting.render`},
Please refer to the documentation of the previously listed definition.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> from colour import SpectralPowerDistribution
>>> data_1 = {
... 500: 0.004900,
... 510: 0.009300,
... 520: 0.063270,
... 530: 0.165500,
... 540: 0.290400,
... 550: 0.433450,
... 560: 0.594500
... }
>>> data_2 = {
... 500: 0.323000,
... 510: 0.503000,
... 520: 0.710000,
... 530: 0.862000,
... 540: 0.954000,
... 550: 0.994950,
... 560: 0.995000
... }
>>> spd1 = SpectralPowerDistribution(data_1, name='Custom 1')
>>> spd2 = SpectralPowerDistribution(data_2, name='Custom 2')
>>> multi_spd_plot([spd1, spd2]) # doctest: +SKIP
"""
canvas(**kwargs)
cmfs = get_cmfs(cmfs)
if use_spds_colours:
illuminant = ILLUMINANTS_RELATIVE_SPDS['D65']
x_limit_min, x_limit_max, y_limit_min, y_limit_max = [], [], [], []
for spd in spds:
wavelengths, values = spd.wavelengths, spd.values
shape = spd.shape
x_limit_min.append(shape.start)
x_limit_max.append(shape.end)
y_limit_min.append(min(values))
y_limit_max.append(max(values))
if use_spds_colours:
XYZ = spectral_to_XYZ(spd, cmfs, illuminant) / 100
if normalise_spds_colours:
XYZ = normalise_maximum(XYZ, clip=False)
RGB = np.clip(XYZ_to_sRGB(XYZ), 0, 1)
pylab.plot(wavelengths,
values,
color=RGB,
label=spd.strict_name,
linewidth=1)
else:
pylab.plot(wavelengths, values, label=spd.strict_name, linewidth=1)
settings = {
'x_label':
'Wavelength $\\lambda$ (nm)',
'y_label':
'Spectral Power Distribution',
'x_tighten':
True,
'y_tighten':
True,
'legend':
True,
'legend_location':
'upper left',
'limits': (min(x_limit_min), max(x_limit_max), min(y_limit_min),
max(y_limit_max))
}
settings.update(kwargs)
return render(**settings)
def multi_spd_plot(spds,
cmfs='CIE 1931 2 Degree Standard Observer',
use_spds_colours=False,
normalise_spds_colours=False,
**kwargs):
"""
Plots given spectral power distributions.
Parameters
----------
spds : list
Spectral power distributions to plot.
cmfs : unicode, optional
Standard observer colour matching functions used for spectrum creation.
use_spds_colours : bool, optional
Use spectral power distributions colours.
normalise_spds_colours : bool
Should spectral power distributions colours normalised.
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> from colour import SpectralPowerDistribution
>>> data1 = {400: 0.0641, 420: 0.0645, 440: 0.0562}
>>> data2 = {400: 0.134, 420: 0.789, 440: 1.289}
>>> spd1 = SpectralPowerDistribution('Custom1', data1)
>>> spd2 = SpectralPowerDistribution('Custom2', data2)
>>> multi_spd_plot([spd1, spd2]) # doctest: +SKIP
"""
canvas(**kwargs)
cmfs = get_cmfs(cmfs)
if use_spds_colours:
illuminant = ILLUMINANTS_RELATIVE_SPDS.get('D65')
x_limit_min, x_limit_max, y_limit_min, y_limit_max = [], [], [], []
for spd in spds:
wavelengths, values = tuple(zip(*spd.items))
shape = spd.shape
x_limit_min.append(shape.start)
x_limit_max.append(shape.end)
y_limit_min.append(min(values))
y_limit_max.append(max(values))
if use_spds_colours:
XYZ = spectral_to_XYZ(spd, cmfs, illuminant) / 100
if normalise_spds_colours:
XYZ = normalise_maximum(XYZ, clip=False)
RGB = np.clip(XYZ_to_sRGB(XYZ), 0, 1)
pylab.plot(wavelengths,
values,
color=RGB,
label=spd.title,
linewidth=2)
else:
pylab.plot(wavelengths, values, label=spd.title, linewidth=2)
settings = {
'x_label':
'Wavelength $\\lambda$ (nm)',
'y_label':
'Spectral Power Distribution',
'x_tighten':
True,
'legend':
True,
'legend_location':
'upper left',
'limits': (min(x_limit_min), max(x_limit_max), min(y_limit_min),
max(y_limit_max))
}
settings.update(kwargs)
boundaries(**settings)
decorate(**settings)
return display(**settings)
def CIE_1931_chromaticity_diagram_colours_plot(
surface=1,
samples=4096,
cmfs='CIE 1931 2 Degree Standard Observer',
**kwargs):
"""
Plots the *CIE 1931 Chromaticity Diagram* colours.
Parameters
----------
surface : numeric, optional
Generated markers surface.
samples : numeric, optional
Samples count on one axis.
cmfs : unicode, optional
Standard observer colour matching functions used for diagram bounds.
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> CIE_1931_chromaticity_diagram_colours_plot() # doctest: +SKIP
True
"""
if is_scipy_installed(raise_exception=True):
from scipy.spatial import Delaunay
settings = {'figure_size': (64, 64)}
settings.update(kwargs)
canvas(**settings)
cmfs = get_cmfs(cmfs)
illuminant = DEFAULT_PLOTTING_ILLUMINANT
triangulation = Delaunay(XYZ_to_xy(cmfs.values, illuminant),
qhull_options='QJ')
xx, yy = np.meshgrid(np.linspace(0, 1, samples),
np.linspace(0, 1, samples))
xy = tstack((xx, yy))
xy = xy[triangulation.find_simplex(xy) > 0]
XYZ = xy_to_XYZ(xy)
RGB = normalise(XYZ_to_sRGB(XYZ, illuminant), axis=-1)
x_dot, y_dot = tsplit(xy)
pylab.scatter(x_dot, y_dot, color=RGB, s=surface)
settings.update({
'x_ticker': False,
'y_ticker': False,
'bounding_box': (0, 1, 0, 1),
'bbox_inches': 'tight',
'pad_inches': 0
})
settings.update(kwargs)
ax = matplotlib.pyplot.gca()
matplotlib.pyplot.setp(ax, frame_on=False)
boundaries(**settings)
decorate(**settings)
return display(**settings)
def colour_checker_plot(colour_checker='ColorChecker 2005', **kwargs):
"""
Plots given colour checker.
Parameters
----------
colour_checker : unicode, optional
Color checker name.
\**kwargs : dict, optional
Keywords arguments.
Returns
-------
bool
Definition success.
Raises
------
KeyError
If the given colour rendition chart is not found in the factory colour
rendition charts.
Examples
--------
>>> colour_checker_plot() # doctest: +SKIP
True
"""
canvas(**kwargs)
colour_checker, name = COLOURCHECKERS.get(colour_checker), colour_checker
if colour_checker is None:
raise KeyError(
('Colour checker "{0}" not found in '
'factory colour checkers: "{1}".').format(
name, sorted(COLOURCHECKERS.keys())))
_name, data, illuminant = colour_checker
colour_parameters = []
for _index, label, x, y, Y in data:
XYZ = xyY_to_XYZ((x, y, Y))
RGB = XYZ_to_sRGB(XYZ, illuminant)
colour_parameters.append(
ColourParameter(label.title(), np.clip(np.ravel(RGB), 0, 1)))
background_colour = '0.1'
matplotlib.pyplot.gca().patch.set_facecolor(background_colour)
width = height = 1.0
spacing = 0.25
across = 6
settings = {
'standalone': False,
'width': width,
'height': height,
'spacing': spacing,
'across': across,
'text_size': 8,
'margins': (-0.125, 0.125, -0.5, 0.125)}
settings.update(kwargs)
multi_colour_plot(colour_parameters, **settings)
text_x = width * (across / 2) + (across * (spacing / 2)) - spacing / 2
text_y = -(len(colour_parameters) / across + spacing / 2)
pylab.text(text_x,
text_y,
'{0} - {1} - Colour Rendition Chart'.format(
name, RGB_COLOURSPACES.get('sRGB').name),
color='0.95',
clip_on=True,
ha='center')
settings.update({
'title': name,
'facecolor': background_colour,
'edgecolor': None,
'standalone': True})
boundaries(**settings)
decorate(**settings)
return display(**settings)
def single_spd_plot(spd,
cmfs='CIE 1931 2 Degree Standard Observer',
out_of_gamut_clipping=True,
**kwargs):
"""
Plots given spectral power distribution.
Parameters
----------
spd : SpectralPowerDistribution
Spectral power distribution to plot.
out_of_gamut_clipping : bool, optional
Whether to clip out of gamut colours otherwise, the colours will be
offset by the absolute minimal colour leading to a rendering on
gray background, less saturated and smoother.
cmfs : unicode
Standard observer colour matching functions used for spectrum creation.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`colour.plotting.render`},
Please refer to the documentation of the previously listed definition.
Returns
-------
Figure
Current figure or None.
References
----------
- :cite:`Spiker2015a`
Examples
--------
>>> from colour import SpectralPowerDistribution
>>> data = {
... 500: 0.0651,
... 520: 0.0705,
... 540: 0.0772,
... 560: 0.0870,
... 580: 0.1128,
... 600: 0.1360
... }
>>> spd = SpectralPowerDistribution(data, name='Custom')
>>> single_spd_plot(spd) # doctest: +SKIP
"""
axes = canvas(**kwargs).gca()
cmfs = get_cmfs(cmfs)
spd = spd.copy()
spd.interpolator = LinearInterpolator
wavelengths = cmfs.wavelengths[np.logical_and(
cmfs.wavelengths >= max(min(cmfs.wavelengths), min(spd.wavelengths)),
cmfs.wavelengths <= min(max(cmfs.wavelengths), max(spd.wavelengths)),
)]
values = spd[wavelengths]
colours = XYZ_to_sRGB(
wavelength_to_XYZ(wavelengths, cmfs),
ILLUMINANTS['CIE 1931 2 Degree Standard Observer']['E'],
apply_encoding_cctf=False)
if not out_of_gamut_clipping:
colours += np.abs(np.min(colours))
colours = DEFAULT_PLOTTING_ENCODING_CCTF(normalise_maximum(colours))
x_min, x_max = min(wavelengths), max(wavelengths)
y_min, y_max = 0, max(values)
polygon = Polygon(np.vstack([
(x_min, 0),
tstack((wavelengths, values)),
(x_max, 0),
]),
facecolor='none',
edgecolor='none')
axes.add_patch(polygon)
axes.bar(x=wavelengths,
height=max(values),
width=1,
color=colours,
align='edge',
clip_path=polygon)
axes.plot(wavelengths, values, color='black', linewidth=1)
settings = {
'title': '{0} - {1}'.format(spd.strict_name, cmfs.strict_name),
'x_label': 'Wavelength $\\lambda$ (nm)',
'y_label': 'Spectral Power Distribution',
'limits': (x_min, x_max, y_min, y_max),
'x_tighten': True,
'y_tighten': True
}
settings.update(kwargs)
return render(**settings)
def CIE_1931_chromaticity_diagram_colours_plot(
surface=1,
samples=4096,
cmfs='CIE 1931 2 Degree Standard Observer',
**kwargs):
"""
Plots the *CIE 1931 Chromaticity Diagram* colours.
Parameters
----------
surface : numeric, optional
Generated markers surface.
samples : numeric, optional
Samples count on one axis.
cmfs : unicode, optional
Standard observer colour matching functions used for diagram bounds.
Other Parameters
----------------
\**kwargs : dict, optional
{:func:`boundaries`, :func:`canvas`, :func:`decorate`,
:func:`display`},
Please refer to the documentation of the previously listed definitions.
Returns
-------
Figure
Current figure or None.
Examples
--------
>>> CIE_1931_chromaticity_diagram_colours_plot() # doctest: +SKIP
"""
settings = {'figure_size': (64, 64)}
settings.update(kwargs)
canvas(**settings)
cmfs = get_cmfs(cmfs)
illuminant = DEFAULT_PLOTTING_ILLUMINANT
triangulation = Delaunay(XYZ_to_xy(cmfs.values, illuminant),
qhull_options='QJ')
xx, yy = np.meshgrid(np.linspace(0, 1, samples),
np.linspace(0, 1, samples))
xy = tstack((xx, yy))
xy = xy[triangulation.find_simplex(xy) > 0]
XYZ = xy_to_XYZ(xy)
RGB = normalise_maximum(XYZ_to_sRGB(XYZ, illuminant), axis=-1)
x_dot, y_dot = tsplit(xy)
pylab.scatter(x_dot, y_dot, color=RGB, s=surface)
settings.update({
'x_ticker': False,
'y_ticker': False,
'bounding_box': (0, 1, 0, 1)
})
settings.update(kwargs)
ax = matplotlib.pyplot.gca()
matplotlib.pyplot.setp(ax, frame_on=False)
boundaries(**settings)
decorate(**settings)
return display(**settings)
def CIE_1931_chromaticity_diagram_colours_plot(
surface=1.25,
spacing=0.00075,
cmfs='CIE 1931 2 Degree Standard Observer',
**kwargs):
"""
Plots the *CIE 1931 Chromaticity Diagram* colours.
Parameters
----------
surface : numeric, optional
Generated markers surface.
spacing : numeric, optional
Spacing between markers.
cmfs : unicode, optional
Standard observer colour matching functions used for diagram bounds.
\*\*kwargs : \*\*
Keywords arguments.
Returns
-------
bool
Definition success.
Examples
--------
>>> CIE_1931_chromaticity_diagram_colours_plot() # doctest: +SKIP
True
"""
cmfs, name = get_cmfs(cmfs), cmfs
illuminant = ILLUMINANTS.get('CIE 1931 2 Degree Standard Observer').get(
'E')
XYZs = [value for key, value in cmfs]
x, y = tuple(zip(*([XYZ_to_xy(x) for x in XYZs])))
path = matplotlib.path.Path(tuple(zip(x, y)))
x_dot, y_dot, colours = [], [], []
for i in np.arange(0, 1, spacing):
for j in np.arange(0, 1, spacing):
if path.contains_path(matplotlib.path.Path([[i, j], [i, j]])):
x_dot.append(i)
y_dot.append(j)
XYZ = xy_to_XYZ((i, j))
RGB = normalise(XYZ_to_sRGB(XYZ, illuminant))
colours.append(RGB)
pylab.scatter(x_dot, y_dot, color=colours, s=surface)
settings = {
'no_ticks': True,
'bounding_box': [0, 1, 0, 1],
'bbox_inches': 'tight',
'pad_inches': 0
}
settings.update(kwargs)
bounding_box(**settings)
aspect(**settings)
return display(**settings)
