On Mon, May 16, 2022 at 11:04 AM Robert Kern <robert.kern@gmail.com> wrote:
I don't know of any work that apply's something like Welch's method per se to the LS periodogram, but if you squint, you can kind of see Welch's method as an ad hoc instance of a multitaper method, and there is a more rigorous application of multitaper methods to the LS periodogram. Not implemented in Python so far as I know, but there is an R implementation.

  https://iopscience.iop.org/article/10.3847/1538-3881/ab7fa1 

Looking through this paper more carefully, I see it cites this paper that implements multitaper spectral estimate on regularly-sampled but gappy time series like yours.

  https://core.ac.uk/download/pdf/237709478.pdf

The corresponding MATLAB code is BSD-licensed:

  https://www.mathworks.com/matlabcentral/fileexchange/71909-mdmwps

--
Robert Kern