US20080243017A1 - Breathing sound analysis for estimation of airlow rate - Google Patents

Breathing sound analysis for estimation of airlow rate Download PDF

Info

Publication number
US20080243017A1
US20080243017A1 US11692745 US69274507A US2008243017A1 US 20080243017 A1 US20080243017 A1 US 20080243017A1 US 11692745 US11692745 US 11692745 US 69274507 A US69274507 A US 69274507A US 2008243017 A1 US2008243017 A1 US 2008243017A1
Authority
US
Grant status
Application
Patent type
Prior art keywords
breathing
sounds
apparatus according
patient
arranged
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US11692745
Inventor
Zahra Moussavi
Azadeh Yadollahi
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TR TECHNOLOGIES Inc
Original Assignee
Telecommunications Res Labs
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B7/00Instruments for auscultation
    • A61B7/003Detecting lung or respiration noise
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/087Measuring breath flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • A61B5/4818Sleep apnoea
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6814Head
    • A61B5/6815Ear
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6822Neck
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6813Specially adapted to be attached to a specific body part
    • A61B5/6825Hand
    • A61B5/6826Finger
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/683Means for maintaining contact with the body
    • A61B5/6838Clamps or clips
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/726Details of waveform analysis characterised by using transforms using Wavelet transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/0002Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0015Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
    • A61B5/002Monitoring the patient using a local or closed circuit, e.g. in a room or building
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/0002Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0015Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
    • A61B5/0022Monitoring a patient using a global network, e.g. telephone networks, internet
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Detecting, measuring or recording for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7232Signal processing specially adapted for physiological signals or for diagnostic purposes involving compression of the physiological signal, e.g. to extend the signal recording period

Abstract

Apparatus for use detection of apnea includes a microphone mounted in the ear of the patient for detecting breathing sounds and a second external microphone together with an oximetric sensor. A transmitter at the patient compresses and transmits the signals to a remote location where there is provided a detector module for receiving and analyzing the signals to extract data relating to the breathing. The detector uses the entropy or range of the signal to generate an estimate of air flow while extracting extraneous snoring and heart sounds and to analyze the estimate of air flow using Otsu's threshold to detect periods of apnea and/or hypopnea. A display provides data of the detected apnea/hypopnea episodes and related information for a clinician.

Description

  • This invention relates to an apparatus for use in breathing sound analysis for estimation of airflow rate.
  • This application is related to a co-pending Application filed on the same day as this application under Attorney Docket No. 84201-1402 and entitled BREATHING SOUND ANALYSIS FOR DETECTION OF SLEEP APNEA/HYPOPNEA EVENTS.
  • BACKGROUND OF THE INVENTION
  • Acoustical respiratory flow estimation has drawn much attention in recent years due to difficulties in airflow measurement. In clinical respiratory and/or swallowing assessment, flow is usually measured by spirometry devices, such as pneumotachograph, nasal cannulae connected to a pressure transducer, heated thermistor or anemometry. Airflow is also measured by indirect means, i.e., detection of chest and/or abdominal movements using respiratory inductance plethysmography (RIP), strain gauges, or magnetometers. The most reliable measurement of airflow is achieved by a mouth piece or facemask connected to a pneumotachograph. However, this device cannot be used during the swallowing assessment. Therefore, when recording sound during a swallow, flow is usually measured by nasal cannulae connected to a pressure transducer. Potentially, this method could be an inaccurate measure of airflow because the air leaks around the nasal cannulae. In addition, if the subject breathes through the mouth, flow is not registered at all.
  • For these reasons, the combined use of nasal cannulae connected to a pressure transducer and the measurement of respiratory inductance plethomogoraphy to monitor volume changes has been recommended as the best approach in recording flow to assess respiratory and swallowing patterns. However, application of these techniques has some disadvantages, especially when studying young children or patients with neurological impairments, where the study of swallowing is clinically important. Although the application of nasal cannulae may seem a minor intrusion, it can produce agitation in children and patients with neurological impairment. In addition, applying the RIP devices is difficult in children with neurological impairment as their poor postural control and physical deformities can make it challenging to ensure stable positioning.
  • In one of the first attempts at acoustical flow estimation, researchers attempted to estimate flow from tracheal sound by investigating eight different methods in the two categories of “reference curve” and “hierarchical clustering analysis”. The results showed a mean error between 13-15% of the measured flow for seven of the methods, with 31% for the eighth method. In the works by another group, flow estimation using either tracheal or lung sounds was achieved by investigating different models with about 90% overall accuracy over different flow rates from low to high flow rates. In these studies the exponential model between flow and average power of tracheal sound was found to be superior to other models.
  • In another study, the tracheal sound envelope was investigated for flow estimation. The tracheal sound was band-pass filtered in the range of 200-1000 Hz and then a Hilbert transform was applied to the filtered signal. The transformed signal was used to calculate the tracheal sound envelope and to estimate the flow from the calculated envelope by a linear model. The estimated flow was then used to measure ventilation, but the flow estimation error was not reported. The flow rate in that study was constant at tidal flow and half of the recorded flow signal was used to calibrate the model.
  • All of the above mentioned methods assume that at least some samples of breath sound with known flow at each flow rate were available to derive the model coefficients for flow estimation. Capturing respiratory sounds at different flow rates for calibration may not always be possible prior to assessment especially when assessing young children, patients with neurological impairments and/or patients in emergency conditions.
  • Analysis of breathing sounds from a patient for determination of sleep apnea and/or hypopnea is proposed in a paper entitled “Validation of a New System of Tracheal sound Analysis for the diagnosis of Sleep Apnea-Hypopnea Syndrome” by Nakano et al in “SLEEP” Vol27 No. 5 published in 2004. This constitutes a research paper postulating that sleep apnea can be detected by breathing sound analysis but providing no practical details for a system which may be used in practise. It is believed that no further work has been published and no commercial machine has arisen from this paper.
  • U.S. Pat. No. 6,241,683 (Macklem) issued Jun. 5th 2001 discloses a method for estimating air flow from breathing sounds where the system determines times when sounds are too low to make an accurate determination and uses an interpolation method to fill in the information in these times. Such an arrangement is of course of no value in detecting apnea or hypopnea since it accepts that the information in these times is inaccurate.
  • SUMMARY OF THE INVENTION
  • It is one object of the invention to provide an apparatus for use in analysis of breathing of a patient.
  • According to a first aspect of the invention there is provided an apparatus comprising:
  • a microphone arranged to be located on the patient for detecting breathing sounds;
  • a detector module for receiving and analyzing the signals to extract data relating to the breathing;
  • the detector module being arranged to analyze the signals to generate an estimate of air flow;
  • and a display for displaying for a clinician the estimated air flow rate relative time;
  • wherein the detector module is arranged to calculate a function representing the range of the signal or the entropy of the signal providing an estimate of air flow during breathing.
  • Preferably the detector module is arranged to cancel heart sounds from the function.
  • In one preferred method, the function is the range of the signal which is defined as the log of the difference between minimum and maximum of the signal within each short window (i.e. 100 ms) of data.
  • In another preferred method, the function is the entropy of the signal which is defined by the following formula:
  • H ( p ) = t = 1 N p i log p i ,
  • where pi is the probability distribution function of the ith event.
  • Preferably the extraneous sounds are removed by the detector module prior to flow estimation.
  • Preferably the display is arranged to display the estimated air flow versus time in any desired time length being chosen by the user.
  • Preferably the display is capable of zoom-in and zoom-out functions in the same window.
  • Preferably the display is capable of playing the breathing sounds in any data window.
  • Preferably the microphone is arranged to be located in the ear of the patient.
  • Preferably the microphone in the ear includes a transmitter arranged for wireless transmission to a receiver.
  • Preferably the apparatus is arranged such that inspiration and expiration are monitored by an initial calibration wherein the patient is instructed to initialize the system by taking a deep breath, hold it, start up of the monitoring system, then exhale and continue breathing.
  • Preferably an estimate of flow rate is calibrated using a look-up table of previously measured flow-sound relationship data that is sorted based on characteristics of the subjects.
  • Preferably the characteristics in the look-up table include BMI, gender, height, neck circumference, and smoking history of the subject.
  • Preferably the detector module is arranged to cancel heart sounds.
  • According to a second aspect of the invention there is provided an apparatus for use in use in analysis of breathing of a patient during sleep comprising:
  • a microphone arranged to be located on the patient for detecting breathing sounds;
  • a detector module for receiving and analyzing the signals to extract data relating to the breathing;
  • the detector module being arranged to analyze the signals to generate an estimate of air flow during inspiration and expiration;
  • and a display for displaying the data for a clinician;
  • wherein the apparatus is arranged such that inspiration and expiration are monitored by an initial calibration wherein the patient is instructed to initialize the system by taking a deep breath, hold it, start up of the monitoring system, then exhale and continue breathing.
  • According to a third aspect of the invention there is provided an apparatus for use in use in analysis of breathing of a patient during sleep comprising:
  • a microphone arranged to be located on the patient for detecting breathing sounds;
  • a detector module for receiving and analyzing the signals to extract data relating to the breathing;
  • the detector module being arranged to analyze the signals to generate an estimate of air flow during inspiration and expiration;
  • and a display for displaying the data for a clinician;
  • wherein an estimate of flow rate is calibrated using a look-up table of previously measured flow-sound relationship data that is sorted based on characteristics of the subjects.
  • This proposal aims to develop a prototype of an integrated system to acquire, de-noise, analyze the tracheal respiratory sounds, estimate airflow acoustically.
  • Long distance monitoring and diagnostic aid tools provide large financial saving to both the health care system and families. This proposal will provide a novel system to both developing a new and yet simple diagnostic tool for sleep apnea disorder, and also a new way to connect the specialists and physicians with patients either in remote areas or even at their homes. Aside from its obvious benefit for covering the remote areas with equal opportunity for health care, it also reduces the long waiting list for sleep studies. From a public health perspective, non-invasive and inexpensive methods to determine airway responses across all ages and conditions would present a major step forward in the management of sleep apnea disorders.
  • Long distance monitoring and diagnostic aid tools provide large financial savings to both the health care system and the patient's families. From a public health perspective, non-invasive and inexpensive methods to determine airway responses across all ages and conditions would present a major step forward in the management of sleep apnea disorders.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • One embodiment of the invention will now be described in conjunction with the accompanying drawings in which:
  • FIG. 1 is a schematic illustration of a sleep apnea detection apparatus according to the present invention.
  • FIG. 2 is an illustration of a typical screen displaying the data to the physician.
  • FIG. 3( a) is a graphical representation of Tracheal sound entropy;
  • FIG. 3( b) is a graphical representation of entropy after applying nonlinear median filter (star marks represents the estimated apnea segments)
  • FIG. 3( c) is a graphical representation of flow signal (solid line) along with the estimated (dotted line) and real (dashed line) apnea segments for a typical subject.
  • FIG. 4 is a graphical representation of Mean and standard deviation values of errors in estimating apnea periods for different subjects
  • FIG. 5 is a block diagram illustrating the adaptive filtering scheme for removing the snoring sounds from the signal using the signal recorded by the auxiliary microphone in the vicinity of the patient.
  • DETAILED DESCRIPTION
  • One of the reasons to record many signals in a sleep study is the inaccuracy of those recorded signals in sleep apnea detection when they are used as a single measure. For example, nasal cannulae are used to measure airflow; however, when the patient breathes through the mouth, the nasal cannulae register nothing and hence give a false positive detection error for apnea. Therefore, combination of nasal pressure plus thermistor and End-tidal carbon dioxide concentration in the expired air (ETCO2) is used to have a qualitative measure of respiratory airflow. The abdominal movement recordings are mainly used to detect respiratory effort and hence to distinguish between central and obstructive sleep apnea. The ECG signals are also being used for detecting heart rate variability and another measure for apnea detection as well as monitoring patient's heart condition during the night. The combination of EOG (Electrooculogram), EEG (Electroencephalogram), and EMG (Electromyogram) signals are used for assessing the rapid eye movement (REM) sleep stage that is characterized by desynchronization of the EEG and loss of muscle tone. Recording these signals are necessary if insight in sleep quality is saught for diagnosis of certain sleep disorders. The most important information that doctors seek from a complete sleep study is the duration and frequency of apnea and/or hypopnea and the blood's Oxygen saturation (SaO2) level of the patient during the apnea. Oxygen level usually drops during the apnea and will rise quickly with awakenings. However, oximetry alone does not detect all cases of sleep apnea.
  • As the first and most important information of a sleep study is an accurate measure of duration and frequency of apnea during sleep, the present arrangement provides a fully automated system to detect apnea with only one single sensor that can also easily be applied by the patient at home and detect apnea acoustically; hence reducing the need for a complete laboratory sleep study.
  • The apparatus provides an integrated system for remote and local monitoring and assessment of sleep apnea as a diagnostic aid for physicians and allows the following:
  • To record the SaO2 data simultaneously with respiratory sound signals through either a neck band with a microphone mounted in a chamber placed over the supra-sternal notch, or by a microphone from inside the ear, followed by a signal conditioning unit.
  • To screen the raw data, separate snoring and other adventitious sounds from breath sounds, estimate flow from the sounds and detect apnea and/or hypopnea episodes, determine the duration and frequency of the apnea episodes and finally display the raw data, estimated airflow and display the estimated airflow with marked detected apnea/hypopnea along with related information (duration, frequency and the corresponded SaO2 data).
  • FIG. 1 shows the apparatus for sleep apnea detection that can also be used as a home-care device while being connected to a clinical diagnostic center for online monitoring.
  • From a public health perspective, non-invasive and inexpensive methods to determine airway responses across all ages and conditions would present a major step forward in the management of sleep apnea disorders
  • The apparatus consists of six modules that permit sleep apnea detection diagnosis. The clinical diagnosis can be performed either locally (e.g. at a clinical diagnostic center) and/or remotely (e.g. at clinician's office/home). The apparatus will support several clinicians simultaneously carrying out clinical work on different patients. Likewise, patients can be monitored either locally (e.g. at a clinical diagnostic center) and/or remotely (e.g. at patient's home). The apparatus will also support many patients being concurrently monitored.
  • The apparatus has the following modules
  • Collector module 10,
  • Transmitter module 11,
  • Organizer module 12,
  • Detector module 13,
  • Interface module 14, and
  • Manager module 15.
  • Collector module 10 captures physiological signals from different body parts. The body parts include a microphone and transmitter 20 at the ear or over the neck by a wireless microphone mounted in chamber with a neckband for recording sounds, a sensor 22 at the fingers of the patient for recording oximetry data and an external microphone 21 for recording sound from the environment around the patient. Other signals can be detected in some cases from other body parts if the physicians request other biological signals, but this is not generally intended herein. The collector module locally transfers wirelessly the signals to the Transmitter module 11.
  • Transmitter module 11 receives biological signals from the Collector module 10, securely transmits those signals and receives the signals at the diagnostic center for its delivery to the Organizer module 12.
  • The Transmitter module 11 consists of two components; The Transmitter Sender (S) and the Transmitter Receiver (R). The Transmitter Sender together with the Collector module resides at the patient location. The Transmitter Sender receives and store temporally signals from the Collector, and securely and reliably transfers the signals to the Transmitter Receiver. The Transmitter Receiver resides at the diagnostic center location. The Transmitter Receiver securely and reliably accepts the signals from the Transmitter Sender, and forwards the signals to the Organizer for the signal management and processing. There is one pair of collector—transmitter modules per patient being monitored.
  • Inter-Transmitter components signal transmission can occurred locally for those cases when the Collector-Transmitter Sender resides in the same center (e.g. at a diagnostic facility) or remotely for those cases when the Collector-Transmitter Sender resides externally (e.g. at a patient home). The transmission can be wireless or wired (e.g. through the internet/intranet).
  • Organizer module 12 receives all captured signals from the Transmitter module, organizes and classifies received signals per patient/physician and prepares the signals for its processing by the Detector module. The Organizer module simultaneously supports receiving many signals from different patients that is signals from collector—transmitter module pairs.
  • Detector module 13 pre-processes and analyzes the patient biological signals, and performs the sleep apnea detection. The Detector performs snoring sound detection and separation prior to the apnea/hypopnea detection. The Detector has self-calibrated acoustical respiratory airflow estimation and phase detection utilized in respiratory and sleep apnea assessments
  • Interface module 14 provides the graphical user interface to the clinicians. The Interface module gives a secure, reliable, user-friendly, interactive access to the analysis performed by the Detector and it is organized per patient physician. The Interface module consists of two main components: the Interface Master (M) and the Interface Client (C). The Interface Master serves the information to the Interface Client(s), while the Interface Client provides the access to the clinicians, Several Interface Clients can run concurrently giving out results to several clinicians. The Interface Client can be executed locally (e.g. intranet) or externally (e.g. internet).
  • Manager module 15 provides the application management functions. It provides the graphical user interface to the application administrator at the diagnostic center location.
  • All system/application parameters are setup at the Manager module. The system/application parameters configure the apparatus for its proper operation.
  • The collector module microphone may comprise a neck band with a microphone mounted in a chamber placed over the supra-sternal notch. However the preferred arrangement as shown in FIG. 1 schematically comprises a wireless microphone inside the ear or by a microphone mounted in a chamber with a neck band to record respiratory sounds followed by a suitable signal conditioning unit depending on the type of the used sensor. The second sensor 21 collects sound from the environment around. The third sensor 22 collects the conventional SaO2 data or other oximetry data. The three sensors allow from the patient simultaneous data acquisition of the sound signals and the SaO2 data.
  • There are two options for recording respiratory sounds: using the ear microphone or the neck microphone. The very small miniature ear microphone is inserted into a piece of foam which has open ends and inserted to inside the microphone. The small preamplifier of the microphone is placed behind the ear similar to a hearing aid device, The ear microphone includes a wireless transmitter which is placed behind the ear, the miniature microphone and the foam for securing the microphone inside the ear. In case of neck microphone, it is inserted in a chamber (with the size of a loony) which allows about 2 mm distance between the microphone and the skin when the chamber is placed over supra-sternal notch of the trachea of the patient with double sided adhesive ring tapes. The neck microphone will come with a neck band mainly for the comfort of the patient and also to keep the wire of the microphone free of touching the skin. In either case, the preamplifier and transmitter of the wireless microphone can be placed in the pocket of the subject. Alternatively the whole element mounted in the ear canal includes the pre-amplifier and transmitter for complete wireless operation.
  • The detector module pre-processes and analyzes the recorded signal in order to provide a user friendly, smart and interactive interface for the physician as a monitoring and diagnostic aid tool. The software in this part will de-noise the recorded sound, separate snoring sounds, estimate the flow acoustically, detect apnea and/or hypopnea episodes, count the duration and the frequency of their occurrence, display the estimated flow with marked apnea episodes as shown in FIG. 3 along with the related information.
  • The respiratory sounds either from the ear or from the neck of the patient will be recorded by a small wireless microphone. A Transmitter—Sender Module DSP board is designed to receive the analog signal, amplify and filter the signal, digitize it with a minimum of 5120 Hz sampling rate and store it as a binary file.
  • The SaO2 data simultaneously with the respiratory sounds is digitized with 5120 Hz sampling rate and stored in a binary file for the entire duration of the sleep at the collector module.
  • The detector module signal processing of the sound signals has three stages. First an automated algorithm finds the artifacts (that normally appear as impulses in the signal) and removes them from further analyses. Secondly, the snoring sounds, if they exist, are identified and separated from the respiratory sounds. Finally, from the cleaned respiratory sounds the entropy of the signal is calculated, the effect of heart sounds is removed, and apnea episodes are detected by the technique as described hereinafter. The average duration of the apnea episodes, their frequency of occurrence and whether they are associated with snoring, is presented as part of the information in the GUI interface for the physician.
  • Artifacts (usually due to movement) appear as very short duration pulses in the recorded signal. Wavelet analysis is a highly reliable method with high accuracy to automatically detect these artifacts. On the other hand, snoring sounds are musical sounds which appear with harmonic components in the spectrogram of the recorded signal. Detection of snoring sounds is similar to detection of crackle sounds in the lung sounds. Multi-scale product of the wavelet coefficients is used to detect and separate the snoring sounds. Techniques for the application of digital signal processing techniques on biological signals including noise and adventitious sounds separation are known.
  • Once the respiratory sound signal is pre-processed and cleaned of extra sounds, the entropy of the signal is calculated. As heart sounds have overlap with respiratory sounds at low frequencies and this is more pronounced at very low flow rate (the case of hypopnea), the effect of heart sounds has to be cancelled from the entropy or the range parameter of the signals prior to apnea detection. This is described in more detail hereinafter.
  • Then, from the entropy or the range parameter of the signal, the apnea episodes are identified using Otsu's thresholding method as described hereinafter.
  • The flow estimation method as described hereinafter is enhanced to make the method self-calibrated. That enables the apparatus to estimate the actual amount of flow. Finally, the episodes of hypopnea and apnea are marked; their duration and frequency of occurrence during the entire sleep is presented on the interface module GUI display as a diagnostic aid to the physician.
  • Depending on the type of microphone used, both sounds result in the same apnea detection episodes and flow estimation while the tuning of the algorithm for each sound signal requires slight modification, i.e. the threshold or the parameters of the flow estimation model are different.
  • The apnea detection algorithm requires a snoring separation algorithm. This can use one or more of the following principles:
  • Applying Wavelet analysis to detect and mark the snoring sounds in the time-frequency domain.
  • Applying the adaptive filter cancellation technique to remove the snoring sounds from the signal using the signal recorded by the auxiliary microphone in the vicinity of the patient.
  • An automated algorithm can be provided to clean the recorded breath sound signal from all extra plausible noises such as cough sounds, swallowing sounds, vocal noise (in case the patient talks while dreaming), and artifacts due to movements following the apnea detection algorithm on the cleaned signal and validate the results. These extraneous sounds will be removed using wavelet analysis for localization and several different filter banks to remove each type of noises either automatically or at the users command.
  • Display
  • The interface module 14 provides a display of the detected apnea/hypopnea episodes and related information for a clinician. The display includes a display 30 of airflow versus time is plotted with apnea and hypopnea episodes marked on the screen.
  • The display includes oximetry data 31 plotted in association with the estimated airflow.
  • The display has touch screen controls 32, 33 providing zoom-in and zoom-out functions in the same window for both airflow and oximetry data simultaneously.
  • The display is capable of playing the breathing and snoring sounds in any zoomed-in or zoomed-out data window, that is the sounds are stored to allow an actual rendition of those sounds to the clinician to study the sounds at or around an apnea event, The display is capable of displaying the extracted information about the frequency and duration of apnea/hypopnea episodes, and their association with the level of oximetry data in a separate window for the clinician.
  • Apnea Detection
  • Referring now to FIGS. 4( a), 4(b) and 4(c), further detail of the Sleep Apnea detection components is now described.
  • In order to smooth the calculated entropy or range parameter, it is segmented into windows of 200 ms with 50% overlap between adjacent windows. Each window was then presented by its median value which is not sensitive to jerky fluctuation of the signal.
  • Next, the smoothed entropy or the smoothed range signal is classified into two groups of breathing and apnea using a nonparametric and unsupervised method for automatic threshold selection using the principles of OTSU.
  • In Otsu's method the threshold is chosen such that the variance between classes is maximized. The between-class variance is defined as the sum of variances of all classes respect to the total mean value of all classes:
  • σ B 2 = w 0 ( μ 0 - μ T ) 2 + w 1 ( μ 1 - μ T ) 2 , ( 1 )
  • where w1i (i=1,2) are the probability and mean values of the classes, respectively and μT is the average of total values.
  • and the optimum threshold k″ is selected so as:
  • σ B 2 ( k * ) = max 1 k < L σ B 2 ( k ) . ( 2 )
  • The average of entropy or range values is another statistical measure that can be used to detect apnea segments. In this study both the Otsu and the average value of entropy or range value were used to define the classification threshold as:

  • Thr=min{k″,m},  (3)
  • where k″ is the Otsu threshold and m is the average of the entropy or range values.
  • FIG. 4 presents (a) Tracheal sound entropy, (b) entropy after applying nonlinear median filter (star marks represents the estimated apnea segments) and c) flow signal (solid line) along with the estimated (doffed line) and real (dashed line) apnea segments for a typical subject. Comparing the results depicted in FIG. 4( a) and FIG. 4( b), the effect of applying median filter is evident. The star marks in FIG. 4( b) show the estimated apnea segments. Investigating the results depicted in FIG. 4( c) it is clear that the proposed method detects all the apnea segments and classifies them correctly from the breath segments.
  • In this arrangement a new acoustical method for apnea detection is proposed which is based on tracheal sound entropy or range value. The method is fast and easy to be implemented, which makes it suitable for on-line applications.
  • Removal of snoring sounds by time-frequency filtering techniques may have some problems due to the fact that snoring sounds also have strong low frequency components, in which the acoustical apnea detection is based on. As an alternative, the snoring sounds can be recorded by another auxiliary microphone in the vicinity of the subject. This signal will not have breathing sounds and can be used as a noise reference. The apparatus then uses adaptive filtering for noise (snore) cancellation.
  • Snoring sounds are musical sounds which appear with harmonic components in the spectrogram of the recorded signal. We record the snoring sounds by an auxiliary microphone in the vicinity of the patient. Using the source of noise (recorded by the auxiliary microphone) adaptive filtering will cancel the snoring sounds from the breath and snoring sounds recorded over the neck or inside the ear of the patient.
  • FIG. 5 illustrates the block diagram of the adaptive filtering scheme. The filter has two inputs, the primary input and the reference signal. The primary input, x(t), (the microphone over the neck or inside the ear) contains an interference, m(n), (snoring sounds) along with the information bearing signal, b(n), (tracheal sound). The reference input, r(n), (the auxiliary microphone) represents a version of interference with undetectable information bearing (tracheal sounds) signal. The output of the RLS FIR filter, y(n), is close to the interference component of the primary signal. Therefore, the output of the adaptive filter, e(n), is the minimum mean square error estimate of the information bearing signal, {circumflex over (b)}(n).
  • Computational demand of the smart, automated algorithm to run 8 hours of sleep data can be high. As an alternative, the algorithms are written in C++ code that increases the speed of the algorithms compared to a high level signal processing software such as MATLAB. With fast, state-of-art new computers, this will not be a problem considering that this system will replace the 4 hours labor work of an sleep lab technician (the usual time to analyze one PSG patient's data) with a few minutes of processing time.
  • Flow Estimation Flow Calibration
  • In order to provide an effective flow estimation method it desirable to provide another sound channel recorded over the lung to be used for respiratory phase detection and second it is desirable to provide one breath with known flow from the patient to calibrate (tune) the model to that patient. This calibration (tuning) is necessary because there is a wide variation of flow-sound relationship between the subjects due to their different chest size, lung capacity, gender, age, etc.
  • Respiratory Phase Detection
  • Thus the present arrangement provides a method of respiratory phase detection with only one channel breath sound (Tracheal sound signal).
  • In this method the patient is required to have a deep breath, hold it, start the program and then exhale and keep breathing normally but with different flow rates from low to high for 30 seconds, This 30 second data that starts with expiration phase is used by the program to derive the necessary information for phase detection of the rest of breath sounds. The phase detection algorithm is:
  • 1. Sequester the 30 second initialization data into 100 ms segments with 50% overlap between the successive segments.
  • 2. Calculate the average power (in dB) of each segment over the range of 150-450 Hz. The valleys of the resultant signal, which looks like a rectified sinusoid, determine the onsets of the breaths.
  • 3 Knowing that the first phase is the exhalation, label the initialization data as inspiration/expiration phases. Also by comparing the max power in each phase, label them as low, tidal and high flow rates.
  • 4. Calculate the mean value of the average power (this time calculated over the range of 500-1200 Hz of each segment) of the top 20% of each phase and store it for inspiration and expiration phases separately.
  • 5. Calculate the ratio of the mean of the average power calculated in Step 4 between the inspiration and expiration phases.
  • 6. Apply this ratio as a threshold to the rest of the data to determine respiratory phases. For example, if the ratio of inspiration and expiration is calculated as 1.2, and the ratio of any known phase respect to the adjacent phase (calculated with the same method) is equal to 0.8, it means that the first known phase is expiration and the second one is inspiration.
  • Automatic Self Calibration
  • Since having one breath with known flow defeats the purpose of eliminating the flow measurement, in this arrangement is provided a method of automatic self calibration using a data bank. The concept includes a very large data bank of breathing sounds (tracheal sound) of people. This data bank is sorted based on body-mass-index (BMI), age, gender, and smoking history of the subjects. This data is used to match the patient's BMI and other information to suggest the known flow-sound relationship required for calibration.
  • De-Noising and Adventitious Sound Removal
  • Since the patient might have some respiratory diseases that may cause some adventitious sounds, i.e., crackle sounds or wheezes, an algorithm is required to be run by the choice of the user (the clinician) to remove all adventitious sounds prior to flow estimation.
  • This algorithm has two parts: adventitious sound localization and removal. For adventitious localization the arrangement herein uses multi-scale (level 3) product of wavelet coefficients and applies a running threshold of mean plus three times of standard deviation to detect and localize the adventitious sounds. Then, the segments including artefacts will be removed in time-frequency domain, the signal will be interpolated by spline interpolation and the breath sound signal will be reconstructed in time domain by taking the inverse of the spectrogram.
  • Flow Estimation Using Entropy or Range Parameter
  • 1. Band-pass filter the tracheal sounds in the frequency range of 75 to 600 Hz and normalize the signal.
  • 2. Sequester the band-pass filtered signal into segments of 50 ms (512 samples) with 75% overlap between the successive segments.
  • 3. Let x(t) be the signal in each segment. The range value in each segment can be defined as:

  • L r=log[mean(x|x>[max(x)*(1−r/100])−mean(x|x<[max(x)*r/100])]  (1)
  • where x is the tracheal sound signal in each window and mean( ) is the average value, and r=1,
    or:

  • Lr=log[var(x)],  (2)
  • where var(x) is the variance of the signal in each segment.
  • 4. The other feature that can be used for flow estimation is the entropy of the signal in each segment. Let {X1, . . . , XN} represent the values of signal x in each segment. Estimate the probability density function (pdf) of signal x(t), p(x), in each window using the Normal kernel estimator:
  • p ^ k ( x ) = 1 N t = 1 N 1 h K ( x - X i h ) , ( 3 )
  • where N is the number of samples (205), K is the Gaussian kernel function (K(x)=(2π)−1/2 exp(−x2/2)) and h is the kernel bandwidth. For Gaussian kernel the optimum h is approximated as:

  • h opt=1.06{circumflex over (σ)}(x)N −0.2  (4)
  • where {circumflex over (σ)}(x) is the estimated standard deviation of the signal x(t) in each window. Calculate the Shannon entropy in each segment:
  • L r = - t = 1 N p t log ( p i ) . ( 5 )
  • 5. Use the modified linear model (6) to estimate flow from tracheal sounds entropy or range (Eq. 1, 2 or 5) feature;
  • F est = C 1 ( mean ( L ph ) mean ( L base ) ) L ph + C 2 , ( 6 )
  • where C1 and C2 are the model coefficients derived from the one breath with known flow, Lph=[L1, . . . , Lw] is a vector representing the entropy or range value of the signal in the upper 40% values of each respiratory phase (inspiration or expiration), w is the number of segments in the upper 40% values of each respiratory phase and Li is the entropy or range values of tracheal sound in each segment (Eq. 1, 2 or 5). Similarly, Lbase is the same vector that is calculated using the base respiratory phase signal. Base respiratory phase is the one breath that is assumed to be available with known flow to calibrate the model.
  • Heart Sounds Localization
  • 1. Band-pass filter the tracheal sound records in the range of 75-2500 Hz to remove motion artifacts and high-frequency noises.
  • 2. Divide the filtered signal into segments of 20 ms (205 samples) with 50% overlap between successive segments.
  • 3. Let x(t) be the signal in each segment. The range value in each segment can be defined as:

  • L r=log[mean(x|x>[max(x)*(1−r/100])−mean(x|x<[max(x)*r/100])]  (7)
  • where x is the tracheal sound signal in each window and mean( ) is the average value, and r=1,
  • or:

  • Lr=log[var(x)],  (8)
  • where var(x) is the variance of the signal in each segment.
  • The other feature that can be used for heart sounds localization is entropy of the signal. Let {X1, . . . , XN} represent the values of signal x in each segment. Estimate the probability density function (pdf) of signal x(t), p(x), in each window using the Normal kernel estimator:
  • p ^ k ( x ) = 1 N i = 1 N 1 h K ( x - X i h ) , ( 9 )
  • where N is the number of samples (205) K is the Gaussian kernel function (K(x)=(2π)−1/2 exp(−x2/2)) and h is the kernel bandwidth. For Gaussian kernel the optimum h is approximated as:

  • h opt=1.06{circumflex over (σ)}(x)N −0.2.  (10)
  • where {circumflex over (σ)}(x) is the estimated standard deviation of the signal x(t) in each window.
  • 4. Calculate the Shannon entropy in each segment:
  • H ( p ) = - i = 1 N p i log ( p i ) . ( 11 )
  • 5. Define average plus standard deviation value (μ+σ) of the calculated entropy or range value as the threshold for heart sounds localization.
  • 6. Mark the segments with entropy or range values of higher than this threshold as heart sounds-included segments.
  • Removing the Effects of Heart Sounds
  • 1. Localize heart sounds with the method mentioned above.
  • 2. Calculate the range or entropy values for the segments void of heart sounds.
  • 3. Apply spline interpolation to estimate the values of the entropy or range value in the segments including heart sounds. This technique effectively cancels the effect of heart sounds on the entropy or range values of the tracheal sound.
  • EXAMPLE 1
  • Eight healthy subjects (3 males) aged 33.1±6.6 years with body mass index of 23.3±3.5 participated in this study. Tracheal sound was recorded using Siemens accelerometer (EMT25C) placed over supra-sternal notch using double adhesive tapes. Respiratory flow signal was measured by a mouth piece pneumotachograph (Fleisch No. 3) connected to a differential pressure transducer (Validyne, Northridge, Calif.). The subjects were instructed to breathe at very shallow flow rates with different periods of breath hold (2, 4, 6 sec) to simulate apnea. Tracheal sound and flow signals were recorded and digitized simultaneously at a 10240 Hz sampling rate.
  • Feature Extraction
  • Among several features of tracheal sound such as the sound's mean amplitude, average power and entropy used for flow estimation, entropy and the range of signal have been shown to be the best features following flow variation. Therefore, in this study tracheal sounds entropy was used to detect apnea (breath hold in the experiments of this study) without the use of the measured flow signal. However, the recorded flow signal was used for validation of the acoustically detected apnea.
  • Tracheal sound signal was band-pass filtered in the range of [75-600] Hz, and then segmented into 50 ms (512 samples) windows with 50% overlap between the adjacent windows. In each window the tracheal sound probability density function (pdf) was estimated based on kernel methods. Then, using the method described earlier in this document Shannon entropy was calculated in each window that represents the changes in the signal's pdf. The effect of heart sounds which is most evident in the frequency range below 200 Hz was removed by the method introduced earlier in this document.
  • FIG. 3 shows the calculated entropy and its corresponding flow signal for a typical subject. By comparing the signals depicted in FIG. 1( a) and FIG. 1( c) (solid line), it is evident that the values of entropy in the breath hold segments are smaller than those during breathing.
  • It should be noted that when localizing the segments including heart sounds, it is nearly impossible to find out the exact boundaries of heart sounds segments. Therefore, there is always a trade off between the amount of heart sounds interference in respiratory sounds and the amount of respiratory sounds information missing during heart sounds cancellation. The high peaks in the calculated entropy (FIG. 3 a) are related to the heart sounds components remained in the tracheal sound. FIG. 4 displays the mean and standard deviation values of length and lag errors in estimating apnea periods for different subjects.
  • EXAMPLE 2
  • In this study 10 healthy subjects of the previous participated. Subjects were in two age groups: 5 adults (all female) 29±8 years old and 5 children (3 female) 9.6±1.7 years old. Respiratory sounds were recorded using Siemens accelerometers (EMT25C) placed over supra-sternal notch and the upper right lobe lung. Respiratory flow was measured by a pneumotachograph (Fleisch No. 3) connected to a differential pressure transducer (Validyne, Northridge, Calif.). Subjects were instructed to breathe at 5 different flow rates with 5 breaths at each target flow followed by a 10s of breath hold at the end of experiment. In this study the shallow (<6 ml/s/kg), low (6-9 ml/s/kg), medium (12-18 ml/s/kg), high (18-27 ml/s/kg) and very high (>27 ml/s/kg) target flow rates were investigated. Tracheal sound signals were used for flow estimation while the lung sound signal in correspondence with tracheal sound signals were used for respiratory phase detection. The onsets of breaths were detected by running a threshold on the average power of the tracheal sounds and detecting the valleys of the signal. Since lung sounds are much louder during inspiration as opposed to expiration, then by comparing the average power of the lung and tracheal breath sounds it can easily and accurately be determined which phases are inspiration or expiration.
  • As described above, the best performance for estimating flow from tracheal sound entropy was achieved in the frequency range of 75-600 Hz. Tracheal sound signals were used for flow estimation while the lung sound signal in correspondence with tracheal sound signal were used for respiratory phase detection as mentioned above.
  • As described above, the best performance for estimating flow from tracheal sound entropy was achieved in the frequency range of [75-600] Hz. This is in accordance with the fact that the main energy components of tracheal sound exists in the frequency range below 600-800 Hz. Thus, tracheal sound was band-pass filtered in this range followed by segmenting the band-pass filtered signal into segments of 50 ms (512 samples) with 75% overlap between the successive segments.
  • When studying tracheal sound in the frequency range below 300 Hz, heart sounds are the main source of interference that changes the time and frequency characteristics of the tracheal sound. Therefore, the presence of heart sounds will cause an error which can become significant in flow estimation in very shallow breathing, when most of the signal's energy is concentrated at low frequencies. Hence, in this study the effect of heart sounds on the extracted parameters was cancelled by using the same method as described above.
  • Since various modifications can be made in my invention as herein above described, and many apparently widely different embodiments of same made within the spirit and scope of the claims without department from such spirit and scope, it is intended that all matter contained in the accompanying specification shall be interpreted as illustrative only and not in a limiting sense.

Claims (19)

  1. 1. Apparatus for use in analysis of breathing of a patient comprising:
    a microphone arranged to be located on the patient for detecting breathing sounds;
    a detector module for receiving and analyzing the signals to extract data relating to the breathing;
    the detector module being arranged to analyze the signals to generate an estimate of air flow;
    and a display for displaying for a clinician the estimated air flow rate relative time;
    wherein the detector module is arranged to calculate a function representing the range of the signal or the entropy of the signal providing an estimate of air flow during breathing.
  2. 2. The apparatus according to claim 1 wherein the detector module is arranged to cancel heart sounds from the function.
  3. 3. The apparatus according to claim 1 wherein the function using entropy is defined by the following formula.
    H ( p ) = i = 1 N p t log p i ,
    where p, is the probability distribution function of the ith event.
  4. 4. The apparatus according to claim 1 wherein the function is the range of signal which is defined as the log of the difference between minimum and maximum of the signal within each short window (i.e. 100 ms) of data.
  5. 5. The apparatus according to claim 1 wherein the extraneous sounds are removed by the detector module prior to flow estimation.
  6. 6. The apparatus according to claim 1 wherein the display is arranged to display the estimated air flow versus time in any desired time length being chosen by the user.
  7. 7. The apparatus according to claim 1 wherein the display is capable of zoom-in and zoom-out functions in the same window.
  8. 8. The apparatus according to claim 16 wherein the display is capable of playing the breathing sounds in any data window.
  9. 9. The apparatus according to claim 1 wherein the microphone is arranged to be located in the ear of the patient.
  10. 10. The apparatus according to claim 9 wherein the microphone in the ear includes a transmitter arranged for wireless transmission to a receiver.
  11. 11. The apparatus according to claim 1 wherein the apparatus is arranged such that inspiration and expiration are monitored by an initial calibration wherein the patient is instructed to initialize the system by taking a deep breath, hold it, start up of the monitoring system, then exhale and continue breathing.
  12. 12. The apparatus according to claim 1 wherein an estimate of flow rate is calibrated using a look-up table of previously measured flow-sound relationship data that is sorted based on characteristics of the subjects.
  13. 13. The apparatus according to claim 1 wherein the characteristics in the look-up table include BMI, gender, height, neck circumference, and smoking history of the subject.
  14. 14. The apparatus according to claim 1 wherein the detector module is arranged to cancel heart sounds.
  15. 15. Apparatus for use in use in analysis of breathing of a patient during sleep comprising:
    a microphone arranged to be located on the patient for detecting breathing sounds;
    a detector module for receiving and analyzing the signals to extract data relating to the breathing;
    the detector module being arranged to analyze the signals to generate an estimate of air flow during inspiration and expiration:
    and a display for displaying the data for a clinician;
    wherein the apparatus is arranged such that inspiration and expiration are monitored by an initial calibration wherein the patient is instructed to initialize the system by taking a deep breath, hold it, start up of the monitoring system, then exhale and continue breathing.
  16. 16. The apparatus according to claim 15 wherein an estimate of flow rate is calibrated using a look-up table of previously measured flow-sound relationship data that is sorted based on characteristics of the subjects.
  17. 17. The apparatus according to claim 16 wherein the characteristics in the look-up table include BMI, gender, height, neck circumference, and smoking history of the subject.
  18. 18. Apparatus for use in use in analysis of breathing of a patient during sleep comprising:
    a microphone arranged to be located on the patient for detecting breathing sounds;
    a detector module for receiving and analyzing the signals to extract data relating to the breathing;
    the detector module being arranged to analyze the signals to generate an estimate of air flow during inspiration and expiration;
    and a display for displaying the data for a clinician;
    wherein an estimate of flow rate is calibrated using a look-up table of previously measured flow-sound relationship data that is sorted based on characteristics of the subjects.
  19. 19. The apparatus according to claim 18 wherein the characteristics in the look-up table include BMI, gender, height, neck circumference, and smoking history of the subject.
US11692745 2007-03-28 2007-03-28 Breathing sound analysis for estimation of airlow rate Abandoned US20080243017A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US11692729 US7559903B2 (en) 2007-03-28 2007-03-28 Breathing sound analysis for detection of sleep apnea/popnea events
US11692745 US20080243017A1 (en) 2007-03-28 2007-03-28 Breathing sound analysis for estimation of airlow rate

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US11692729 US7559903B2 (en) 2007-03-28 2007-03-28 Breathing sound analysis for detection of sleep apnea/popnea events
US11692745 US20080243017A1 (en) 2007-03-28 2007-03-28 Breathing sound analysis for estimation of airlow rate
PCT/CA2008/000581 WO2008116318A8 (en) 2007-03-28 2008-03-27 Breathing sound analysis for estimation of air flow rate

Publications (1)

Publication Number Publication Date
US20080243017A1 true true US20080243017A1 (en) 2008-10-02

Family

ID=39788005

Family Applications (2)

Application Number Title Priority Date Filing Date
US11692729 Active 2027-05-01 US7559903B2 (en) 2007-03-28 2007-03-28 Breathing sound analysis for detection of sleep apnea/popnea events
US11692745 Abandoned US20080243017A1 (en) 2007-03-28 2007-03-28 Breathing sound analysis for estimation of airlow rate

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US11692729 Active 2027-05-01 US7559903B2 (en) 2007-03-28 2007-03-28 Breathing sound analysis for detection of sleep apnea/popnea events

Country Status (2)

Country Link
US (2) US7559903B2 (en)
WO (1) WO2008116318A8 (en)

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080310662A1 (en) * 2007-06-15 2008-12-18 Davidson Terence M Earpiece snoring sound transmitter
US20100210962A1 (en) * 2009-02-13 2010-08-19 Jingping Xu Respiratory signal detection and time domain signal processing method and system
US20110015702A1 (en) * 2009-07-15 2011-01-20 Ternes David J Remote sensing in an implantable medical device
US20110015704A1 (en) * 2009-07-15 2011-01-20 Ternes David J Physiological vibration detection in an implanted medical device
US20110015703A1 (en) * 2009-07-15 2011-01-20 Ternes David J Remote pace detection in an implantable medical device
US20110071378A1 (en) * 2009-09-24 2011-03-24 Nellcor Puritan Bennett Llc Signal Processing Warping Technique
US20110092839A1 (en) * 2008-11-17 2011-04-21 Toronto Rehabilitation Institute Mask and method for use in respiratory monitoring and diagnostics
US20110301485A1 (en) * 2010-06-04 2011-12-08 Yongji Fu Data binning method and system for estimating respiratory airflow from body sound signal
US20120271872A1 (en) * 2009-06-21 2012-10-25 Tom Chau Procedure for denoising dual-axis swallowing accelerometry signals
US20130102908A1 (en) * 2010-03-31 2013-04-25 Nanyang Technological University Air Conduction Sensor and a System and a Method for Monitoring a Health Condition
US20130253357A1 (en) * 2010-05-24 2013-09-26 University Of Manitoba System and methods of acoustical screening for obstructive sleep apnea during wakefulness
US8579794B2 (en) 2008-05-02 2013-11-12 Dymedix Corporation Agitator to stimulate the central nervous system
US20130324877A1 (en) * 2012-05-30 2013-12-05 National University Corporation Chiba University Apparatus for determining respiratory condition
US20140066724A1 (en) * 2011-02-18 2014-03-06 Matias Zanartu System and Methods for Evaluating Vocal Function Using an Impedance-Based Inverse Filtering of Neck Surface Acceleration
US20140249852A1 (en) * 2013-03-04 2014-09-04 Hello Inc. Methods using patient monitoring devices with unique patient IDs and a telemetry system
US8834347B2 (en) 2008-08-22 2014-09-16 Dymedix Corporation Anti-habituating sleep therapy for a closed loop neuromodulator
US20140276174A1 (en) * 2013-03-15 2014-09-18 Elwha Llc Correlating a condition of a subject with a degree of sleep apnea being experienced by the subject
WO2014153246A2 (en) * 2013-03-14 2014-09-25 Aliphcom Sleep management implementing a wearable data-capable device for snoring-related conditions and other sleep disturbances
US20140330095A1 (en) * 2011-09-07 2014-11-06 University Of Manitoba System and methods for estimating respiratory airflow
US20150051501A1 (en) * 2011-06-03 2015-02-19 Brian M. Dugan Bands for measuring biometric information
US20150126888A1 (en) * 2012-05-10 2015-05-07 University Of Washington Through Its Center For Commercialization Sound-based spirometric devices, systems and methods
US9232910B2 (en) 2008-11-17 2016-01-12 University Health Network Method and apparatus for monitoring breathing cycle by frequency analysis of an acoustic data stream
WO2016049285A1 (en) * 2014-09-25 2016-03-31 Aedio, Inc. Systems and methods for digital predictive disease exacerbation and pre-emptive treatment
US9649087B2 (en) 2011-05-17 2017-05-16 University Health Network Method and device for apnea and hypopnea detection
US9801590B2 (en) 2011-05-17 2017-10-31 University Health Network Breathing disorder identification, characterization and diagnosis methods, devices and systems
US9949667B2 (en) 2008-11-17 2018-04-24 University Health Network Mask and method for use in respiratory monitoring and diagnostics

Families Citing this family (46)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5915380A (en) 1997-03-14 1999-06-29 Nellcor Puritan Bennett Incorporated System and method for controlling the start up of a patient ventilator
US20060155202A1 (en) * 2005-01-12 2006-07-13 Arand Patricia A Hemodynamic assessment/adjustment
US20060243280A1 (en) * 2005-04-27 2006-11-02 Caro Richard G Method of determining lung condition indicators
WO2006117780A3 (en) 2005-04-29 2007-09-07 Oren Gavriely Cough detector
US9779751B2 (en) * 2005-12-28 2017-10-03 Breath Research, Inc. Respiratory biofeedback devices, systems, and methods
WO2007079068A3 (en) * 2005-12-28 2007-12-13 Nirinjan Bikko Breathing biofeedback device
US20120179066A1 (en) * 2006-10-18 2012-07-12 Yuan Ze University Sleeping quality monitor system and a method for monitoring a physiological signal
US8021310B2 (en) 2006-04-21 2011-09-20 Nellcor Puritan Bennett Llc Work of breathing display for a ventilation system
US7784461B2 (en) 2006-09-26 2010-08-31 Nellcor Puritan Bennett Llc Three-dimensional waveform display for a breathing assistance system
US9986926B2 (en) * 2007-10-26 2018-06-05 Inovise Medical, Inc. Q-onset ventricular depolarization detection in the presence of a pacemaker
US8185190B2 (en) * 2008-01-29 2012-05-22 Inovise Medical, Inc. Assessment of ischemia, and risk of sudden cardiac death, via heart-functionality parameter and acoustic cardiographic monitoring
US8348852B2 (en) * 2008-03-06 2013-01-08 Inovise Medical, Inc. Heart-activity sound monitoring
DE102008016855A1 (en) * 2008-04-02 2009-10-15 Pfeifer, Thomas, Dr. Apparatus and method for influencing the sleep of a person, the control unit and computer program product
US8577448B2 (en) * 2008-10-14 2013-11-05 Inovise Medical, Inc. Differential apneic detection in aid of diagnosis and treatment
FR2937873B1 (en) * 2008-11-03 2012-08-24 Assistance Publique Hopitaux Paris disagreement quantification system between a patient on a respirator and a corresponding support device.
US20100256460A1 (en) * 2009-04-03 2010-10-07 The General Electric Company Wearable Monitoring System
US20100262031A1 (en) * 2009-04-14 2010-10-14 Yongji Fu Method and system for respiratory phase classification using explicit labeling with label verification
WO2011010384A1 (en) * 2009-07-24 2011-01-27 富士通株式会社 Sleep apnea syndrome examination device and program
EP2477544A4 (en) * 2009-09-14 2017-11-01 Sleep Methods System and method for anticipating the onset of an obstructive sleep apnea event
US9409013B2 (en) 2009-10-20 2016-08-09 Nyxoah SA Method for controlling energy delivery as a function of degree of coupling
US9415216B2 (en) 2009-10-20 2016-08-16 Nyxoah SA Devices for treatment of sleep apnea
US8409108B2 (en) * 2009-11-05 2013-04-02 Inovise Medical, Inc. Multi-axial heart sounds and murmur detection for hemodynamic-condition assessment
US8335992B2 (en) 2009-12-04 2012-12-18 Nellcor Puritan Bennett Llc Visual indication of settings changes on a ventilator graphical user interface
US8924878B2 (en) 2009-12-04 2014-12-30 Covidien Lp Display and access to settings on a ventilator graphical user interface
US9119925B2 (en) 2009-12-04 2015-09-01 Covidien Lp Quick initiation of respiratory support via a ventilator user interface
US9262588B2 (en) 2009-12-18 2016-02-16 Covidien Lp Display of respiratory data graphs on a ventilator graphical user interface
US8499252B2 (en) 2009-12-18 2013-07-30 Covidien Lp Display of respiratory data graphs on a ventilator graphical user interface
EP2533693B1 (en) * 2010-02-12 2018-08-22 Koninklijke Philips N.V. Method and apparatus for processing a cyclic physiological signal
US20110230778A1 (en) * 2010-03-18 2011-09-22 Yungkai Kyle Lai Methods and devices for continual respiratory monitoring using adaptive windowing
US20110295138A1 (en) * 2010-05-26 2011-12-01 Yungkai Kyle Lai Method and system for reliable inspiration-to-expiration ratio extraction from acoustic physiological signal
US9844336B2 (en) 2010-08-26 2017-12-19 Ben Gurion University Of The Negev Research And Development Authority Apparatus and method for diagnosing obstructive sleep apnea
CN103167828B (en) 2010-10-20 2016-02-24 皇家飞利浦电子股份有限公司 Method and apparatus for patients with obstructive sleep apnea diagnosis conscious of suspended
JP5605204B2 (en) * 2010-12-15 2014-10-15 ソニー株式会社 Respiratory signal processing apparatus and the processing method, and program
US8482418B1 (en) 2011-02-18 2013-07-09 Pursuit Enterprises Method and apparatus for monitoring and treatment of sleep-related conditions
CN104066476A (en) 2011-09-30 2014-09-24 尼科索亚股份有限公司 Apparatus and method to control an implant
WO2013057637A1 (en) * 2011-10-20 2013-04-25 Koninklijke Philips Electronics N.V. A system and method for characterizing an upper airway using speech characteristics
JP6310401B2 (en) * 2012-03-01 2018-04-11 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method of processing a signal representative of the physiological rhythm, system and computer program
JP6019659B2 (en) * 2012-03-27 2016-11-02 富士通株式会社 Apnea episode determination device, apnea episode determination method, and apnea episode determination program
US9814438B2 (en) 2012-06-18 2017-11-14 Breath Research, Inc. Methods and apparatus for performing dynamic respiratory classification and tracking
US8548588B1 (en) 2012-09-21 2013-10-01 Inovise Medical, Inc. CRM-device ventricular-pacing blanking control
RU2015129091A (en) 2012-12-17 2017-01-23 Конинклейке Филипс Н.В. diagnosis of sleep apnea and method for generating information using a non-intrusive audio analysis
CN104994906B (en) 2012-12-19 2017-05-24 维斯卡迪亚公司 Improve hemodynamic function by stimulating the diaphragm asymptomatic
EP3005189A1 (en) * 2013-05-31 2016-04-13 Koninklijke Philips N.V. System and method for transferring a group of related files as one logical unit
CN103735267B (en) * 2014-01-02 2016-03-30 上海大学 An apparatus snore based on screening osahs
US9950129B2 (en) 2014-10-27 2018-04-24 Covidien Lp Ventilation triggering using change-point detection
WO2017135899A1 (en) * 2016-02-03 2017-08-10 Nanyang Technological University Methods for detecting a sleep disorder and sleep disorder detection devices

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5797852A (en) * 1993-02-04 1998-08-25 Local Silence, Inc. Sleep apnea screening and/or detecting apparatus and method
US6241683B1 (en) * 1998-02-20 2001-06-05 INSTITUT DE RECHERCHES CLINIQUES DE MONTRéAL (IRCM) Phonospirometry for non-invasive monitoring of respiration
US20020165462A1 (en) * 2000-12-29 2002-11-07 Westbrook Philip R. Sleep apnea risk evaluation
US6852084B1 (en) * 2000-04-28 2005-02-08 Peter V. Boesen Wireless physiological pressure sensor and transmitter with capability of short range radio frequency transmissions
US20060064037A1 (en) * 2004-09-22 2006-03-23 Shalon Ventures Research, Llc Systems and methods for monitoring and modifying behavior

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5265624A (en) * 1990-09-06 1993-11-30 Edentec Stimulation collar
US5123425A (en) * 1990-09-06 1992-06-23 Edentec Obstructive sleep apnea collar
DE4138702A1 (en) * 1991-03-22 1992-09-24 Madaus Medizin Elektronik Process and apparatus for diagnosis and quantitative analysis of apnea, and for simultaneous determination of other diseases
US5844996A (en) 1993-02-04 1998-12-01 Sleep Solutions, Inc. Active electronic noise suppression system and method for reducing snoring noise
DE19538473A1 (en) * 1995-10-16 1997-04-17 Map Gmbh Device and method for the quantitative analysis of sleep disorders
US6213955B1 (en) * 1998-10-08 2001-04-10 Sleep Solutions, Inc. Apparatus and method for breath monitoring
US6290654B1 (en) * 1998-10-08 2001-09-18 Sleep Solutions, Inc. Obstructive sleep apnea detection apparatus and method using pattern recognition
US6171258B1 (en) * 1998-10-08 2001-01-09 Sleep Solutions, Inc. Multi-channel self-contained apparatus and method for diagnosis of sleep disorders
US6415033B1 (en) * 1999-09-15 2002-07-02 Ilife Systems, Inc. Physiological condition monitors utilizing very low frequency acoustic signals
US6575916B2 (en) * 2000-03-24 2003-06-10 Ilife Solutions, Inc. Apparatus and method for detecting very low frequency acoustic signals
US6706002B1 (en) * 2000-03-24 2004-03-16 Ilife Systems, Inc. System and method for remotely monitoring at least one physiological characteristic of a child
US6666830B1 (en) * 2000-08-17 2003-12-23 East River Ventures, Lp System and method for detecting the onset of an obstructive sleep apnea event
US6935335B1 (en) * 2000-08-17 2005-08-30 Ilife Systems, Inc. System and method for treating obstructive sleep apnea
US6517497B2 (en) * 2000-12-13 2003-02-11 Ge Medical Systems Information Technologies, Inc. Method and apparatus for monitoring respiration using signals from a piezoelectric sensor mounted on a substrate
US20050197588A1 (en) * 2004-03-04 2005-09-08 Scott Freeberg Sleep disordered breathing alert system
US7747323B2 (en) * 2004-06-08 2010-06-29 Cardiac Pacemakers, Inc. Adaptive baroreflex stimulation therapy for disordered breathing
US20070118054A1 (en) * 2005-11-01 2007-05-24 Earlysense Ltd. Methods and systems for monitoring patients for clinical episodes

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5797852A (en) * 1993-02-04 1998-08-25 Local Silence, Inc. Sleep apnea screening and/or detecting apparatus and method
US6241683B1 (en) * 1998-02-20 2001-06-05 INSTITUT DE RECHERCHES CLINIQUES DE MONTRéAL (IRCM) Phonospirometry for non-invasive monitoring of respiration
US6852084B1 (en) * 2000-04-28 2005-02-08 Peter V. Boesen Wireless physiological pressure sensor and transmitter with capability of short range radio frequency transmissions
US20020165462A1 (en) * 2000-12-29 2002-11-07 Westbrook Philip R. Sleep apnea risk evaluation
US20060064037A1 (en) * 2004-09-22 2006-03-23 Shalon Ventures Research, Llc Systems and methods for monitoring and modifying behavior

Cited By (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080310662A1 (en) * 2007-06-15 2008-12-18 Davidson Terence M Earpiece snoring sound transmitter
US8579794B2 (en) 2008-05-02 2013-11-12 Dymedix Corporation Agitator to stimulate the central nervous system
US8834346B2 (en) 2008-08-22 2014-09-16 Dymedix Corporation Stimulus sequencer for a closed loop neuromodulator
US8834347B2 (en) 2008-08-22 2014-09-16 Dymedix Corporation Anti-habituating sleep therapy for a closed loop neuromodulator
US9949667B2 (en) 2008-11-17 2018-04-24 University Health Network Mask and method for use in respiratory monitoring and diagnostics
US9232910B2 (en) 2008-11-17 2016-01-12 University Health Network Method and apparatus for monitoring breathing cycle by frequency analysis of an acoustic data stream
US20110092839A1 (en) * 2008-11-17 2011-04-21 Toronto Rehabilitation Institute Mask and method for use in respiratory monitoring and diagnostics
US20100210962A1 (en) * 2009-02-13 2010-08-19 Jingping Xu Respiratory signal detection and time domain signal processing method and system
US8992446B2 (en) * 2009-06-21 2015-03-31 Holland Bloorview Kids Rehabilitation Hospital Procedure for denoising dual-axis swallowing accelerometry signals
US20120271872A1 (en) * 2009-06-21 2012-10-25 Tom Chau Procedure for denoising dual-axis swallowing accelerometry signals
US8515534B2 (en) 2009-07-15 2013-08-20 Cardiac Pacemakers, Inc. Remote pace detection in an implantable medical device
US9867987B2 (en) 2009-07-15 2018-01-16 Cardiac Pacemakers, Inc. Remote detection of physiological parameter using an implantable medical device
US8301241B2 (en) 2009-07-15 2012-10-30 Cardiac Pacemakers, Inc. Remote pace detection in an implantable medical device
US9381358B2 (en) 2009-07-15 2016-07-05 Cardiac Pacemakers, Inc. Remote detection of physiological parameter using an implantable medical device
US8285373B2 (en) 2009-07-15 2012-10-09 Cardiac Pacemakers, Inc. Remote sensing in an implantable medical device
US8588906B2 (en) 2009-07-15 2013-11-19 Cardiac Pacemakers, Inc. Physiological vibration detection in an implanted medical device
US20110015703A1 (en) * 2009-07-15 2011-01-20 Ternes David J Remote pace detection in an implantable medical device
US20110015704A1 (en) * 2009-07-15 2011-01-20 Ternes David J Physiological vibration detection in an implanted medical device
US8626281B2 (en) 2009-07-15 2014-01-07 Cardiac Pacemakers, Inc. Remote sensing in an implantable medical device
US8942799B2 (en) 2009-07-15 2015-01-27 Cardiac Pacemakers, Inc. Remote pace detection in an implantable medical device
US8700146B2 (en) 2009-07-15 2014-04-15 Cardiac Pacemakers, Inc. Remote pace detection in an implantable medical device
US8934970B2 (en) 2009-07-15 2015-01-13 Cardiac Pacemakers, Inc. Physiological vibration detection in an implanted medical device
US20110015702A1 (en) * 2009-07-15 2011-01-20 Ternes David J Remote sensing in an implantable medical device
US8862226B2 (en) 2009-07-15 2014-10-14 Cardiac Pacemakers, Inc. Remote sensing in an implantable medical device
US20110071378A1 (en) * 2009-09-24 2011-03-24 Nellcor Puritan Bennett Llc Signal Processing Warping Technique
US8840562B2 (en) * 2009-09-24 2014-09-23 Covidien Lp Signal processing warping technique
US9826955B2 (en) * 2010-03-31 2017-11-28 Nanyang Technological University Air conduction sensor and a system and a method for monitoring a health condition
US20130102908A1 (en) * 2010-03-31 2013-04-25 Nanyang Technological University Air Conduction Sensor and a System and a Method for Monitoring a Health Condition
US20130253357A1 (en) * 2010-05-24 2013-09-26 University Of Manitoba System and methods of acoustical screening for obstructive sleep apnea during wakefulness
US9931073B2 (en) * 2010-05-24 2018-04-03 University Of Manitoba System and methods of acoustical screening for obstructive sleep apnea during wakefulness
US8603007B2 (en) * 2010-06-04 2013-12-10 Sharp Laboratories Of America, Inc. Data binning method and system for estimating respiratory airflow from body sound signal
US20110301485A1 (en) * 2010-06-04 2011-12-08 Yongji Fu Data binning method and system for estimating respiratory airflow from body sound signal
US20140066724A1 (en) * 2011-02-18 2014-03-06 Matias Zanartu System and Methods for Evaluating Vocal Function Using an Impedance-Based Inverse Filtering of Neck Surface Acceleration
US9649087B2 (en) 2011-05-17 2017-05-16 University Health Network Method and device for apnea and hypopnea detection
US9801590B2 (en) 2011-05-17 2017-10-31 University Health Network Breathing disorder identification, characterization and diagnosis methods, devices and systems
US20150051501A1 (en) * 2011-06-03 2015-02-19 Brian M. Dugan Bands for measuring biometric information
US9974481B2 (en) * 2011-06-03 2018-05-22 Brian M. Dugan Bands for measuring biometric information
US20140330095A1 (en) * 2011-09-07 2014-11-06 University Of Manitoba System and methods for estimating respiratory airflow
US10004452B2 (en) * 2011-09-07 2018-06-26 University Of Manitoba System and methods for estimating respiratory airflow
US20150126888A1 (en) * 2012-05-10 2015-05-07 University Of Washington Through Its Center For Commercialization Sound-based spirometric devices, systems and methods
US10028675B2 (en) * 2012-05-10 2018-07-24 University Of Washington Through Its Center For Commercialization Sound-based spirometric devices, systems and methods
US20130324877A1 (en) * 2012-05-30 2013-12-05 National University Corporation Chiba University Apparatus for determining respiratory condition
US9848831B2 (en) * 2012-05-30 2017-12-26 Nihon Kohden Corporation Apparatus for determining respiratory condition
US9298882B2 (en) * 2013-03-04 2016-03-29 Hello Inc. Methods using patient monitoring devices with unique patient IDs and a telemetry system
US20140249852A1 (en) * 2013-03-04 2014-09-04 Hello Inc. Methods using patient monitoring devices with unique patient IDs and a telemetry system
WO2014153246A2 (en) * 2013-03-14 2014-09-25 Aliphcom Sleep management implementing a wearable data-capable device for snoring-related conditions and other sleep disturbances
WO2014153246A3 (en) * 2013-03-14 2014-11-13 Aliphcom Sleep management implementing a wearable data-capable device for snoring-related conditions and other sleep disturbances
US20140276174A1 (en) * 2013-03-15 2014-09-18 Elwha Llc Correlating a condition of a subject with a degree of sleep apnea being experienced by the subject
US9707121B2 (en) 2013-03-15 2017-07-18 Elwha Llc Treating sleep apnea with negative pressure
WO2016049285A1 (en) * 2014-09-25 2016-03-31 Aedio, Inc. Systems and methods for digital predictive disease exacerbation and pre-emptive treatment

Also Published As

Publication number Publication date Type
WO2008116318A8 (en) 2009-07-16 application
US7559903B2 (en) 2009-07-14 grant
WO2008116318A1 (en) 2008-10-02 application
US20080243014A1 (en) 2008-10-02 application

Similar Documents

Publication Publication Date Title
US7314451B2 (en) Techniques for prediction and monitoring of clinical episodes
Abeyratne et al. Pitch jump probability measures for the analysis of snoring sounds in apnea
Hilton et al. Evaluation of frequency and time-frequency spectral analysis of heart rate variability as a diagnostic marker of the sleep apnoea syndrome
US20020169384A1 (en) Method and device for sleep monitoring
US20080275349A1 (en) Monitoring, predicting and treating clinical episodes
Taplidou et al. Wheeze detection based on time-frequency analysis of breath sounds
US20090024044A1 (en) Data recording for patient status analysis
US6142952A (en) Method and apparatus for detection and diagnosis of airway obstruction degree
US6918878B2 (en) Methods and systems for monitoring respiration
US7640055B2 (en) Self-adaptive system for the analysis of biomedical signals of a patient
US20130144190A1 (en) Sleep apnea detection system
US6415174B1 (en) ECG derived respiratory rhythms for improved diagnosis of sleep apnea
US6091973A (en) Monitoring the occurrence of apneic and hypopneic arousals
US5301679A (en) Method and system for analysis of body sounds
US7785262B2 (en) Method and apparatus for diagnosing respiratory disorders and determining the degree of exacerbations
US20050119586A1 (en) Systems and methods for respiratory event detection
US7207948B2 (en) Systems and methods for monitoring cough
US5058600A (en) Graphical readout of laryngotracheal spectra and airway monitor
US20110288431A1 (en) Method and apparatus for monitoring breathing cycle by frequency analysis of an acoustic data stream
US20060198533A1 (en) Method and system for continuous monitoring and diagnosis of body sounds
US20140180154A1 (en) Non-invasive monitoring of respiratory rate, heart rate and apnea
Fenton et al. Automated spectral characterization of wheezing in asthmatic children
US6241683B1 (en) Phonospirometry for non-invasive monitoring of respiration
Behar et al. A review of current sleep screening applications for smartphones
US20060111635A1 (en) Sleep staging based on cardio-respiratory signals

Legal Events

Date Code Title Description
AS Assignment

Owner name: TELECOMMUNICATIONS RESEARCH LABORATORY, CANADA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:MOUSSAVI, ZAHRA;YADOLLAHI, AZADEH;REEL/FRAME:019402/0016

Effective date: 20070518

AS Assignment

Owner name: TR TECHNOLOGIES INC., CANADA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:TELECOMMUNICATIONS RESEARCH LABORATORIES;REEL/FRAME:022841/0374

Effective date: 20090610