y_sky, x_sky = np.where(np.logical_and(bad_solution_mask == 0, np.logical_and(skyMask != 0, deadMask == 1)))
# tic = time()
for sec in np.arange(0, nSecInFile, integrationTime):
jd = (
startJD + sec / (24.0 * 3600.0) + integrationTime / 2.0 / (24.0 * 3600.0)
) # add seconds offset to julian date, move jd to center of bin
print count, jd
count += 1
times.append(jd)
titles.append("%.6f" % jd)
frameData = ob.getAperturePixelCountImage(
firstSec=sec,
integrationTime=integrationTime,
y_values=y_values,
x_values=x_values,
y_sky=y_sky,
x_sky=x_sky,
apertureMask=apertureMask,
skyMask=skyMask,
weighted=True,
)
# frameData = ob.getPixelCountImage(firstSec=sec,integrationTime=integrationTime,weighted=True)
backgroundPerPixel.append(frameData["SkyCountSubtractedPerPixel"])
lightcurve = frameData["lightcurve"]
frame = frameData["image"]
showFrame = np.array(frame)
showframes.append(showFrame)
frame[deadMask == 0] = np.nan
frame[bad_solution_mask == 1] = np.nan
frame[apertureMask == 0] = np.nan
frames.append(frame)
