Next Article in Journal
An Improved Time-Series Model Considering Rheological Parameters for Surface Deformation Monitoring of Soft Clay Subgrade
Previous Article in Journal
Hyperspectral Imaging in Environmental Monitoring: A Review of Recent Developments and Technological Advances in Compact Field Deployable Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Different Electrocardiography with Vectorcardiography Transformations

Department of Cybernetics and Biomedical Engineering, Faculty of Electrical Engineering and Computer Science, VSB–Technical University of Ostrava, 17. listopadu 15, 708 33 Ostrava, Czech Republic
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(14), 3072; https://doi.org/10.3390/s19143072
Submission received: 11 June 2019 / Revised: 4 July 2019 / Accepted: 9 July 2019 / Published: 11 July 2019
(This article belongs to the Section Biosensors)

Abstract

:
This paper deals with transformations from electrocardiographic (ECG) to vectorcardiographic (VCG) leads. VCG provides better sensitivity, for example for the detection of myocardial infarction, ischemia, and hypertrophy. However, in clinical practice, measurement of VCG is not usually used because it requires additional electrodes placed on the patient’s body. Instead, mathematical transformations are used for deriving VCG from 12-leads ECG. In this work, Kors quasi-orthogonal transformation, inverse Dower transformation, Kors regression transformation, and linear regression-based transformations for deriving P wave (PLSV) and QRS complex (QLSV) are implemented and compared. These transformation methods were not yet compared before, so we have selected them for this paper. Transformation methods were compared for the data from the Physikalisch-Technische Bundesanstalt (PTB) database and their accuracy was evaluated using a mean squared error (MSE) and a correlation coefficient (R) between the derived and directly measured Frank’s leads. Based on the statistical analysis, Kors regression transformation was significantly more accurate for the derivation of the X and Y leads than the others. For the Z lead, there were no statistically significant differences in the medians between Kors regression transformation and the PLSV and QLSV methods. This paper thoroughly compared multiple VCG transformation methods to conventional VCG Frank’s orthogonal lead system, used in clinical practice.

1. Introduction

Electrocardiography (ECG) is currently the most widely used method of sensing electrical activity in the human heart. Vectorcardiography (VCG) is an investigation method that was previously compared to ECG several times and was evaluated as a useful investigation method [1,2,3]. It achieved a higher sensitivity for the detection of hypertrophy as well as ischaemic heart disease [4,5,6,7]. Three leads scanned in three mutually orthogonal axes (vertical, transversal, and sagittal) are used for VCG measurements, but this is not commonly used in clinical practice due to the more complex interpretation of the recordings. Mathematical transformations can be made between ECG and VCG leads. Previously, ECG transformations from VCG leads were used, but later only 12-lead ECG spread and VCG leads were derived from ECG leads instead of measurements. These transformations are based on transformation coefficients that have been defined by various authors [8,9,10]. For example, the authors used thoracic models to obtain the coefficients, or, based on the similarity of one of the ECG leads, they determined the orthogonal lead. Transformation methods have also been developed based on a regression approach for multiple recordings of simultaneously measured ECG and VCG.
VCG is a diagnostic and less specific method than the standard 12-lead ECG measurement [1]. This method can be measured using several different lead systems. In clinical practice, Frank’s orthogonal lead system is the most widely used. It is measured using three orthogonal leads described by Frank, which measure cardiac activity with the same sensitivity. The problem of needing to connect multiple electrodes to measure both ECG and VCG was solved by synthesis from the currently recorded ECG leads. Compared to ECG, VCG has been recognized as a very useful investigative clinical method for cardiac activity sensing [1,3]. In addition, it shows the balance of cardiac cycle forces at any moment and contains more information in comparison to ECG that may have diagnostic significance, such as the display of temporal or phase relationships between complexes. The importance of VCG is, for example, in the detection and localization of acute myocardial infarction, right ventricular hypertrophy, and Tawar arm blockade. Despite all the benefits of VCG, ECG is increasingly being used.
VCG recognizes three planes, sagittal, transversal, and frontal. The frontal one is located between the XY leads, the transverse one between the XZ leads, and the sagittal one between the YZ leads. The electrical activity of the heart is described by three loops that represent the individual phases, and their contours, rotation, or the direction of the cardiac axis are taken into account. The first loop corresponds to the P wave, the second one, which is also the largest loop, corresponds to the QRS complex, and the third one corresponds to the T wave. The loops can be displayed in a 1-D image as three scalar recordings (Figure 1), in a 2-D image using three planes (Figure 2), or in one 3-D image (Figure 3). The axes are represented as time and voltage dependence in the 1-D image. In the case of 2-D and 3-D images, the dependence between the individual leads in mV is displayed. In the analysis, the greatest attention was paid to the QRS complex, which is oval in shape and faces in the same direction as the cardiac axis of the heart [2].
The Frank’s orthogonal lead system is most widely used for its simplicity and is almost orthogonal. This is a bipolar circuit. It consists of seven electrodes whose position is marked with capital letters: I, E, C, A, M, F, and H.Each electrode has its specific location, see Figure 4. The E electrode is located on the front, in the middle of the chest, and the M electrode is located opposite, in the rear, on the back. The I electrode is located in the right central axillary line and the A electrode is in the left central axillary line. The position of the C electrode is obtained by halving the distance between the A electrode and E electrode. The F Electrode is positioned on the left leg and the H electrode on the neck. Furthermore, one electrode attached to the right leg is used as a ground electrode. On the X-axis, the left side of the body is positive with respect to the right side of the body, on the Y-axis, the left foot is a positive electrode with respect to the head, and on the Z-axis, the back of the chest is positive with respect to the front [2]. The individual electrodes are routed to a resistor network. Malmivuo et al., 1995 [11] explained the correct values of each resistors in their book. The signals on the leads are derived using mathematical Equations (1)–(3), where P represents the potentials on the individual electrode clips and I, E, C, A, M, F, H represent individual electrodes [8]. Macfarlane et al., 2010 [12] explained calculations of Equations (1)–(3), which are based on the resistor network in Figure 4. In addition to Frank’s orthogonal leads, other leads, such as McFee–Parungao leads [13] and Schmitt–Simonson [14] leads called Svec-III, are used for sensing VCG [15,16,17].
P X = 0.610 · A + 0.171 · C 0.781 · I .
P Y = 0.655 · F + 0.345 · M 1.000 · H .
P Z = 0.133 · A + 0.736 · M 0.264 · I 0.374 · E 0.231 · C .
There are many different transformation methods of 3-lead VCG to 8-lead ECG, or 3-lead ECG to 8-lead VCG. Kors et al., 1990 [18] compared the Kors quasi-orthogonal transform to the inverse Dower transform, the Kors regression transform, and the Frank’s orthogonal lead system. A total of 90 recordings were used. The regression method achieved the best results in this study, while the quasi-orthogonal method achieved the worst results. Other quasi-orthogonal transforms were described in articles [19,20,21,22]. The authors use different coefficients in equations to perform ECG to VCG transformations. Burger et al., 1952 [23] introduced first linear transformation method. They showed that it would be possible to transform ECG to VCG. Levkov, 1987 [3] introduced five different transformation methods in his work. One transformation method is based on one-dipole model and four transformation methods are derived by regression. His methods provide high accuracy, particularly th the methods based on regression. Edenbrandt et al., 1994 [4] tested the inverse Dower transformation. Authors compared this method with the directly measured Frank’s leads, the quasi-orthogonal method described by Bjerl, and the quasi-orthogonal method implemented by Marquette Electronics Inc. control systems. The comparison was performed on 80 people with different types of heart attacks. The inverse Dower transformation achieved the best QRS complex image with respect to other methods and its recording was the most similar to the Frank’s. Guillem et al., 2006 [9] tested linear regression-based transformations for deriving the P wave (PLSV) and QRS complex (QLSV). They compared these methods with the inverse Dower transformation methods during the QRS complex and P wave detection on 247 recordings. In the P wave transformation, the PLSV method achieved the best results and the QLSV method achieved the best results for QRS complex detection. Both methods achieved better results than the inverse Dower method. Acar and Koymen, 1999 [24] optimized the singular value decomposition (SVD) transformation to derive VCG indications. They used data from 23 patients. Their approach does not perform transformation of three VCG leads but provides three non-correlated orthogonal leads. Dawson et al., 2009 [25] used the Dower transform and statistical method called affine transform to derive 8-lead ECG from 3-lead VCG and their inverse forms to derive 3-lead VCG from 8-lead ECG. They concluded that for subjects with myocardial infarction and healthy control subjects, affine transform provides better accuracy than Dower transform. Maheshwari et al., 2016 [26] showed that principal component analysis (PCA) could also be used for deriving the VCG from ECG. They compared PCA with the inverse Dower transform and the Kors regression transform. Their results support the theoretical basis behind the PCA-based reconstruction methodology and show better performance then the compared transforms. Vozda et al., 2014 [27] used transformation methods based on artificial neural networks. They showed that these transformation methods offer even higher accuracy then other reviewed transforms.
In this paper, Kors quasi-orthogonal transform, inverse Dower transform, and Kors regression transform were chosen, because they are the most commonly used transformation methods. A very large number of transformation methods and matrices are optimizing the transformations of the QRS complex. Methods based on linear regression transformations are focused not only on the QRS complex, but also on optimizing the transformations of the P wave. That is the main reason why this transform was chosen.

2. Methods

The first transformations were designed to derive ECG from VCG leads. This was due to a reduction in the number of electrodes attached compared to the current ECG and VCG measurements. The first one was the Dower transformation, which was derived using Frank’s thoracic model and by calculating the linear coefficients [3]. Later, only 12-lead ECG spread in clinical practice while VCG ceased to be measured. Therefore, VCG leads were more often derived from ECG leads. 3 VCG leads X, Y, and Z are derived by transforming the 12-lead ECG. To transform the 12-lead ECG into 3 lead VCG, only the values of 6 thoracic leads, lead I and lead II, i.e., 8 linearly independent leads are used. The remaining 4 leads are linearly dependent and can therefore be calculated. At the same time, linearly dependent leads are not used because they could exacerbate the resulting transformation [3,8,9]. The mathematical transformation is based on the multiplication of matrices according to Equation (4), where V is a VCG matrix whose lines correspond to 3 VCG leads, M is a transformation matrix of the selected method, and E is a matrix whose lines are the individual ECG leads.
V = M · E .

2.1. Kors Quasi-Orthogonal Transformation

This method is based on the assumption that the X lead is similar to the V6 lead to a certain extent, and the Y lead is similar to lead II. Z Lead partially corresponds to the similarity of the negative half of the V2 lead. Hence, all input leads have coefficients equal to zero, except for three selected leads similar to orthogonal leads. The resulting VCG leads can be obtained using Equations (5)–(7) [10,18].
X = V 6 .
Y = I I .
Z = 0.5 · V 2 .

2.2. Inverse Dower Transformation

This method is based on the mathematical pseudoinversion of the Dower’s method of deriving ECG from VCG. Orthogonal leads X, Y, and Z are expressed by a linear combination of 8 leads. Transformation coefficients are represented in Table 1 [4,8,28,29].

2.3. Kors Regression Transformation

This is a statistical method for which Kors used mathematical regression [8,18]. The regression method is used when recordings from a large number of patients are available. Kors, while measuring ECG and VCG, derived mathematically the coefficients for the regression method, see Table 2. These coefficients allow ECG to be transformed into VCG. The reconstruction coefficients are obtained by minimizing the mean squared error (MSE) according to Equation (9) [8,18,29,30].

2.4. Linear Regression-Based Transformations

Least-square value (LSV) is based on a regression approach. By minimizing the MSE, the coefficients for two methods focused on another ECG interval were derived, see Table 3 and Table 4. Using the interval from the beginning to the end of the P wave, a matrix for transformation with emphasis on P wave (PLSV transform) was obtained. Using the interval from the beginning to the end of the QRS complex, a matrix for transformation with emphasis on the QRS complex (QLSV transform) was obtained [9,10].

3. Evaluation Parameters

The Physikalisch-Technische Bundesanstalt (PTB) [31,32,33] diagnostic database was used to test the ECG to VCG transformation. The database contains recordings measured by the Benjamin Franklin Cardiology Department, which were measured in healthy patients and patients with various heart diseases. Of the total of 289 entities, 549 recordings are recorded in the database. A total of 1 to 5 recordings were measured for each patient. This database has mixed recordings that consists of healthy entities and patients with pathological heart diseases. The individual recordings contain 15 simultaneously measured leads. The leads measured comprise of 12 ECG leads (Eithoven, Goldberg, and Wilson leads) and 3 VCG leads (Frank’s leads). The sampling frequency of each signal is 1 kHz. For the evaluation, 50 recordings of healthy patients were chosen from this database. This database is publicly accessible on Physionet [31,32,33].
The correlation coefficient (R) and the mean squared error (MSE) were chosen to test the accuracy of the transformation by the individual methods in relation to the reference, Frank’s leads. Subsequently, a statistical analysis using the Mann–Whitney nonparametric median test was selected.

3.1. Correlation Coefficient

Correlation indicates the relationship between two recordings or variables. It is often used in statistics, where a linear dependency between recordings is sought. The correlation coefficient does not have a unit and its value is in the interval of <−1,1>. If the correlation value is 0, there is no correlation between the two variables. If the value is −1, it means indirect linearity between the recordings. The value of a correlation equaling 1 then means direct linearity, i.e., similarity between the recordings evaluated. The equation for calculating the correlation coefficient is given by Equation (8), where V is the original value of the VCG measured, D V is the value of the VCG derived from the ECG, and n is the number of samples [34].
R = i = 1 n ( V · D V ) i = 1 n V 2 i = 1 n D V 2 .

3.2. Mean Squared Error

The mean squared error compares two recordings based on similarities or differences. To evaluate the transformation’s accuracy, we compare MSE calculated between the transformed signal and the directly measured VCG in each lead. The VCG uses MSE to evaluate the amplitude error of the recording. The unit is the squared unit of the variable measured. The closer the result is to zero, the more similar or identical the recordings are [35]. The equation for calculating MSE is given by Equation (9), where V is the original value of the VCG measured, D V is the value of the VCG derived from the ECG, and n is the number of samples.
M S E = 1 n i = 1 n ( V i D V i ) 2 .

3.3. Statistical Analysis

When testing the hypotheses, the significance level of α is determined. 95% is considered a reliable estimate and the significance level is the remaining 5% (0.05). If the difference between the two files tested is less than α , then the recordings are different at the significance level [36,37]. In statistics, the mean and the median are most often determined for performing a statistical test. Statistical analysis is based on hypotheses testing. The zero hypothesis (H 0 ) represents the equilibrium (e.g., equality between the median results) and the alternative hypothesis (H A ) then expresses the equilibrium violation (e.g., the difference between the median results). Based on the selected statistical test, either the zero hypothesis is rejected and the alternative hypothesis is accepted, or the zero hypothesis is not rejected and the alternative hypothesis is rejected. When an alternative hypothesis is reached between two compared transformations, the testing will be assessed by the significance level. For statistical analysis, we chose from statistical tests, which are divided into parametric and nonparametric ones. Parametric tests work with normal input data distribution and have a higher test power than nonparametric tests. The data used in this article do not originate from the normal distribution, therefore, the Mann–Whitney nonparametric median test is used [36,37].

4. Results

Before conducting individual transformations, signal pre-processing was performed using a frequency selective filter with a finite impulse response (FIR filter). Based on recommendations [38,39,40], filter cut-off frequencies should be selected from 0.05 to 100 Hz. Since we only analyzed the recordings of healthy volunteers, the limit values were set from 0.2 to 100 Hz with a filter order of 500, which is in line with the recommendation stated in the articles [38,39,40].
Subsequently, R-oscillations were detected using the Pan–Tompkins detector. The individual positions of the R-oscillations determine the beginning of the heart rate as R 300 ms and the end of the heart rate as R + 400 ms. Thus, all heart rates of the individual leads were cut out. The resulting transformations could be influenced by extrasystoles, so they were removed using the library function in MATLAB. Once the extrasystoles were removed, the individual heart rates were averaged, so that the average heart rate for each ECG and VCG lead was created to perform the transformation and the transformation accuracy check.
The test was performed on the recordings from the PTB database and the MSE and R result vector in relation to the reference Frank’s leads was obtained for each method. Since the mean value of the results was different from the median of the results, it was determined that the data did not come from a normal distribution. Statistical testing of data normality was performed at work [27]. Therefore, the Mann–Whitney nonparametric median test was selected for statistical analysis. If the resulting p-value was less than the significance level of α = 0.05 , the zero hypothesis would be rejected in favor of the alternative hypothesis. This would mean that the difference between the medians compared was statistically significant. The determination of the zero and alternative hypotheses can be seen in Equations (10) and (11), where the zero hypothesis indicates that the median results were statistically identical, and the alternative hypothesis indicates that the difference between the median results was statistically significant. The test was bilateral, so the medians for Kors regression transformation and the method compared must be determined.
H 0 : x ˜ 1 = x ˜ 2 .
H A : x ˜ 1 x ˜ 2 .

4.1. Visual Evaluation

The testing focused on the regression method described by Kors because, according to studies, this transformation has achieved the best results in comparison to other transformations. Statistical analysis was performed for recordings measured on healthy volunteers that were received from the PTB database. For comparison, box plots are created, see Figure 5 and Figure 6. In the figures and tables, the result of the Kors quasi-orthogonal transformation will be referred to as Quasi, the result of inverse Dower transformation will be referred to as Dower, the result of Kors regression transformation will be referred to as Regress, the result of the PLSV method will be referred to as PLSV, and the result of the QLSV method will be referred to as QLSV. It can be concluded from the figure that Kors regression transformation was the most accurate because the median MSE was lower than the other methods for X and Y leads, see Figure 5. Similarly, it can be seen that the median R was higher than the other methods for the X and Y leads, see Figure 6. In the box plots, the median is represented by the horizontal line and the mean by the circle. Red crosses represent outliers, which are commonly used in statistics to represent samples which are removed from evaluation. Statistical verification that this method is better must be proved by a statistical test.

4.2. Statistical Analysis of MSE Results

Figure 5 indicates that the MSE medians of Kors regression transformation for X and Y leads are lower than the medians of the other transformation methods and, for Z lead, it could not be assessed. According to the statistical test, all p-values for the X and Y leads were less than the significance level of α = 0.05 , see Table 5, and, thus, it was impossible to reject the zero hypothesis. This means that Kors regression transformation has been recognized as more accurate for X and Y leads than all other transformation methods, based on the medians of the MSE results.
For Z lead, the p-values between Kors regression transformation and Kors quasi-orthogonal transformation, and between Kors regression transformation and inverse Dower transformation were less than the significance level of α = 0.05 , therefore, the zero hypothesis could not be rejected here. Thus, Kors regression transformation was recognized as more accurate for Z lead than Kors quasi-orthogonal transformation and inverse Dower transformation, based on the median of the MSE results. For Z lead, the p-values between Kors regression transformation and the PLSV method, and Kors regression transformation and the QLSV method were greater than the significance level of α = 0.05 , therefore, the zero hypothesis was rejected in favor of the alternative hypothesis. This indicates that, for Z lead, it cannot be asserted that Kors regression transformation would achieve a more accurate transformation than the PLSV and QLSV methods, hence, the difference between the methods was not statistically significant.

4.3. Statistical Analysis of R Results

Figure 6 indicates that the R medians of Kors regression transformation for X and Y leads were higher than the medians of the other transformation methods and, for Z lead, it could not be assessed with certainty. According to the statistical test, all p-values for the X and Y leads were higher than the significance level of α = 0.05 , see Table 6, and, thus, it was impossible to reject the zero hypothesis. This means that Kors regression transformation has been recognized as more accurate for X and Y leads than all other transformation methods, based on the medians of the R results. For Z lead, the p-values between Kors regression transformation and Kors quasi-orthogonal transformation, and between Kors regression transformation and inverse Dower transformation were less than the significance level of α = 0.05 , therefore, it was impossible to reject the zero hypothesis here, either. Thus, Kors regression transformation was recognized as more accurate for Z lead than Kors quasi-orthogonal transformation and inverse Dower transformation, based on the medians of the R results. For Z lead, the p-valuess between Kors regression transformation and the PLSV method, and Kors regression transformation and the QLSV method were greater than the significance level of α = 0.05 , therefore, the zero hypothesis was rejected in favor of the alternative hypothesis. This indicates that, for Z lead, it cannot be asserted that Kors regression transformation would achieve a more accurate transformation than the PLSV and QLSV methods based on the medians of the R results, hence, the difference between the methods was not statistically significant.

5. Discussion

In Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, you can see 2-D and 3-D comparisons of tested tranformations. For visual comparison, we chose four random recordings out of the whole tested data set. 2-D figures were selected, as they are ideal for comparisons to the reference Frank’s orthogonal lead system. According to the MSE and R parameters, Kors regression transformation provided the best results, which also corresponds to the mentioned figures.
Kors regression transformation has advantages over Frank’s orthogonal lead system. It is much easier to identify myocardial infarction [41,42,43,44,45,46], hypertrophy [47,48,49], and ischaemic heart disease [50,51] from the VCG compared to the ECG. However it is important to properly estimate the VCG from the ECG. Transformations need a lower number of required leads in comparison to Frank’s orthogonal lead system combined with the conventional ECG system. To better diagnose heart diseases, doctors often require both signals at once. In the case of a successful extraction of the VCG from the ECG, it is possible to identify multiple symptoms, such as the QRS-T angle, QRS-T spatial angle, heart axis slope, and QRS loop vector size [52]. A high accuracy of Kors regression transformation could lead to a complete deprecation of Frank’s orthogonal lead system, which also provides higher comfort for patients.
Apart from the presented solutions, a number of different VCG transformations were tested by other teams, with varying results. Other teams carried out experiments based on affine transform, PCA, or transformation methods based on artificial neural networks that offer even higher accuracy [25,26,27]. These transformations will be a basis for our future research.
VCG systems based on transformations are much cheaper to conventional Frank’s orthogonal lead system, since it does not need a dedicated hardware platform. Newly graduated doctors could potentially choose this system for their clinical practice, since it has significantly lower initial costs.

6. Conclusions

The accuracy of the ECG transformation to the VCG is very important for possibly determining symptoms of hypertrophy, ischaemic heart disease, etc. A total of five different transformation methods, which were tested on PTB database recordings, were selected for the testing. The results of the transformations were compared with the original Frank’s leads measured on the basis of the MSE and R. From the results, the medians of the individual results were determined and statistical analysis was performed with emphasis on Kors regression transformation in relation to other transformation methods. The results were determined separately for each lead. A nonparametric statistical test was used to calculate the p-values by means of which the significance of the difference of values between the transformations was determined.
During testing, Kors regression transformation was statistically the most accurate for X and Y leads compared to other transformation methods. The smallest accuracy of transformation was provided by Kors quasi-orthogonal transformation. The difference between Kors regression transformation and the PLSV and QLSV methods was not statistically significant for Z lead. Overall, the results showed that Kors regression transformation was the most accurate of the methods tested and could be significant for transforming the ECG into the VCG for automatic computer analysis of ECG recordings based on VCG symptoms. Evaluating the accuracy of transformation methods for patients suffering from some heart disease and comparing the effects of transformations on pathological recordings could be the subject of the next work. Furthermore, newer types of transformation methods could be added. In the future, ECG transformations into VCG could be used in clinical practice, which could facilitate diagnosing heart diseases.
2-D/3-D graphs (Figure 7, Figure 8, Figure 9 and Figure 10) were reconstructed after processing tested transformations via statistical analysis. 2-D/3-D graphs consisting of VCG loops are very important in clinical practice, as they are provide more information than the ECG and therefore help doctors identify heart diseases. Comparisons of visual signals correlates with statistical analysis. Kors regression transformation tend to be the most similar to Frank’s orthogonal lead system, so its safe to say that it could be used for morphological analysis of 2-D/3-D graphs (for example pathological phenomena).

Author Contributions

Conceptualization, R.J., and R.M.; Methodology, R.J., R.M., and L.D.; Software R.J.; Validation, R.J., R.M., and L.D.; Formal Analysis, R.J., and R.M.; Investigation, R.J., R.M., and L.D.; Resources, R.J.; Data Curation, R.J.; Writing—Original Draft Preparation, R.J., R.M., and L.D.; Writing—Review & Editing, R.J., R.M., and L.D.; Visualization, R.J.; Supervision, R.M.; Project Administration, R.M.; Funding Acquisition, R.M.

Funding

This article was supported by the Ministry of Education of the Czech Republic (Project No. SP2019/85 and SP2019/118). This work was supported by the European Regional Development Fund in the Research Centre of Advanced Mechatronic Systems project, project number CZ.02.1.01/0.0/0.0/16_019/0000867 within the Operational Programme Research, Development and Education. This work was supported by the European Regional Development Fund in A Research Platform focused on Industry 4.0 and Robotics in Ostrava project, CZ.02.1.01/0.0/0.0/17_049/0008425 within the Operational Programme Research, Development and Education.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Rubel, P.; Benhadid, I.; Fayn, J. Quantitative assessment of eight different methods for synthesizing Frank VCGs from simultaneously recorded standard ECG leads. J. Electrocardiol. 1991, 24, 197–202. [Google Scholar] [CrossRef]
  2. Winsor, T. Primer of Vectorcardiography; Lea & Febiger: Philadelphia, PA, USA, 1972. [Google Scholar]
  3. Levkov, C.L. Orthogonal electrocardiogram derived from the limb and chest electrodes of the conventional 12-lead system. Med. Biol. Eng. Comput. 1987, 25, 155. [Google Scholar] [CrossRef]
  4. Edenbrandt, L.; Houston, A.; Macfarlane, P.W. Vectorcardiograms synthesized from 12-lead ECGs: A new method applied in 1792 healthy children. Pediatr. Cardiol. 1994, 15, 21–26. [Google Scholar] [CrossRef] [PubMed]
  5. Lingman, M.; Hartford, M.; Karlsson, T.; Herlitz, J.; Rubulis, A.; Caidahl, K.; Bergfeldt, L. Transient repolarization alterations dominate the initial phase of an acute anterior infarction—A vectorcardiography study. J. Electrocardiol. 2014, 47, 478–485. [Google Scholar] [CrossRef]
  6. Sederholm, M. The origin of monitoring of acute myocardial infarction with continuous vectorcardiography. J. Electrocardiol. 2014, 47, 418–424. [Google Scholar] [CrossRef]
  7. Cortez, D.; Bos, J.M.; Ackerman, M.J. Vectorcardiography identifies patients with electrocardiographically concealed long QT syndrome. Heart Rhythm 2017, 14, 894–899. [Google Scholar] [CrossRef] [PubMed]
  8. Macfarlane, P.W.; Edenbrandt, L.; Pahlm, O. 12-Lead Vectorcardiography; Butterworth-Heinemann: Oxford, UK, 1995. [Google Scholar]
  9. Guillem, M.S.; Sahakian, A.V.; Swiryn, S. Derivation of orthogonal leads from the 12-lead ECG. Accuracy of a single transform for the derivation of atrial and ventricular waves. In Proceedings of the 2006 Computers in Cardiology, Valencia, Spain, 17–20 September 2006; pp. 249–252, ISBN 978-1-4244-2532-7. [Google Scholar]
  10. Hyttinen, J.A.; Viik, J.J.; Eskola, H.; Malmivuo, J.A. Optimization and comparison of derived Frank VECG lead systems employing an accurate thorax model. In Proceedings of the Computers in Cardiology 1995, Vienna, Austria, 10–13 September 1995; pp. 385–388. [Google Scholar]
  11. Malmivuo, P.; Malmivuo, J.A.; Plonsey, R. Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  12. Macfarlane, P.W.; Van Oosterom, A.; Pahlm, O.; Kligfield, P.; Janse, M.; Camm, J. Comprehensive Electrocardiology; Springer Science Business Media: Berlin, Germany, 2010. [Google Scholar]
  13. McFee, R.; Parungao, A. An orthogonal lead system for clinical electrocardiography. Am. Heart J. 1961, 62, 93–100. [Google Scholar] [CrossRef]
  14. Simonson, E.; Schmitt, O.H. Present status of vectorcardiography. AMA Arch. Int. Med. 1955, 96, 574–590. [Google Scholar]
  15. Malmivuo, J.A. The SVEC III vectorcardiographic lead system. IEEE Eng. Med. Biol. Mag. 2004, 23, 47–51. [Google Scholar] [CrossRef] [PubMed]
  16. Kudaiberdiev, Z. Vectorcardiographic assessment of acute hypoxia effects in pulmonary hypertension due to chronic bronchitis. J. Electrocardiol. 2007, 40, S64. [Google Scholar] [CrossRef]
  17. Fischmann, E.J.; Elliott, B.J. Experimental comparison of “parallel grid leads” with simple bipolar, and the SVEC-III, Frank, and McFee-Parungao systems. I. Sagittal leads. Am. Heart J. 1964, 67, 792–803. [Google Scholar] [CrossRef]
  18. Kors, J.A.; Van Herpen, G.; Sittig, A.C.; Van Bemmel, J.H. Reconstruction of the Frank vectorcardiogram from standard electrocardiographic leads: Diagnostic comparison of different methods. Eur. Heart J. 1990, 11, 1083–1092. [Google Scholar] [CrossRef] [PubMed]
  19. Bjerle, P.; Arvedson, O. Comparison of Frank vectorcardiogram with two different vectorcardiograms derived from conventional ECG-leads. Proc. Eng. Found. 1986, 11, 13–26. [Google Scholar]
  20. van Bemmel, J.H.; Kors, J.A.; van Herpen, G. Combination of diagnostic classifications from ECG and VCG computer interpretations. J. Electrocardiol. 1992, 25, 126–130. [Google Scholar] [CrossRef]
  21. Edenbrandt, L.; Pahlm, O. Vectorcardiogram synthesized from a 12-lead ECG: Superiority of the inverse Dower matrix. J. Electrocardiol. 1988, 21, 361–367. [Google Scholar] [CrossRef]
  22. Schreck, D.M.; Fishberg, R.D. Derivation of the 12-lead electrocardiogram and 3-lead vectorcardiogram. Am. J. Emerg. Med. 2013, 31, 1183–1190. [Google Scholar] [CrossRef]
  23. Burger, H.C.; Van Milaan, J.B.; Den Boer, W. Comparison of different systems of vectorcardiography. Br. Heart J. 1952, 14, 401. [Google Scholar] [CrossRef]
  24. Acar, B.; Koymen, H. SVD-based on-line exercise ECG signal orthogonalization. IEEE Trans. Biomed. Eng. 1999, 46, 311–321. [Google Scholar] [CrossRef] [Green Version]
  25. Dawson, D.; Yang, H.; Malshe, M.; Bukkapatnam, S.T.; Benjamin, B.; Komanduri, R. Linear affine transformations between 3-lead (Frank XYZ leads) vectorcardiogram and 12-lead electrocardiogram signals. J. Electrocardiol. 2009, 42, 622–630. [Google Scholar] [CrossRef]
  26. Maheshwari, S.; Acharyya, A.; Schiariti, M.; Puddu, P.E. Frank vectorcardiographic system from standard 12 lead ECG: An effort to enhance cardiovascular diagnosis. J. Electrocardiol. 2016, 49, 231–242. [Google Scholar] [CrossRef]
  27. Vozda, M.; Peterek, T.; Cerny, M. Novel Method for Deriving Vectorcardiographic Leads Based on Artificial Neural Networks. In Proceedings of the 41st International Congress on Electrocardiol, Bratislava, Slovakia, 4–7 June 2014. [Google Scholar]
  28. Guillem, M.S.; Climent, A.M.; Bollmann, A.; Husser, D.; Millet, J.; Castells, F. Limitations of Dower’s inverse transform for the study of atrial loops during atrial fibrillation. Pacing Clin. Electrophysiol. 2009, 32, 972–980. [Google Scholar] [CrossRef] [PubMed]
  29. Cortez, D.L.; Schlegel, T.T. When deriving the spatial QRS-T angle from the 12-lead electrocardiogram, which transform is more Frank: Regression or inverse Dower? J. Electrocardiol. 2010, 43, 302–309. [Google Scholar] [CrossRef] [PubMed]
  30. Cortez, D.; Sharma, N.; Devers, C.; Devers, E.; Schlegel, T.T. Visual transform applications for estimating the spatial QRS–T angle from the conventional 12-lead ECG: Kors is still most Frank. J. Electrocardiol. 2014, 47, 12–19. [Google Scholar] [CrossRef] [PubMed]
  31. Bousseljot, R.; Kreiseler, D.; Schnabel, A. Nutzung der EKG-Signaldatenbank CARDIODAT der PTB uber das Internet. Biomed. Tech. Biomed. Eng. 1995, 40, 317–318. [Google Scholar] [CrossRef]
  32. Laguna, P.; Mark, R.G.; Goldberg, A.; Moody, G.B. A database for evaluation of algorithms for measurement of QT and other waveform intervals in the ECG. In Proceedings of the Computers in Cardiology 1997, Lund, Sweden, 7–10 September 1997; pp. 673–676. [Google Scholar]
  33. Goldberger, A.L.; Amaral, L.A.; Glass, L.; Hausdorff, J.M.; Ivanov, P.C.; Mark, R.G.; Mietus, J.E.; Moody, G.B.; Peng, C.K.; Stanley, H.E. PhysioBank, PhysioToolkit, and PhysioNet: Components of a new research resource for complex physiologic signals. Circulation 2000, 101, e215–e220. [Google Scholar] [CrossRef] [PubMed]
  34. Ahlgren, P.; Jarneving, B.; Rousseau, R. Requirements for a cocitation similarity measure, with special reference to Pearson’s correlation coefficient. J. Am. Soc. Inf. Sci. Technol. 2003, 54, 550–560. [Google Scholar] [CrossRef]
  35. Wang, Z.; Bovik, A.C. Mean squared error: Love it or leave it? A new look at signal fidelity measures. IEEE Signal Process. Mag. 2009, 26, 98–117. [Google Scholar] [CrossRef]
  36. Andel, J. Fundamentals of Mathematical Statistics; SHTL: Prague, Czech Republic, 1978. [Google Scholar]
  37. Bris, R.; Litschmannova, M. Statistics I. For Combined and Distance Study; VSB–TU Ostrava: Ostrava, Czech Republic, 2004. [Google Scholar]
  38. Buendia-Fuentes, F.; Arnau-Vives, M.A.; Arnau-Vives, A.; Jimenez-Jimenez, Y.; Rueda-Soriano, J.; Zorio-Grima, E.; Osa-Sáez, A.; Martínez-Dolz, L.V.; Almenar-Bonet, L.; Palencia-Perez, M.A. High-bandpass filters in electrocardiography: Source of error in the interpretation of the ST segment. ISRN Cardiol. 2012, 2012, 706217. [Google Scholar] [CrossRef]
  39. Venkatachalam, K.L.; Herbrandson, J.E.; Asirvatham, S.J. Signals and signal processing for the electrophysiologist: Part I: Electrogram acquisition. Circ. Arrhythm. Electrophysiol. 2011, 4, 965–973. [Google Scholar] [CrossRef]
  40. Venkatachalam, K.L.; Herbrandson, J.E.; Asirvatham, S.J. Signals and signal processing for the electrophysiologist: Part II: Signal processing and artifact. Circ. Arrhythm. Electrophysiol. 2011, 4, 974–981. [Google Scholar] [CrossRef]
  41. Choudhuri, S.; Ghosal, T.; Goswami, D.P.; Sengupta, A. Planarity of the spatial qrs loop of vectorcardiogram is a crucial diagnostic and prognostic parameter in acute myocardial infarction. Med. Hypotheses 2019, 130, 109251. [Google Scholar] [CrossRef]
  42. Perez-Riera, A.R.; Barbosa-Barros, R.; Daminello-Raimundo, R.; de Abreu, L.C.; de Almeida, M.C.; Rankinen, J.; Baeub Soler, F.; Nikus, K. Re-evaluating the electro-vectorcardiographic criteria for left bundle branch block. Ann. Noninvasive Electrocardiol. 2019, e12644. [Google Scholar] [CrossRef] [PubMed]
  43. Prabhakararao, E.; Dandapat, S. A Weighted SVM Based Approach for Automatic Detection of Posterior Myocardial Infarction Using VCG Signals. In Proceedings of the 2019 National Conference on Communications (NCC), Bangalore, India, 20–23 February 2019; pp. 1–6. [Google Scholar]
  44. Costa, C.M.; Silva, I.S.; de Sousa, R.D.; Hortegal, R.A.; Regis, C.D.M. The association between reconstructed phase space and Artificial Neural Networks for vectorcardiographic recognition of myocardial infarction. J. Electrocardiol. 2018, 51, 443–449. [Google Scholar] [CrossRef] [PubMed]
  45. Sadhukhan, D.; Datta, J.; Pal, S.; Mitra, M. Automated Identification of Myocardial Infarction Using a Single Vectorcardiographic Feature. In Proceedings of the International Conference on Modelling and Simulation, London, UK, 28–29 September 2017; pp. 641–651. [Google Scholar]
  46. Tripathy, R.K.; Dandapat, S. Detection of myocardial infarction from vectorcardiogram using relevance vector machine. Signal Image Video Process. 2017, 11, 1139–1146. [Google Scholar] [CrossRef]
  47. Bonomini, M.P.; Ingallina, F.J.; Barone, V.; Valentinuzzi, M.E.; Arini, P.D. Comparison of electrocardiographic and vectorcardiographic planes on a set of left ventricular hypertrophy patients. In Proceedings of the VI Latin American Congress on Biomedical Engineering CLAIB 2014, Parana, Argentina, 29–31 October 2014; pp. 564–567. [Google Scholar]
  48. Cortez, D.; Schlegel, T.T.; Ackerman, M.J.; Bos, J.M. ECG-derived spatial QRS-T angle is strongly associated with hypertrophic cardiomyopathy. J. Electrocardiol. 2017, 50, 195–202. [Google Scholar] [CrossRef] [PubMed]
  49. Nikus, K.; Perez-Riera, A.R.; Konttila, K.; Barbosa-Barros, R. Electrocardiographic recognition of right ventricular hypertrophy. J. Electrocardiol. 2018, 51, 46–49. [Google Scholar] [CrossRef] [PubMed]
  50. Mather, A.N.; Maredia, N.; Plein, S. Ischemic Heart Disease. In Cardiovascular MR Manual; Springer: Cham, Switzerland, 2015; pp. 371–410. [Google Scholar]
  51. Cruces, P.D.; Arini, P.D. Quaternion-based study of angular velocity of the cardiac vector during myocardial ischaemia. Int. J. Cardiol. 2017, 248, 57–63. [Google Scholar] [CrossRef] [PubMed]
  52. Vozda, M.; Cerny, M. Methods for derivation of orthogonal leads from 12-lead electrocardiogram: A review. Biomed. Signal Process. Control 2015, 19, 23–34. [Google Scholar] [CrossRef]
Figure 1. 1-D image of three scalar vectorcardiographic (VCG) recordings. (a) X lead, (b) Y lead, and (c) Z lead.
Figure 1. 1-D image of three scalar vectorcardiographic (VCG) recordings. (a) X lead, (b) Y lead, and (c) Z lead.
Sensors 19 03072 g001
Figure 2. 2-D image using three VCG planes. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Figure 2. 2-D image using three VCG planes. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Sensors 19 03072 g002
Figure 3. 3-D image of VCG.
Figure 3. 3-D image of VCG.
Sensors 19 03072 g003
Figure 4. Frank’s orthogonal lead system.
Figure 4. Frank’s orthogonal lead system.
Sensors 19 03072 g004
Figure 5. Box plot for comparison of the methods based on the mean squared error (MSE).
Figure 5. Box plot for comparison of the methods based on the mean squared error (MSE).
Sensors 19 03072 g005
Figure 6. Box plot for comparison of the methods based on the correlation coefficient.
Figure 6. Box plot for comparison of the methods based on the correlation coefficient.
Sensors 19 03072 g006
Figure 7. Comparison of transformation methods in 2-D, random patient 1. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Figure 7. Comparison of transformation methods in 2-D, random patient 1. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Sensors 19 03072 g007
Figure 8. Comparison of transformation methods in 3-D, random patient 1.
Figure 8. Comparison of transformation methods in 3-D, random patient 1.
Sensors 19 03072 g008
Figure 9. Comparison of transformation methods in 2-D, random patient 2. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Figure 9. Comparison of transformation methods in 2-D, random patient 2. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Sensors 19 03072 g009
Figure 10. Comparison of transformation methods in 3-D, random patient 2.
Figure 10. Comparison of transformation methods in 3-D, random patient 2.
Sensors 19 03072 g010
Figure 11. Comparison of transformation methods in 2-D, random patient 3. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Figure 11. Comparison of transformation methods in 2-D, random patient 3. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Sensors 19 03072 g011
Figure 12. Comparison of transformation methods in 3-D, random patient 3.
Figure 12. Comparison of transformation methods in 3-D, random patient 3.
Sensors 19 03072 g012
Figure 13. Comparison of transformation methods in 2-D, random patient 4. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Figure 13. Comparison of transformation methods in 2-D, random patient 4. (a) X lead to Y lead, (b) X lead to Z lead, and (c) Y lead to Z lead.
Sensors 19 03072 g013
Figure 14. Comparison of transformation methods in 3-D, random patient 4.
Figure 14. Comparison of transformation methods in 3-D, random patient 4.
Sensors 19 03072 g014
Table 1. Transformation matrix coefficients for inverse Dower transformation.
Table 1. Transformation matrix coefficients for inverse Dower transformation.
LeadXYZ
V1−0.1720.057−0.229
V2−0.074−0.019−0.310
V30.122−0.106−0.246
V40.231−0.022−0.063
V50.2390.0410.055
V60.1940.0480.108
I0.156−0.2270.022
II−0.0100.8870.102
Table 2. Coefficients for Kors regression transformation.
Table 2. Coefficients for Kors regression transformation.
LeadXYZ
V1−0.1300.060−0.430
V20.050−0.020−0.060
V3−0.010−0.050−0.140
V40.1400.060−0.200
V50.060−0.170−0.110
V60.5400.1300.310
I0.380−0.0700.110
II−0.0700.930−0.230
Table 3. Coefficients for the P wave (PLSV) method.
Table 3. Coefficients for the P wave (PLSV) method.
LeadXYZ
V10.1470.0230.184
V20.0580.0850.163
V30.0370.0030.190
V40.1390.0330.119
V50.2320.0600.023
V60.2260.1040.043
I0.1990.1460.085
II0.0180.5030.130
Table 4. Coefficients for the QRS complex (QLSV) method.
Table 4. Coefficients for the QRS complex (QLSV) method.
LeadXYZ
V1−0.2660.088−0.319
V20.027−0.088−0.198
V30.0650.003−0.167
V40.1310.042−0.099
V50.2030.047−0.009
V60.220.0670.06
I0.37−0.1310.184
II−0.1540.717−0.114
Table 5. Table of the medians of the MSE results for individual transformation methods and individual leads and the resulting p-values between the medians of Kors regression transformation and the medians of other transformation methods for individual leads.
Table 5. Table of the medians of the MSE results for individual transformation methods and individual leads and the resulting p-values between the medians of Kors regression transformation and the medians of other transformation methods for individual leads.
Mean Squared ErrorXYZ
Regress Median2.2024 · 10 4 3.8445 · 10 4 1.9622 · 10 3
Quasi Median1.5799 · 10 3 1.9417 · 10 3 3.8956 · 10 3
Dower Median1.4880 · 10 3 8.7119 · 10 4 6.5558 · 10 3
PLSV Median1.1881 · 10 3 7.2627 · 10 4 2.3116 · 10 3
QLSV Median1.1887 · 10 3 1.0313 · 10 3 2.5286 · 10 3
p-value (Regress × Quasi)2.4982 · 10 11 3.5755 · 10 12 1.9268 · 10 4
p-value (Regress × Dower)1.2882 · 10 9 5.7956 · 10 5 1.1888 · 10 7
p-value (Regress × PLSV)2.7719 · 10 8 3.6261 · 10 4 0.2878
p-value (Regress × QLSV)4.9807 · 10 9 3.2827 · 10 7 0.3576
Table 6. Table of the medians of the R results for individual transformation methods and individual leads and the resulting p-valuess between the medians of Kors regression transformation and the medians of other transformation methods for the individual leads.
Table 6. Table of the medians of the R results for individual transformation methods and individual leads and the resulting p-valuess between the medians of Kors regression transformation and the medians of other transformation methods for the individual leads.
Correlation CoefficientXYZ
Regress Median0.9983970.9943430.981839
Quasi Median0.9867840.9862230.948075
Dower Median0.9903450.9828410.944712
PLSV Median0.9903500.9820360.966346
QLSV Median0.9946950.9767520.971982
p-values (Regress × Quasi)2.5482 · 10 9 3.2827 · 10 4 3.0635 · 10 3
p-values (Regress × Dower)6.0021 · 10 8 2.4887 · 10 4 3.6261 · 10 4
p-values (Regress × PLSV)2.1785 · 10 6 9.7206 · 10 5 0.2485
p-values (Regress × QLSV)7.1535 · 10 6 1.4784 · 10 6 0.1612

Share and Cite

MDPI and ACS Style

Jaros, R.; Martinek, R.; Danys, L. Comparison of Different Electrocardiography with Vectorcardiography Transformations. Sensors 2019, 19, 3072. https://doi.org/10.3390/s19143072

AMA Style

Jaros R, Martinek R, Danys L. Comparison of Different Electrocardiography with Vectorcardiography Transformations. Sensors. 2019; 19(14):3072. https://doi.org/10.3390/s19143072

Chicago/Turabian Style

Jaros, Rene, Radek Martinek, and Lukas Danys. 2019. "Comparison of Different Electrocardiography with Vectorcardiography Transformations" Sensors 19, no. 14: 3072. https://doi.org/10.3390/s19143072

APA Style

Jaros, R., Martinek, R., & Danys, L. (2019). Comparison of Different Electrocardiography with Vectorcardiography Transformations. Sensors, 19(14), 3072. https://doi.org/10.3390/s19143072

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop