A New Spike Detection Algorithm for Extracellular Neural Recordings

Signals from extracellular electrodes in neural systems record voltages resulting from activity in many neurons. Detecting action potentials (spikes) in a small number of specific (target) neurons is difficult because many neurons, both near and more distant, contribute to the signal at the electrode. We consider some nearby neurons as target neurons (providing a signal) and all the other contributions to the signal as noise. A new algorithm for spike detection has been developed: this applies a cepstrum of bispectrum (CoB) estimated inverse filter to provide blind equalization. This technique is based on higher order statistics, and seeks to find a sequence of event times or delta sequence. We show that the CoB-based technique can achieve a 98% hit rate on an extracellular signal containing three spike trains at up to 0 dB SNR. Threshold setting for this technique is discussed, and we show the application of the technique to some real signals. We compare performance with four established techniques and report that the CoB-based algorithm performs best.

[1]  J. Letelier,et al.  Spike sorting based on discrete wavelet transform coefficients , 2000, Journal of Neuroscience Methods.

[2]  I. Miller Probability, Random Variables, and Stochastic Processes , 1966 .

[3]  W. Marsden I and J , 2012 .

[4]  R. Quian Quiroga,et al.  Unsupervised Spike Detection and Sorting with Wavelets and Superparamagnetic Clustering , 2004, Neural Computation.

[5]  A. S. Spinelli,et al.  A simple method for efficient spike detection in multiunit recordings , 2007, Journal of Neuroscience Methods.

[6]  Jing Wang,et al.  A spike detection method in EEG based on improved morphological filter , 2007, Comput. Biol. Medicine.

[7]  Leslie S. Smith,et al.  A tool for synthesizing spike trains with realistic interference , 2007, Journal of Neuroscience Methods.

[8]  Bruno O. Shubert,et al.  Random variables and stochastic processes , 1979 .

[9]  Aarnout Brombacher,et al.  Probability... , 2009, Qual. Reliab. Eng. Int..

[10]  Sung June Kim,et al.  Neural spike sorting under nearly 0-dB signal-to-noise ratio using nonlinear energy operator and artificial neural-network classifier , 2000, IEEE Transactions on Biomedical Engineering.

[11]  Christof Koch,et al.  Electrical Interactions via the Extracellular Potential Near Cell Bodies , 1999, Journal of Computational Neuroscience.

[12]  M.R. Raghuveer,et al.  Bispectrum estimation: A digital signal processing framework , 1987, Proceedings of the IEEE.

[13]  Leslie S. Smith,et al.  Smoothing and thresholding in neuronal spike detection , 2006, Neurocomputing.

[14]  J. Okada,et al.  Multineuronal spike classification based on multisite electrode recording, whole-waveform analysis, and hierarchical clustering , 1999, IEEE Transactions on Biomedical Engineering.

[15]  A.F. Atiya,et al.  Recognition of multiunit neural signals , 1992, IEEE Transactions on Biomedical Engineering.

[16]  J. F. Kaiser,et al.  On a simple algorithm to calculate the 'energy' of a signal , 1990, International Conference on Acoustics, Speech, and Signal Processing.

[17]  C. L. Nikias,et al.  Signal processing with higher-order spectra , 1993, IEEE Signal Processing Magazine.

[18]  Manuel Rodrı́guez,et al.  Unsupervised classification of neural spikes with a hybrid multilayer artificial neural network , 1998, Journal of Neuroscience Methods.

[19]  Joel W. Burdick,et al.  Spike detection using the continuous wavelet transform , 2005, IEEE Transactions on Biomedical Engineering.

[20]  Sung June Kim,et al.  A wavelet-based method for action potential detection from extracellular neural signal recording with low signal-to-noise ratio , 2003, IEEE Transactions on Biomedical Engineering.

[21]  Paul Wintz,et al.  Digital image processing (2nd ed.) , 1987 .

[22]  X. Yang,et al.  A totally automated system for the detection and classification of neural spikes , 1988, IEEE Transactions on Biomedical Engineering.

[23]  Athanasios Papoulis,et al.  Probability, Random Variables and Stochastic Processes , 1965 .

[24]  Leslie S. Smith,et al.  A novel technique for spike detection in extracellular neurophysiological recordings using Cepstrum of Bispectrum , 2008, 2008 16th European Signal Processing Conference.

[25]  Jacqueline Walker,et al.  Application of higher order statistics techniques to EMG signals to characterize the motor unit action potential , 2005, IEEE Transactions on Biomedical Engineering.

[26]  Petros Maragos,et al.  On amplitude and frequency demodulation using energy operators , 1993, IEEE Trans. Signal Process..

[27]  Clark F. Olson Maximum-Likelihood Image Matching , 2002, IEEE Trans. Pattern Anal. Mach. Intell..

[28]  Reid R. Harrison,et al.  A low-power integrated circuit for adaptive detection of action potentials in noisy signals , 2003, Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE Cat. No.03CH37439).

[29]  Jacqueline Walker,et al.  Cepstrum of bispectrum - A new approach to blind system reconstruction , 2008, Signal Process..

[30]  R. Chandra,et al.  Detection, classification, and superposition resolution of action potentials in multiunit single-channel recordings by an on-line real-time neural network , 1997, IEEE Transactions on Biomedical Engineering.

[31]  Hagai Bergman,et al.  A personal computer-based spike detector and sorter: implementation and evaluation , 1992, Journal of Neuroscience Methods.

[32]  Edward M. Schmidt,et al.  Computer separation of multi-unit neuroelectric data: a review , 1984, Journal of Neuroscience Methods.

[33]  K. Horch,et al.  Separation of action potentials in multiunit intrafascicular recordings , 1992, IEEE Transactions on Biomedical Engineering.

[34]  WangJing,et al.  A spike detection method in EEG based on improved morphological filter , 2007 .

[35]  Steven M. Bierer,et al.  Multi-channel spike detection and sorting using an array processing technique , 1999, Neurocomputing.

[36]  Patrick D. Wolf,et al.  Evaluation of spike-detection algorithms fora brain-machine interface application , 2004, IEEE Transactions on Biomedical Engineering.

[37]  S. Mukhopadhyay,et al.  A new interpretation of nonlinear energy operator and its efficacy in spike detection , 1998, IEEE Transactions on Biomedical Engineering.

[38]  Clark F. Olson,et al.  Maximum-likelihood template matching , 2000, Proceedings IEEE Conference on Computer Vision and Pattern Recognition. CVPR 2000 (Cat. No.PR00662).

[39]  M. Abeles,et al.  Multispike train analysis , 1977, Proceedings of the IEEE.

[40]  J. Csicsvari,et al.  Intracellular features predicted by extracellular recordings in the hippocampus in vivo. , 2000, Journal of neurophysiology.

[41]  Tom Fawcett,et al.  An introduction to ROC analysis , 2006, Pattern Recognit. Lett..

[42]  D. C. Tam,et al.  Multi-unit spike discrimination using wavelet transforms , 1997, Comput. Biol. Medicine.

[43]  Paul Wintz,et al.  Instructor's manual for digital image processing , 1987 .

[44]  Leslie S. Smith,et al.  Finding events in noisy signals , 2007 .