code: def rayleigh_spectrogram(self, stride, fftlength=None, overlap=0,
                               nproc=1, **kwargs):
          
          specgram = spectral.average_spectrogram(
              self,
              spectral.rayleigh,
              stride,
              fftlength=fftlength,
              overlap=overlap,
              nproc=nproc,
              **kwargs
          )
          specgram.override_unit('')
          return specgram
docs: Calculate the Rayleigh statistic spectrogram of this `TimeSeries`
  
          Parameters
          ----------
          stride : `float`
              number of seconds in single PSD (column of spectrogram).
  
          fftlength : `float`
              number of seconds in single FFT.
  
          overlap : `float`, optional
              number of seconds of overlap between FFTs, default: ``0``
  
          nproc : `int`, optional
              maximum number of independent frame reading processes, default
              default: ``1``
  
          Returns
          -------
          spectrogram : `~gwpy.spectrogram.Spectrogram`
              time-frequency Rayleigh spectrogram as generated from the
              input time-series.
  
          See Also
          --------
          TimeSeries.rayleigh
              for details of the statistic calculation