exclude the last order

hypergravity · Jan 17, 2021 · bf9daec · bf9daec
1 parent f010aa7
commit bf9daec
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/E9G10/bfosc_pipeline.py b/E9G10/bfosc_pipeline.py
@@ -164,11 +164,13 @@ def read_star(fp_star):
     if np.any(np.abs(intersect[i] - intersect[i + 1:]) < 5):
         ind_valid_ap[i] = False
 ap_interp = ap_illusion[ind_valid_ap]
+ap_interp = ap_interp[ap_interp[:, -1] < 1070]
 try:
     assert ap_interp.shape[0] == 12
 except:
-    print("invalid number of apertures, please check images")
+    print("invalid number of apertures, please check images!")
     print("ap_interp", ap_interp.shape)
+    print("Note that the FEAT template has 11 orders!")
 ap_interp = ap_interp[1:]
 print("  |- final aperture shape", ap_interp.shape)