Local wavelet-based filtering of electromyographic signals to eliminate the electrocardiographic-induced artifacts in patients with spinal cord injury

Matthew Nitzken, Nihit Bajaj, Sevda Aslan, Georgy Gimel’farb, Ayman El-Baz, Alexander Ovechkin

BioImaging Laboratory, Bioengineering Department, University of Louisville, Louisville, USA.
Department of Computer Science, The University of Auckland, Auckland, New Zealand.
Department of Neurological Surgery, University of Louisville, Louisville, USA Frazier Rehab Institute, Louisville, USA.
Department of Neurological Surgery, University of Louisville, Louisville, USA;Frazier Rehab Institute, Louisville, USA.

DOI: 10.4236/jbise.2013.67A2001                

Abstract

Surface Electromyography (EMG) is a standard method used in clinical practice and research to assess motor function in order to help with the diagnosis of neuromuscular pathology in human and animal models. EMG recorded from trunk muscles involved in the activity of breathing can be used as a direct measure of respiratory motor function in patients with spinal cord injury (SCI) or other disorders associated with motor control deficits. However, EMG potentials recorded from these muscles are often contaminated with heart-induced electrocardiographic (ECG) signals. Elimination of these artifacts plays a critical role in the precise measure of the respiratory muscle electrical activity. This study was undertaken to find an optimal approach to eliminate the ECG artifacts from EMG recordings. Conventional global filtering can be used to decrease the ECG-induced artifact. However, this method can alter the EMG signal and changes physiologically relevant information. We hypothesize that, unlike global filtering, localized removal of ECG artifacts will not change the original EMG signals. We develop an approach to remove the ECG artifacts without altering the amplitude and frequency components of the EMG signal by using an externally recorded ECG signal as a mask to locate areas of the ECG spikes within EMG data. These segments containing ECG spikes were decomposed into 128 sub-wavelets by a custom-scaled Morlet Wavelet Transform. The ECG-related subwavelets at the ECG spike location were removed and a de-noised EMG signal was reconstructed. Validity of the proposed method was proven using mathematical simulated synthetic signals and EMG obtained from SCI patients. We compare the Rootmean Square Error and the Relative Change in Variance between this method, global, notch and adaptive filters. The results show that the localized wavelet-based filtering has the benefit of not introducing error in the native EMG signal and accurately removing ECG artifacts from EMG signals.

Keywords

EMG; De-Noising; ECG; Local Wavelet Filtering

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Devivo, M.J. (2012) Epidemiology of traumatic spinal cord injury: Trends and future implications. Spinal Cord, 50, 365-372.
[2]	Thomas, C.K., Zaidner, E.Y., Calancie, B., Broton, J.G. and Bigland-Ritchie, B.R. (1997) Muscle weakness, paralysis, and atrophy after human cervical spinal cord injury. Experimental Neurology, 148, 414-423. doi:10.1006/exnr.1997.6690
[3]	El Masry, W.S., Tsubo, M., Katoh, S., El Miligui, Y.H.S. and Khan, A. (1996) Validation of the American Spinal Injury Association (ASIA) motor score and the National Acute Spinal Cord Injury Study (NASCIS) motor score. Spine, 21, 614. doi:10.1097/00007632-199603010-00015
[4]	von Tscharner, V. (2000) Intensity analysis in time-frequency space of surface myoelectric signals by wavelets of specified resolution. Journal of Electromyography and Kinesiology, 10, 433-445. doi:10.1016/S1050-6411(00)00030-4
[5]	Clancy, E., Bouchard, S., et al. (2001) Estimation and application of EMG amplitude during dynamic contractions. IEEE Engineering in Medicine and Biology Magazine, 20, 47-54.
[6]	Redfern, M., Hughes, R., et al. (1993) High-pass filtering to remove electrocardiographic interference from torso EMG recordings. Clinical Biomechanics, 8, 44-48.
[7]	Drake, J.D. and Callaghan, J.P. (2006) Elimination of electrocardiogram contamination from electromyogram signals: An evaluation of currently used removal techniques. Journal of Electromyography and Kinesiology, 16, 175-187. doi:10.1016/j.jelekin.2005.07.003
[8]	Rodríguez-Carreno, I., Malanda-Trigueros, A., Gila-Useros, L., Navallas-Irujo, J. and Rodríguez-Falces, J. (2006) Filter design for cancellation of baseline-fluctuation in needle EMG recordings. Computer Methods and Programs in Biomedicine, 81, 79-93. doi:10.1016/j.cmpb.2005.11.002
[9]	Schweitzer, T.W., Fitzgerald, J.W., Bowden, J.A. and Lynne-Davies, P. (1979) Spectral analysis of human inspiratory diaphragmatic electromyograms. Journal of Applied Physiology, 46, 152-165.
[10]	van Boxtel, A. (2001) Optimal signal bandwidth for the recording of surface EMG activity of facial, jaw, oral, and neck muscles. Psychophysiology, 38, 22-34. doi:10.1111/1469-8986.3810022
[11]	Zschorlich, V.R. (1989) Digital filtering of EMG-signals. Electromyography and Clinical Neurophysiology, 29, 8186.
[12]	Hof, A.L. (2009) A simple method to remove ECG artifacts from trunk muscle EMG signals. Journal of Electromyography & Kinesiology, 19, 554-555. doi:10.1016/j.jelekin.2008.11.007
[13]	Lu, G.H., Brittain, J.-S., Holland, P., Yianni J., Green, A.L., Stein, J.F., Aziz, T.Z. and Wang, S.Y. (2009) Removing ECG noise from surface EMG signals using adaptive filtering. Neuroscience Letters, 462, 14-19. doi:10.1016/j.neulet.2009.06.063
[14]	Marque, C., Bisch, C., Dantas, R., Elayoubi, S., Brosse, V. and Pérot, C. (2005) Adaptive filtering for ECG rejection from surface EMG recordings. Journal of Electromyography and Kinesiology, 15, 310-315. doi:10.1016/j.jelekin.2004.10.001
[15]	Bloch, R. (1983) Subtraction of electrocardiographic signal from respiratory electromyogram. Journal of Applied Physiology, 55, 619-623.
[16]	Deng, Y.C., Wolf, W., Schnell, R. and Appel, U. (2000) New aspects to event-synchronous cancellation of ECG interference: An application of the method in diaphragmatic EMG signals. IEEE Transactions on Biomedical Engineering, 47, 1177-1184. doi:10.1109/10.867924
[17]	Hyvarinen, A., Karhunen, J. and Oja, E. (2001) Independent component analysis. Wiley-Interscience, Hoboken. doi:10.1002/0471221317
[18]	Costa Junior, J.D., Ferreira, D.D., Nadal, J. and Miranda de Sa, A.M.F.L. (2010) Reducing electrocardiographic artifacts from electromyogram signals with independent component analysis. 2010 Annual International Conference of the IEEE on Engineering in Medicine and Biology Society (EMBC), Buenos Aires, 31 August-4 September 2010, 4598-4601.
[19]	Carre, P., Leman, H., Fernandez, C. and Marque, C. (1998) Denoising of the uterine EHG by an undecimated wavelet transform. IEEE Transactions on Biomedical Engineering, 45, 1104-1113. doi:10.1109/10.709554
[20]	Daubechies, I. (1992) Ten lectures on wavelets. Society for Industrial and Applied Mathematics (SIAM), Philadelphia.
[21]	Mallat, S.G. (1989) A theory for multiresolution signal decomposition: The wavelet representation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 11, 674-693. doi:10.1109/34.192463
[22]	Zhou, W.D. and Gotman, J. (2004) Removal of EMG and ECG artifacts from EEG based on wavelet transform and ICA. 26th Annual International Conference of the IEEE on Engineering in Medicine and Biology Society, IEMBS’04, San Francisco, 1-5 September 2004, 392-395.
[23]	Taelman, J., Van Huffel, S. and Spaepen, A. (2007) Wavelet-independent component analysis to remove electrocardiography contamination in surface electromyography. 29th Annual International Conference of the IEEE on Engineering in Medicine and Biology Society, EMBS 2007, Lyon, 22-26 August 2007, 682-685.
[24]	von Tscharner, V., Eskofier, B. and Federolf, P. (2011) Removal of the electrocardiogram signal from surface EMG recordings using non-linearly scaled wavelets. Journal of Electromyography and Kinesiology, 21, 683688. doi:10.1016/j.jelekin.2011.03.004
[25]	Swain, D.P., Abernathy, K.S., Smith, C.S., Lee, S.J. and Bunn, S.A. (1994) Target heart rates for the development of cardiorespiratory fitness. Medicine & Science in Sports & Exercise, 26, 112-116.
[26]	Grossmann, A. and Morlet, J. (1984) Decomposition of Hardy functions into square integrable wavelets of constant shape. SIAM Journal on Mathematical Analysis, 15, 723-736. doi:10.1137/0515056
[27]	Torrence, C. and Compo, G.P. (1998) A practical guide to wavelet analysis. Bulletin of American Meteorological Society, 79, 61-78. doi:10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2
[28]	McGill, K.C., Lateva, Z.C. and Marateb, H.R. (2005) EMGLAB: An interactive EMG decomposition program. Journal of Neuroscience Methods, 149, 121-133. doi:10.1016/j.jneumeth.2005.05.015
[29]	Karthik, R. (2006) ECG simulation using MATLAB. College of Engineering, Guindy, Anna University, Chennai
[30]	Nehorai, A. and Porat, B. (1986) Adaptive comb filtering for harmonic signal enhancement. IEEE Transactions on Acoustics, Speech and Signal Processing, 34, 1124-1138.