Authors

1 MS.c Student, Electrical Engineering Department, Iran University of Science and Technology, Tehran, Iran

2 Assistant Professor, Electrical Engineering Department, Iran University of Science and Technology, Tehran, Iran

Abstract

The extracellular recording from the brain's single neurons is known as a popular method in neuroscience and neuro-rehabilitation engineering. These recordings include the activity of all neurons around the electrode, for better use of which, spike sorting methods should be utilized to obtain the activity of single neurons. Based on the structural properties of the neuron, such as its dendritic tree, and the distance and direction of it relative to the electrode, it can be claimed that the form of its spike waveform is unique and constant. However, spike sorting under low signal-to-noise ratio (SNR) conditions is always accompanied with challenges. A spike sorting algorithm usually consists of three sections including the spike detection, feature extraction, and classification. In this paper, a method based on optimization of continuous wavelet coefficients is presented which is effective in low SNR values. In the proposed method, after the calculation of the parameterized wavelet coefficients, using the Euclidean distance and the area under the receiver operator characteristic curve, the best parameters were chosen to increase the separation of the features, so that a suitable scale was first found with the Euclidean distance criterion and then the translation parameter was obtained with the second criterion. In this research k-means algorithm was used for the clustering as a simple but efficient method. For evaluation, three simulated data sets were made in 9 different SNRs with a modeled background noise. The obtained results from simulated data showed that the optimization of parameters in continuous wavelet transform using the proposed algorithm could effectively improve the spike sorting performance compared to principal component analysis method.

Keywords

Main Subjects

[1]    Thomas Jessell, Siegelbaum, S., & Hudspeth, A. J., Principles of neural science, New york: McGraw-hill, 2000.
[2]    W. Huang, "A Framework for Real Time and Automatic Spike Sorting of Multichannel Neuronal Activity," International Journal of Enhanced Research in Science, Technology & Engineering, 2017.
[3]    Kandel, E.R., Schwartz, J.H. and Jessell, T.M. eds., Principles of neural science, New York: McGraw-hill, 2000.
[4]    M. S. Lewicki, "A review of methods for spike sorting: the detection and classification of neural action potentials. Network," Computation in Neural Systems, vol. 9, no. 4, pp. 53-78, 1998.
[5]    Gibson, S., Judy, J. W., & Marković, D. , "Spike sorting: The first step in decoding the brain: The first step in decoding the brain," IEEE Signal processing magazine, vol. 29, no. 1, pp. 124-143, 2012.
[6]    Mukhopadhyay, S., & Ray, G. C., "A new interpretation of nonlinear energy operator and its efficacy in spike detection," IEEE Transactions on biomedical engineering, vol. 45, no. 2, pp. 180-187, 1998.
[7]    Maragos, P., Kaiser, J. F., & Quatieri, T. F. , "On amplitude and frequency demodulation using energy operators," IEEE Transactions on signal processing, vol. 41, no. 4, p. 1532–1550, 1993.
[8]    Chandra, R., & Optican, L. M., "Detection, classification, and superposition resolution of action potentials in multiunit single-channel recordings by an on-line real-time neural network," IEEE Transactions on Biomedical Engineering, vol. 44, no. 5, p. 403–412, 1997.
[9]    T. Sapatinas, "Discriminant Analysis and Statistical Pattern Recognition," Journal of the Royal Statistical Society, vol. 168, no. Series A (Statistics in Society), pp. 635-636, 2005.
[10]W. Cai, J. Xu, and L. Huang,, "Low Power SI Class E Power Amplifier And RF Switch for HealthCare," Informatics Engineering, an International Journal, vol. 4, no. 4, pp. 7-13, 2016.
[11]Dempster, A. P., Laird, N. M., & Rubin, D. B., "Maximum likelihood from incomplete data via the EM algorithm," Journal of the royal statistical society, no. Series B (Methodological), pp. 1-38, 1997.
[12]Hazan, L., Zugaro, M., & Buzsáki, G., "Klusters, NeuroScope, NDManager: a free software suite for neurophysiological data processing and visualization," Journal of neuroscience methods, vol. 155, no. 2, pp. 207-216, 2006.
[13]Zhang, C., Zhang, X., Zhang, M. Q., & Li, Y., "Neighbor number, valley seeking and clustering," Pattern Recognition Letters, vol. 28, no. 2, p. 173–180, 2007.
[14]Blatt, M., Wiseman, S., & Domany, E. , "Data clustering using a model granular magnet," Neural Comput., vol. 9, no. 8, p. 1805–1842, 1997.
[15]Quiroga, R. Q., Nadasdy, Z., & Ben-Shaul, Y., "Unsupervised spike detection and sorting with wavelets and superparamagnetic clustering," Neural computation, vol. 16, no. 8, pp. 1661-1687, 2004.
[16]Einevoll, G. T., Franke, F., Hagen, E., Pouzat, C., & Harris, K. D., "Towards reliable spike-train recordings from thousands of neurons with multielectrodes," Current opinion in neurobiology, vol. 22, no. 1, pp. 11-17, 2012.
[17]Harris, K. D., Hirase, H., Leinekugel, X., Henze, D. A., & Buzsáki, G. , "Temporal interaction between single spikes and complex spike bursts in hippocampal pyramidal cells," Neuron, vol. 32, no. 1, pp. 141-149, 2001.
[18]Chung, J. E., Magland, J. F., Barnett, A. H., Tolosa, V. M., Tooker, A. C., Lee, K. Y., ... & Greengard, L. F., "A fully automated approach to spike sorting," Neuron , vol. 95, no. 6, pp. 1381-1394, 2017.
[19]Abeles, M., & Goldstein, M. H., "Multispike train analysis," Proceedings of the IEEE, vol. 65, no. 5, pp. 762-773, 1997.
[20]Glasen, E. M., & Marks, W. B., "On-line Separation of Interleaved Neuronal Pulse Sequences, Data Acquisition Process," Bio. Med, vol. 5, pp. 137-156, 1968.
[21]Shalchyan, V., Jensen, W., & Farina, D., "Spike detection and clustering with unsupervised wavelet optimization in extracellular neural recordings," IEEE Transactions on Biomedical Engineering, vol. 59, no. 9, pp. 2576-2585, 2012.
[22]Yang, Y., & Mason, A. J., "Frequency band separability feature extraction method with weighted haar wavelet implementation for implantable spike sorting," IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 25, no. 6, pp. 530-538, 2017.
[23]Letelier, J. C., & Weber, P. P., "Spike sorting based on discrete wavelet transform coefficients," Journal of neuroscience methods, vol. 101, no. 2, p. 93–106, 2000.
[24]Vargas-Irwin, C., & Donoghue, J. P., "Automated spike sorting using density grid contour clustering and subtractive waveform decomposition," Journal of neuroscience methods, vol. 164, no. 1, pp. 1-18, 2007.
[25]"Center for BioMedical Signal Processing and Computation (CBMSPC)," University of California, Irvine, [Online]. Available: http://cbmspc.eng.uci.edu/SOFTWARE/SPIKEDETECTION/templates.mat.
[26]Khorasani, A., Beni, N. H., Shalchyan, V., & Daliri, M. R. , "Continuous Force Decoding from Local Field Potentials of the Primary Motor Cortex in Freely Moving Rats," Scientific Reports, vol. 6, p. 35238, 2016.
[27]           Box, G. E. P., G. M. Jenkins, and G. C. Reinsel., Time Series Analysis: Forecasting and Control, NJ: Prentice Hall: Englewood Cliffs, 1994.
[28]K. P. Soman, Insight into wavelets: From theory to practice, PHI Learning Pvt. Ltd, 2010.
[29]Pavlov, A., Makarov, V. A., Makarova, I., & Panetsos, F. , "Sorting of neural spikes: when wavelet based methods outperform principal component analysis," Natural Computing, vol. 6, no. 3, pp. 269-281, 2007.
[30]Hulata, E., Segev, R., & Ben-Jacob, E., "A method for spike sorting and detection based on wavelet packets and Shannon's mutual information," Journal of Neuroscience Methods, vol. 117, no. 1, pp. 1-12, 2002.