From 7b94b07c13c6757546dd416073305650938d9820 Mon Sep 17 00:00:00 2001 From: Kristoffer Gryte Date: Fri, 26 Sep 2025 12:23:54 +0200 Subject: [PATCH] Correct interpolation and dark pixel offset --- ramses_calibrate.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/ramses_calibrate.py b/ramses_calibrate.py index 5ae6856..e08be7e 100644 --- a/ramses_calibrate.py +++ b/ramses_calibrate.py @@ -199,7 +199,7 @@ def raw2cal_Air(spec, msdate, serialn, if len(calsensind) > 1: calsensdates = [SAMDateTime_Air[i] for i in calsensind] tdeltas = [(msdate-x) for x in calsensdates - if (msdate-x).days > 0] + if (msdate-x).days > 0] calsensind = [calsensind[tdeltas.index(min(tdeltas))]] Cal = CalData[calsensind[0]] @@ -209,8 +209,8 @@ def raw2cal_Air(spec, msdate, serialn, dp1 = Cal.ini.DarkPixelStart dp2 = Cal.ini.DarkPixelStop wave = [0.0] * 256 - for i in range(1, len(wave)+1, 1): - wave[i-1] = (Cal.ini.c0s) + (Cal.ini.c1s*(i+1)) +\ + for i in range(0, len(wave), 1): + wave[i] = (Cal.ini.c0s) + (Cal.ini.c1s*(i+1)) +\ (Cal.ini.c2s*(i+1)**2) + (Cal.ini.c3s*(i+1)**3) t0 = 8192 @@ -222,7 +222,7 @@ def raw2cal_Air(spec, msdate, serialn, B = B0 + (t1/t0*B1) C = M - B # subtract dark offset, - Offset = np.mean(C[dp1-1:dp2]) # dark pixels + Offset = np.mean(C[dp1:dp2+1]) # dark pixels D = C-Offset E = D*(t0/t1) # Scale the spectrum to the Air calibration