# An effective way of J wave separation based on multilayer NMF.

1. IntroductionIn recent years, as a new index of the ECG of ventricular bipolar, J wave draws more and more attention in the clinical diagnosis. We can consider that J wave is located between QRS wave of electrocardiogram and ST segment [1]. It represents the depolarization end and bipolarization start, under normal circumstances of the heart, ventricular finally depolarization and ventricular earliest bipolarization have the mutual transitional zone which is overlapped about 10 ms in humans, if J point deviates from the baseline, that is, the J point shift, in early repolarization syndrome, acute myocardial ischemia, pericarditis and bundle branch block, and so forth [2]. If J point is shifting and forming the special dome type or rush type, it means that J wave happened.

As the ECG signals we observed are compounded by the baseline drift, power frequency interference and myoelectricity interference noise, and so on, in most cases the original J wave signal is mixed in the normal ST segment [3]. So we cannot observe the characteristics of J wave from mixed signals. In clinic, the ECG and vectorcardiogram provide limited diagnostic information. These results are caused by some complex symptoms or clinical characteristics, and it is hard to decide the relative treatment. From the point view of signal processing, it is difficult to get accurate J wave from ECG with noise interference. It is obvious that the noise must be removed so as to get the desired information from the ECG. Then the mixed signals including J wave and normal ECG signals will be separated effectively.

It is necessary for us to study digital characteristics of J wave and we should extract J wave from ECG signal exclusively. Although the method of processing ECG signals has been very mature, methods of extracting J wave signals from ECG signal are rare to find and the extraction precision is undesirable. It will affect not only the J wave numerical characteristics directly, but also the analysis of the late classification accuracy [4]. In such situation of the clinical diagnosis, it is impossible to diagnose the healthy condition for patients with J wave. This paper puts forward a set of related methods to achieve the desired result. Commonly, the steps of signal processing can be realized as follows:

(1) preprocessing the mixed signals with the method of the principal component analysis (PCA) in order to filter out noises and initialize,

(2) separating the original J wave from mixed signals based on multilayer fixed-point NMF decomposition,

(3) analyzing the separated J wave from point view of the spectrum, power spectrum density, and spectra in digital domain.

The simulation results show that the methods in this paper have better precision and realize the characteristics extraction of J wave. It can be used as the judgment basis for doctors.

2. Nonnegative Matrix Decomposition

Although the independent component analysis (ICA) and sparse component analysis (SCA) algorithm based on blind source separation (BSS) are the main methods, they are no longer applicable in practical on the condition that source signals are dependent or nonsparse. In particular, for complex ECG signal, the separation is poor. It is obvious that the nonnegative matrix factorization is a new way of blind source separation to solve the dependent source signals [5]. Therefore, this project adopts the improved NMF algorithm to extract J wave. Because the NMF has no limits of statistical independence source or sparse source, it can realize the source separation better and the algorithm has such features as fast convergence rate, sparse, nonnegative, and dimensionality reduction [6]. The optimized NMF algorithm can meet the practical requirements. The extraction results are desirable, such as extracting J wave from ECG signal accurately with small error and simple procedures. We can get more ideal extraction results.

Mathematical model of BSS problem can be expressed as follows:

X = AS, (1)

where X = [[x.sub.1], [x.sub.2],...,[x.sub.m]] [member of] [R.sup.mxN] is the observed mixed signal matrix; S = [[s.sub.1],[s.sub.2],...,[s.sub.n]] [member of] [R.sup.nxN] is unknown source signal matrix; and A [member of] [R.sup.mxn] is mixed matrix [7]. Suppose N [greater than or equal to] m [greater than or equal to] n, and n is known or can be estimated. The goal of BSS is to isolate unknown A and S from the mixed signal matrix X.

The basic model of NMF is to decompose the matrix V = [[v.sub.1],[v.sub.2],...,[v.sub.m]] [member of] [R.sup.mxN] into the matrix of W = [[w.sub.i], [w.sub.2],...,[w.sub.m]] [member of] [R.sup.mxn] and H = [[h.sub.1],[h.sub.2],...,[h.sub.n]] [member of] [R.sup.mxN], where [W.sub.ij] and [H.sub.ij] are nonnegative under the limiting condition and the formula is defined as follows:

V = WH. (2)

Comparing (1) and (2), we can take the observation data X in BSS as the matrix V to be decomposed in NMF. If W and H are equal to A and S, respectively, then the NMF can be used to solve the problem of BSS and even the J wave separation.

3. The Preprocessing before Extracting

Taking the mean square error (MSE) into consideration, PCA method with the property of optimal orthogonal decomposition is the most representative in the field of signal processing and could compress ECG signal to reduce the dimension and remove the noises to get the maximum signalto-noise ratio. It saves the subsequent decomposition time and improves the separation precision. As the final result of NMF is always the local optimum, the initialization of the objective functions W and H will affect the entire process of the NMF algorithm directly. That is to say, initializing W and H differently will produce different separation results [8]. Comparing to the traditional methods, eventually PCA is chosen for initializing. PCA can find the projection direction with the most concentrated energy in the sample space. The projection direction of the observation signal V will produce reconstruction signals and reconstruction error defined as [E.sub.mse] = E[[parallel][bar.V] - WH[parallel].sup.2], where [bar.V] is a collection standing for the mean of ECG signal; W is the transformation matrix; and H is the projection coefficient matrix. We could calculate H after setting a minimum range for [E.sub.mse] and taking the appropriate W. Then we can obtain the initial transformed matrixes [W.sub.0] and [H.sub.0] of W and H by nonnegative processing. The procedures are [w.sub.ij] = max(0, [w.sub.ij]) and [h.sub.ij] = max(0, [h.sub.ij]). Therefore, PCA has great advantages of preprocessing in the NMF algorithm.

4. Fixed-Point NMF

Compared to the different algorithm of NMF, fixed-point NMF is suitable to extract J wave from the ECG signal. The advantage of this algorithm is that the divergence deviation is zero when V = WH [9]. We can judge the advantages and disadvantages of the algorithm by calculating the final value of corresponding objective function and the algorithm can simplify the evaluation of separation accuracy. The objective function is shown in the following formula:

D(V[parallel]WH) = [summation over (i,j)][[V.sub.ij] [log.sub.2] [V.sub.ij]/[(WH).sub.ij] - [V.sub.ij] + [(WH).sub.ij]], (3)

where W and H are nonnegative and when and only when V = WH is satisfied, formula (3) can obtain the minimum value 0.

The iteration rules of W and H are shown in the following formula:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)

where [W.sub.0] and [H.sub.0] are the initial nonnegative matrix. The divergence deviation D(V[parallel]WH) is monotonically increasing under the rules of iteration and the value of D(V[parallel]WH) is minimum and the algorithm is convergent when W and H are local optimum [10].

5. Selection of the Extraction Scheme

After lots of surveys, we found the constraint conditions of NMF algorithm are that the matrix elements to be decomposed are nonnegative [11]. ECG signal data acquired from a lead usually form a one-dimensional matrix. For example, we can check the specific value of ECG from Massachusetts Institute of Technology and Beth Israel Hospital (MIT-BIH) [12]. The ECG signal matrix V is usually composed of positive and negative elements, so the signal matrix must be improved to extract J wave from ECG by the method of NMF. In order to satisfy the limiting conditions of NMF, we classify the ECG data into positive and negative elements. The positive elements of ECG signal matrix V should be selected from signal matrix to form a new matrix [V.sub.a] and the rest of the absolute value of the negative elements to form matrix [V.sub.b]. The positions of elements in the new matrices are the same as original, and other remained elements are considered as zero. [V.sub.a] and [V.sub.b] are decomposed, respectively, based on fixed-point NMF. The whole process can be described as [V.sub.a] = [W.sub.a][H.sub.a] and [V.sub.b] = [W.sub.b][H.sub.b]. The [W.sub.a] and [W.sub.b] are together a mixed matrix W. The [H.sub.a] and [H.sub.b] are together a J wave matrix H, where the [H.sub.a] after decomposing [V.sub.a] is J wave signals at the positive axis and, similarly, [H.sub.b] is the absolute value of J wave signal at the negative axis. Therefore, the whole J wave source signal is extracted and then we can further analyze the features of hump type J wave. The J wave extraction flowchart is shown in Figure 1.

6. The Design of Multilayer NMF

It is difficult to get the higher separation precision with single layer NMF factorization. In particular, for some complex data and the function with local minimum, the result is not always ideal. In order to improve the performance of the NMF, this paper adopted a simple approach with hierarchy and multistage execution so that continuous nonnegative matrix factorization is performed [13].

In the first step, we use NMF algorithm to perform V = [W.sub.1][H.sub.1] as a basic factorization. Then in second phase, we can use the same or different update rules and the results of the first stage to perform a similar factorization [H.sub.1] = [W.sub.2][H.sub.2]. Carrying out the multilevel factorization continuously, we could get the required result at last. This process can be repeated several times until satisfying the standard. In each step, we usually make the decomposition performance gradually improve [14]. This model is V = [W.sub.1] [W.sub.2]...[W.sub.1][H.sub.1] and the fundamental matrix is defined as V = [V.sub.1] [V.sub.2]...[V.sub.L]. It means that we build a system with many layers or L hybrid subsystem to improve the separation accuracy.

Comparing to the different layer procession of fixed-point NMF, we find that the more the number of layers is, the slower the processing speed and the extraction precision are not necessarily better. Considering the aspects of processing speed and accuracy comprehensively, we choose the 7-layer fixed-point NMF decomposition.

7. Experiment Results and Analysis

Taking the hump J wave, for example, the figures of source signals and mixed signals are shown in Figures 2 and 3. The source signals have been preprocessed and mixed signal is generated by sparse mixed randomly. In order to make the mixed signal closer to the real signal, we need to mix the preprocessed source signals randomly. By this way, we can get a hump shape J wave that is close to the real signal and the research will be meaningful.

7.1. The Result of 7-Layer NMF Separation. The figures of separation results are shown in Figures 4 and 5. An estimated signal by 7-layer fixed-point NMF decomposing shown in Figure 4 is very similar to the source signals in Figure 2. Figure 5 shows the signal interference ratios (SIRs) of estimated signal and the mixing matrix A after separation. We can find that the average SIR of estimated signals from four channels is 156.5343 dB and the SIR of the mixing matrix A is 178.8577 dB. In a word, the experiment is successful by analyzing the separation results of Figures 4 and 5.

7.2. The Result of Single Layer NMF Separation. Asacomparison, the figures of single layer NMF separation results are shown in Figures 6 and 7. Figure 6 shows estimated signals by single layer fixed-point NMF decomposing. Figure 7 shows the SIRs of four-channel estimated signals and mixing matrix A which are 109.5713 dB and 91.935 dB. It verifies the feasibility of 7-layer separating method once again and can separate the J wave from the ECG with high accuracy. In the separated signals, x1 and x2 form a normal ECG and x3 and x4 form a hump shape J wave signal.

Compared to single layer NMF, 7-layer NMF separation is more ideal. The reason is that we cannot obtain final W and H by single layer NMF in decomposition to make the objective function converge completely, which leads to the separated signals in some areas which are different to the real signal to a great extent. The W and H obtained by one-layer decomposition of multilayer NMF can make the objective function convergent and the separated signals close to the real signals. The disadvantage is that the more the layers, the longer the time spent in the separation process. After taking comprehensive consideration, we adopt the 7-layer NMF decomposition and receive the ideal results.

7.3. The Result by Using ICA Separation. ICA is the most common method in BSS problem [15], of which fast independent component analysis (Fast-ICA) is the most respective. It has the advantage of fast processing speed and stabilization performance [16]. The following is the result of using FastICA and the experiment result of the proposed scheme and the Fast-ICA are compared.

Compared to the results of our method, we can conclude that, in Figure 8, J wave and ECG signals by Fast-ICA are different from the source signals to a great extent. In Figure 9, the SIR of separating signals is 18.7589 dB, and the mixing matrix A is only 30.8253 dB. Through the separation effect of J wave, the proposed methods are superior to the classic FastICA.

8. J Wave Transformation in the Field of Digital Domain

The spectrum, power spectrum density, and spectra of the extracted hump shape J wave are shown in Figures 10,11, and 12. Comparing to Figures 10-12, we can determine that the initial signal has the characteristics of hump-shaped J wave when they are similar to the J wave. Therefore, when the J wave is not obvious and even hidden in the ST segment, this method can be used to compare the figures of the text with the patient sECG feature in the digital domain. Combing with the clinical experiences of doctors, they can provide a basis for doctors to judge whether the patients have an unconspicuous J wave.

9. Conclusion

This paper introduced a decomposing approach of multilayer fixed-point NMF and J wave can be successfully extracted from the abnormal ECG. The best initial method and decomposing layers are proposed and work well. At the same time, we conclude the spectrum, power spectrum density, and spectra of J wave in the digital domain. It can be considered as a basis for the future study of J wave characters and is helpful for doctors to diagnose the heart disease.

http://dx.doi.org/10.1155/2014/217067

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant no. 61371062); International Cooperation Project of Shanxi Province (Grant no. 20140810291); Scientific Research Project for the Returned Overseas Chinese Scholars of Shanxi Province (Grant no. 2013-032); International Cooperation Project of Shanxi Province (Grant no. 2012081031).

References

[1] X. Yang and X. Li, "J wave research history and status quo," China Cardiac Pacemaker and Electrophysiology, vol. 13, no. 1, pp. 50-51, 1999.

[2] G. X. Yan, "Abnormal J wave syndrome," Journal of Clinical Electrocardiology, vol. 1, pp. 3-9, 2007

[3] G. X. Yan, Q. H. Yao, D. Q. Wang, and C. C. Cui, "Electrocardiographic J wave and J wave syndromes," Chinese Journal of Cardiac Arrhythmias, vol. 8, no. 6, pp. 360-365, 2004.

[4] Z. Zhao, H. Lu, and C. Xu, "An blind source separation algorithm based on constrained NMF," Piezoelectrics and Acoustooptics, vol. 32, no. 6, pp. 1049-1052, 2010.

[5] F. Potet, P. Mabo, G. le Coq et al., "Novel Brugada SCN5A mutation leading to ST segment elevation in the inferior or the right precordial leads," Journal of Cardiovascular Electrophysiology, vol. 14, no. 2, pp. 200-203, 2003.

[6] M. D. Plumbley, "Algorithms for nonnegative independent component analysis," IEEE Transactions on Neural Networks, vol. 14, no. 3, pp. 534-543, 2003.

[7] A. Cichocki, R. Zdunek, and S.-I. Amari, "New algorithms for non-negative matrix factorization in applications to blind source separation," in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '06), vol. 5, pp. V621-V624, Toulouse, France, May 2006.

[8] A. Pascual-Montano, J. M. Carazo, K. Kochi, D. Lehmann, and R. D. Pascual-Marqui, "Nonsmooth nonnegative matrix factorization (nsNMF)," IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 28, no. 3, pp. 403-415, 2006.

[9] F.-Y. Wang, Y. Wang, T.-H. Chan, and C.-Y. Chi, "Blind separation of multichannel biomedical image patterns by nonnegative least-correlated component analysis," in Pattern Recognition in Bioinformatics, vol. 4146 of Lecture Notes in Computer Science, pp. 151-162, 2006.

[10] W. S. Ouedraogo, A. Souloumiac, M. Ja'idane, and C. Jutten, "Non-negative blind source separation algorithm based on minimum aperture simplicial cone," IEEE Transactions on Signal Processing, vol. 62, no. 2, pp. 376-389, 2014.

[11] P. O. Hoyer, "Non-negative matrix factorization with sparseness constraints," Journal of Machine Learning Research, vol. 5, pp. 1457-1469, 2004.

[12] Y. Luo, R. H. Hargraves, A. Belle et al., "A hierarchical method for removal of baseline drift from biomedical signals: application in ECG analysis," The Scientific World Journal, vol. 2013, Article ID 896056,10 pages, 2013.

[13] A. Cichocki and R. Zdunek, "Multilayer nonnegative matrix factorisation," Electronics Letters, vol. 42, no. 16, pp. 947-948, 2006.

[14] A. Cichocki and R. Zdunek, "Multilayer nonnegative matrix factorization using projected gradient approaches," International Journal ofNeural Systems, vol. 17, no. 6, pp. 431-446, 2007.

[15] Y. Matsuda and K. Yamaguchi, "Linear multilayer ica using adaptive pca," Neural Processing Letters, vol. 30, no. 2, pp. 133-144, 2009.

[16] D. Li, J. Zhao, H. Liu, and D. Hao, "The application ofFastICA combined with related function in blind signal separation," Mathematical Problems in Engineering, vol. 2014, Article ID 953745, 9 pages, 2014.

Deng-ao Li, Jing-ang Lv, Ju-min Zhao, and Jin Zhang

College of Information Engineering, Taiyuan University of Technology, Taiyuan 030024, China

Correspondence should be addressed to Deng-ao Li; lidengao@tyut.edu.cn

Received 20 June 2014; Revised 25 August 2014; Accepted 1 September 2014; Published 12 October 2014

Academic Editor: Irena Cosic

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article; nonnegative matrix factorization |
---|---|

Author: | Li, Deng-ao; Lv, Jing-ang; Zhao, Ju-min; Zhang, Jin |

Publication: | Computational and Mathematical Methods in Medicine |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 3222 |

Previous Article: | Privacy-preserving self-helped medical diagnosis scheme based on secure two-party computation in wireless sensor networks. |

Next Article: | Fetal ECG extraction from abdominal signals: a review on suppression of fundamental power line interference component and its harmonics. |

Topics: |