alternate FFT calculation using lomb.py

nksaunders / giants

Finding Giant Planets in TESS FFIs

2 stars 3 forks source link

alternate FFT calculation using lomb.py #10

Closed skgrunblatt closed 4 years ago

skgrunblatt commented 4 years ago

osample=5
nyq=283

freq, amp, nout, jmax, prob = lomb.fasper(lc.time,lc.flux, osample, 3.)
freq = 1000.*freq/86.4
binn = freq[1]-freq[0]
fts = 2.*amp*np.var(lc.flux*1e6)/(np.sum(amp)*binn)
fts = scipy.ndimage.filters.gaussian_filter(fts, 4)
use=np.where(freq < nyq+150)
freq=freq[use]
fts=fts[use]

skgrunblatt commented 4 years ago

lomb.py.zip