def print_new_dead_new_alive():
start_orbit = 21847
end_orbit = 21890
# start_orbit = 22820
# end_orbit = 22890
# start_orbit = 24023
# end_orbit = 24051
# start_orbit = 27229
# end_orbit = 27258
# start_orbit = 32726
# end_orbit = 32870
# start_orbit = 49230
# end_orbit = 49259
m1 = Mask()
m1.apply_pixel_window(394, 620)
m1.load_ascii(start_orbit)
m2 = Mask()
for orbit in range(start_orbit, end_orbit):
m2.load_ascii(orbit)
print(orbit, "new dead:", m1.get_new_dead(m2), "new alive:", m1.get_new_alive(m2))
return
def compare_combined_30_vs_32(orbit):
m30 = Mask()
m30.load_sdmf30_crit(orbit, "combined")
m32 = Mask()
m32.load_sdmf32_mask(orbit, "combinedFlag")
# for i in range(1024):
# print(i, m30.mask[i], m32.mask[i])
newdead = m30.get_new_dead(m32)
newalive = m30.get_new_alive(m32)
print("3.0->3.2 dead:", newdead.size, newdead)
print("3.0->3.2 alive:", newalive.size, newalive)
return
def compare_error_30_vs_32(orbit):
m30 = Mask()
m30_combined = Mask()
m30_combined.load_ascii(orbit)
m30.load_sdmf30_crit(orbit, "residual")
m32 = Mask()
m32.load_sdmf32_figure(orbit, "darkError")
newdead = m30.get_new_dead(m32)
newalive = m30.get_new_alive(m32)
print("new errory:", newdead.size, newdead)
print("new errory but already bad", np.sum(np.in1d(newdead, np.where(m30_combined.mask))))
print("real new errory", newdead[np.in1d(newdead, np.where(m30_combined.mask), invert=True)])
print("new not-errory:", newalive.size, newalive)
print("new not-errory but still bad", np.sum(np.in1d(newalive, np.where(m30_combined.mask))))
print("real new not-errory:", newalive[np.in1d(newalive, np.where(m30_combined.mask), invert=True)])
return
