Keywords: Magnetocardiography, MCG, QRS complex, selective averaging, spline interpolation 

Introduction

In a previous paper [11] we discussed QRS amplitude and shape variability in magnetocardiograms which were mainly caused by respiration. In this abstract we review these results in regard to selective averaging of QRS complexes.

Methods

The magnetocardiagram of a healthy volunteer was measured in a magnetically shielded room (Vacuumschmelze Hanau, AK3b) at the Biomagnetic Center Jena using the 62 channel biomagnetometer (Philips, Hamburg) [12]. The healthy state of this volunteer was proven by a history free of cardiac symptoms, normal physical examination, normal 12-lead electrocardiogram, and normal findings in echocardiography. The subject was lying in a supine position and the two dewars were positioned above the thorax so that they covered the field maxima. The subject was instructed to breathe normal but very uniform. We selected the channel providing the highest signal amplitude at R peak. The position of the pick-up coil of this channel was approximately 12 mm right and 170 mm below the jugulum. We recorded  a time period of 100 s at a sampling rate of 1000 Hz. An analog first order highpass filter (0.036 Hz, 3 dB) was applied to the analog signals.  An example of two QRS complexes measured is given in Figure 1.
 

Figure 1. Example of a measured MCG-signal.
 

Modeling and Computation

A scalable QRS function SQRS is used to model the QRS complexes:
 

      S_QRS(t) = A ^. N_QRS (L  ^.  t - t_beat) + S₀ + S₁^.(t - t_beat)
 

The function N_QRS (t) is the normalized QRS signal. We use the parameters  A (amplitude), t_beat (event time), and (S₀ , S₁) (describing the linear baseline) for every QRS complex. The parameter L describes changes in the length of QRS complex (shortened and prolonged QRS complexes).
We use a standard MCM method to compute a first average of the QRS complex. Next, the averaged QRS complex is used as a new template for a second MCM computation. The normalized signal N_QRS (t) is then estimated by performing  a third order spline interpolation on the averaged QRS signal calculated with the two step MCM method. Figure 2 shows the averaged signal of a one heart cycle and the spline interpolation of the averaged QRS complex.
 
 

Figure 2. Averaged MCG-signal with QRS starting at t=0 ms and spline interpolation of the averaged  signal in the time interval from -20 ms to 169 ms.
 

Subsequently, the spline function N_QRS (t) is employed to fit the function
 

      S_QRS(t) = S_QRS (t, A, S₀ , S₁, t_beat, N_QRS (t))
 

to every QRS complex through a nonlinear regression (Levenberg-Marquardt-Method [13]). The regression parameters are the amplitude A, the time event t_beat, and linear baseline parameters (S₀, S₁). In order to avoid spline interpolation errors introduced by the fitting interval edges, we use a shorter time interval (10 ms on both sides) for the nonlinear regression. The regression for L is applied in the time interval 30 ms to 89 ms (QRS complex duration), where this parameter is relevant.
 

Results

Figure 3. Regression parameters over t_beat, their empirical distributions, and power spectra.

The modulation of both A and L fluctuate up to  ± 3 % and have a typical cycle duration of the respiration. The histograms on the right side in Figure 3 illustrate  that  the distribution basically consists  of two peaks, the inhaled and exhaled state peak, respectively. The power spectrum of the amplitude in Figure 3 exhibit a peak at approximately 0.30 Hz, which equals 18 respiration cycles. The time series of the linear baseline parameters (S₀, S₁) shows an influence of environmental noise on the linear baseline. However, the power spectra have a small peak at 0.30 Hz, that means a weak influence of respiration on the linear baseline.
There was no significant cross correlation between all parameters.
MCG based source reconstructions require the use of multi-channel systems. During the MCG recording it is important that the arrangement of signal sources and the volume conductor is fixed with respect to the sensor positions. Therefore, we quantify the influence of the alternating arrangement caused by respiration on the signal preprocessing in sensor arrays. As a first step we analyze the cross correlation of the estimated regression parameters for two adjacent pick-up coil positions.
In Figure 4 the regression parameters of both sensors for 90 heart cycles are plotted against each other, where the statistical regression errors are shown as bars. The solid lines illustrate the linear regressions. The gradient of the linear regression function is 1.5 for the amplitude A. This means, respiration causes a greater absolute fluctuation of the QRS amplitude at the position of sensor 2 than at the position of sensor 1 although the absolute signal is larger in sensor 1.
 

Figure 4. Averaged signals (a) and linear cross correlation for the regression parameters of two adjacent measurement channels (b-f) at the MCG field maximum.

Conclusions

References

[1] Achenbach, S., Moshage, W., Flüg, M., Ropers, D., Schibgilla, V. and Bachmann, K., "Vergleich der Time-Domain-Analyse von Late Potentials im
  signalgemittelten Elektrokardiogramm und Magnetokardiogramm", Biomedizinische Technik, 41 (suppl. 1), pp. 298-299, 1996.

[2] Mäkijärvi, M., "Magnetocardiography and risk analysis", in Hoke, M., Erné, S., Okada, Y. and Romani, G. (eds) "Biomagnetism: Clinical aspects", Elsevier Science Publishers  B.V., pp. 523-529, 1992.

[3] Mäkijärvi, M., Montonen, J., Toivonen, L., Leiniö, M., Siltanen, P. and Katila, T., "High-resolution magnetocardiography can identify ventricular tachycardia patients after myocardial infarction", in Hoke, M., Erné, S., Okada, Y. and Romani, G. (eds) "Biomagnetism: Clinical aspects", Elsevier Science Publishers  B.V., pp. 483-486, 1992.

[4] Mäkijärvi, M., Montonen, J., Toivonen, L., Leiniö, M., Siltanen, P. and Katila, T., "High-resolution and signal-averaged electrocardiography to separate post-myocardial infarction patients with and without ventricular tachycardia", European Heart Journal, 15, pp. 189-199, 1994.

[5] Simson, M., "Use of signals in the terminal QRS complex of identify patients with ventricular tachycardia after myocardial infarction", Circulation, 64(2), pp. 235-242, 1981.

[6] Leder, U., Haueisen, J., Huck, M. and Nowak, H., "Non-invasive imaging of arrhythmogenic left-ventricular myocardium after infarction", Lancet, 352, p. 1825, 1998.

[7] Moshage, W., Achenbach, S., Weikl, A., Göhl, K., Abraham-Fuchs, K., Schneider, S. and Bachmann, K., "Progress in Biomagnetic Imaging of Heart Arrhythmias", in Baert and Heuck (eds) "Frontiers in European Radiology", vol. 8, Berlin, Springer Verlag, pp.1-19, 1991.

[8] Wischmann, H.-A., Dössel, O. and Fuchs, M., "Effect of the signal-to-noise ratio on the quality of linear estimation reconstruction of distributed current sources", Brain Topography, 5(2), pp. 189-194, 1992.

[9] Escalona, O., Mitchell, R., Balderson, D. and Harron, D., "Fast and reliable QRS alignment technique for high-frequency analysis of signal-averaged ECG", Med. Biol. Eng. Comput., Suppl., pp. 137-146, 1993.

[10] Horacek, B., Purcell, C., Lamothe, R., Merritt, R., Kafer, C., Paryalwar, S., Dey, S., Leon,  L. and Stroink, G., "The effect of torso geometry on magnetocardiographic isofield maps", Physics in Medicine and Biology, 32(1), pp. 121-124, 1987.

[11] Huck, M., Haueisen, J., Hönecke, O., Fritschi, T., Leder, U. and Nowak, H., "QRS Amplitude and Shape Variability in Magnetocardiograms", PACE, vol. 23, no. 2, pp. 234-242, 2000.

[12] Dössel, O., David, B., Fuchs, M., Krüger, J., Lüdeke, K.-M. and Wischmann,  H.-A., "A 31-Channel SQUID-System for Biomagnetic Imaging", Applied Superconductivity, 1(10-12), pp. 1813-1825, 1993.

[13] Press, W., Teukolsky, S., Vetterling, W. and Flannery, B., "Numerical Recipes in C", Cambridge University Press, second edition, 1992.