def rollPhaseTask(clip):
centColRow = clip['centroids.cent_colrow']
flags = clip['extract.flags']
rot = arclen.computeArcLength(centColRow, flags > 0)
rollPhase = rot[:, 0]
rollPhase[flags > 0] = -9999 #A bad value
clip['rollPhase'] = {'rollPhase': rollPhase}
return clip
def rollPhaseTask(clip):
centColRow = clip['centroids.cent_colrow']
flags = clip['extract.flags']
rot = arclen.computeArcLength(centColRow, flags>0)
rollPhase = rot[:,0]
rollPhase[flags>0] = -9999 #A bad value
clip['rollPhase'] = {'rollPhase':rollPhase}
return clip
def exampleDiffImgCentroiding():
k2id = 206103150
campaign = 3
ar = mastio.K2Archive()
fits, hdr = ar.getLongTpf(k2id, campaign, header=True)
hdr0 = ar.getLongTpf(k2id, campaign, ext=0)
cube = tpf.getTargetPixelArrayFromFits(fits, hdr)
idx = np.isfinite(cube)
cube[~idx] = 0 #Remove Nans
flags = fits['QUALITY']
ccdMod = hdr0['module']
ccdOut = hdr0['output']
#Compute roll phase
llc = ar.getLongCadence(k2id, campaign)
time= llc['TIME']
cent1 = llc['MOM_CENTR1']
cent2 = llc['MOM_CENTR2']
centColRow = np.vstack((cent1, cent2)).transpose()
rot = arclen.computeArcLength(centColRow, flags>0)
rollPhase = rot[:,0]
rollPhase[flags>0] = -9999 #A bad value
prfObj = prf.KeplerPrf("/home/fergal/data/keplerprf")
bbox = getBoundingBoxForImage(cube[0], hdr)
period = 4.1591409
epoch = fits['time'][491]
dur = 3.0
out, log = measureDiffOffset(period, epoch, dur, time, prfObj, \
ccdMod, ccdOut, cube, bbox, rollPhase, flags)
idx = out[:,1] > 0
mp.clf()
mp.plot(out[:,3]-out[:,1], out[:,4]- out[:,2], 'ro')
return out
def exampleDiffImgCentroiding():
k2id = 206103150
campaign = 3
ar = mastio.K2Archive()
fits, hdr = ar.getLongTpf(k2id, campaign, header=True)
hdr0 = ar.getLongTpf(k2id, campaign, ext=0)
cube = tpf.getTargetPixelArrayFromFits(fits, hdr)
idx = np.isfinite(cube)
cube[~idx] = 0 #Remove Nans
flags = fits['QUALITY']
ccdMod = hdr0['module']
ccdOut = hdr0['output']
#Compute roll phase
llc = ar.getLongCadence(k2id, campaign)
time = llc['TIME']
cent1 = llc['MOM_CENTR1']
cent2 = llc['MOM_CENTR2']
centColRow = np.vstack((cent1, cent2)).transpose()
rot = arclen.computeArcLength(centColRow, flags > 0)
rollPhase = rot[:, 0]
rollPhase[flags > 0] = -9999 #A bad value
prfObj = prf.KeplerPrf("/home/fergal/data/keplerprf")
bbox = getBoundingBoxForImage(cube[0], hdr)
period = 4.1591409
epoch = fits['time'][491]
dur = 3.0
out, log = measureDiffOffset(period, epoch, dur, time, prfObj, \
ccdMod, ccdOut, cube, bbox, rollPhase, flags)
idx = out[:, 1] > 0
mp.clf()
mp.plot(out[:, 3] - out[:, 1], out[:, 4] - out[:, 2], 'ro')
return out
def example():
ar = mastio.K2Archive()
kepid = 206103150 #A wasp planet
fits = ar.getLongCadence(kepid, 3)
flags = fits['SAP_QUALITY']
cent1 = fits['MOM_CENTR1']
cent2 = fits['MOM_CENTR2']
fits, hdr = ar.getLongTpf(kepid, 3, header=True)
cube = tpf.getTargetPixelArrayFromFits(fits, hdr)
# cube *= gain
#Compute roll phase
centColRow = np.vstack((cent1, cent2)).transpose()
rot = arclen.computeArcLength(centColRow, flags > 0)
rollPhase = rot[:, 0]
rollPhase[flags > 0] = -9999 #A bad value
cadenceInTransit = 490
diff, oot = constructK2DifferenceImage(cube, cadenceInTransit, rollPhase,
flags)
return diff, oot
