#!/usr/bin/env python3
# Copyright © 2014 Mattias Andrée (maandree@member.fsf.org)
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU Affero General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Affero General Public License for more details.
#
# You should have received a copy of the GNU Affero General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
import math
from colour import *
DATADIR = '.'
i_size = 2 ** 8
o_size = 2 ** 16
r_curve = [i / (i_size - 1) for i in range(i_size)]
g_curve = [i / (i_size - 1) for i in range(i_size)]
b_curve = [i / (i_size - 1) for i in range(i_size)]
clip_result = True
cmf_2deg_cache = None
cmf_10deg_cache = None
def curves(r, g, b):
'''
Generate a tuple of curve–parameter pairs
@param r The red parameter
@param g The green parameter
@param b The blue parameter
@return ((r_curve, r), (g_curve, g), (b_curve, b))
'''
return ((r_curve, r), (g_curve, g), (b_curve, b))
def series_d(temperature):
'''
Calculate the colour for a blackbody temperature
@param temperature:float The blackbody temperature in kelvins, must be inside [4000, 7000]
@return :(float, float, float) The red, green and blue components of the white point
'''
x = 0
ks = ((0.244063, 0), (0.09911, 1), (2.9678, 2), (-4.6070, 3))
if temperature > 7000:
ks = ((0.237040, 0), (0.24748, 1), (1.9018, 2), (-2.0064, 3))
for (k, d) in ks:
x += k * 10 ** (d * 3) / temperature ** d
y = 2.870 * x - 3.000 * x ** 2 - 0.275
return to_srgb(x, y, 1.0)
def simple_whitepoint(temperature):
'''
Calculate the colour for a blackbody temperature using a simple, but inaccurate, algorithm
@param temperature:float The blackbody temperature in kelvins, not guaranteed for values outside [1000, 40000]
@return :(float, float, float) The red, green and blue components of the white point
'''
r, g, b = 1, 1, 1
temp = temperature / 100
if temp > 66:
temp -= 60
r = 1.292936186 * temp ** 0.1332047592
g = 1.129890861 * temp ** -0.0755148492
else:
g = 0.390081579 * math.log(temp) - 0.631841444
if temp <= 19:
b = 0
elif temp < 66:
b = 0.543206789 * math.log(temp - 10) - 1.196254089
return (r, g, b)
def cmf_2deg(temperature):
'''
Calculate the colour for a blackbody temperature using raw CIE 1931 2 degree CMF data with interpolation
@param temperature:float The blackbody temperature in kelvins, clipped to [1000, 40000]
@return :(float, float, float) The red, green and blue components of the white point
'''
if cmf_2deg_cache is None:
with open(DATADIR + '/2deg', 'rb') as file:
cmf_2deg_cache = file.read()
cmf_2deg_cache.decode('utf-8', 'error').split('\n')
cmf_2deg_cache = [[float(x) for x in x_y.split(' ')] for x_y in cmf_2deg_cache]
temperature = min(max(0, temperature), 1000)
x, y = 0, 0
if (temp % 100) == 0:
(x, y) = temperature[(temp - 1000) // 100]
else:
temp -= 1000
(x1, y1) = temperature[temp // 100]
(x2, y2) = temperature[temp // 100 + 1]
temp = (temp % 100) / 100
x = x1 * temp + x2 * (1 - temp)
y = y1 * temp + y2 * (1 - temp)
return to_srgb(x, y, 1.0)
def cmf_10deg(temperature):
'''
Calculate the colour for a blackbody temperature using raw CIE 1964 10 degree CMF data with interpolation
@param temperature:float The blackbody temperature in kelvins, clipped to [1000, 40000]
@return :(float, float, float) The red, green and blue components of the white point
'''
if cmf_2deg_cache is None:
with open(DATADIR + '/10deg', 'rb') as file:
cmf_2deg_cache = file.read()
cmf_2deg_cache.decode('utf-8', 'error').split('\n')
cmf_2deg_cache = [[float(x) for x in x_y.split(' ')] for x_y in cmf_2deg_cache]
temperature = min(max(0, temperature), 1000)
x, y = 0, 0
if (temp % 100) == 0:
(x, y) = temperature[(temp - 1000) // 100]
else:
temp -= 1000
(x1, y1) = temperature[temp // 100]
(x2, y2) = temperature[temp // 100 + 1]
temp = (temp % 100) / 100
x = x1 * temp + x2 * (1 - temp)
y = y1 * temp + y2 * (1 - temp)
return to_srgb(x, y, 1.0)
def temperature(temperature, algorithm, linear_rgb = True):
'''
Change colour temperature according to the CIE illuminant series D
@param temperature:float The blackbody temperature in kelvins
@param algorithm:(float)→(float, float, float) Algorithm for calculating a white point, for example `series_d` or `simple_whitepoint`
@param linear_rgb:[bool] Whether to use linear RGB, otherwise sRG is used
'''
if temperature == 6500:
return
(r, g, b) = algorithm(temperature)
if linear_rgb:
for curve in (r_curve, g_curve, b_curve):
for i in range(i_size):
R, G, B = r_curve[i], g_curve[i], b_curve[i]
(R, G, B) = standard_to_linear(R, G, B)
r_curve[i], g_curve[i], b_curve[i] = R, G, B
rgb_brightness(r, g, b)
if linear_rgb:
for curve in (r_curve, g_curve, b_curve):
for i in range(i_size):
R, G, B = r_curve[i], g_curve[i], b_curve[i]
(R, G, B) = linear_to_standard(R, G, B)
r_curve[i], g_curve[i], b_curve[i] = R, G, B
def divide_by_maximum():
'''
Divide all colour components by the value of the most prominent colour component for each colour
'''
for i in range(i_size):
m = max([abs(x) for x in (r_curve[i], g_curve[i], b_curve[i])])
if m != 0:
for curve in (r_curve, g_curve, b_curve):
curve[i] /= m
def rgb_contrast(r, g, b):
'''
Apply contrast correction on the colour curves using sRGB
@param r:float The contrast parameter for the red curve
@param g:float The contrast parameter for the green curve
@param b:float The contrast parameter for the blue curve
'''
for (curve, level) in curves(r, g, b):
if not level == 1.0:
for i in range(i_size):
curve[i] = (curve[i] - 0.5) * level + 0.5
def cie_contrast(level):
'''
Apply contrast correction on the colour curves using CIE XYZ
@param level:float The brightness parameter
'''
if not level == 1.0:
for i in range(i_size):
(x, y, Y) = to_ciexyy(r_curve[i], g_curve[i], b_curve[i])
(r_curve[i], g_curve[i], b_curve[i]) = to_rgb(x, y, Y * level)
def rgb_brightness(r, g, b):
'''
Apply brightness correction on the colour curves using sRGB
@param r:float The brightness parameter for the red curve
@param g:float The brightness parameter for the green curve
@param b:float The brightness parameter for the blue curve
'''
for (curve, level) in curves(r, g, b):
if not level == 1.0:
for i in range(i_size):
curve[i] *= level
def cie_brightness(level):
'''
Apply brightness correction on the colour curves using CIE XYZ
@param level:float The brightness parameter
'''
if not level == 1.0:
for i in range(i_size):
(x, y, Y) = to_ciexyy(r_curve[i], g_curve[i], b_curve[i])
(r_curve[i], g_curve[i], b_curve[i]) = to_rgb(x, y, Y * level)
def gamma(r, g, b):
'''
Apply gamma correction on the colour curves
@param r:float The gamma parameter for the red curve
@param g:float The gamma parameter for the green curve
@param b:float The gamma parameter for the blue curve
'''
for (curve, level) in curves(r, g, b):
if not level == 1.0:
for i in range(i_size):
curve[i] **= level
def sigmoid(r, g, b):
'''
Apply S-curve correction on the colour curves
@param r:float? The sigmoid parameter for the red curve
@param g:float? The sigmoid parameter for the green curve
@param b:float? The sigmoid parameter for the blue curve
'''
for (curve, level) in curves(r, g, b):
if level is not None:
for i in range(i_size):
try:
curve[i] = 0.5 - math.log(1 / curve[i] - 1) / level
except:
curve[i] = 0;
def clip():
'''
Clip all values belowed the actual minimum and above actual maximums
'''
for curve in (r_curve, g_curve, b_curve):
for i in range(i_size):
curve[i] = min(max(0.0, curve[i]), 1.0)
temperature(6500, series_d, True)
divide_by_maximum()
temperature(6500, simple_whitepoint, True)
clip()
rgb_contrast(1.0, 1.0, 1.0)
cie_contrast(1.0)
rgb_brightness(1.0, 1.0, 1.0)
cie_brightness(1.0)
gamma(1.0, 1.0, 1.0)
sigmoid(None, None, None)
clip()
for curve in (r_curve, g_curve, b_curve):
for i in range(i_size):
curve[i] = int(curve[i] * (o_size - 1) + 0.5)
if clip_result:
curve[i] = min(max(0, curve[i]), (o_size - 1))
print(r_curve)
print(g_curve)
print(b_curve)
