 Research
 Open Access
 Published:
Spatiotemporal representation of cardiac vectorcardiogram (VCG) signals
BioMedical Engineering OnLine volume 11, Article number: 16 (2012)
Abstract
Background
Vectorcardiogram (VCG) signals monitor both spatial and temporal cardiac electrical activities along three orthogonal planes of the body. However, the absence of spatiotemporal resolution in conventional VCG representations is a major impediment for medical interpretation and clinical usage of VCG. This is especially so because timedomain features of 12lead ECG, instead of both spatial and temporal characteristics of VCG, are widely used for the automatic assessment of cardiac pathological patterns.
Materials and methods
We present a novel representation approach that captures critical spatiotemporal heart dynamics by displaying the real time motion of VCG cardiac vectors in a 3D space. Such a dynamic display can also be realized with only one lead ECG signal (e.g., ambulatory ECG) through an alternative lagreconstructed ECG representation from nonlinear dynamics principles. Furthermore, the trajectories are color coded with additional dynamical properties of spacetime VCG signals, e.g., the curvature, speed, octant and phase angles to enhance the information visibility.
Results
In this investigation, spatiotemporal VCG signal representation is used to characterize various spatiotemporal pathological patterns for healthy control (HC), myocardial infarction (MI), atrial fibrillation (AF) and bundle branch block (BBB). The proposed color coding scheme revealed that the spatial locations of the peak of T waves are in the Octant 6 for the majority (i.e., 74 out of 80) of healthy recordings in the PhysioNet PTB database. In contrast, the peak of T waves from 31.79% (117/368) of MI subjects are found to remain in Octant 6 and the rest (68.21%) spread over all other octants. The spatiotemporal VCG signal representation is shown to capture the same important heart characteristics as the 12lead ECG plots and more.
Conclusions
Spatiotemporal VCG signal representation is shown to facilitate the characterization of spacetime cardiac pathological patterns and enhance the automatic assessment of cardiovascular diseases.
Background
The electrocardiogram (ECG) signals are recorded on the body surface to track the continuous dynamic details of cardiac functioning. Such valuable realtime information is usually unavailable in static and discrete clinical laboratory tests, for e.g., computer imaging, chest xray, and blood enzyme test. Even if routine laboratory examinations are performed multiple times per day, discontinuity often fails to prevent the lethal consequences from acute cardiac disorders. There is an increasing awareness that realtime ECG monitoring is an essential tool for the early identification of cardiac pathological patterns because it tracks cardiac dynamic behaviors as opposed to static screenshots.
However, one lead ECG signals only capture one perspective temporal view of the spacetime excitation and propagation of cardiac electrical activities. Multiple lead ECG systems, for e.g., 12lead ECG and 3lead vectorcardiogram (VCG), are designed to capture the multidirectional view of spacetime cardiac electrical activities. The 12lead ECG is more commonly used than the 3lead VCG because medical doctors are accustomed to using them in clinical applications. It has thus proven its value, time tested, and considered as the Gold Standard. However, much of that information is redundant and even in that, only a small fraction of the data is used in the analysis by the physicians based on experience, expertise, and oftentimes on the memorization of ECG signals for different cardiac disorders. This is a difficult task and cardiologists are constantly looking for more accurate and effective alternatives. In addition, when it comes to the automated computer analysis of heart electrical activity, 12lead ECG signals have higher dimensionality than 3lead VCG and can potentially introduce the “curse of dimensionality” problem. VCG surmounts not only the information loss from only one or two ECG signals but also the dimensionality problems induced by the 12lead ECG signals. In this paper, we propose a new dynamic display of VCG to characterize both spatial and temporal patterns of cardiac electrical activities as an alternative for the conventional 12lead ECG.
VCG signals monitor the cardiac electrical activities along three orthogonal planes of the body, namely, frontal, transverse, and sagittal. As shown in Figure 1, the VCG vector loops contain 3 dimensional recurring, nearperiodic patterns of heart dynamics [1]. Each heart cycle consists of three loops corresponding to P, QRS, and T wave activities. The VCG signals are traditionally projected onto different planes (XY plane, XZ plane and YZ plane) to capture the time correlations, or plotted as a static attractor in a 3D space that provides the topological relationships. The absence of combined spatial and temporal information in the VCG representations reported earlier expects the interpreters to have knowledge on both spatial and temporal patterns of cardiac events.
With rapid advancements in information technology and the availability of computing hardware at reasonable costs, representation of 3D VCG loops are not constrained by computational resources anymore and this has generated renewed interest in VCG since the 1990’s. Dower and his colleagues [2–4] conducted pioneering research based on Frank’s tank torso model studies and introduced a linear transformation matrix to derive 12lead ECG signals from 3lead VCG signals. Such a transform was shown to preserve clinically useful information pertinent to heart dynamics. It may be noted that the 3 orthogonal leads, e.g. the corrected Frank leads, contain the all necessary information but very often the signal amplitudes are low and superimposed by noise. In such cases the redundancy is extremely useful. Besides, the orthogonal leads are acquired via uncomfortable electrode positions (i.e., on the patient back). Therefore, the inverse transform is often used to derive the 3lead VCG from the 12lead ECG [5]. Some even consider Dower transformation matrix as “generalized” or “universal” transformation matrix [2–4, 6] although there is a need to develop different transformation matrices for healthy subjects and patients with cardiac disorders [7].
For a given cluster of subjects, e.g., healthy control subjects in certain age and gender group, the transformation studies show statistical equivalence between the 12lead ECG and 3lead VCG signals [2–4, 6, 7]. However, the spatial information of cardiac vectors is hardly recognized in the conventional 12lead electrocardiogram (ECG) plots, and temporal aspects are less discernible in the traditional visualization of 3D static VCG vector trajectory. In addition, it is difficult for human beings to visually project a spatial VCG vector into any specified angle in the 12lead measurement system, which is the traditional way for the interpretation of ECG signals. It may be noted that the 12lead system present the necessary information in more than 90% of the cases. The VCG leads are useful specifically for deriving details that are pertinent to myocardial infarctions. Therefore, the physicians very often interpret in parallel the 12 leads and the VCG loops.
It is well known that the 12lead ECG is heavily depended on temporal information, for e.g., intervals and durations, but such temporal aspects are absent in the static VCG representation. Hence, this poses extra difficulties to closely relate the 12lead ECG characteristics (e.g., QT interval, ST elevation) to the patterns of 3lead Frank X, Y, Z VCG signals. Although the medical interpretation and clinical usage of VCG have been investigated by many researchers [8–13], the dynamic representation of both spatial and temporal aspects in VCG signals are rarely, if any, to be found. Most previous investigations studied the pathological patterns in the static VCG signals and exemplified the advantages of VCG in the automatic assessment of certain cardiovascular disorders. Spatial VCG signals were shown to not only facilitate in the basic understanding of the electrical phenomena associated with the heart but also disclose pathological characteristics unknown, or not feasible, from the ECG signals [14, 15]. This was experimentally validated for different cardiovascular disorders including the right and left ventricular hypertrophy, singly or combined, intraventricular block, and myocardial infarction [14, 15]. This paper presents a dynamic VCG signal representation approach to capture critical spatiotemporal heart dynamics by displaying the real time motion of 3dimensional VCG cardiac vectors on a computer screen.
Materials and methods
In this investigation, we have used 3lead VCG signals gathered from the PTB database available in the PhysioNet [16, 17]. Each of the recordings in the PTB database contains 15 simultaneously recorded signals, namely, the conventional 12lead ECGs and the three orthogonal Frank XYZ lead VCG signals digitized at 1 kHz, with a 16 bit resolution over a range of ±16.384 mV. The database used consists of 80 HC recordings acquired from 54 healthy volunteers and 368 MI recordings from 148 patients. The recordings were collected at the Department of Cardiology of University Clinic Benjamin Franklin in Berlin, Germany.
Description of the method
As shown in Figure 2, the VCG vector loops contain 3 dimensional recurring, nearperiodic patterns of heart dynamics, which can be visualized in the X, Y and Z space domain with time entering implicitly. The dynamic VCG signal representation provides an easier way to understand, interpret, and use spacetime information of cardiac electrical activities. This approach includes presenting the real time motion of cardiac vectors in 3D space and color coding of cardiac vector movements with some dynamical properties, for e.g., curvature, spatial octant [18], velocity, and phase angle. For one dimensional heart monitoring signals, an alternative lagreconstructed representation from nonlinear dynamics principles is provided. In addition, the Poincaré sectioning of the 3D VCG vector loops extracts the homogeneous ECG ensembles for selective averaging. The cycletocycle selfsimilarities and variations provide valuable information about the majority clustering of ECG morphology and heart rate variability. The proposed dynamic representation of cardiac vectors can be rotated freely on the computer monitor screen for the spatiotemporal analysis.

A.
Spatiotemporal VCG signal representationIn the Frank XYZ lead system, VCG is represented as three orthogonal scalar measurements with respect to time as given in Eq.1. The dynamic VCG signal representation embeds the cardiac vector, composed of three scalar measurements, in real time. As shown in Figure 2, three scalar x, y, z components are plotted in the top and the simultaneous 3D movement of cardiac vectors in the bottom.
Additional file 1: Animation video for the spatiotemporal VCG representation of a patient with myocardial infarction, anterior. The top plot displays VCG signals in three vector components as a function of time, and the bottom part shows the real time cardiac vector movement in the 3D space. Head (green) gives the current position of cardiac vector. Tail (red) indicates the direction and rotation of cardiac vector movements. (AVI 1 MB)
Therefore, this explicitly real time spatiotemporal VCG representation makes it easier to integrate with prior knowledge and experiences of timebased ECG. As shown in Figure 2, this representation consists of three components, namely, head (green), body (red), and tail (blue). Head gives the current position of the cardiac vector. Body records a short history of the cardiac vector movements which clearly indicates where the current vector is from. It avoids the confusion regarding which group of heart activity the current cardiac vector belongs to as they usually intersect at the isoelectric points. The tail provides full history pertinent to the complete topological shape of VCG state space. By following the cardiac vector movement with respect to time, the P, QRS, and T waves will be easily located in the VCG state space (see Figure 2).

B.
Lagreconstructed spatiotemporal ECG representation
In some realworld cases, for e.g., ambulatory ECG, only a single channel ECG is available in lieu of the complete measurements of three cardiac vector components. Takens embedding theorem [19] states that the individual measurements contain adequate information to reconstruct pseudo state space because of the high couplings in realworld complex systems. Although a single channel ECG is hardly presented in the spacetime form as cardiac vectors, a pseudo state space (attractor) can be reconstructed from delayed coordinates of the single measurement y(t) as
$$\underset{\u035e}{v}\left({t}_{i}\right)=\left[y\left({\text{t}}_{i}\right),y\left({\text{t}}_{i}+\tau \right)\text{,}\phantom{\rule{0.25em}{0ex}}y\left({\text{t}}_{i}+2\tau \right)\right]$$(2)
where τ is the time delay. The optimal time delay τ is selected to minimize mutual information function M(τ), defined as
where p(t,τ) is the joint density function, and p(t) and p(t + τ) are marginal density functions of y(t) and y(t + τ), respectively [19]. As shown in Figure 3, the animation of lagreconstructed ECG attractor will provide spacetime information when the multidimensional ECG signals are not readily available. It may be noted that the lagreconstruction approach was widely used in physics to explore the nonlinear dynamics underlying the ECG signals [20–22]. Examples of the nonlinear dynamical quantifiers may include Lyapunov exponents, recurrence statistics, and correlation dimension, which describes the characteristics of the signals. This paper presents the lag reconstruction approach to facilitate the realtime spatiotemporal representation when 3lead VCG, i.e., actual state space, is not available.

C.
Color coding of spatiotemporal VCG signal representationThe hardware apparatus for recording “color vectorcardiogram” was previously designed using camera, filter, lens and oscilloscope to facilitate the understanding and increase the diagnostic scope of vectorcardiography [23]. This paper presents a software colorcoding scheme to incorporate additional dynamical attributes of spatiotemporal VCG signals as opposed to use hardware devices. It will not only overcome the color resolution drawbacks from the limitations of hardware but also reduce the economic costs of VCG monitoring systems. The dynamical properties for the coloring of VCG vector loops may include, but not limited to, speed (v′(t) = Δv / Δt, and Δv =  v(t) v(t+ Δt)), phase angles, octant numbers, and curvature. The phase angles provide similar information as octant numbers, but in finegrained scales. The phase angles of cardiac vectors can be determined as the following equation:
The eight octants delimit the cardiac vectors in the scale of 90^{0}. The VCG trajectory in the 3D space is color coded with eight different colors with respect to the octant numbers. As can be noted from Table 1, a binary number with three bits (binary digits) is used to designate each octant. Binary code 0 is used to represent the negative directions in X, Y, and Z axes and 1 for the positive directions [18]. For example, if the octant lies in the (−, +, −) XYZ directions, binary coding for this will be (010) and the resulting octant number is 2 (binary coding: 0 × 2^{2} + 1 × 2^{1} + 0 × 2^{0} = 2).If the VCG is denoted as three orthogonal scalar measurements with respect to time as $v\left(t\right)=<{v}_{x},{v}_{y},{v}_{z}>=<f\left(t\right),g\left(t\right),h\left(t\right)>,a\le t\le b$ and the cardiac trajectory is traversed as t increases from a to b, then it can be shown that the length of such a space curve is $L={\int}_{a}^{b}\left{v}^{\prime}\left(t\right)\rightdt={\int}_{a}^{b}\sqrt{{\left[{f}^{\prime}\left(t\right)\right]}^{2}+{\left[{g}^{\prime}\left(t\right)\right]}^{2}+{\left[{h}^{\prime}\left(t\right)\right]}^{2}}dt$. The arc length s as a function of t is $s\left(t\right)={\int}_{a}^{t}\left{v}^{\prime}\left(u\right)\rightdu={\int}_{a}^{t}\sqrt{{\left[{f}^{\prime}\left(u\right)\right]}^{2}+{\left[{g}^{\prime}\left(u\right)\right]}^{2}+{\left[{h}^{\prime}\left(u\right)\right]}^{2}}du$. The unit tangent vector giving the speed direction at a particular point is calculated by the following formula: $T\left(t\right)={v}^{\prime}\left(t\right)/\left{v}^{\prime}\left(t\right)\right$. The curvature, at a given point, of the cardiac VCG trajectory is defined as the rate of change of the unit tangent vector with respect to arc length. In other words, the curvature at a point indicates how fast the trajectory is bending at that point. Therefore, the curvature $\kappa \left(t\right)$ of 3lead VCG trajectory is calculated as follows:
The magnitudes of VCG dynamical properties (speed, phase angle, curvature etc.) are mapped into a color scale. Thus, the color coded VCG signals can be plotted in real time to provide extra information besides X, Y, Z, and time scale. As shown in Figure 4, the color coded spatiotemporal VCG representation incorporates extra cardiac vector movement information so as to facilitate doctor’s interpretation of valuable spatiotemporal patterns. Figure 4 displays the speed of each vector movement in a specific color corresponding to the color bar included in the right side of the figure. As illustrated by the color variations in Figure 4 for a patient with the bundle branch block disorder, cardiac vectors in QRS loop move the fastest, T wave the second, P wave the third, and isoelectric points the slowest. It may also be noted that there are two curly R peaks (see Figure 4) which are designated as “M waves”, i.e., one of the typical pathological patterns for bundle branch block patients. The “M wave” pattern shows that both ventricles are not depolarized simultaneously. The delay in the blocked bundle branch allows the unblocked ventricle to begin depolarizing before the blocked ventricle. This kind of slightly lag effect in one ventricle produces pathological ‘M’ wave and widens the QRS loop.

D.
ECG Ensemble CharacterizationThe 3D display of cardiac vector loops from several cycles commonly show the nearperiodic patterns but with hidden temporal variations between heart cycles. The aligned ECG ensembles will contribute to the gathering of homogeneous majority clusters of ECG signals for selective averaging. The Poincaré section (see Figure 5 (a)) is utilized to characterize the spatiotemporal cardiac patterns in the form of ensembles as well as detailed beattobeat variations [18, 24]. Here, Poincaré section is a 2 dimensional hyperplane intersecting with the state space trajectories. The recurrence property of VCG trajectory shows that for every ϵ < 0 and almost every cardiac vector $v\left(i\right),\exists j>0$ such that $v\left(i\right)v\left(j\right)<\u03f5$, in effect, the trajectories of VCG vectors remain bounded. Those points at which the trajectory intersects the Poincaré section follow a return map. Figure 5 (bd) shows the aligned heart beats along the X, Y, Z axes. It may be noted that heart cycles are sharing similar morphologies in any of three orthogonal directions but there exist beattobeat variations due to heart rate variability. The time elapsed for the completion of one VCG cycle (P, QRS, T) provides the RR intervals [24]. Heart rate variability will drive some ensembles move faster, i.e., some ECG strands have shorter intervals, compared to the others.
Results
This present investigation is aimed at a new dynamic display of VCG to characterize both spatial and temporal patterns of cardiac electrical activities. The temporal information in the spatiotemporal VCG representation is essential for the identification of cardiac arrhythmia (abnormal heart rhythm such as Bradycardia and Tachycardia). In the 12lead ECG, Bradycardia is identified as a resting heart rate of <60 beats per minute in the time domain. Similarly, cardiac vectors from Bradycardia patients will rotate for less than 60 cardiac cycles per minute in the dynamic VCG signal representation. The implementation and advantages of dynamic VCG representation are shown as follows:

A.
Dynamic VCG Implementation DetailsMost previous investigations utilized the lagreconstructed ECG representation to quantify nonlinear dynamical patterns, e.g., recurrence statistics and study their correlations with cardiac disorders. The mathematical formulation of recurrence quantifiers is detailed in our previous investigation [25, 26] and other references [20, 27]. The recurrence plot is defined as: ${T}_{i,j}:=\Theta \left(\u03f5\overrightarrow{x}\left(i\right)\overrightarrow{x}\left(j\right)\right)$, where ϵ is a cutoff distance and Θ is the Heaviside function. As shown in Figure 6, the recurrence plot of 3lead VCG (see Figure 6 a) yields similar recurrence patterns as the lagreconstructed ECG state space (see Figure 6 b) for the recording of patient104/s0306lre. The black dot indicates that the distance between the states $\overrightarrow{x}\left(j\right)$and $\overrightarrow{x}\left(j\right)$ is below the cutoff distance ε. The texture patterns in the recurrence plots reveal information of the ECG signal, e.g., the diagonal structures represent the nearperiodic patterns. Furthermore, six quantitative features are extracted from recurrence plots to analyze the underlying processes and hidden cardiac rhythms. The recurrence quantifiers include recurrence rate (RR), determinism (DET), maximal length of diagonal structures (LMAX), entropy (ENT), laminarity (LAM) and trapping time (TT) [25–27].The values of recurrence quantifiers are RR:94.63, DET:99.54, LMAX:404, ENT:6.55, LAM:99.64, TT:55.87 for 3lead VCG (see Figure 6 a), and RR:99.01, DET:97.44, LMAX:410, ENT:5.65, LAM:98.40, TT:37.69 for the lagreconstructed ECG state space (see Figure 6 b). In addition, we extracted the recurrence quantification statistics for 80 healthy controls and calculated the relative errors. The relative error is calculated as $\left{x}_{\mathit{true}}{x}_{\mathit{pseudo}}\right/{x}_{\mathit{true}})$, where ${x}_{\mathit{true}}$ is the quantifier computed from the VCG state space, and ${x}_{\mathit{pseudo}}$ is from the lagreconstructed pseudo state space. The relative errors provide a good measure of how good the recurrence statistics from the lagreconstructed ECG state space is relative to the 3lead VCG. As shown in Figure 6 (c), the box plot is used to visualize the distribution of relative errors. The red line in the middle of boxplot represents the median, the blue box shows the lower quartile and upper quartile of data distributions, and the black dash lines represent the most extreme values within 1.5 times the interquartile range. Figure 6 (c) shows that the relative errors of DET, LMAX, and LAM are less than 6.6% for all 80 subjects. The relative errors of RR are less than 18.9% for half of the subjects, and the median relative errors of ENT are less than 13.4%. However, the median relative errors of TT are around 43.2% and higher than the other five recurrence quantifiers. It may be noted that the lagreconstructed ECG state space approximates the measures of some dynamical properties with small relative errors when the 3lead VCG is not readily available.Figure 7 shows the colorcoded spatiotemporal VCG representation for a healthy control subject. It may be noted that the realtime movements of cardiac vectors can be examined using animations of the dynamic VCG signals (view Additional file 4 for animations enclosed separately). As shown in Figure 7, the peak of T wave loop is in the Octant 6 (XYZ: ++−) with the yellow color. This is true for most of the healthy control cases. Statistical analysis showed that only 6 out of 80 healthy recordings in PTB database are away from Octant 6. In contrast, the peak of T waves from 31.79% (117/368) of MI subjects are found to remain in Octant 6 and the rest (68.21%) spread over all other octants. This indicates that the spatial directions of ventricle repolarization are deviated from octant 6 for the majority (i.e., 68.21%) of MI subjects. Although this is analogous to the T wave inversion patterns in the time domain, it may be noted that a variety of T wave inversions is shown in different leads of the 12lead ECG. In the spatiotemporal domain, this present investigation revealed that the spatial directions of ventricle repolarization belong to the octant 6 for the majority (i.e., 74/80) of health controls.Therefore, the spatiotemporal VCG representation makes irregular ventricle repolarization behaviors more clearcut and easier for computer implementation. In addition, the initial portion of QRS loops is ‘Q wave’ when the initial vectors are in opposite direction from the maximal R vector (i.e., scalar product negative), and similarly for terminal vectors as ‘S wave’. A different color scheme can also be added by identifying within the initial 3D QRS loop with colored dots or other suitable labels 20 ms and 30 ms time points to add temporal information about Q wave duration. The Q wave and S wave can be located by combining the temporal XYZ display in the top and the real time motion of cardiac vectors in the bottom of spatiotemporal VCG representation.

B.
Benefits of Dynamic VCG representationThe dynamic VCG representation further enhances our previous investigations [18, 25] to characterize the spatiotemporal VCG patterns and automatically assess cardiovascular conditions. We showed that recurrence quantifiers extracted from 3lead VCG are good indicators of myocardial infarction [18, 25]. Dynamic VCG representation possesses advantages to disclose spacetime pathological characteristics unknown, or not feasible, from the timedomain ECG signals. Figure 8 shows color coded spatiotemporal VCG representation for a dysrhythmia and atrial fibrillation patient. Atrial fibrillation [28], due to continuous rapidfiring of many atrial automaticity foci, typically produces tiny and wavy ECG spikes instead of identifiable P waves. In the dynamic VCG representation, P loops is not clearly distinguishable and the VCG trajectory appears to be chaotic before entering the QRS loops. Such erratic cardiac vector movements during atrial electrical activity closely correlate with the rapidfiring of multiple irritable atrial foci. Multidimensional spatial view provides a complete picture of abnormal atrial electrical activities and facilitates the computer based multivariate analysis. If the spatiotemporal cardiac electrical activity is projected along the direction of V1V2 leads, it is easier to visualize the prominent fibrillation waves in the time domain. However, the projection of spatiotemporal cardiac electrical activity diminishes important multidimensional information of cardiac pathological behaviors into one dimension. As such, there is an information loss in spite of the fact of easy visualization. In other words, if a spatial vector is projected into an axis that does not contain this vector, there will be information loss.It may also be noted in Figure 8 that the RR interval is longer than the normal case mainly because of the long time taken for the excitation of atrium (isopotential duration). In addition, the color changes of the VCG trajectory with respect to octant numbers are vastly different from the normal case shown in Figure 7. The disappearance of red and the appearance of green color instead in the QRS loop also indicate the abnormal locations of cardiac electrical axis, which are pathological patterns for some diseases. The cardiac electrical axis refers to the depolarization direction of ventricular myocardium, which is conventionally defined based on the QRS complexes of leads I, II, III, AVL, AVF, AVR in the frontal plane and V1V6 in the horizontal plane. This is similar to determine the direction of ventricular repolarization in 3D space. Therefore, the distribution of the QRS loop along various octants, i.e., colors of the QRS loop, provides an analogous measure of the heart electrical axis.As shown in Figures 4, 7 and 8, pathological patterns of cardiac disorders were characterized from the spatiotemporal VCG representation. Spatiotemporal VCG signal representation enhances the automatic assessment of cardiovascular diseases. The software colorcoding scheme facilitates in the quantification of spacetime cardiac pathological patterns.
Conclusions
In this investigation, we have made an attempt to capture the spatiotemporal characteristics of the VCG signals by viewing the cardiac vectors in real time on a computer monitor screen instead of a static signal output on a paper which at best can record 3 dimensions. This approach overcomes the drawbacks of conventional static VCG representation and provides concurrent spatial and temporal resolutions. The alternative lagreconstructed approach from nonlinear dynamics principles addresses the sometime difficult situations for the study of cardiac state space when there is only 1lead ECG signal available. The spatiotemporal VCG representation incorporates additional dynamical properties of cardiac vector movements (curvature, velocity, octant, and phase angle etc.) with the color coding scheme. Furthermore, Poincaré sectioning of the 3D VCG vector loops extracts the homogeneous ECG ensembles for the information on cardiac cycletocycle selfsimilarities and variations. It is shown that the proposed dynamic VCG approach surmounts some drawbacks of both 12lead ECG and static VCG representation, and provides critical spatial as well as temporal information of the heart dynamics. The cardiovascular pathological patterns are found to be effectively captured by this new 3D dynamic representation approach. The presence of both spatial and temporal characteristics in dynamic representation improves the automatic assessment of cardiovascular diseases with the use of VCG signals.
Dedication
This paper is respectfully dedicated to the memory of Dr. Ranga Komanduri (19422011) who was the originator of this manuscript.
References
 1.
Malmivuo J, Plonsey R: Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields. Oxford University Press, USA; 1995.
 2.
Dower GE, Yakush A, Nazzal SB, Jutzy RV, Ruiz CE: Deriving the 12lead electrocardiogram from four (EASI) electrodes. J Electrocardiol 1988, 21(Supplement 1):S182S187.
 3.
Dower GE, Machado HB, Osborne JA: On deriving the electrocardiogram from vectorcardiographic leads. Clin Cardiol 1980, 3(2):87–95.
 4.
Dower GE, Machado HB: XYZ data interpreted by a 12lead computer program using the derived electrocardiogram. J Electrocardiol 1979, 12(3):249–261. 10.1016/S00220736(79)800588
 5.
Levkov CL: Orthogonal electrocardiogram derived from the limb and chest electrodes of the conventional 12lead system. Med Biol Eng Comput 1987, 25(2):155–164. 10.1007/BF02442844
 6.
Hubert VP: Advantages of three lead cardiographic recordings. Ann N Y Acad Sci 1965, 126: 873–881. Computation for Cardiovascular Research
 7.
Dawson D, Yang H, Malshe M, Bukkapatnam STS, Benjamin B, Komanduri R: Linear affine transformations between 3lead (Frank XYZ leads) vectorcardiogram (VCG) and 12lead electrocardiogram (ECG) signals. J Electrocardiol 2009, 42(6):622–630. 10.1016/j.jelectrocard.2009.05.007
 8.
Benchimol A, Reich F, Desser KB: Comparison of the electrocardiogram and vectorcardiogram for the diagnosis of left atrial enlargement. J Electrocardiol 1976, 9(3):215–218. 10.1016/S00220736(76)800489
 9.
Ng J, Sahakian A, Swiryn S: Vector analysis of atrial activity from surface ECGs recorded during atrial fibrillation. Comput Cardiol 2002, 29: 21–24.
 10.
Olson C, Wamer RA: The quantitative 3diemnsional vectorcardiogram. J Electrocardiol 2000, 33(Supplement 1):176S74.
 11.
Talwar KK, Radhakrishnan S, Hariharan V, Bhatia ML: Spatial vectorcardiogram in acute inferior wall myocardial infarction: its utility in identification of patients prone to complete heart block. Int J Cardiol 1989, 24(3):289–292. 10.1016/01675273(89)900065
 12.
van Oosterom A, Ihara Z, Jacquemet V, Hoekema R: Vectorcardiographic lead systems for the characterization of atrial fibrillation. J Electrocardiol 2006, 40(4):343.e1–343.e11.
 13.
Willems JL, Lesaffre E, Pardaens J: Comparison of the classification ability of the electrocardiogram and vectorcardiogram. Am J Cardiol 1987, 59(1):119–124. 10.1016/S00029149(87)800838
 14.
Wolff L: The vectorcardiographic diagnosis of myocardial infarction. Chest 1955, 27(3):263–281.
 15.
Wolff L, Richman J, Soffe A: Spatial vectorcardiography, review and critique. N Engl J Med 1953, 248(19):810–819. 10.1056/NEJM195305072481905
 16.
Proceedings of the Computers in Cardiology, 2006 2006.
 17.
Goldberger A, Amaral L, Glass L, Hausdorff J, Ivanov P, Mark R, Mietus J, Moody G, Peng C, Stanley H: PhysioBank, PhysioToolkit, and PhysioNet: components of a new research resource for complex physiologic signals. Circulation 2000, 23: e215e220.
 18.
Yang H, Bukkapatnam STS, Le Te, Komanduri R: Identification of myocardial infarction (MI) using spatiotemporal heart dynamics. Med Eng Phys 2002. Available online 20 September 2011, ISSN 1350–4533 10.1016/j.medengphy.2011.08.009
 19.
Kantz H, Schreiber T: Nonlinear Time Series Analysis. Cambridge University Press, Cambridge; 1997.
 20.
Zbilut JP, Thomasson N, Webber CL: Recurrence quantification analysis as a tool for nonlinear exploration of nonstationary cardiac signals. Med Eng Phys 2002, 24(1):53–60. 10.1016/S13504533(01)001126
 21.
Stefanovska A, Bracic M: Physics of the human cardiovascular system. Contemp Phys 1999, 40(1):31–55. 10.1080/001075199181693
 22.
Stanley HE, Amaral LA, Goldberger AL, Havlin S, Ivanov PC, Peng CK: Statistical physics and physiology: monofractal and multifractal approaches. Physica A 1999, 270(1–2):309–324.
 23.
Kinoshita S, Kobayashi T, Takahashi K: Color vectorcardiograms recorded in two unusual cases. Chest 1974, 65(6):673–674. 10.1378/chest.65.6.673
 24.
Yang H, Bukkapatnam STS, Komanduri R: Nonlinear adaptive wavelet analysis of electrocardiogram signals. Phys Rev E 2007, 76(2):026214.
 25.
Yang H: Multiscale recurrence quantification analysis of spatial cardiac vectorcardiogram (VCG) signals. Biomed Eng, IEEE Trans 2011, 58(2):339–347.
 26.
Yang H, Bukkapatnam STS, Barajas LG: Local recurrence based performance prediction and prognostics in the nonlinear and nonstationary systems. Pattern Recognit 2011, 44(8):1834–1840. 10.1016/j.patcog.2011.01.010
 27.
Marwan N, Carmen Romano M, Thiel M, Kurths J: Recurrence plots for the analysis of complex systems. Physics Reports 2007, 438(5–6):237–329.
 28.
Bukkapatnam STS, Komanduri R, Yang H, Rao P, Lih W, Malshe M, Raff LM, Benjamin B, Rockley M: Classification of atrial fibrillation episodes from sparse electrocardiogram data. J Electrocardiol 2008, 41(4):292–299. 10.1016/j.jelectrocard.2008.01.004
Acknowledgments
This project is funded by grants from the USF new faculty startup fund, National Science Foundation (NSF) (CMMI1000978, CMMI0700680), and the generous support from OSU Tulsa Center for Health Sciences. The authors also thank, A.H. Nelson, Jr. Endowed Chair in Engineering (RK), and AT&T Endowed Professorship (SB) for additional financial support.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contribution
HY contributed to the development of spatiotemporal methods and tools, evaluated the data, performed the data analysis, and drafted the manuscript as part of his PhD work under the advisement of the other two coauthors (SB and RK). SB contributed to the design of the study and revised the manuscript. RK conceived the study and revised the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Animation video for Lagreconstructed dynamic ECG representation for a normal subject.
Additional file 2: Time delay coordinates of the original one dimensional ECG time series are used to embed the 3D manifold. (AVI 1 MB)
Animation video for dynamic VCG representation of a patient with bundle branch block.
Additional file 3: The 3D VCG trajectory is color coded with the potential variation rate listed in the right side. Two curly R peaks designated as “M waves” are shown in the QRS loop. The ‘M’ wave pattern shows that both ventricles are not depolarized simultaneously. (AVI 730 KB)
Animation video for color coded dynamic VCG representation of a healthy control (HC) subject. T wave is in yellow color and at the normal position – Octant 6.
Additional file 4: The appearance of red color in QRS loop indicates the usual position of cardiac electrical axis. (AVI 1 MB)
Animation video for Color coded dynamic VCG representation of a patient with dysrhythmia and atrial fibrillation.
Additional file 5: The irregular P wave morphology shows abnormal excitations from SA node to Atrial. The disappearance of red and appearance of green color in QRS loop also indicate abnormal locations of cardiac electrical axis. (AVI 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Yang, H., Bukkapatnam, S.T. & Komanduri, R. Spatiotemporal representation of cardiac vectorcardiogram (VCG) signals. BioMed Eng OnLine 11, 16 (2012). https://doi.org/10.1186/1475925X1116
Received:
Accepted:
Published:
Keywords
 Vectorcardiogram (VCG)
 Electrocardiogram (ECG)
 Spatiotemporal representation
 Colorcoding scheme