from omg_dosimetry import LUT
from omg_dosimetry.i_o import retrieve_demo_file
import matplotlib.pyplot as plt

img_path_1 = retrieve_demo_file("C14_calib-18h-1_001.tif")
img_path_2 = retrieve_demo_file("C14_calib-18h-2_001.tif")
my_path = img_path_1.parent
doses = [0.0, 100.0, 200.0, 400.0, 650.0, 950.0]
lut = LUT(my_path, doses, crop_top_bottom = 650) # Crop needed because a glass on the scanner
lut.plot_fit()
plt.show()