netpyne.support.morlet

Module with support for Morlet wavelets

Time-frequency representation using Morlet wavelets

Original version written by Shane Lee (Brown University) Modified by Sam Neymotin (NKI; added phase calculations, saving/passing in Morlet, specifying different frequency steps and frequencies used for calculations (e.g. logarithmic frequency instead of linear) don’t detrend the signal > 1 time; don’t waste RAM by copying/storing time-series in MorletSpec class; subtract mean from time series within the wavelet class)

Functions:

index2ms(idx, sampr)
ms2index(ms, sampr)
Morlet(sampr, freq[, width])	Morlet's wavelet for frequency f and time t Wavelet normalized so total energy is 1 freq: specific frequency (Hz) width: number of cycles for the Morlet wavelet output returned: final units are 1/s
MorletVec(sig, sampr, freq, width[, m, getphase])	Final units of y: signal units squared.

Classes:

MorletSpec(tsvec, sampr[, freqmin, freqmax, ...])

netpyne.support.morlet.index2ms(idx, sampr)

netpyne.support.morlet.ms2index(ms, sampr)

netpyne.support.morlet.Morlet(sampr, freq, width=7.0)

Morlet’s wavelet for frequency f and time t Wavelet normalized so total energy is 1 freq: specific frequency (Hz) width: number of cycles for the Morlet wavelet output returned: final units are 1/s

netpyne.support.morlet.MorletVec(sig, sampr, freq, width, m=None, getphase=False)

Final units of y: signal units squared. For instance, a signal of Am would have Am^2 The energy is calculated using Morlet’s wavelets; also returns phase when getphase==True sig: input signal sampr: sampling rate (Hz) freq: frequency m: Morlet wavelet width: number of cycles of Morlet wavelet

class netpyne.support.morlet.MorletSpec(tsvec, sampr, freqmin=1.0, freqmax=250.0, freqstep=1.0, width=7.0, getphase=False, lfreq=None)

Bases: object

Methods:

plot_to_ax(ax_spec, dt)
transform(tsvec[, getphase])

plot_to_ax(ax_spec, dt)

transform(tsvec, getphase=False)