NeuroDSP: A package for neural digital signal processing

Abstract

Neural oscillations are widely studied using methods based on the Fourier transform, which models data as sums of sinusoids. For decades these Fourier-based approaches have successfully uncovered links between oscillations and cognition or disease. However, because of the fundamental sinusoidal basis, these methods might not fully capture neural oscillatory dynamics, because neural data are both nonsinusoidal and non-stationary. Here, we present a new analysis framework, complementary to Fourier analysis, that quantifies cycle-by-cycle time-domain features. For each cycle, the amplitude, period, and waveform symmetry are measured, the latter of which is missed using conventional approaches. Additionally, oscillatory bursts are algorithmically identified, allowing us to investigate the variability of oscillatory features within and between bursts. This approach is validated on simulated noisy signals with oscillatory bursts and outperforms conventional metrics. Further, these methods are applied to real data—including hippocampal theta, motor cortical beta, and visual cortical alpha—and can differentiate behavioral conditions.