Systems and methods for filtering autocorrelation peaks and detecting harmonics

• US 10,537,289 B2
• Filed: 02/23/2018
• Issued: 01/21/2020
• Est. Priority Date: 02/20/2014
• Status: Active Grant

• Current Assignee: Covidien LP (Medtronic Plc)
• Original Assignee: Covidien LP (Medtronic Plc)
• Inventors: Ochs, James, McGonigle, Scott, Addison, Paul
• Primary Examiner(s): Flory, Christopher A
• Application Number:
  US15/903,821
• Publication Number:
  US 20180177460A1
• Time in Patent Office: 697 Days
• Field of Search: 600476, 600529
• CPC Class Codes: A61B 5/02416 using photoplethysmograph signals, e.g. generated by infrared radiation A61B5/14552 takes precedence
• CPC Class Codes: A61B 5/0816 Measuring devices for examining respiratory frequency measuring frequency of electric signals G01R23/00
• CPC Class Codes: A61B 5/14551 for measuring blood gases
• CPC Class Codes: A61B 5/7203 for noise prevention, reduction or removal
• CPC Class Codes: A61B 5/7246 using correlation, e.g. template matching or determination of similarity
• CPC Class Codes: A61B 5/7278 Artificial waveform generation or derivation, e.g. synthesising signals from measured signals A61B5/327 takes precedence

Abstract

Systems and methods are provided for determining respiration information from physiological signals such as PPG signals. A physiological signal is processed to generate at least one respiration information signal and an autocorrelation sequence is generated based on the at least one respiration information signal. In some embodiments, a respiration peak is identified based on the autocorrelation sequence and a composite peak is generated based on the identified peak and at least one previous respiration peak. Respiration information is calculated based on the composite peak. In some embodiments, a determination is made whether the autocorrelation sequence includes an undesired harmonic. When the autocorrelation sequence includes an undesired harmonic, the autocorrelation sequence may not be used in the calculation of respiration information.

Claims

1. A method comprising:

receiving a photoplethysmograph (PPG) signal;

processing, with processing equipment, plural windows of the PPG signal to generate plural respiration information signals;

generating, with the processing equipment, plural autocorrelation sequences based on the plural respiration information signals;

identifying first and second respiration peak segments of respective first and second autocorrelation sequences where the first and second autocorrelation sequences are derived from different data windows of said PPG signal;

determining, with the processing equipment, whether the first and second autocorrelation sequences include an undesired harmonic; and

calculating, with the processing equipment, respiration information based at least in part on the first and second autocorrelation sequences when it is determined that the autocorrelation sequences do not include an undesired harmonic.

2. The method of claim 1, further comprising calculating, with the processing equipment, a metric based on the autocorrelation sequences, wherein the metric is indicative of whether a particular pattern of peaks is present in the autocorrelation sequences and wherein the determination of whether the autocorrelation sequences includes an undesired harmonic is based on the metric.

3. The method of claim 2, wherein the particular pattern comprises a repeating sequence of small and large peaks.

4. The method of claim 1, further comprising:

identifying, with the processing equipment, a plurality of local maxima in the autocorrelation sequences;

calculating, with the processing equipment, a plurality of amplitudes associated with the plurality of local maxima; and

comparing, with the processing equipment, the plurality of amplitudes, wherein determining whether the autocorrelation sequences includes an undesired harmonic is based on the comparison of the plurality of amplitudes.

5. The method of claim 4, wherein the plurality of amplitudes comprises a sequence of amplitudes, and wherein comparing the plurality of amplitudes comprises determining whether subsequent amplitudes in the sequence are greater than the multiplication of the previous amplitudes in the sequence and a coefficient.

6. The method of claim 1, wherein determining whether the autocorrelation sequence includes an undesired harmonic is based on a threshold, the method further comprising modifying the threshold based on whether a previous autocorrelation sequence was determined to include an undesired harmonic.

7. The method of claim 1, further comprising not using an autocorrelation sequence in the calculation of the respiration information when it is determined that the autocorrelation sequence includes an undesired harmonic.

8. The method of claim 1, wherein calculating the respiration information based at least in part on autocorrelation sequences comprises:

identifying a respiration peaks of the autocorrelation sequences;

generating a composite peak based on a respiration peak and at least one previous respiration peak; and

calculating respiration information based at least in part on the composite peak.

9. The method of claim 8, wherein calculating the respiration information based at least in part on the composite peak comprises:

identifying a local maxima of the composite peak that exceeds a threshold; and

calculating the respiration information based at least in part on a lag index associated with the local maxima.

10. The method of claim 1, wherein the respiration information comprises respiration rate.

11. A system comprising:

an input for receiving a photoplethysmograph (PPG) signal; and

processing equipment configured for;

receiving a photoplethysmograph (PPG) signal;

processing plural windows of the PPG signal to generate to generate plural respiration information signals;

generating, with the processing equipment, plural autocorrelation sequences based on the plural respiration information signals;

identifying first and second respiration peak segments of respective first and second autocorrelation sequences where the first and second autocorrelation sequences are derived from different data windows of said PPG signal;

determining, with the processing equipment, whether the first and second autocorrelation sequences include an undesired harmonic; and

calculating, with the processing equipment, respiration information based at least in part on the first and second autocorrelation sequences when it is determined that the autocorrelation sequences do not include an undesired harmonic.

12. The system of claim 11, wherein the processing equipment is further configured for calculating a metric based on the autocorrelation sequences, wherein the metric is indicative of whether a particular pattern of peaks is present in the autocorrelation sequences and wherein the determination of whether the autocorrelation sequences includes an undesired harmonic is based on the metric.

13. The system of claim 12, wherein the particular pattern comprises a repeating sequence of small and large peaks.

14. The method of claim 11, wherein the processing equipment is further configured for:

identifying a plurality of local maxima in the autocorrelation sequences;

calculating a plurality of amplitudes associated with the plurality of local maxima; and

comparing the plurality of amplitudes, wherein determining whether the autocorrelation sequences include an undesired harmonic is based on the comparison of the plurality of amplitudes.

15. The system of claim 14, wherein the plurality of amplitudes comprises a sequence of amplitudes, and wherein comparing the plurality of amplitudes comprises determining whether subsequent amplitudes in the sequence are greater than the multiplication of the previous amplitudes in the sequence and a coefficient.

16. The system of claim 11, wherein determining whether the autocorrelation sequences include an undesired harmonic is based on a threshold, and wherein the processing equipment is further configured for modifying the threshold based on whether a previous autocorrelation sequence was determined to include an undesired harmonic.

17. The system of claim 11, wherein the processing equipment is further configured for not using the autocorrelation sequences in the calculation of the respiration information when it is determined that an autocorrelation sequence includes an undesired harmonic.

18. The system of claim 11, wherein the respiration information based at least in part on the autocorrelation sequences comprises:

identifying plural respiration peaks of the autocorrelation sequence;

generating a composite peak based on a respiration peak and at least one previous respiration peak; and

calculating respiration information based at least in part on the composite peak.

19. The system of claim 18, wherein calculating the respiration information based at least in part on the composite peak comprises:

identifying a local maxima of the composite peak that exceeds a threshold; and

calculating the respiration information based at least in part on a lag index associated with the local maxima.

20. The system of claim 11, wherein the respiration information comprises respiration rate.