[SciPy-Dev] Lomb-Scargle periodogram
Wed Dec 22 12:38:13 CST 2010
On Wed, Dec 22, 2010 at 1:20 PM, Warren Weckesser
> On Mon, Dec 20, 2010 at 3:27 PM, Pim Schellart <P.Schellart@astro.ru.nl>
>> Dear SciPy developers,
>> in my research field of Astronomy we often encounter timeseries with
>> uneven temporal sampling.
>> Such timeseries cannot be analyzed with the standard Fast Fourier
>> The most frequently used tool to analyze such timeseries is the
>> Lomb-Scargle periodogram, developed by Lomb (1976) and further extended by
>> Scargle (1982).
>> Currently, to my knowledge, SciPy does not include a routine to calculate
>> such periodograms.
>> Therefore I would like to contribute the following code which efficiently
>> calculates the periodogram.
>> It is licensed under the terms of the BSD open source license but I would
>> be willing to release it under any other language if this is necessary for
>> possible inclusion with SciPy.
>> The routine is written in Fortran 95 and wrapped with f2py, but I would be
>> willing to rewrite it in C if required.
>> Also included is a Python test script.
>> Please provide me with any feedback in case you would consider including
> I think an implementation of Lomb-Scargle and related algorithms would be a
> great addition to SciPy. With the 0.9 release coming soon, and also being
> holiday season, I'm not sure how soon anyone will be able to look at your
> code. Could you add this as an enhancement request to the Trac system?
> Here's the link:
> Create an account, and then create a new ticket. You can attach the code to
> the ticket.
> Best regards,
> P.S. There is a python/numpy implementation here:
I would also be interested, but I am, like scipy, still on fortran 77,
g77 on older MingW
> SciPy-Dev mailing list
More information about the SciPy-Dev