search

joshbriegal / roto

This package provides tools for extracting rotation periods from astrophysical time series.
MIT License
3 stars 1 forks source link

Add minimum period in LombScargle #20

Closed joshbriegal closed 3 years ago

joshbriegal commented 3 years ago

e.g. from other code base. To prevent spuriously short periods based on sampling.

def calculate_periodogram(self) -> Periodogram:
        min_timestamp_diff = np.min(np.diff(self.data.timeseries))
        maximum_frequency = 1.0 / (self.params.nyquist_factor * min_timestamp_diff)
        return Periodogram(
            *self.lomb_scargle.autopower(maximum_frequency=maximum_frequency)
        )