1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
import numpy as np
from colour import *
from colour.difference import delta_E_CIE1976
from colour.models import eotf_inverse_sRGB
from gsoc_common import D65, D65_xy, jakob_hanika, model_sd, plot_comparison
# These numbers are taken from Jakob and Hanika's Jupyter notebook.
RGB_ref = np.array([0.79264853, 0.4, 0.63703843]) # *linear* sRGB
cc_ref = np.array([ 18.70184886, -13.19804478, 2.12180137])
if __name__ == "__main__":
#XYZ = sRGB_to_XYZ(eotf_inverse_sRGB(RGB_ref), D65_xy)
XYZ = np.array([1/3, 1/4, 1/3])
Lab = XYZ_to_Lab(XYZ, D65_xy)
print("Target: X=%g, Y=%g, Z=%g, L=%g, a=%g, b=%g" % (*XYZ, *Lab))
ccp, error = jakob_hanika(XYZ, D65, D65_xy)
reference_sd = model_sd(cc_ref)
reference_XYZ = sd_to_XYZ(reference_sd, illuminant=D65)
reference_Lab = XYZ_to_Lab(reference_XYZ, D65_xy)
found_sd = model_sd(ccp)
found_XYZ = sd_to_XYZ(found_sd, illuminant=D65)
found_Lab = XYZ_to_Lab(found_XYZ, D65_xy)
print("Our results differ from the reference by ΔE = %g" \
% delta_E_CIE1976(reference_Lab, found_Lab))
plot_comparison(XYZ, found_sd, "Reference", error, D65)
|