How to remove common systematics using basis vectors (CBVs)¶
Cotrending Basis Vectors are used to remove long term trends in light curves. They are built from the most common trends observed in each channel. You can read more about CBVs in Demystifying Kepler Data. They can be used to clean lightcurves of common trends experience by all targets.
Below is an example workflow to remove CBVs from an light curve in Kepler; Tabby’s star.
from lightkurve import search_lightcurvefile lcf = search_lightcurvefile(8462852, quarter=8).download() lc = lcf.SAP_FLUX #Remove nans and outliers lc = lc.remove_nans() %matplotlib inline lc.plot(linestyle='-');
Downloading URL https://mast.stsci.edu/api/v0/download/file?uri=mast:Kepler/url/missions/kepler/lightcurves/0084/008462852/kplr008462852-2011073133259_llc.fits to /home/gb/.lightkurve-cache/mastDownload/Kepler/kplr008462852_lc_Q111111111111111111/kplr008462852-2011073133259_llc.fits ... [Done]
The light curve above has a significant dip due to the astrophysical
variability from Tabby’s star. However there are other long term trends
and small ramps due to heating in the spacecraft, which we would like to
remove. We can use CBV’s to do this. Let’s import the CBV corrector from
from lightkurve.correctors import KeplerCBVCorrector
We can use the
plot_cbvs function to take a look at the trends we
are going to remove. Let’s look at just the first 8 trends.
We probably don’t want to correct all these trends out of Tabby’s star. Let’s pick just the first 5.
corr_lc = KeplerCBVCorrector(lcf).correct(cbvs=[1,2,3,4,5])
We can plot the two light curves on the same plot window to see the difference.
ax = lc.plot(color='C3', label='SAP Flux', linestyle='-') corr_lc.plot(ax=ax, color='C2', label='CBV Corrected SAP Flux', linestyle='-');
Now the CBVs have been corrected only the astrophysical variability is left, the variability from spacecraft heating has been removed.