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)