Pulsation artifact removal from intra-operatively recorded local field potentials using sparse signal processing and data-specific dictionary

Neural recordings frequently get contaminated by ECG or pulsation artifacts. These large amplitude components can mask the neural patterns of interest and make the visual inspection process difficult. The current study describes a sparse signal representation strategy that targets to denoise pulsation artifacts in local field potentials (LFPs) recorded intraoperatively. To estimate the morphology of the artifact, we first detect the QRS-peaks from the simultaneously recorded ECG trace as an anchor point. After the LFP data has been epoched with respect to each beat, a pool of raw data segments of a specific length is generated. Using the K-singular value decomposition (K-SVD) algorithm, we constructed a data-specific dictionary to represent each contaminated LFP epoch in a sparse fashion. Since LFP is aligned to each QRS complex and the background neural activity is uncorrelated to the anchor points, we assumed that constructed dictionary will be formed to mainly represent the pulsation artifact. In this scheme, we performed an orthogonal matching pursuit to represent each LFP epoch as a linear combination of the dictionary atoms. The denoised LFP data is thus obtained by calculating the residual between the raw LFP and its approximation. We discuss and demonstrate the improvements in denoised data and compare the results with respect to principal component analysis (PCA). We noted that there is a comparable change in the signal for visual inspection to observe various oscillating patterns in the alpha and beta bands. We also see a noticeable compression of signal strength in the lower frequency band (<13 Hz), which was masked by the pulsation artifact, and a strong increase in the signal-to-noise ratio (SNR) in the denoised data.Clinical Relevance- Pulsation artifact can mask relevant neural activity patterns and make their visual inspection difficult. Using sparse signal representation, we established a new approach to reconstruct the quasiperiodic pulsation template and computed the residue signal to achieve noise-free neural activity.