# Patent

## Title
Respiratory rate measurement using a combination of respiration signals

## Abstract

A method and system for measuring respiratory rate are disclosed. In a first aspect, the method comprises measuring at least one respiration signal and filtering the respiration signal using a lowpass filter. The method includes peak-picking the respiration signal to determine the respiratory rate and determining a quality metric of the respiratory rate. In a second aspect, the system comprises a wireless sensor device coupled to a user via at least one electrode, wherein the wireless sensor device includes a processor and a memory device coupled to the processor, wherein the memory device stores an application which, when executed by the processor, causes the processor to carry out the steps of the method.

## Background
FIELD OF THE INVENTION
The present invention relates to sensors, and more particularly, to a sensor device utilized to measure respiratory rate using a combination of respiration signals.
BACKGROUND
Respiratory rate (RR) is an important vital sign that is useful for examining periods of distress. High and low respiratory rates are often symptoms of serious diseases including heart failure (Cheyne-Stokes breathing), obstructive sleep apnea (cessation of breathing), and metabolic acidosis (hyperventilation). The accurate measurement of respiratory rates using non-intrusive sensors enables the continuous monitoring of a person's respiratory rate.
Conventional methods of measuring a person's respiratory rate include measuring each breath when the person is at rest over a predetermined time period by counting how many times the person's chest rises. Devices such as stethoscopes are often utilized by doctors to assist as well. However, these conventional methods are inefficient, intrusive, and require time consuming analysis to measure the respiratory rate. Therefore, there is a strong need for a cost-effective solution that overcomes the above issue by utilizing a combination of sensors and signals to measure RR thereby providing increased accuracy over the measurement of RR using any single sensor/signal. The present invention addresses such a need.
SUMMARY OF THE INVENTION
A method and system for measuring respiratory rate are disclosed. In a first aspect, the method comprises measuring at least one respiration signal and filtering the respiration signal using a lowpass filter. The method includes peak-picking the respiration signal to determine the respiratory rate and determining a quality metric of the respiratory rate.
In a second aspect, the system comprises a wireless sensor device coupled to a user via at least one electrode, wherein the wireless sensor device includes a processor and a memory device coupled to the processor, wherein the memory device stores an application which, when executed by the processor, causes the processor to measure at least one respiration signal and filter the respiration signal using a lowpass filter. The system further causes the processor to peak-pick the respiration signal to determine the respiratory rate and to determine a quality metric of the respiratory rate.



BRIEF DESCRIPTION OF THE DRAWINGS
The accompanying figures illustrate several embodiments of the invention and, together with the description, serve to explain the principles of the invention. One of ordinary skill in the art will recognize that the embodiments illustrated in the figures are merely exemplary, and are not intended to limit the scope of the present invention.
FIG. 1 illustrates a wireless sensor device in accordance with an embodiment.
FIG. 2 illustrates a flow chart of a method for measuring respiratory rate in accordance with an embodiment.
FIG. 3 illustrates a flow chart of an accelerometer respiration process in accordance with an embodiment.
FIG. 4 illustrates a flow chart of a peak detection process in accordance with an embodiment.
FIG. 5 illustrates a flow chart of a quality metric process in accordance with an embodiment.
FIG. 6 illustrates a diagram that compares the quality metric to the respiration signal in accordance with an embodiment.
FIG. 7 illustrates a flow chart of a QRS respiration process 700 in accordance with an embodiment.
FIG. 8 illustrates a diagram of respiration signal computations in accordance with an embodiment.



DETAILED DESCRIPTION
The present invention relates to sensors, and more particularly, to a sensor device utilized to measure respiratory rate using a combination of respiration signals. The following description is presented to enable one of ordinary skill in the art to make and use the invention and is provided in the context of a patent application and its requirements. Various modifications to the preferred embodiment and the generic principles and features described herein will be readily apparent to those skilled in the art. Thus, the present invention is not intended to be limited to the embodiments shown but is to be accorded the widest scope consistent with the principles and features described herein.
Utilizing a combination of sensors and signals to estimate respiratory rate (RR) provides increased accuracy over the conventional utilization of any single sensor/signal. A method and system in accordance with the present invention's combination requires knowledge of the quality of the respiratory signal obtained from each individual sensor. The respiratory signals from a variety of sources can be utilized as inputs to the method and system in accordance with the present invention. These sources include but are not limited to mechanical signals from a chest-mounted accelerometer, EKG-derived respiration signals, bend sensors, stretch sensors (respiratory inductance plethysomography), impedance pneumography, flow sensors, and CO2 sensors.
A method and system in accordance with the present invention determines a person's respiratory rate using respiration signals derived from a plurality of sources. In one embodiment, the plurality of sources include a chest-mounted accelerometer where the correct axis is perpendicular to gravity, a QRS area/amplitude where the area/amplitude of the QRS complex of an EKG signal is computed, a respiratory sinus arrhythmia (RSA) where times between subsequent R peaks of the EKG signal (R-R intervals) are computed, and any other respiration signals of interest.
After a wireless sensor device attached to a person and the wireless sensor device detects a plurality of respiration signals, each signal is filtered by the wireless sensor device to get activity in a respiratory frequency range. In one embodiment, the filtering is done by a lowpass filter at 0.7 Hertz (Hz). An application stored on the wireless sensor device conducts peak-picking of each respiration signal to detect a number of breaths that the person takes and computes the respiratory rate by using a changing threshold based on a standard deviation of the signal at that time point.
The wireless sensor device analyzes the respiratory rate and determines a quality metric between 0 and 1 based on regularity of minima times and amplitudes and a noisiness of the respiration signal. The quality metric enables the wireless sensor device to determine whether each derived respiratory rate can be trusted. The quality metric is inversely correlated to expected error with 0 denoting high error and 1 denoting low error. Additionally, smoothing is performed on the respiratory rate to reduce effects of missed/extraneous peaks using a trimmed mean.
One of ordinary skill in the art readily recognizes that a variety of wireless sensor devices can be utilized to measure the person's respiratory rate including portable wireless sensor devices with embedded circuitry in a patch form factor and that would be within the spirit and scope of the present invention.
To describe the features of the present invention in more detail, refer now to the following description in conjunction with the accompanying Figures.
FIG. 1 illustrates a wireless sensor device 100 in accordance with an embodiment. The wireless sensor device 100 includes a sensor 102, a processor 104 coupled to the sensor 102, a memory 106 coupled to the processor 104, an application 108 coupled to the memory 106, and a transmitter 110 coupled to the application 108. The sensor 102 obtains data from the user and transmits the data to the memory 106 and in turn to the application 108. The processor 104 executes the application 108 to process ECG signal information of the user. The information is transmitted to the transmitter 110 and in turn relayed to another user or device.
In one embodiment, the sensor 102 comprises two electrodes to measure cardiac activity and an accelerometer to record physical activity and posture and the processor 104 comprises a microprocessor. One of ordinary skill in the art readily recognizes that a variety of devices can be utilized for the processor 104, the memory 106, the application 108, and the transmitter 110 and that would be within the spirit and scope of the present invention.
FIG. 2 illustrates a flow chart of a method 200 for measuring respiratory rate in accordance with an embodiment. Referring to FIGS. 1 and 2 together, the method 200 comprises the wireless sensor device 100 measuring at least one respiration signal, via 202, and filtering the respiration signal using a lowpass filter, via 204. The method 200 includes peak-picking the respiration signal to determine a number of breaths and the respiratory rate, via 206, and determining a quality metric of the respiratory rate, via 208. In one embodiment, at least one respiration signal is derived from a variety of sources including but not limited to a chest-mounted accelerometer, a QRS-area and a QRS-amplitude of an EKG signal, and a respiratory sinus arrhythmia (RSA).
With a chest-mounted accelerometer embedded within a wireless sensor device that has been attached to a person, breathing can be detected during periods of low activity. During inspiration, the chest expands and the angle of the accelerometer to gravity changes slightly. An angle change of e yields a magnitude change of sine approximately equal to θwhen θis close to 0. When the person is upright or lateral, the best axis to examine is the anterior-posterior axis and when the person is supine, the best axis to examiner is the inferior-superior axis (head to toe).
The respiration signal derived from a chest-mounted accelerometer has minima and maxima that correspond to inspiration and expiration. The minima and maxima can correspond to either the inspiration or the expiration. The magnitudes of the respiration signal's peaks are not uniform so a method and system in accordance with the present invention utilizes an accelerometer respiration process with a time-varying threshold based on a standard deviation of the respiration signal in a small predetermined time period window. By determining differences between the peak times of the respiration signal, an instantaneous respiratory rate is computed.
FIG. 3 illustrates a flow chart of an accelerometer respiration process 300 in accordance with an embodiment. The accelerometer respiration process 300 obtains rotated acceleration signals Ax, Ay, Az, via 302, and computes respiratory signal (R[n]) based upon the correct axis that is perpendicular to gravity, via 304. A rolling standard deviation (σ[n]) of the respiratory signal R[n] is computed within a predetermined 45 second window (n) per the equation: σ[n]=√{square root over (Σ(R[i]−R)2)}, via 306. Using the rolling standard deviation, a peak threshold is computed per the equation: th[n]=α*σ[n], via 308, where α is the number of standard deviations (e.g. 1.3) for the peak threshold th[n].
After computing the peak threshold, a peak detection process is utilized to determine the minima and maxima of the respiratory signal R[n], via 310. A respiratory rate quality metric is determined using a number of features that judge the regularity of the respiration peaks, via 312, and an instantaneous respiratory rate (RR) is computed as the first difference between minima times per the equation: RR=60/(tmin[i]−tmin[i−1]), via 314. Respiratory rates outside of an interquartile range (25-75%) within the predetermined 45 second window are discarded, via 316, and a trimmed mean respiratory rate is computed, via 318, to smooth the computed instantaneous RR.
In one embodiment, the respiratory signal (R[n]) computation via 304 determines an axis that is perpendicular to a gravity vector (G=[gx, gy, gz]. When a person is upright, the axis that is perpendicular to gravity is the Z axis (anterior-posterior axis) where G×[1;0;0]=[0;0;1]. When a person is supine, the axis that is perpendicular to gravity is the Y axis (inferior-superior axis) where G x [1;0;0]=[0;1;0]. When a person is lateral decubitus, the axis that is perpendicular to gravity is empirically found as the Z axis where G x [0;1;0]=[0;0;1]. Thus, the more that gravity is pointing in the X direction, the greater the contribution of the Y axis to the cross product.
In this embodiment, low-pass filter accelerations below 0.02 Hz for gravity G and low-pass filter accelerations below 0.7 Hz for respiration ALP are determined. The angle θ between the gravity and the X axis is computed per the equation: θ=cos−1({right arrow over (G)}·{right arrow over (i)})=cos−1gx. Using θ, the perpendicular-to-gravity direction for the respiration signal is computed per the following set of equations, where P is the direction for computing the respiration signal, and C=∥{right arrow over (g)}×[sin θ cos θ 0]T∥:




  
    
      
        
          P
          =
          
            
              G
              →
            
            ×
            
              
                [
                
                  
                    
                      
                        sin
                        ⁢
                        
                          
                        
                        ⁢
                        θ
                      
                    
                  
                  
                    
                      
                        cos
                        ⁢
                        
                          
                        
                        ⁢
                        θ
                      
                    
                  
                  
                    
                      0
                    
                  
                
                ]
              
              /
              C
            
          
        
        ⁢
        
          

        
        ⁢
        
          P
          =
          
            
              G
              →
            
            ×
            
              
                [
                
                  
                    
                      
                        
                          1
                          -
                          
                            g
                            x
                            2
                          
                        
                      
                    
                  
                  
                    
                      
                        g
                        x
                      
                    
                  
                  
                    
                      0
                    
                  
                
                ]
              
              /
              C
            
          
        
        ⁢
        
          

        
        ⁢
        
          P
          =
          
            
              [
              
                
                  
                    
                      
                        -
                        
                          g
                          x
                        
                      
                      ⁢
                      
                        g
                        z
                      
                    
                  
                
                
                  
                    
                      
                        g
                        z
                      
                      ⁢
                      
                        
                          1
                          -
                          
                            g
                            x
                            2
                          
                        
                      
                    
                  
                
                
                  
                    
                      
                        g
                        x
                        2
                      
                      -
                      
                        
                          g
                          y
                        
                        ⁢
                        
                          
                            1
                            -
                            
                              g
                              x
                              2
                            
                          
                        
                      
                    
                  
                
              
              ]
            
            /
            
              C
              .
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




Using P, an accelerometer-based respiratory signal is computed per the following set of equations, where R is the respiratory signal:




  
    
      
        
          R
          =
          
            P
            ·
            
              A
              LP
            
          
        
        ⁢
        
          

        
        ⁢
        
          R
          =
          
            
              [
              
                
                  
                    
                      
                        -
                        
                          g
                          x
                        
                      
                      ⁢
                      
                        g
                        z
                      
                    
                  
                
                
                  
                    
                      
                        g
                        z
                      
                      ⁢
                      
                        
                          1
                          -
                          
                            g
                            x
                            2
                          
                        
                      
                    
                  
                
                
                  
                    
                      
                        g
                        x
                        2
                      
                      -
                      
                        
                          g
                          y
                        
                        ⁢
                        
                          
                            1
                            -
                            
                              g
                              x
                              2
                            
                          
                        
                      
                    
                  
                
              
              ]
            
            ·
            
              
                A
                LP
              
              /
              
                C
                .
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




In one embodiment, the peak detection/peak-picking process that is utilized to determine the minima and maxima of the respiratory signal via 310 utilizes a changing threshold based on a standard deviation of the respiration signal at a certain time point. The peak detection process ensures that a peak-to-peak amplitude of the respiration signal is greater than a certain threshold (th[n]) that depends on the variability of the respiration signal.
FIG. 4 illustrates a flow chart of a peak detection process 400 in accordance with an embodiment. The peak detection process 400 can be utilized for respiration signals from a variety of sources including but not limited to the chest-mounted accelerometer, QRS-area, and RSA. The peak detection process 400 sets left_min=inf, temp_max=inf, and peak_count=0 via 402 and obtains a sample of a respiration signal R[n] via 404. A sign of the first difference of the respiration signal is computed per the equation: d[n]=sign(R[n]−R[n−1]), via 406.
The peak detection process 400 analyzes whether d[n]*d[n−1]=−1 via 408. If no (d[n]*d[n−1] does not equal −1), the process 400 returns back to 404 to get another sample of the respiration signal R[n]. If yes (d[n]*d[n−1]=−1), the process 400 determines a peak-to-peak threshold th[n−1] via 410 and analyzes whether d[n]=−1 via 412. If no (d[n] does not equal −1), the process 400 analyzes whether temp_max>R[n−1]+th[n−1] via 414 and if yes (d[n]=−1), the process 400 analyzes whether R[n−1]>temp_max and R[n−1]>left_min +th[n] via 416. If no (R[n−1] is not greater than temp_max or R[n−1] is not greater than left_min +th[n]), then the process 400 returns back to 404 to get another sample of the respiration signal R[n]. If yes (R[n−1]>temp_max and R[n−1]>left_min+th[n]), then the process 400 sets temp_max=R[n−1]and temp_max_time=n−1 and returns back to 404 to get another sample of the respiration signal R[n].
If temp_max>R[n−1]+th[n−1], then the process 400 determines that a maxima is found at temp_max_time with a magnitude of temp_max and also determines that a minima is found at left_min_time with a magnitude of left_min via 420. If temp_max is not greater than R[n−1]+th[n−1], then the process 400 analyzes whether R[n−1]<left_min via 418. If no (R[n−1] is not less than left_min), then the process 400 returns back to 404 to get another sample of the respiration signal R[n]. If yes (R[n−1]<left_min), then the process sets left_min=R[n−1] and left_min_time=n−1 and then returns back to 404 to get another sample of the respiration signal R[n].
Therefore, if the peak detection process 400 determines that the first difference of the respiration signal (R[n]−R[n−1]) moves from positive to negative, the process 400 analyzes whether a peak of the respiration signal is greater than left_min+th[n] and whether the peak is greater than the last temp_max. If both of these conditions are met, the peak is set as the new temp_max. If the peak detection process 400 determines a derivative of the respiration signal moves from negative to positive, the process 400 analyzes whether the valley is less than th[n] of the temp_max and if so, the process 400 sets the temp_max as an actual maxima and the left_min as an actual minima. Otherwise, if the valley is less than the previous left_min, then the process 400 sets the valley as the new left_min. Every time a minima/maxima pair is found, the process 400 sets the current temp_max to −∞ and sets the current left_min to ∞.
In one embodiment, the respiratory rate quality metric determination via 312 is derived within a predetermined time period (e.g. 45 seconds) and is derived based on features including but not limited to a coefficient of variation (standard deviation over the mean) of a difference in minima times (ct), a coefficient of variation (standard deviation over the mean) of peak-to-peak values (cp), a mean peak-to-peak value (mp), and a ratio of a number of picked extrema to a total number of all local extrema (rext). In this embodiment, the quality metric is weighted using previously learned weights and transformed using an exponential to get the quality metric to a value between 0 and 1.
The quality metric for each input signal is estimated in 45 second windows using the aforementioned features that judge the regularity of the respiration peaks. A smaller ct equates to more regular peak times and thus a higher respiration signal quality. A smaller cp equates to more regular peak heights and thus a higher respiration signal quality. A larger mp indicates motion artifacts/noise and thus a lower respiration signal quality. A rext closer to 1 indicates the respiration signal has a minimal spurious peaks and thus a higher respiration signal quality.
The features are weighted using previously learned weights to get an estimation of the error (E) in the respiration signal. Weights are learned by least-squares regression on a collected dataset so that the learned weights allow for the best prediction of the error. In one embodiment, E is computed as E=w1*ct+w2*cp+w3*mp+w4*rext+w5 where typical weights include but are not limited to w1=3.8, w2=7.9, w3=9.8, w4=8.4, and w5=8.3. The error (E) is transformed using an exponential to get a quality metric (Q) between 0 and 1 per the equation Q=e−E/γ, where γ=5. The γ parameter is variable depending on the desired correspondence between estimated error and quality. The smaller the value of γ, the more the quality is reduced for a given error. When the estimated error (E) is high, Q gets closer to 0 which indicates a lower respiration signal quality. When the estimated error (E) is low and closer to 0, then Q gets closer to 1 which indicates a higher respiration signal quality.
If the mean peak-to-peak value mp is below a threshold min_p2p, the quality of the respiration signal (Q) is reduced by multiplying it by mp/min_p2p. Because mp/min_p2p is less than 1, the quality will become smaller. The small peak-to-peak values occur when the respiration signal is too small, and near noise, thereby making it more difficult to accurately estimate a respiratory rate from this respiration signal. Therefore, by reducing the quality metric in these situations, the algorithm is prevented from relying upon these signals.
FIG. 5 illustrates a flow chart of a quality metric process 500 in accordance with an embodiment. The quality metric process 500 obtains minima and maxima within a 45 second time window of the respiration signal as Rmin[i], Rmax[i], tmin[i], and tmax[i] via 502. A coefficient of variation of the difference in minima times (ct) is computed by the quality metric process 500 per the equation: ct=std(tmin[i+1]−tmin[i])/mean(tmin[i+1]−tmin[i]), via 504. For each minima, the process 500 averages the maxima on each side and subtracts the minima to get an average peak-to-peak (P[i]) per the equation: P[i]=0.5*(Rmax[i]+Rmax[i+1])−Rmin[i], via 506. The process 500 also obtains the minima and maxima from the last 45 second time window (Rhalfmin[i] and Rhalfmax[i]) using half the original threshold (th[n]/2) via 508.
Once the average peak-to-peak P[i] is computed, the process 500 computes a coefficient of variation of the mean peak-to-peak amplitudes (cp) per the equation: cp=std(P)/mean(P), via 510, and computes a mean of the peak-to-peak amplitudes in the 45 second time window per the equation: mp=mean(P), via 512. Once the minima and maxima from the last 45 second time window are obtained, the process 500 computes a ratio of the number of picked extrema to the total number of local extrema per the equation: rext=Npicked/Nall, via 514.
Following the computation of ct, cp, mp, and rext, the process 500 estimates an error (E) per the equation: E=w1*ct+w2*cp+w3*mp+w4*rext+w5, via 516, and sets E=0 if E is less than 0 via 518. The process 500 transforms the estimated error (E) to a quality metric (Q) for the current 45 second time window per the equation: Q=e−E/γ, via 520. The process 500 reduces the quality of the respiration signal when peak-to-peak is small if mp<min_p2p to result in the equation Q=Q*mp/min_p2p via 522.
The quality metric process 500 can be utilized for respiration signals from a variety of sources including but not limited to the chest-mounted accelerometer, a QRS-area, a QRS-amplitude, and an RSA. For the varying respiration signals, the weights for estimating the error (E) are calibrated for the particular respiration signal that is utilized. One of ordinary skill in the art readily recognizes that large motion artifacts in the accelerometer (or other respiration signal sourcing mechanism) data reduces the quality. FIG. 6 illustrates a diagram 600 that compares the quality metric to the respiration signal in accordance with an embodiment. As the respiration signal experiences noise during the start and end of the time window, the quality metric is closer to 0 but as the respiration signal experiences minimal noise during the middle of the time window, the quality metric is closer to 1.
In addition to deriving a respiration signal from a chest-mounted accelerometer, the respiration signal can also be derived via a QRS wave of an EKG signal. In one embodiment, the sensor 102 comprises electrodes that enable the wireless sensor device 100 to measure an EKG signal of a person when attached. The QRS wave of the EKG is modulated with respiration due to the mechanical effect of breathing on the heart. Breathing shifts the heart within the thorax and changes the cardiac axis which results in a change in QRS-amplitude during inspiration/expiration.
In one embodiment, the QRS-area or QRS-amplitude is computed as a measure of modulation of the QRS wave due to respiration using a QRS respiration process. In this embodiment, an EKG signal is detected by the wireless sensor device 100 and the EKG signal is filtered using 20 Hz of high-pass filtering. After the filtering, R-wave peaks are detected and an R-wave peak is analyzed utilizing a predetermined time window (e.g. 100 milliseconds) that is centered around the R-wave peak.
The absolute value of the respiration signal in the predetermined time window is summed which results in a single QRS-area number for each QRS complex. These numbers are resampled at 4 Hz to result in a smoother and more accurate respiratory rate. For the QRS-amplitude, the minimum value within the 100 millisecond (ms) predetermined time window is subtracted from the maximum value and the resulting output is resampled at 4 Hz to result in a smoother and more accurate respiratory rate.
The QRS respiration process is similar to the accelerometer respiration process 300 except for the computation of the respiration signal. A quality metric process of the QRS based respiration signal is generated using the same computations as the quality metric process 500. FIG. 7 illustrates a flow chart of a QRS respiration process 700 in accordance with an embodiment. The QRS respiration process 700 obtains an EKG signal via 702 and high-pass filters the EKG signal at 20 Hz via 704. The QRS respiration process 700 determines R-peaks via 706 and computes a QRS-area/QRS-amplitude respiration signal via 708. The steps 710-722 of the QRS respiration process 700 resemble steps 306-318 of the accelerometer respiration process 300.
In one embodiment, the QRS-area computation via 708 comprises obtaining the latest R-peak time tR and extracting a 100 ms time period window around the current R-peak to determine a current QRS wave. An absolute value of the QRS wave is determined per the equation and an area of the QRS wave is computed by summing the absolute values of the EKG signal within the 100 ms time period window per the equation




  
    ∑
    
      i
      =
      
        
          t
          R
        
        -
        
          50
          ⁢
          
            
          
          ⁢
          ms
        
      
    
    
      
        t
        R
      
      +
      
        50
        ⁢
        
          
        
        ⁢
        ms
      
    
  
  ⁢
  
    
  
  ⁢
  
    
      
      
        EKG
        ⁡
        
          [
          i
          ]
        
      
      
    
    .
  




The QRS-area derived respiration signal is resampled at 4 Hz to determine the QRS-area derived respiratory rate.

In one embodiment, the QRS-amplitude computation via 708 comprises the same process except that the amplitude is computed as a maximum/minimum instead of summing the values in the 100 ms time period window. Therefore, after obtaining the latest R-peak time and extracting a 100 ms time period window around the current R-peak, the maximum and minimum values in the 100 ms time period window are determined as QRSMAX and QRSMIN. The QRS-amplitude is computed as the maximum minus the minimum or QRSMAX−QRSMIN and the QRS-amplitude derived respiration signal is resample at 4 Hz to determine the QRS-amplitude derived respiratory rate.
In addition to deriving a respiration signal from a chest-mounted accelerometer and the QRS wave of an EKG signal, the respiration signal can also be derived from heart rates computed via an EKG signal. Besides changing the QRS-area/QRS-amplitude, breathing also changes the heart rate such that the heart rate increases during inspiration and decreases during expiration. This is known as the respiratory sinus arrhythmia (RSA). The R-R intervals of the EKG signal are calculated to determine another signal that is correlated to breathing as R-R intervals decrease with inspiration and increase with expiration.
In one embodiment, an RSA respiration process is utilized to determine a respiration signal from heart rates computed from an EKG signal. The RSA respiration process is similar to the QRS-area computation via 708 except that instead of computing the area of each QRS complex, the time between the current and last QRS complex is computed and then resampled at 4 Hz to determine the RSA derived respiratory rate. RSA decreases with age thereby making the RSA respiration process derived respiration signal less reliable and increasing the importance of the quality metric determination.
Respiration signals can be determined via the aforementioned accelerometer respiration process 300 (resulting in RRACCEL), the QRS respiration process 700 (resulting in RRQRS), and the RSA respiration process (resulting in RRRSA). Using the outputted quality metric for each of these respiration signals, a combined and weighted respiratory rate (RRweighted) is computed per the equation:




  
    RR
    weighed
  
  =
  
    
      
        RR
        ACCEL
      
      ⁢
      
        
          Q
          ACCEL
        
        
          
            Q
            ACCEL
          
          +
          
            Q
            QRS
          
          +
          
            Q
            RSA
          
        
      
    
    +
    
      
        RR
        QRS
      
      ⁢
      
        
          Q
          QRS
        
        
          
            Q
            ACCEL
          
          +
          
            Q
            QRS
          
          +
          
            Q
            RSA
          
        
      
    
    +
    
      
        RR
        RSA
      
      ⁢
      
        
          
            Q
            RSA
          
          
            
              Q
              ACCEL
            
            +
            
              Q
              QRS
            
            +
            
              Q
              RSA
            
          
        
        .
      
    
  




Each respiration rate is weighted by a value defined as each associated quality metric divided by the sum of all of the quality metrics. For example, if RRACCEL is unreliable due to heavy motion artifacts and noise, QACCEL will be close to 0 and that component of the RRweighted equation will not contribute much to the final respiratory rate. FIG. 8 illustrates a diagram 800 of respiration rate computations in accordance with an embodiment. In addition to the combined and weighted respiratory rate (RRweighted), the diagram 800 includes respiratory rates derived from an accelerometer, QRS-area, RSA, and a NOX RIP-band reference.
Weighting the respiration rates that are derived by the accelerometer (RRACCEL) and the EKG signals (RRQRS and RRRSA) by their respective quality metrics overcomes issues such as motion artifacts, increased respiration rates, and patient ages that can diminish respiration signal quality. For example, utilizing the RRweighted combination results in an output that weighs the RRACCEL less during the presence of motion artifacts, weighs RRQRS less during high respiration rates, and weighs RRRSA less if the patient is elderly.
QRS-area and RSA derived respiration signals cannot detect breathing rates of higher than half the heart rate (HR) because they are EKG-based. This is known as cardiac aliasing. A fundamental property of all sampled signals is that it is not possible to distinguish any frequency content that is larger than 1⁄2 the sampling rate (Fs) which is the Nyquist Frequency (NS). Any frequency content above NS cannot be reliably distinguished, and this is called aliasing. The heart rate serves as the sampling rate of the EKG-derived respiration signal and respiration rates above 1⁄2HR cannot be detected. As a result, if the respiration rate is above 1⁄2HR, then only the respiratory rate derived from the accelerometer (RRACCEL) is utilized per the equation:

IF (RRACCEL>HR/2) and QACCEL≧τACCEL. then RRweighted=RRACCEL.

As above described, the method and system allow for measuring a respiratory rate using a combination of respiration signals. A wireless sensor device determines a plurality of respiration signals from a variety of sources including but not limited to a chest-mounted accelerometer, a QRS-area of an EKG signal, a QRS-amplitude of an EKG signal, and an RSA. The wireless sensor device then computes a quality metric for each of the plurality of respiration signals, and combines the plurality of quality metrics to compute a weighted respiratory rate. The weighted respiratory rate is an accurate measurement of a person's respiratory rate that reduces the effects of quality issues including but not limited to motion artifacts and the person's age.
A method and system for measuring respiratory rate has been disclosed. Embodiments described herein can take the form of an entirely hardware implementation, an entirely software implementation, or an implementation containing both hardware and software elements. Embodiments may be implemented in software, which includes, but is not limited to, application software, firmware, resident software, microcode, etc.
The steps described herein may be implemented using any suitable controller or processor, and software application, which may be stored on any suitable storage location or computer-readable medium. The software application provides instructions that enable the processor to cause the receiver to perform the functions described herein.
Furthermore, embodiments may take the form of a computer program product accessible from a computer-usable or computer-readable storage medium providing program code or program instructions for use by or in connection with a computer or any instruction execution system. For the purposes of this description, a computer-usable or computer-readable storage medium can be any apparatus that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.
The computer-readable storage medium may be an electronic, magnetic, optical, electromagnetic, infrared, semiconductor system (or apparatus or device), or a propagation medium. Examples of a computer-readable storage medium include a semiconductor or solid state memory, magnetic tape, a removable computer diskette, a random access memory (RAM), a read-only memory (ROM), a rigid magnetic disk, and an optical disk. Current examples of optical disks include DVD, compact disk-read-only memory (CD-ROM), and compact disk-read/write (CD-RAN).
Although the present invention has been described in accordance with the embodiments shown, one of ordinary skill in the art will readily recognize that there could be variations to the embodiments and those variations would be within the spirit and scope of the present invention. Accordingly, many modifications may be made by one of ordinary skill in the art without departing from the spirit and scope of the appended claims.

## Claims

1. A system for determining respiratory rates of a user, the system comprising:
at least one electrode adapted to be coupled to the user, wherein the at least one electrode measures at least one respiration signal,
wherein to measure at least one respiration signal is to measure a plurality of respiration signals to determine a plurality of respiration rates (RR1 . . . N), wherein the application further causes the processor to:
determine a quality metric (Q1 . . . N) for each of the plurality of respiration rates; and
combine the plurality of respiration rates by weighting the quality metric for each of the plurality of respiration rates based on a quality process to calculate a weighted respiration rate (RRweighted);


a processor coupled to the at least one electrode; and
a memory device coupled to the processor, wherein the memory device stores an application which, when executed by the processor, causes the processor to:
filter the at least one respiration signal measured by the at least one electrode using a lowpass filter;
peak-pick the filtered respiration signal to determine a respiratory rate by using a standard deviation of the filtered respiration signal over a changing threshold; and
determine a quality metric of the respiratory rate by using previously learned weights of features of the respiratory rate, wherein the features of the respiratory rate include a coefficient of variation of a difference in minima times (ct), a coefficient of variation of peak-to-peak values (cp), a mean peak-to-peak value (mp), and a ratio of a number of picked extrema to a total number of all local extrema (rext).

2. The system of claim 1, wherein the at least one respiration signal is derived using any of a chest-mounted accelerometer, a QRS-area of an EKG signal, a QRS-amplitude of the EKG signal, and a respiratory sinus arrhythmia (RSA).

3. The system of claim 1, wherein the quality metric is weighted using the previously learned weights and transformed using an exponential to get the quality metric to a value between 0 and 1.

4. The system of claim 1, wherein the quality metric is derived within a predetermined time period.

5. The system of claim 4, wherein the predetermined time period is 45 seconds.

6. The system of claim 1, wherein the application further causes the processor to:
perform a smoothing of the respiratory rate to reduce effects of missed and extraneous peaks using a trimmed mean.

7. The system of claim 1, wherein the weighted respiration rate (RRweighted) comprises:


RRweighted=RR1*(Q1/(Q1+Q2+ . . . +QN))+RR2*(Q2/(Q1+Q2+ . . . +QN))+RRN*(QN/(Q1+Q2+ . . . QN)).