IP Library Granted Patent US 11,304,624
Granted Patent B2
US 11,304,624 · App. 16/196,987 · Granted Apr 19, 2022

Method and apparatus for performing dynamic respiratory classification and analysis for detecting wheeze particles and sources

Inventors: Charalampos-Christos Stamatopoulos (Athens, GR); Nirinjan Bikko Yee (Walnut Creek, CA)
Assignee: AireHealth Inc.
A61B5/0823A61B5/024A61B5/0205A61B5/0816A61B5/0826A61B5/6823A61B5/725A61B5/7278A61B5/7465A61B7/003A61B7/04G06N3/08G10L25/66A61B5/6819A61B5/6831A61B5/742G10L25/06G10L25/18G10L25/30G16H50/20
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 11,304,624
App. No.
16/196,987
Granted
Apr 19, 2022
Kind
B2
Abstract

A method for analyzing an audio respiratory signal comprises capturing the audio respiratory signal from a subject using a microphone and partitioning the audio respiratory signal into a plurality of overlapping frames. The method further comprises calculating a fourier transform for each frame and determining a magnitude spectrum using the fourier transform of the plurality of overlapping frames. Additionally, the method comprises extracting a spectrogram using the magnitude spectrum and analyzing the spectrogram to determine characteristics pertaining to wheeze sounds in the audio respiratory signal.

Claims (67)

1. A method for analyzing an audio respiratory signal, the method comprising:

capturing the audio respiratory signal from a subject using a microphone;

partitioning the audio respiratory signal into a plurality of overlapping frames and a plurality of non-overlapping frames;

calculating a fourier transform for each frame;

determining a magnitude spectrum using the fourier transform of the plurality of overlapping frames;

extracting a spectrogram using the magnitude spectrum;

analyzing the spectrogram to determine characteristics pertaining to wheeze sounds in the audio respiratory signal;

creating an artificial crackle comprising a filtered impulse response frame created by filtering a delta function;

calculating a cross correlation function between each of the non-overlapping frames and the artificial crackle;

determining if the cross correlation function exceeds a predetermined threshold; and

responsive to a determination that the cross correlation function exceeds the predetermined threshold, identifying a respective frame for which the cross correlation function exceeds the predetermined threshold as associated with crackling.

2. The method of claim 1 , wherein the artificial crackle is created by filtering the delta function in the filtered impulse response frame with a narrow IIR bandpass filter.

3. The method of claim 1 , further comprising:

storing a timestamp and intensity corresponding to each frame identified as associated with crackling; and

extracting a plurality of characteristics pertaining to crackling using frames identified as associated with crackling.

4. The method of claim 1 , further comprising:

extracting a magnified spectrogram using the spectrogram, wherein the magnified spectrogram is operable to determine a number of wheeze particles.

5. The method of claim 1 , further comprising:

extracting an interpolated spectrogram using the spectrogram, wherein the interpolated spectrogram is operable to provide a clearer view of wheeze particles by subtracting out residual airflow energies from the spectrogram.

6. The method of claim 1 , wherein the analyzing further comprises:

executing an edge-detection procedure on the spectrogram to highlight wheeze particles in the spectrogram.

7. The method of claim 6 , wherein the analyzing further comprises:

for each column in the spectrogram, performing the following:

(a) store locations of elements with high values resulting from the edge-detection procedure into a vector;

(b) compare a location of the first element in the vector with locations of remaining elements in the vector;

(c) designate any of the remaining elements where a respective location is a multiple of the location of the first element as a harmonic of the first element; and

(d) repeat steps a), b) and c) for each element in the vector.

8. A non-transitory computer-readable storage medium having stored thereon, computer executable instructions that, if executed by a computer system cause the computer system to perform a method for analyzing an audio respiratory signal, the method comprising:

capturing the audio respiratory signal from a subject using a microphone;

partitioning the audio respiratory signal into a plurality of overlapping frames and a plurality of non-overlapping frames;

calculating a fourier transform for each frame;

extracting a spectrogram using the fourier transform of the plurality of overlapping frames;

analyzing the spectrogram to determine characteristics pertaining to wheeze sounds in the audio respiratory signal;

creating an artificial crackle comprising a filtered impulse response frame created by filtering a delta function;

calculating a cross correlation function between each of the non-overlapping frames and the artificial crackle;

determining if the cross correlation function exceeds a predetermined threshold; and

responsive to a determination that the cross correlation function exceeds the predetermined threshold, identifying a respective frame for which the cross correlation function exceeds the predetermined threshold as associated with crackling.

9. The non-transitory computer-readable storage medium of claim 8 , wherein the artificial crackle is created by filtering the delta function in the filtered impulse response frame with a narrow IIR bandpass filter.

10. The non-transitory computer-readable storage medium of claim 8 , wherein the method further comprises:

storing a timestamp and intensity corresponding to each frame identified as associated with crackling; and

extracting a plurality of characteristics pertaining to crackling using frames identified as associated with crackling.

11. The non-transitory computer-readable storage medium of claim 8 , wherein the method further comprises:

extracting a magnified spectrogram using the spectrogram, wherein the magnified spectrogram is operable to determine a number of wheeze particles.

12. The non-transitory computer-readable storage medium of claim 8 , wherein the method further comprises:

extracting an interpolated spectrogram using the spectrogram, wherein the interpolated spectrogram is operable to provide a clearer view of wheeze particles by subtracting out residual airflow energies from the spectrogram.

13. A system for detecting wheeze and crackle from an audio respiratory signal, the system comprising:

a spirometer comprising a first microphone, wherein the first microphone is operable to capture the audio respiratory signal from a subject;

a memory coupled to the spirometer and operable to store the audio respiratory signal, wherein the memory further comprises an application for detecting wheeze and crackle from a breathing session stored therein; and

a processor coupled to said memory and said spirometer, the processor being configured to operate in accordance with said application to:

capture the audio respiratory signal from a subject using a microphone;

partition the audio respiratory signal into a plurality of overlapping frames and a plurality of non-overlapping frames;

calculate a fourier transform for each frame;

determine a magnitude spectrum using the fourier transform of the plurality of overlapping frames;

extract a spectrogram using the magnitude spectrum;

analyze the spectrogram to determine characteristics pertaining to wheeze sounds in the audio respiratory signal;

create an artificial crackle comprising a filtered impulse response frame created by filtering a delta function;

calculate a cross correlation function between each of the non-overlapping frames and the artificial crackle;

determine if the cross correlation function exceeds a predetermined threshold; and

responsive to a determination that the cross correlation function exceeds the predetermined threshold, identify a respective frame for which the cross correlation function exceeds the predetermined threshold as associated with crackling.

14. The system of claim 13 , wherein the artificial crackle is created by filtering the delta function in the filtered impulse response frame with a narrow IIR bandpass filter.

15. The system of claim 13 , wherein the processor is further configured to:

store a timestamp and intensity corresponding to each frame identified as associated with crackling; and

extract a plurality of characteristics pertaining to crackling using frames identified as associated with crackling.

16. The system of claim 13 , wherein the processor is further configured to:

extract a magnified spectrogram using the spectrogram, wherein the magnified spectrogram is operable to determine a number of wheeze particles.

17. The system of claim 13 , wherein the processor is further configured to:

extract an interpolated spectrogram using the spectrogram, wherein the interpolated spectrogram is operable to provide a clearer view of wheeze particles by subtracting out residual airflow energies from the spectrogram.

Assignments (5)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jun 5, 2025
From: VUAANT, INC. D/B/A CARE.AI
To: STRYKER CORPORATION
Reel/Frame 071329/0108 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Oct 25, 2023
From: AIREHEALTH, INC.
To: VUAANT, INC.
Reel/Frame 065339/0296 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Oct 9, 2023
From: BREATHRESEARCH, INC.
To: AIREHEALTH INC.
Reel/Frame 065159/0684 →
MERGER Recorded Oct 4, 2023
From: BREATHRESEARCH, INC.
To: AIREHEALTH, INC.
Reel/Frame 065118/0550 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Nov 20, 2018
From: STAMATOPOULOS, CHARALAMPOS-CHRISTOS; YEE, NIRINJAN BIKKO
To: BREATH RESEARCH, INC.
Reel/Frame 047557/0534 →
Continuity (4)
Continuation In Part 15641262 · Jul 4, 2017
Continuation In Part 13920655 · Jun 18, 2013
Provisional Application 61661267 · Jun 18, 2012
Related Publication 20190192047A1 · Jun 27, 2019
Cited By (1)
US 12,527,493