Skip to content
BY-NC-ND 4.0 license Open Access Published by De Gruyter September 30, 2016

Derivation of the respiratory rate from directly and indirectly measured respiratory signals using autocorrelation

  • Fabian Schrumpf EMAIL logo , Matthias Sturm , Gerold Bausch and Mirco Fuchs

Abstract

The estimation of respiratory rates from contineous respiratory signals is commonly done using either fourier transformation or the zero-crossing method. This paper introduces another method which is based on the autocorrelation function of the respiratory signal. The respiratory signals can be measured either directly using a flow sensor or chest strap or indirectly on the basis of the electrocardiogram (ECG). We compare our method against other established methods on the basis of real-world ECG signals and use a respiration-based breathing frequency as a reference. Our method achieved the best agreement between respiration rates derived from directly and indirectly measured respiratory signals.

1 Introduction

The continuous measurement of the respiratory rate is particularly important in the medical environment as well as in the area of sports applications. Here, the most widely used approach is the measurement of the movements due to breathing using a chest strap or an impedance-based pneumographic device [1], [2]. However, there are also algorithms that derive the respiratory frequency noninvasively from other biosignals. It has been shown that the amplitude of the QRS-complex is modulated by chest movements during breathing, which results in a changing R-peak amplitude that correlates highly with the actual respiratory signal (taken from a chest strap or a flow sensor) [3]. From this signal one can estimate the breathing rate using a great variety of algorithms. The most commonly used methods are the Fast-Fourier Transform (FFT-method) and the zero-crossing method. The former estimates the rate from the strongest frequency component in the respiration signal, whereas the latter is based on determining the average duration of each breath [3], [4].

Interferences in the estimated respiration signal can lead to deviations in the estimated respiration rate. In this paper we propose a new method that allows a less interference-prone estimation of the respiratory rate. We applied the method to respiration signals derived both from ECG and a chest strap and evaluated the results in comparison to FFT and zero-crossing based respiration rate estimation methods.

2 Methods

2.1 Estimation of the respiratory rate

As a first step the respiratory signal is extracted from the ECG-signal sampled with the sampling frequency fs. For this purpose, a Hamilton-Tompkins algorithm [5] is used to detect all QRS-complexes in the ECG signal (comp. Figure 1). Next, the amplitude of the R-to-S peak is computed, which provides a robust estimation of the signal amplitude with respect to the baseline [6]. These extracted amplitudes are then interpolated by means of cubic splines to obtain the respiratory signal with the same sampling frequency as the ECG signal. An example for an extracted respiratory signal can be seen in Figure 1. Aside from a phase shift there is a high correlation between the extracted and the reference signal.

Figure 1 Top: baseline-corrected ECG signal with marked R- and S-peaks for the estimation of the R-to-S amplitude; bottom: reference respiratory signal (blue) compared to ECG-derived respiratory signal (orange). Amplitudes are normalized, arbitrary units.
Figure 1

Top: baseline-corrected ECG signal with marked R- and S-peaks for the estimation of the R-to-S amplitude; bottom: reference respiratory signal (blue) compared to ECG-derived respiratory signal (orange). Amplitudes are normalized, arbitrary units.

In order to estimate the respiration rate, the fundamental frequency of the extracted respiratory signal has to be calculated. As mentioned above, this can be done using methods like the FFT in the frequency domain [4] or the zero-crossing method in the time domain [3]. In this paper we present a different approach which is based on the autocorrelation function of the estimated respiratory signal. The autocorrelation function Ψ(i) of the respiratory signal xresp is defined in Eq. 1 [7]:

(1)Ψ(i)=n=0Nxresp(n)xresp(ni)

Here, N and i denote the number of samples and the displacement of the signal towards itself, respectively. Generally, the autocorrelation function describes the similarity of a signal to a shifted version of itself. In the case of periodic signals, this function shows periodic maxima. The displacements i at which those maxima occur correspond to integer multiples of the fundamental period of the signal (here, i.e. the respiratory signal). Employing the Wiener-Khinchin theorem [8] allows an efficient calculation of Eq. 1 on the basis of an FFT. Next, the intervals ti of all neighbouring local maxima of the autocorrelation function are determined. Disturbances in the extracted respiratory signal or its non-stationarity may cause additional local maxima. Figure 2 shows one such example. Here it can be seen that high-frequency disturbances of a respiratory signal lead to local maxima (denoted in red) in the autocorrelation function that are unrelated to the fundamental period of the respiratory signal. To suppress the influence of these additional maxima an empirically estimated threshold is used so that only maxima exceeding this threshold were included in subsequent calculations. To further weaken the influence of highly deviating time intervals their median is usesd as the final estimate for the duration of a breath cycle. The respiratory frequency in breaths per minute (bpm) can then be calculated using Eq. 2.

Figure 2 Top: ECG-derived respiratory signal (example); bottom: local maxima belonging to the fundamental period (green) and local maxima caused by artifacts (red).
Figure 2

Top: ECG-derived respiratory signal (example); bottom: local maxima belonging to the fundamental period (green) and local maxima caused by artifacts (red).

(2)fresp=60fs/Tresp

The first reference method to estimate respiratory rates is the FFT. It is used to calculate the spectrum of the respiratory signal. The frequency component fmax with the highest amplitude is thaken as the respiratory frequency. This frequencies can be converted to the respiratory rate in bpm using Eq. 3.

(3)fresp=60fmax

For the zero-crossing method, i.e. the second reference method, the first step is to center the signal at zero by subtracting the mean respiratory signal from each sample. After that, all zero-crossings ti are detected. The respiratory rate can finally be estimated by calculating the mean of all detected zero-crossings using Eq. 4.

(4)fresp=602Ni=1Nti

2.2 Simulation setup

Simultaneous recordings of ECG and respiration were used to evaluate the method. The recordings of 13 subjects with a length of 30 min each were taken from the Combined Measurement of ECG, Breathing and Seismocardiogram (CEBS) database [9], [10]. All datasets contained an ECG-signal as well as a respiratory signal measured by a chest strap. The datasets were converted to a MATLAB-readable format using the WFDB toolbox [11]. From each dataset ECG and respiratory signal were obtained and resampled to a sampling frequency of fs = 250 Hz. Both signals were subdivided into sliding windows of length Tmeas = 30 s and an overlap of 15 s. Baseline removal was applied to the ECG-signal [12]. Furthermore, the 50 Hz powerline interference was removed using an IIR-Notch-Filter. The respiratory signal was smoothed using an FIR-lowpass-filter (fcut = 2 Hz). The method described in section 2.1 was used for the estimation of the respiratory rate by means of directly measured respiratory signals (chest strap) and ECG signals. Subsequently, the results are then compared to the ones obtained for the FFT-method and the zero-crossing method.

3 Results

Figure 3 shows the raw results of all evaluated methods. The horizontal axis represents estimated respiratory rate based on the ECG and the vertical axis represents the respiratory rate based on the chest strap signal. All methods show a high correlation between respiratory rates estimated from ECG and chest strap. Both the autocorrelation method (r = 0.86, p < 10−3) and the zero-crossing method (r = 0.88, p < 10−3) exhibit a strong linear relationship between both estimated respiratory rates whereas this relation is much weaker for the FFT-method (r = 0,6, p < 10−3). The statistical significance of the correlation was tested using a t-test. It can be seen that the reference-based respiratory rate estimated by means of the FFT-method shows discrete values. This is due to the finite frequency resolution. The ECG-based respiration frequency does not show this effect. This is due to the spline interpolation of the R-to-S peak amplitudes leading to slightly varying sample-sizes for each time window and thus to varying frequency resolutions depending on the length of the time window Tmeas.

Figure 3 Correlation between respiration-derived respiratory rate and ECG-derived estimated respiratory rate using the autocorellation method (A), FFT-method (B) and zero-crossing method (C).
Figure 3

Correlation between respiration-derived respiratory rate and ECG-derived estimated respiratory rate using the autocorellation method (A), FFT-method (B) and zero-crossing method (C).

The error between ECG- and chest-strap-based respiration rates for each method is displayed in Figure 4. A Kolmogorov-Smirnov Test was used to compare the error distibutions. It turned out that all distributions were statistically significantly different (p < 10−3). While the median of the autocorrelation method is located around zero, the median of the FFT-method lies above zero and shows a higher spread. The median of the zero-crossing method is also close to zero. However, its error distribution exhibits an asymmetric deviation from the median. Values for mean and variance of the error are shown in Table 1. It can be seen that the autocorrelation method has a smaller error in average than any of the reference methods. In contrast to the zero-crossing method there is a huge difference regarding the variance of the error between autocorrelation- and FFT-method.

Figure 4 Comparison of the error of respiration-based and ecg-based estimation of the respiratory rate using all methods considered. Boxes show interquartile range (IQR = Q0.75 − Q0.25) and median. Outliers are depicted as dots outside of the boxes.
Figure 4

Comparison of the error of respiration-based and ecg-based estimation of the respiratory rate using all methods considered. Boxes show interquartile range (IQR = Q0.75 − Q0.25) and median. Outliers are depicted as dots outside of the boxes.

Table 1

Mean and variance of the error distributions shown in Figure 4.

μ (err)Var (err)
Autocorrelation−0.153.68
FFT−0.817.54
Zero-crossing−0.442.7

An inter-subject analysis emphasizes these findings. It can bee seen in Figure 5 that the error of the zero-crossing method shows a higher spread than the error of the autocorrelation method. We have omitted the FFT-method since we consider it unsuited here due to its discrete-valued results.

Figure 5 Error of the autocorrelation- and zero-crossing method on an inter-subject level. Boxes show interquartile range (IQR = Q0.75 − Q0.25) and median. Outliers are depicted as dots outside of the boxes.
Figure 5

Error of the autocorrelation- and zero-crossing method on an inter-subject level. Boxes show interquartile range (IQR = Q0.75Q0.25) and median. Outliers are depicted as dots outside of the boxes.

4 Discussion

We compared ECG-derived and chest-strap-derived respiratory rates using three different methods, i.e. FFT, zero-crossing and the newly proposed autocorrelation method. It became apparent that there is a high correlation between the results of zero-crossing and the autocorrelation. Yet, the correlation between the FFT-method and the two others is considerably lower. The results also show a low variance of the error for zero-crossing and autocorrelation as opposed to the FFT-method. However, the ECG-based respiratory rate estimates by means of zero-crossing are systematically lower than the chest-based estimation. Apparently, the autocorrelation method is less prone to interferences in the respiratory signal than the other methods. These interferences can be caused by either artifacts in the ECG signal or as a result of natural changes of the breath frequency in the analysed time interval. Another disadvantage of the FFT-method is its finite frequency resolution in the time window under consideration. It is also possible that an insufficient baseline removal leads to low-frequency components in the derived respiratory signal and therefore to a larger error. This effect was reduced by calculating the amplitude of the R-peak with respect to the S-peak but can not be completely eliminated.

Algorithms used for preprocessing can also have a deteriorating effect on the estimation of the respiratory frequency. FIR-filters used to enhance the signal-to-noise ratio almost always also remove important signal components, which in turn can lead to changed R-peak amplitudes. Other techniques (i.e. wavelet filters) might help to mitigiate this effect.

We have shown that the proposed autocorrelation method is generally feasible for the estimation of the respiration rate from both ECG-based and chest strap based respiratory signals. It is applicable to the long-term monitoring of patients in a clinical environment without directly measuring the respiratory signal. This is particularly useful, e.g. during rehabilitation after a surgical intervention. The used data exhibits a rather low variability in terms of respiratory frequency. An evaluation of our method using very low and very high respiratory frequencies is necessary to further substantiate its usefulness. An extension of our method to PPG signals may also be conceivable.

Author’s Statement

Research funding: This work was funded by the Federal Ministry of Education and Research (BMBF) (FKZ 03FH032IX5). Conflict of interest: Authors state no conflict of interest. Material and Methods: Informed consent: Informed consent has been obtained from all individuals included in this study. Ethical approval: The conducted research is not related to either human or animal use.

References

[1] Grenvik A, Ballou S, McGinley E, Millen JE, Cooley WL, Safar P. Impedance pneumography: comparison between chest impedance changes and respiratory volumes in 11 healthy volunteers. Chest J. 1972;62:439–43.10.1378/chest.62.4.439Search in Google Scholar

[2] Gupta AK. Respiration rate measurement basetd on impedance pneumography. Texas Instruments, Tech. Rep. February, 2011.Search in Google Scholar

[3] Moody GB, Mark RG, Zoccola A, Mantero S. Derivation of respiratory signals from multi-lead ECGs. Compu Cardiol. 1985;12:113–6.Search in Google Scholar

[4] Karlen W, Mattiussi C, Florean D. Sleep and wake classification with ECG and respiratory effort signals. IEEE Trans Bio-Med Circuits Syst. 2009;3:71–8.10.1109/TBCAS.2008.2008817Search in Google Scholar

[5] Hamilton P, Tompkins WJ. Quantitative investigation of QRS detection rules using the MIT/BIH arrhythmia database. IEEE Trans Bio-Med Eng. 1986;33:1157–65.10.1109/TBME.1986.325695Search in Google Scholar

[6] Dobrev D, Daskalov I. Two-electrode telemetric instrument for infant heart rate and apnea monitoring. Med Eng Phys. 1999;20:729–34.10.1016/S1350-4533(98)00091-5Search in Google Scholar

[7] McLeod P, Wyvill G. A smarter way to find pitch. In: Proceedings of International Computer Music Conference, ICMC; 2005.Search in Google Scholar

[8] Box GE, Jenkins GM, Reinsel GC. Autocorrelation function and spectrum of stationary processes. In: Time series analysis. 4th ed. John Wiley & Sons, Inc., 2013, ch. Autocorrel, p. 19–47.Search in Google Scholar

[9] García-González M, Argelagós A, Fernández-Chimeno M, Ramos-Castro J. Differences in QRS locations due to ECG lead: relationship with breathing. In: XIII Mediterranean Conference on Medical and Biological Engineering and Computing 2013, 2014;13:962–4.10.1007/978-3-319-00846-2_238Search in Google Scholar

[10] Goldberger AL, Amaral LA, Glass L, Hausdorff JM, Ivanov PC, Mark RG, et al. Physiobank, physiotoolkit, and physionet components of a new research resource for complex physiologic signals. Circulation. 2000;101:215–20.10.1161/01.CIR.101.23.e215Search in Google Scholar

[11] Silva I, Moody GB. An open-source toolbox for analysing and processing physionet databases in MATLAB and Octave. J Open Res Softw. 2014;2:e27.10.5334/jors.biSearch in Google Scholar PubMed PubMed Central

[12] Zhang F, Chen S, Zhang H, Zhang X, Li G. Bioelectric signal detrending using smoothness prior approach. Med Eng Phys. 2014;36:1007–13.10.1016/j.medengphy.2014.05.009Search in Google Scholar PubMed

Published Online: 2016-9-30
Published in Print: 2016-9-1

©2016 Fabian Schrumpf et al., licensee De Gruyter.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Downloaded on 29.3.2024 from https://www.degruyter.com/document/doi/10.1515/cdbme-2016-0054/html
Scroll to top button