Add another hack script.

[dcpomatic.git] / hacks / gamma_graphs.py

diff --git a/hacks/gamma_graphs.py b/hacks/gamma_graphs.py
index 166a0382858d2ccc58fd5a16df2aefd1c6e3a181..ad75ba908cf48f4a91db8784533546ceb7c7267a 100644 (file)
--- a/hacks/gamma_graphs.py
+++ b/hacks/gamma_graphs.py
@@ -6,14 +6,24 @@ import numpy as np
 x = np.logspace(-3, 0, 100)
 
 plt.loglog(x, pow(x, 2.2))
-plt.loglog(x, pow(x, 2.4))
+# plt.loglog(x, pow(x, 2.4))
 
-linearised = []
+srgb_linearised = []
 for xx in x:
     if xx > 0.04045:
-        linearised.append(pow((xx + 0.055) / 1.055, 2.4))
+        srgb_linearised.append(pow((xx + 0.055) / 1.055, 2.4))
     else:
-        linearised.append(xx / 12.92)
+        srgb_linearised.append(xx / 12.92)
+
+# plt.loglog(x, srgb_linearised)
+
+rec_linearised = []
+for xx in x:
+    if xx > 0.081:
+        rec_linearised.append(pow((xx + 0.099) / 1.099, 1 / 0.45))
+    else:
+        rec_linearised.append(xx / 4.5)
+
+plt.loglog(x, rec_linearised)
 
-plt.loglog(x, linearised)
 plt.show()