syncopy.specest.superlet.superlet#

syncopy.specest.superlet.superlet(data_arr, samplerate, scales, order_max, order_min=1, c_1=3, adaptive=False)[source]#

Performs Superlet Transform (SLT) according to Moca et al. [1] Both multiplicative SLT and fractional adaptive SLT are available. The former is recommended for a narrow frequency band of interest, whereas the latter is better suited for the analysis of a broad range of frequencies.

A superlet (SL) is a set of Morlet wavelets with increasing number of cycles within the Gaussian envelope. Hence the bandwith is constrained more and more with more cycles yielding a sharper frequency resolution. Complementary the low cycle numbers will give a high time resolution. The SLT then is the geometric mean of the set of individual wavelet transforms, combining both wide and narrow-bandwidth wavelets into a super-resolution estimate.

Parameters:

  • data_arr (nD numpy.ndarray) – Uniformly sampled time-series data The 1st dimension is interpreted as the time axis

  • samplerate (float) – Samplerate of the time-series in Hz

  • scales (1D numpy.ndarray) – Set of scales to use in wavelet transform. Note that for the SL Morlet the relationship between scale and frequency simply is s(f) = 1/(2*pi*f) Need to be ordered high to low for adaptive=True

  • order_max (int) – Maximal order of the superlet set. Controls the maximum number of cycles within a SL together with the c_1 parameter: c_max = c_1 * order_max

  • order_min (int) – Minimal order of the superlet set. Controls the minimal number of cycles within a SL together with the c_1 parameter: c_min = c_1 * order_min Note that for admissability reasons c_min should be at least 3!

  • c_1 (int) – Number of cycles of the base Morlet wavelet. If set to lower than 3 increase order_min as to never have less than 3 cycles in a wavelet!

  • adaptive (bool) – Whether to perform fractional adaptive SLT or multiplicative SLT. If set to True, the order of the wavelet set will increase linearly with the frequencies of interest from order_min to order_max. If set to False the same SL will be used for all frequencies.

Returns:

gmean_spec – Complex time-frequency representation of the input data. Shape is (len(scales),) + data_arr.shape

Return type:

numpy.ndarray

Notes