The main purpose of this work is the estimation of the respiratory rate (r R ) from the electrocardiogram (ECG). In this study, the r R was estimated from the same Fit directions of maximum projection (FD) used for multi lead ECG automatic delineation (ML). A previously developed and validated methodology for boundaries location was extended to include wave peaks and estimate FD. The median of power spectral density obtained over the directions based on QRS complex main peak, T wave peak and end spatial loops was used for r R estimation. In a control database, the proposed method yielded more accurate r R estimates (mean absolute error (MAE), 2.64 bpm, SD=3.92) than the estimates based on the single-lead ECG R-peak amplitude (MAE values from 3.29 to 5.26 bpm, SD ≥ 5) and RR series (2.89 to 3.66 bpm, SD >4.6), close to results from EDR method (2.89 bpm, SD=3.63).
IntroductionDirect estimation of respiratory rate (r R ) is done with devices that are intrusive, expensive and uncomfortable for the patient [1]. Therefore the development of alternative methods for indirect r R estimation is of great interest. The respiratory and heart activities are linked through physiological processes. The respiration modulates the heart rate such that it increases during inspiration and decreases during expiration [2]. In the same way, during the respiratory cycle, chest movements and changes in the thorax impedance distribution due to filling and emptying of the lungs cause a rotation of the electrical axis of the heart which has an effect over beat morphology [3,4]. Different signal processing techniques have been developed to extract respiratory information from the modulation of the ECG R-peak amplitudes (RPA) and heart rate (HR) [4][5][6][7][8]. However, the single-lead RPA and RR series have been shown not to yield accurate r R estimates as the respiratory influence in a given ECG lead is subject dependent [4]. Often, in ambulatory and clinical applications, multilead ECG recordings are available, allowing to construct vectocardiographic (VCG) loops from orthogonal leads. In [4] the series of rotation angles of the heart's electrical axis as induced by respiration are estimated from least-squares loop alignment, and the r R estimated spectrally (EDR). The main directions on the wavelet transform (WT) of VCG loops, taken as parallel to each ECG fiducial point, have been successfully used for multilead delineation of wave boundaries [9]. Moreover a strong relation between the respiration and the T wave end based Fit directions of maximum projection (FD) used the by multi lead ECG delineation system (ML) was found in [10]. Taking advantage of this relation, in [11] we have been able to estimate r R with good results from the FD used for T wave end delineation. In this study we aim to combine the information based in FD found for the delineation of QRS peak and T wave peak and end to obtain a single r R estimation.
Methods
DataThe data used in this study was a database of ECG signals from healthy subjects recorded at Uni...