+ 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)