{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Kepler Data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## The Data" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "SIMBAD info: http://simbad.u-strasbg.fr/simbad/sim-id?Ident=KIC7198959\n", "\n", "Lightcurve data from: https://archive.stsci.edu/kepler/publiclightcurves.html" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# !curl -O http://archive.stsci.edu/pub/kepler/lightcurves/0071/007198959/kplr007198959-2009259160929_llc.fits" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Filename: kplr007198959-2009259160929_llc.fits\n", "No. Name Type Cards Dimensions Format\n", " 0 PRIMARY PrimaryHDU 58 () \n", " 1 LIGHTCURVE BinTableHDU 155 4354R x 20C [D, E, J, E, E, E, E, E, E, J, D, E, D, E, D, E, D, E, E, E] \n", " 2 APERTURE ImageHDU 48 (12, 41) int32 \n" ] } ], "source": [ "from astropy.io import fits\n", "hdulist = fits.open('kplr007198959-2009259160929_llc.fits')\n", "hdulist.info()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "XTENSION= 'BINTABLE' / marks the beginning of a new HDU \n", "BITPIX = 8 / array data type \n", "NAXIS = 2 / number of array dimensions \n", "NAXIS1 = 100 / length of first array dimension \n", "NAXIS2 = 4354 / length of second array dimension \n", "PCOUNT = 0 / group parameter count (not used) \n", "GCOUNT = 1 / group count (not used) \n", "TFIELDS = 20 / number of table fields \n", "TTYPE1 = 'TIME ' / column title: data time stamps \n", "TFORM1 = 'D ' / column format: 64-bit floating point \n", "TUNIT1 = 'BJD - 2454833' / column units: barycenter corrected JD \n", "TDISP1 = 'D14.7 ' / column display format \n", "TTYPE2 = 'TIMECORR' / column title: barycenter - timeslice correction\n", "TFORM2 = 'E ' / column format: 32-bit floating point \n", "TUNIT2 = 'd ' / column units: day \n", "TDISP2 = 'E13.6 ' / column display format \n", "TTYPE3 = 'CADENCENO' / column title: unique cadence number \n", "TFORM3 = 'J ' / column format: signed 32-bit integer \n", "TDISP3 = 'I10 ' / column display format \n", "TTYPE4 = 'SAP_FLUX' / column title: aperture photometry flux \n", "TFORM4 = 'E ' / column format: 32-bit floating point \n", "TUNIT4 = 'e-/s ' / column units: electrons per second \n", "TDISP4 = 'E14.7 ' / column display format \n", "TTYPE5 = 'SAP_FLUX_ERR' / column title: aperture phot. flux error \n", "TFORM5 = 'E ' / column format: 32-bit floating point \n", "TUNIT5 = 'e-/s ' / column units: electrons per second (1-sigma) \n", "TDISP5 = 'E14.7 ' / column display format \n", "TTYPE6 = 'SAP_BKG ' / column title: aperture phot. background flux \n", "TFORM6 = 'E ' / column format: 32-bit floating point \n", "TUNIT6 = 'e-/s ' / column units: electrons per second \n", "TDISP6 = 'E14.7 ' / column display format \n", "TTYPE7 = 'SAP_BKG_ERR' / column title: ap. phot. background flux error \n", "TFORM7 = 'E ' / column format: 32-bit floating point \n", "TUNIT7 = 'e-/s ' / column units: electrons per second (1-sigma) \n", "TDISP7 = 'E14.7 ' / column display format \n", "TTYPE8 = 'PDCSAP_FLUX' / column title: aperture phot. PDC flux \n", "TFORM8 = 'E ' / column format: 32-bit floating point \n", "TUNIT8 = 'e-/s ' / column units: electrons per second \n", "TDISP8 = 'E14.7 ' / column display format \n", "TTYPE9 = 'PDCSAP_FLUX_ERR' / column title: ap. phot. PDC flux error \n", "TFORM9 = 'E ' / column format: 32-bit floating point \n", "TUNIT9 = 'e-/s ' / column units: electrons per second (1-sigma) \n", "TDISP9 = 'E14.7 ' / column display format \n", "TTYPE10 = 'SAP_QUALITY' / column title: aperture photometry quality flag \n", "TFORM10 = 'J ' / column format: signed 32-bit integer \n", "TDISP10 = 'B16.16 ' / column display format \n", "TTYPE11 = 'PSF_CENTR1' / column title: PSF-fitted column centroid \n", "TFORM11 = 'D ' / column format: 64-bit floating point \n", "TUNIT11 = 'pixel ' / column units: pixel \n", "TDISP11 = 'F10.5 ' / column display format \n", "TTYPE12 = 'PSF_CENTR1_ERR' / column title: PSF-fitted column error \n", "TFORM12 = 'E ' / column format: 32-bit floating point \n", "TUNIT12 = 'pixel ' / column units: pixel (1-sigma) \n", "TDISP12 = 'E14.7 ' / column display format \n", "TTYPE13 = 'PSF_CENTR2' / column title: PSF-fitted row centroid \n", "TFORM13 = 'D ' / column format: 64-bit floating point \n", "TUNIT13 = 'pixel ' / column units: pixel \n", "TDISP13 = 'F10.5 ' / column display format \n", "TTYPE14 = 'PSF_CENTR2_ERR' / column title: PSF-fitted row error \n", "TFORM14 = 'E ' / column format: 32-bit floating point \n", "TUNIT14 = 'pixel ' / column units: pixel (1-sigma) \n", "TDISP14 = 'E14.7 ' / column display format \n", "TTYPE15 = 'MOM_CENTR1' / column title: moment-derived column centroid \n", "TFORM15 = 'D ' / column format: 64-bit floating point \n", "TUNIT15 = 'pixel ' / column units: pixel \n", "TDISP15 = 'F10.5 ' / column display format \n", "TTYPE16 = 'MOM_CENTR1_ERR' / column title: moment-derived column error \n", "TFORM16 = 'E ' / column format: 32-bit floating point \n", "TUNIT16 = 'pixel ' / column units: pixel (1-sigma) \n", "TDISP16 = 'E14.7 ' / column display format \n", "TTYPE17 = 'MOM_CENTR2' / column title: moment-derived row centroid \n", "TFORM17 = 'D ' / column format: 64-bit floating point \n", "TUNIT17 = 'pixel ' / column units: pixel \n", "TDISP17 = 'F10.5 ' / column display format \n", "TTYPE18 = 'MOM_CENTR2_ERR' / column title: moment-derived row error \n", "TFORM18 = 'E ' / column format: 32-bit floating point \n", "TUNIT18 = 'pixel ' / column units: pixel (1-sigma) \n", "TDISP18 = 'E14.7 ' / column display format \n", "TTYPE19 = 'POS_CORR1' / column title: column position correction \n", "TFORM19 = 'E ' / column format: 32-bit floating point \n", "TUNIT19 = 'pixels ' / column units: pixel \n", "TDISP19 = 'E14.7 ' / column display format \n", "TTYPE20 = 'POS_CORR2' / column title: row position correction \n", "TFORM20 = 'E ' / column format: 32-bit floating point \n", "TUNIT20 = 'pixels ' / column units: pixel \n", "TDISP20 = 'E14.7 ' / column display format \n", "INHERIT = T / inherit the primary header \n", "EXTNAME = 'LIGHTCURVE' / name of extension \n", "EXTVER = 1 / extension version number (not format version) \n", "TELESCOP= 'Kepler ' / telescope \n", "INSTRUME= 'Kepler Photometer' / detector type \n", "OBJECT = 'KIC 7198959' / string version of target id \n", "KEPLERID= 7198959 / unique Kepler target identifier \n", "RADESYS = 'ICRS ' / reference frame of celestial coordinates \n", "RA_OBJ = 291.366301 / [deg] right ascension \n", "DEC_OBJ = 42.784367 / [deg] declination \n", "EQUINOX = 2000.0 / equinox of celestial coordinate system \n", "EXPOSURE= 81.90687481 / [d] time on source \n", "TIMEREF = 'SOLARSYSTEM' / barycentric correction applied to times \n", "TASSIGN = 'SPACECRAFT' / where time is assigned \n", "TIMESYS = 'TDB ' / time system is barycentric JD \n", "BJDREFI = 2454833 / integer part of BJD reference date \n", "BJDREFF = 0.00000000 / fraction of the day in BJD reference date \n", "TIMEUNIT= 'd ' / time unit for TIME, TSTART and TSTOP \n", "TELAPSE = 88.96781229 / [d] TSTOP - TSTART \n", "LIVETIME= 81.90687481 / [d] TELAPSE multiplied by DEADC \n", "TSTART = 169.51002764 / observation start time in BJD-BJDREF \n", "TSTOP = 258.47783993 / observation stop time in BJD-BJDREF \n", "LC_START= 55002.01747542 / mid point of first cadence in MJD \n", "LC_END = 55090.96492052 / mid point of last cadence in MJD \n", "DEADC = 0.92063492 / deadtime correction \n", "TIMEPIXR= 0.5 / bin time beginning=0 middle=0.5 end=1 \n", "TIERRELA= 5.78E-07 / [d] relative time error \n", "TIERABSO= / [d] absolute time error \n", "INT_TIME= 6.019802903270 / [s] photon accumulation time per frame \n", "READTIME= 0.518948526144 / [s] readout time per frame \n", "FRAMETIM= 6.538751429414 / [s] frame time (INT_TIME + READTIME) \n", "NUM_FRM = 270 / number of frames per time stamp \n", "TIMEDEL = 0.02043359821692 / [d] time resolution of data \n", "DATE-OBS= '2009-06-20T00:10:27.146Z' / TSTART as UTC calendar date \n", "DATE-END= '2009-09-16T23:24:11.862Z' / TSTOP as UTC calendar date \n", "BACKAPP = T / background is subtracted \n", "DEADAPP = T / deadtime applied \n", "VIGNAPP = T / vignetting or collimator correction applied \n", "GAIN = 111.9 / [electrons/count] channel gain \n", "READNOIS= 80.176350 / [electrons] read noise \n", "NREADOUT= 270 / number of read per cadence \n", "TIMSLICE= 2 / time-slice readout sequence section \n", "MEANBLCK= 728 / [count] FSW mean black level \n", "LCFXDOFF= 419400 / long cadence fixed offset \n", "SCFXDOFF= 219400 / short cadence fixed offset \n", "CDPP3_0 = 4018.13134765625 / [ppm] RMS CDPP on 3.0-hr time scales \n", "CDPP6_0 = 3200.025146484375 / [ppm] RMS CDPP on 6.0-hr time scales \n", "CDPP12_0= 2240.04345703125 / [ppm] RMS CDPP on 12.0-hr time scales \n", "CROWDSAP= 0.9991 / Ratio of target flux to total flux in op. ap. \n", "FLFRCSAP= 0.9981 / Frac. of target flux w/in the op. aperture \n", "NSPSDDET= 0 / Number of SPSDs detected \n", "NSPSDCOR= 0 / Number of SPSDs corrected \n", "PDCVAR = 8497.4873046875 / Target variability \n", "PDCMETHD= 'regularMap' / PDC algorithm used for target \n", "NUMBAND = 1 / Number of scale bands \n", "FITTYPE1= 'prior ' / Fit type used for band 1 \n", "PR_GOOD1= 0.9999999403953552 / Prior goodness for band 1 \n", "PR_WGHT1= 7.222428E7 / Prior weight for band 1 \n", "PDC_TOT = 0.9671307802200317 / PDC total goodness metric for target \n", "PDC_TOTP= 47.21249771118164 / PDC_TOT percentile compared to mod/out \n", "PDC_COR = 0.9740492701530457 / PDC correlation goodness metric for target \n", "PDC_CORP= 16.63991355895996 / PDC_COR percentile compared to mod/out \n", "PDC_VAR = 0.9556516408920288 / PDC variability goodness metric for target \n", "PDC_VARP= 25.864707946777344 / PDC_VAR percentile compared to mod/out \n", "PDC_NOI = 0.9995011687278748 / PDC noise goodness metric for target \n", "PDC_NOIP= 99.8350601196289 / PDC_NOI percentile compared to mod/out \n", "PDC_EPT = 1.0 / PDC earth point goodness metric for target \n", "PDC_EPTP= 76.93851470947266 / PDC_EPT percentile compared to mod/out \n", "CHECKSUM= '5dE25ZB05bB05ZB0' / HDU checksum updated 2015-09-07T00:10:06Z " ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "hdulist[1].header" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "<Table length=4354>\n", "
| TIME | TIMECORR | CADENCENO | SAP_FLUX | SAP_FLUX_ERR | SAP_BKG | SAP_BKG_ERR | PDCSAP_FLUX | PDCSAP_FLUX_ERR | SAP_QUALITY | PSF_CENTR1 | PSF_CENTR1_ERR | PSF_CENTR2 | PSF_CENTR2_ERR | MOM_CENTR1 | MOM_CENTR1_ERR | MOM_CENTR2 | MOM_CENTR2_ERR | POS_CORR1 | POS_CORR2 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| float64 | float32 | int32 | float32 | float32 | float32 | float32 | float32 | float32 | int32 | float64 | float32 | float64 | float32 | float64 | float32 | float64 | float32 | float32 | float32 |
| 169.520244707 | 0.00276929 | 2965 | 1.04089e+07 | 80.6514 | 17445.8 | 6.81535 | nan | nan | 256 | nan | nan | nan | nan | 658.98490534 | 5.95852e-06 | 48.3103679836 | 6.3673e-05 | -0.00062082 | -0.0417439 |
| 169.540678885 | 0.00276986 | 2966 | 1.02251e+07 | 79.9479 | 17458.6 | 6.81904 | nan | nan | 256 | nan | nan | nan | nan | 658.984756055 | 6.02055e-06 | 48.3645953434 | 6.32755e-05 | -0.000859906 | -0.0414842 |
| 169.561113062 | 0.00277044 | 2967 | 1.00993e+07 | 79.4513 | 17447.1 | 6.81976 | nan | nan | 256 | nan | nan | nan | nan | 658.984589033 | 6.06324e-06 | 48.3997750666 | 6.29868e-05 | -0.00104417 | -0.0411753 |
| 169.581547239 | 0.00277102 | 2968 | 1.01336e+07 | 79.5974 | 17440.8 | 6.83005 | nan | nan | 8576 | nan | nan | nan | nan | 658.984422904 | 6.04838e-06 | 48.390495693 | 6.31236e-05 | -0.00107995 | -0.0406672 |
| 169.601981415 | 0.0027716 | 2969 | 1.02778e+07 | 80.1858 | 17444.1 | 6.82776 | nan | nan | 393600 | nan | nan | nan | nan | 658.984122755 | 5.99424e-06 | 48.3358836774 | 6.35298e-05 | -0.00127112 | -0.0405857 |
| 169.622415492 | 0.00277217 | 2970 | 1.02574e+07 | 80.072 | 17459.8 | 6.82085 | nan | nan | 256 | nan | nan | nan | nan | 658.984584229 | 5.99411e-06 | 48.3514164486 | 6.3458e-05 | -0.00144995 | -0.0403933 |
| 169.642849768 | 0.00277275 | 2971 | 9.98968e+06 | 79.0959 | 17483.7 | 6.82168 | nan | nan | 384 | nan | nan | nan | nan | 658.984344926 | 6.09163e-06 | 48.4532447707 | 6.3e-05 | -0.00170547 | -0.0402648 |
| 169.663283944 | 0.00277332 | 2972 | 9.77098e+06 | 78.2334 | 17483.0 | 6.82184 | nan | nan | 384 | nan | nan | nan | nan | 658.98438983 | 6.17017e-06 | 48.4774463806 | 6.24124e-05 | -0.00160952 | -0.0400706 |
| 169.68371802 | 0.0027739 | 2973 | 1.00747e+07 | 79.3498 | 17492.0 | 6.82596 | nan | nan | 8576 | nan | nan | nan | nan | 658.984164904 | 6.05106e-06 | 48.4005135707 | 6.30872e-05 | -0.00186426 | -0.0400735 |
| 169.704152195 | 0.00277448 | 2974 | 1.16034e+07 | 85.2008 | 17466.9 | 6.84612 | nan | nan | 384 | nan | nan | nan | nan | 658.98393699 | 5.57952e-06 | 48.1491895595 | 6.69151e-05 | -0.00199162 | -0.0393474 |
| ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
| 258.283726632 | 0.00270851 | 7309 | 1.31303e+07 | 90.7647 | 17822.3 | 6.88246 | 1.31696e+07 | 90.8242 | 0 | nan | nan | nan | nan | 658.976566253 | 5.28842e-06 | 48.556622252 | 6.65606e-05 | -0.00338648 | 0.0333565 |
| 258.304159514 | 0.00270789 | 7310 | 1.29471e+07 | 90.1064 | 17815.1 | 6.87775 | 1.29858e+07 | 90.159 | 0 | nan | nan | nan | nan | 658.976785138 | 5.28611e-06 | 48.5489829917 | 6.71985e-05 | -0.00362269 | 0.0343312 |
| 258.324592595 | 0.00270728 | 7311 | 1.27119e+07 | 89.2577 | 17829.2 | 6.86093 | 1.27434e+07 | 89.3091 | 0 | nan | nan | nan | nan | 658.97690318 | 5.31089e-06 | 48.4355448483 | 6.73667e-05 | -0.00367731 | 0.0329374 |
| 258.345025577 | 0.00270666 | 7312 | 1.24702e+07 | 88.3857 | 17808.9 | 6.86379 | 1.25071e+07 | 88.4538 | 0 | nan | nan | nan | nan | 658.976988051 | 5.34708e-06 | 48.2808913315 | 6.72824e-05 | -0.00375214 | 0.0333151 |
| 258.365458458 | 0.00270604 | 7313 | 1.22394e+07 | 87.5396 | 17795.8 | 6.86721 | 1.22863e+07 | 87.6056 | 8192 | nan | nan | nan | nan | 658.976720505 | 5.3801e-06 | 48.1381236156 | 6.72635e-05 | -0.00439123 | 0.0352742 |
| 258.385891439 | 0.00270542 | 7314 | nan | nan | nan | nan | nan | nan | 32800 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan |
| 258.406324519 | 0.0027048 | 7315 | 1.16802e+07 | 85.472 | 17810.2 | 6.84704 | 1.17144e+07 | 85.5362 | 0 | nan | nan | nan | nan | 658.978009467 | 5.4965e-06 | 48.0239072802 | 6.65793e-05 | -0.00313238 | 0.0327484 |
| 258.4267574 | 0.00270418 | 7316 | 1.13527e+07 | 84.2527 | 17809.3 | 6.85194 | 1.1391e+07 | 84.311 | 0 | nan | nan | nan | nan | 658.977539358 | 5.58907e-06 | 48.1019007076 | 6.57001e-05 | -0.00367803 | 0.033705 |
| 258.44719038 | 0.00270356 | 7317 | 1.10874e+07 | 83.2584 | 17793.6 | 6.84324 | 1.11284e+07 | 83.3245 | 0 | nan | nan | nan | nan | 658.977533777 | 5.6692e-06 | 48.1753740325 | 6.50243e-05 | -0.00364378 | 0.0347512 |
| 258.46762346 | 0.00270294 | 7318 | 1.09189e+07 | 82.6193 | 17790.1 | 6.85121 | 1.09536e+07 | 82.6836 | 0 | nan | nan | nan | nan | 658.977162406 | 5.72111e-06 | 48.2262843684 | 6.45673e-05 | -0.0039766 | 0.0342977 |