US20100145210A1 - Multi parametric classification of cardiovascular sounds - Google Patents

Multi parametric classification of cardiovascular sounds Download PDF

Info

Publication number
US20100145210A1
US20100145210A1 US12/308,755 US30875509A US2010145210A1 US 20100145210 A1 US20100145210 A1 US 20100145210A1 US 30875509 A US30875509 A US 30875509A US 2010145210 A1 US2010145210 A1 US 2010145210A1
Authority
US
United States
Prior art keywords
sound
cardiovascular sound
frequency
cardiovascular
parameters
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
US12/308,755
Other languages
English (en)
Inventor
Claus Graff
Samuel Schmidt
Johannes Struijk
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.)
Coloplast AS
Original Assignee
Coloplast AS
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
Application filed by Coloplast AS filed Critical Coloplast AS
Assigned to COLOPLAST A/S reassignment COLOPLAST A/S ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: GRAFF, CLAUS, SCHMIDT, SAMUEL EMIL, STRUIJK, JOHANNES
Publication of US20100145210A1 publication Critical patent/US20100145210A1/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B7/00—Instruments for auscultation
    • A61B7/02—Stethoscopes
    • A61B7/04—Electric stethoscopes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7239—Details of waveform analysis using differentiation including higher order derivatives
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Definitions

  • the present invention relates to methods and systems for classification of heart sounds recorded from a living subject into classes describing whether or not murmurs due to coronary artery stenosis is present in the heart sound.
  • Coronary artery disease is the single most common cause of death from cardiovascular disease in the western world.
  • the heart muscle receives its blood supply through the coronary arteries, and atherosclerosis is the most common pathophysiologic process occurring in the coronary arteries giving rise to coronary artery disease (CAD).
  • Atherosclerosis is a process that builds up plaques within the artery, and the blood flow can therefore be is reduced or even blocked by the plaque.
  • the constantly working heart requires a continuous and efficient blood supply in order to work properly. Defects in the blood supply may be very severe and even fatal. Increasing degrees of luminal diameter reduction or stenosis of the coronary artery will first limit reserve flow, then reduce flow at rest and may finally totally occlude the vessel.
  • MRI magnetic resonance imaging
  • IVUS in vivo intravascular ultrasound
  • OCT optical coherence tomography
  • a clinician may utilize a stethoscope to monitor and record heart sounds in order to detect heart valve diseases.
  • the recorded heart sounds could be digitized, saved and stored as data files for later analysis.
  • Devices have been developed that apply algorithms to electronically recorded auscultatory sounds.
  • One example is an automated blood-pressure monitoring device.
  • Other examples include analysis systems that attempt to automatically detect physiological conditions based on the analysis of auscultatory sounds.
  • artificial neural networks have been discussed as one possible mechanism for analyzing auscultatory sounds and providing an automated diagnosis or suggested diagnosis. Using these conventional techniques, it is difficult to provide an automated device for diagnosis of coronary stenosis using auscultatory sounds. Moreover, it is often difficult to implement the conventional techniques in a manner that may be applied in real-time or pseudo real-time to aid the clinician.
  • the object of the present invention is to solve the above-mentioned problems.
  • a method for classifying a cardiovascular sound recorded from a living subject comprises the step of extracting at least two signal parameters from said cardiovascular sound, said at least two signal parameters characterize at least two different properties of at least a part of said cardiovascular sound.
  • the method further comprises the step of classifying said cardiovascular sound using said at least two signal parameters in a multivariate classification method.
  • the cardiovascular sound related to turbulence consists of at least two components: a broad band component caused by turbulent blood flow colliding with the arterial wall and a narrow banded component related to the resonance frequency of the artery wall, therefore different variables describing different properties are needed in order to perform a robust classification.
  • the different properties describe different characteristics of the cardiovascular sound and would therefore be uncorrelated and therefore provide different information of the cardiovascular sound.
  • Different properties could for instance be the time duration of the diastolic segment of cardiovascular sound, the time duration of the systolic cardiovascular sound, the most dominant frequency component of the sound, the bandwidth of different frequency components, the energy in two frequency bands, the mobility of part of the signal, the complexity of the signal, the power ratio between different parts of the signal, e.g. two different segments or two different frequency bands, morphological characteristics such as correlation ratios between different segments or amplitude change over time.
  • the method could easily be implemented in any kind of data processor unit and therefore be e.g. integrated in a software program which clinicians and doctors could use in order to classify the cardiovascular sound.
  • the method could be integrated in a digital stethoscope and the stethoscope could therefore be used in order to classify a patient's cardiovascular sound. Since doctors and other clinicians are familiar with a stethoscope, they could easily be taught to use the stethoscope to classify the cardiovascular sound. The result is that the classification could assist the doctor or other clinicians to diagnose whether or not the patient suffers from CAD.
  • At least one of said at least two signal parameters is a frequency parameter describing a property in the frequency domain of at least a part of said cardiovascular sound.
  • the frequency components of the cardiovascular sound could be used as a parameter in the multivariable classification method.
  • Frequency parameters are very good parameters for classifying whether or not murmurs due to stenosis are present in a cardiovascular sound because the stenosis would change the frequency components of the cardiovascular sound.
  • At least one of said at least two signal parameters describes a property in the time domain of at least a part of said cardiovascular sound.
  • time properties of the cardiovascular sound could be used as a parameter in the multivariable classification method.
  • Time properties like the mobility or number of turning points are good indicators, whether or not murmurs due to stenosis are present in cardiovascular sound.
  • time and frequency parameters a very robust classification of the cardiovascular is achieved since time and frequency properties are often uncorrelated.
  • At least one of said frequency is parameters is a frequency level parameter describing a frequency level property of at least a part of said cardiovascular sound.
  • a frequency level property of the cardiovascular sound is used in the multivariable classification method.
  • the murmurs would typically change the frequency level of the cardiovascular sound, and by using parameters describing the frequency level of the sound a robust classification of the cardiovascular sound could be achieved.
  • At least one of said at least two signal parameters is a frequency bandwidth parameter describing a frequency bandwidth property of at least a part of said cardiovascular sound.
  • the bandwidth of, for instance, dominating frequency components could be used in the multivariable classification method.
  • the advantage of using a frequency bandwidth property of the cardiovascular sound is that murmurs often has a limited frequency bandwidth, and the frequency bandwidth parameter would therefore be a good indicator of whether or not murmurs due to stenosis are present in the cardiovascular sound.
  • At least one of said frequency level properties characterizes the most powerful frequency component of at least a part of said cardiovascular sound.
  • This parameter is a very useful parameter as the murmurs due to stenosis typically have a dominating frequency component between 200-800 Hz. And if the most powerful frequency component is inside this interval, it would be a good indication of the presence of murmurs due to stenosis.
  • At least one of said frequency bandwidth properties characterizes the bandwidth of the most powerful frequency component of at least a part of said cardiovascular sound.
  • the bandwidth of the most powerful frequency component could be used in the multivariable classification method. This bandwidth would most likely depend on whether or not murmurs due to stenosis are present in the cardiovascular sound.
  • At least one of said time parameters is a property characterizing the mobility of at least a part of said cardiovascular sound.
  • the mobility is a good indicator of whether or not murmurs due to stenosis are present in the cardiovascular sound.
  • the mobility describes the variance of the sound, and since murmurs would cause larger variance in the sound the mobility would be a good indicator.
  • the method further comprises the step of dividing said cardiovascular sound into at least one sub-segment and at least one of said signal parameters is extracted from said at least one sub-segment.
  • the cardiovascular sound could be divided into sub-segments, e.g. into a systolic part and a diastolic part.
  • relevant sub-segments could be used to extract the above-described different parameters.
  • the method further comprises the step of modelling at least a part of said cardiovascular sound and at least one of said signal parameters is extracted from said model.
  • time models and frequency models of the cardiovascular sound or sub-segments of the sound could e.g. be used to extract the above-described parameters.
  • the advantage of using models is that the models could enhance the signal properties, e.g. by using an envelope function or an autoregressive model.
  • models would simplify and optimize the calculation process when the method is implemented in a data processor.
  • the multivariate classification method is a discriminant function.
  • a simple and fast implementation of the classification method is achieved.
  • any number of parameters could be used in the discriminant function, and the different parameters could also be weighted differently depending on the parameters' significance.
  • the discriminant function could also be trained using cardiovascular test sounds is recorded from patients suffering from stenosis and healthy patients. Thereby the weights of the different parameters could be optimized to experimental data.
  • the invention further relates to a system for classifying a cardiovascular sound recorded from a living subject, said system comprises processing means for extracting at least two signal parameters from said cardiovascular sound, said at least two signal parameters characterizes at least two different properties of at least a part of said cardiovascular sound; processing means for classifying said cardiovascular sound using said at least two signal parameters using a multivariate classification method.
  • processing means for classifying a cardiovascular sound can be constructed and hereby the same advantages as described above are achieved.
  • said processing means for extracting at least two signal parameters from said cardiovascular sound is adapted to extract at least one frequency parameter describing a property in the frequency domain of at least a part of said cardiovascular sound.
  • said processing means for extracting at least two signal parameters from said cardiovascular sound is adapted to extract at least one time parameter describing a property in the time domain of at least a part of said cardiovascular sound.
  • said processing means adapted to extract at least one of said frequency parameters are further adapted to extract at least one frequency level parameter describing a frequency level property of at least a part of said cardiovascular sound.
  • said processing means adapted to extract at least one frequency parameter is further adapted to extract at least one frequency bandwidth parameter describing a frequency bandwidth property of at least a part of said cardiovascular sound.
  • said processing means adapted to extract at least one frequency level property is further adapted to extract the most powerful frequency component of at least a part of said cardiovascular sound.
  • said processing means adapted to extract at least one of said frequency bandwidth properties are further adapted to extract the bandwidth of the most powerful frequency component of at least a part of said cardiovascular sound.
  • said processing means for extracting at least one time parameter are further adapted to extract the mobility of at least a part of said cardiovascular sound.
  • said system further comprises processing means for dividing said cardiovascular sound into at least one sub-segment and at least one of said signal parameters is extracted from said at least one sub-segment.
  • said system further comprises processing means for modelling at least a part of said cardiovascular sound and in that said processing means for extracting at least two signal parameters from said cardiovascular sound are further adapted to extract at least one of said parameters from said model.
  • said multivariate classification method used by said processing means for classification of said cardiovascular sound is a discriminant function.
  • the invention further relates to a computer-readable medium having stored therein instructions for causing a processing unit to execute a method as described above.
  • a processing unit to execute a method as described above.
  • the invention further relates to a stethoscope comprising recording means adapted to record a cardiovascular sound from a living subject, storing means adapted to store said recorded cardiovascular sound, a computer-readable medium and a processing unit, said computer-readable medium having stored therein instructions for causing said processing unit to execute a method according to claims 1 - 12 and thereby classify said recorded cardiovascular sound.
  • a stethoscope comprising recording means adapted to record a cardiovascular sound from a living subject, storing means adapted to store said recorded cardiovascular sound, a computer-readable medium and a processing unit, said computer-readable medium having stored therein instructions for causing said processing unit to execute a method according to claims 1 - 12 and thereby classify said recorded cardiovascular sound.
  • the invention further relates to a server device connected to a communication network comprising receiving means adapted to receive a cardiovascular sound recorded form a living subject through said communication network, storing means adapted to store said received cardiovascular sound, a computer-readable medium and a processing unit, said computer-readable medium having stored therein instructions for causing said processing unit to execute a method as described above and thereby classify said received cardiovascular sound.
  • a server device connected to a communication network comprising receiving means adapted to receive a cardiovascular sound recorded form a living subject through said communication network, storing means adapted to store said received cardiovascular sound, a computer-readable medium and a processing unit, said computer-readable medium having stored therein instructions for causing said processing unit to execute a method as described above and thereby classify said received cardiovascular sound.
  • said receiving means are further adapted to receive said cardiovascular sound from a client connected to said communication network.
  • a clinician/doctor could send a cardiovascular sound to the server using a client device such as a laptop.
  • the server could thereafter classify the received cardiovascular sound.
  • the server device further comprises means for sending said classification of said cardiovascular sound to at least one client unit connected to said communication network.
  • the result of the classification can be sent back to a client, and the clinician/doctor can therefore receive the result of the classification.
  • FIG. 1 illustrates a graph of a typical heart sound
  • FIG. 2 illustrates a fluid dynamic model of an arterial stenosis
  • FIG. 3 illustrates an overview in form of a flow diagram of the method according to the present invention
  • FIG. 4 illustrates an embodiment of the system according to the present invention
  • FIG. 5 illustrates another embodiment of the method according to the present invention
  • FIG. 6 illustrates a flow diagram of the segmentation method
  • FIG. 7 illustrates for a heart sound the relationship between the envelope autocorrelation of a cardiac cycle and the cardiac cycle
  • FIG. 8 illustrates the implementation of a Bayesian network used to calculate the probability of a sound being an S 1 , S 2 and noise sound.
  • FIG. 1 illustrates a graph of a typical heart sound recorded by a stethoscope and shows the amplitude (A) of the sound pressure at the y-axis and time (t) at the x-axis.
  • the heart sounds reflect events in the cardiac cycle: the deceleration of blood, turbulence of the blood flow and the closing of valves.
  • the closing of the valves is typically represented by two different heart sounds, the first (S 1 ) and the second (S 2 ) heart sound.
  • the first and second heart sounds are illustrated in the figure, and S 1 marks the beginning of systole which is the part of the cardiac cycle in which the heart muscle contracts, forcing the blood into the main blood vessels, and the end of the diastole which is the part of the heart cycle during which the heart muscle relaxes and expands.
  • S 1 marks the beginning of systole which is the part of the cardiac cycle in which the heart muscle contracts, forcing the blood into the main blood vessels, and the end of the diastole which is the part of the heart cycle during which the heart muscle relaxes and expands.
  • blood fills the heart chambers.
  • the duration of systolic segments is nearly constant around 300 ms for healthy subjects. Given a pulse of 60 beats per minute the duration of a cardiac cycle will be one second on average, and the duration of the diastole will be 700 ms. However, the diastolic durations are not constant, but will vary depending on the subject's pulse. In addition, smaller variations of the
  • FIG. 2 illustrates a fluid dynamic model of an arterial stenosis and shows an artery ( 201 ) with a stenotic lesion ( 202 ).
  • the arrows ( 203 ) indicate the blood flow through the artery.
  • Vortices ( 204 ) will occur when high velocity blood exits a stenotic lesion ( 202 ). These vortices collide with the arterial wall ( 205 ) and are transformed into pressure vibrations that cause the arteries to vibrate at their resonance frequencies. The result is that soundwaves in the form of murmers ( 206 ) with a frequency corresponding to the aterial wall's resonance frequencies are created and emitted from the arterial wall.
  • Resonance frequencies in the arterial segment are increased if a stenosis is present and their frequencies depend on the diameter of the stenotic segment compared to the diameter of the artery. As the severity of a stenosis increases, so does the resonance frequency.
  • the resonance frequency of a partial occluded stenotic artery is most likely between 200 Hz to 1100 Hz.
  • the intensity of the vortice fluctuations depends on the blood flow so that murmurs from the left coronary arteries are most intense during diastole, when the blood flow through these arteries is highest.
  • Murmurs from the right coronary arteries are most intense during diastole if there is a stenosis in branches of the right coronary artery supplying the right-sided cavities, whereas the murmur more likely will be systolic from those branches of the right coronary artery giving arterial blood to the left ventricle.
  • the intensity of murmurs not only depends on the blood flow, but also on the frequency content of a murmur. High murmur frequencies are more suppressed by the chest wall compared to low frequencies.
  • the murmurs caused by the arterial vibrations would affect the graph of a heart sound recorded by e.g. a stethoscope.
  • FIG. 3 illustrates an overview in form of a flow diagram of the method according to the present invention.
  • the method could for instance be implemented as a software program running on a computer or on a microcontroller implemented in a stethoscope.
  • the method starts with an initialization ( 301 ), receiving a test signal ( 302 ), dividing the test signal into relevant segments ( 306 ), filtering the relevant segments ( 307 ); calculating/developing a model of the signal ( 308 ) in relevant segments; extracting different parameters from the signal and the model ( 309 ), performing an analysis of the signal ( 310 ) using the extracted parameters and classifying the relevant segments into two groups: one indicating that the signal contains murmurs due to stenosis ( 311 ), and one indicating that the signal does not contain murmurs due to stenosis ( 312 ).
  • the method receives the test signal ( 302 ) as a data file ( 303 ).
  • the test signal would be the heart sound from a person ( 304 ) recorded and digitalized into a data file, e.g. by a digital stethoscope ( 305 ).
  • the test signal would be similar to the heart sound illustrated in FIG. 1 , however, the duration of the test signal would typically be 5-15 times longer than the signal shown in FIG. 1 .
  • segmentation 306 is performed in order to detect and divide the test signal into segments. The segmentation process would typically detect the heart sounds S 1 and S 2 and thereafter divide the test signal into systolic and diastolic parts.
  • the test signal is filtered ( 307 ), and the filtration process includes an autoregressive filter that reduces white noise in the signal and a band pass filter that only lets frequencies between 450-1100 Hz pass.
  • the test signal would thereafter contain the frequencies caused by the vibrations of the arterial wall when stenosis is present in the artery.
  • the autoregressive filter could be implemented as a Kalman filter that is a powerful estimator of past, present and future states and it can do so even when the precise nature of the modelled system is unknown. This is a desirable feature in the present application when reducing the effects of noise since the exact composition of a murmur is unknown.
  • a first order Kalman filter can reduce the effects of white noise and smooth the noisy heart sound recordings for further processing.
  • the band pass filter could be implemented as a wavelet filter.
  • the Kalman filter is omitted in order to simplify the implementation of the method in e.g. a microprocessor and further to reduce the number of calculations performed by the microprocessor.
  • relevant segments are selected for further analysis.
  • a part of the diastolic segment is selected for further analysis as the murmur due to stenosis is most likely to be audible in the diastolic segment.
  • a mathematical model of the signal in the selected segment is hereafter calculated/developed ( 308 ) using the sampled heart sound in the data file.
  • the model is used to extract parameters that characterize the sound in the segment and could be used to categorize whether or not the murmurs due to stenosis exist in the sound segment.
  • an autoregressive all-pole parametric estimation (AR-model) is used to model the signal.
  • AR-model autoregressive all-pole parametric estimation
  • the sampled sound signal, y from the data file is modelled as a linear combination of M past values of the signal and the present input, u, driving the sound generating process.
  • the model can be described by the following equation:
  • M represents the model order
  • a p the AR coefficients and n the sample number.
  • the AR coefficients are determined through an autocorrelation and by minimizing the error associated with the model.
  • the AR model in this embodiment is used to extract frequency parameters describing the heart sound.
  • parameters are extracted ( 309 ) from the sampled signal and the AR model using signal processing techniques.
  • Some parameters could be extracted from the selected segments.
  • Each parameter characterize the heart sound in the selected segments and could therefore be used to categorize the heart sound, e.g. whether or not murmurs due to stenosis are present in the heart sound.
  • the parameters can in this embodiment be the number of turnings points per signal length, TP; the mobility of the signal, MB; pole magnitude, PM; normalized AR-peak frequency, NF; and AR spectral ratio, SR.
  • the number of turning points TP is extracted from the sampled signal in the time domain, and it is found by calculating the number of turns the signal performs in the time domain per unit time. This could be done by determining the amount of local maxima in a time period.
  • T ⁇ ⁇ P number ⁇ ⁇ of ⁇ ⁇ turns signal ⁇ ⁇ length [ 3.2 ]
  • the mobility MB is extracted from the sampled signal in the time domain and found by calculating the variance, ⁇ y , of the signal in the time domain and the variance of the signal's first derivative, ⁇ y′ .
  • the mobility is hereafter found by:
  • the pole magnitude PM is found by transforming the AR-model into the z-domain and calculating the magnitude of the poles in z-domain described by the AR-spectrum.
  • the normalized AR peak frequency NF is based on the assumption that murmurs due to stenosis are more likely to be found in the diastolic segment than in the systolic segment.
  • the NF is found by calculating the angle of the poles in the AR-spectrum in the z-plane and transforming this into a frequency of both a diastolic segment and a systolic segment. If the absolute difference between the two is less than 25 Hz, which is typical in cases where no murmurs due to stenosis are present, then 25 Hz is subtracted from the diastolic peak frequency. If the average diastolic frequency is more than 50 Hz greater than the average systolic peak frequency, which is typical when murmurs due to stenosis are present, then 25 Hz is added to the average peak diastolic frequency.
  • the AR spectral ratio SR is found by calculating the ratio of the energy in the frequency rang 200-500 Hz to the energy in the frequency range 500-1000 Hz of a diastolic segment.
  • the extracted parameters are thereafter used in a multiparametric discriminant function in order to classify whether or not the sound segment contains murmurs due to stenosis ( 310 ).
  • a linear discriminant function is used to classify the sound segments.
  • the linear discriminant function combines weighted features into a discriminant score g(x) and could be described by:
  • x is the feature vector consisting of the extracted parameters
  • k represents the number of features
  • i represents the classes
  • w is a weight vector that holds the discriminant coefficients.
  • a two class classifier is called a dichotomizer.
  • the discriminant score g(x) is the algebraic distance to the hyper-plane.
  • the discriminant function needs to be trained in order to find the weights values, w, and make a safe and robust classification of the sound segments.
  • the discriminant training procedure needs to be performed before using the system, and the purpose of the procedure is to find the optimal weights values of w so that the hyper plane separates the feature vectors optimally.
  • the training procedure is in one embodiment carried out by using 18 test sounds recorded from 18 test persons where nine test persons have coronary stenosis and nine test persons do not have coronary stenosis.
  • the discriminant training procedure is performed by using the statistical software program SPSS v.12.0 for windows (SPSS inc., Chicago Ill., USA).
  • SPSS statistical software program
  • MB is the mobility of the signal
  • NF the AR-peak frequency
  • PM the pole magnitude
  • SR the AR spectral ratio
  • TP the number of turning point
  • the discriminant function is larger than zero (g(x)>0) then the sound segment does not contain murmurs due to stenosis ( 312 ). On the other hand, if the discriminant function is smaller than zero (g(x) ⁇ 0) then the sound segment contains murmurs due to stenosis ( 311 ).
  • the discriminant function could by a person skilled in the art easily be adjusted to include additional or fewer parameters in order to develop a proper discriminant function that can be used to classify the heart sound. Further parameters could for instance be:
  • the Complexity, CP of the sampled signal in the time domain. This parameter is based on the ratio of the mobility of the first derivative of the signal to the mobility of the signal itself where y′′ is the second derivative of the filtered heart sound signal.
  • the complexity measure is relatively sensitive to noisy signals since it is based on the second derivative.
  • the AR-peak frequency could be extracted and used in the discriminant function.
  • the AR-peak frequency could be found by calculating the angle of the AR poles in the z-plane.
  • the parameters used in the discriminant function could be extracted from different segments of the heart sound, e.g. a number of different diastolic segments where a number of parameters is extracted from each diastolic segment. Thereafter an average value of each parameter could be calculated and used as input in the discriminant function.
  • FIG. 4 a illustrates an embodiment of the system according to the present invention where a server ( 401 ) is programmed to execute the method described in FIG. 3 . Furthermore, the server is connected to a network ( 402 ), e.g. the Internet and adapted to on request to receive and analyze heart sound. Clinicians or other medical professionals would record the heart sound from a patient by a digital stetoscope ( 305 ) and thereafter transmit the digitalized heart sound to a personal computer ( 403 ). The clinician can hereafter send a request to the server in order to have the heart sound analyzed. Once the server has analyzed the heart sound the result is automatically sent back to the clinician.
  • FIG. 4 b illustrates a flow diagram of the process and the communication between the personal computer ( 403 ) and the server.
  • the left hand side represents the client side ( 410 ) and the right hand side represents the server side ( 411 ).
  • the client sends a heart sound in digital form to the server ( 412 ).
  • the server performs the method illustrated in FIG. 3 and sends ( 413 ) the result of the analysis back to the client where it is displayed ( 414 ) to the clinician.
  • the clinician could hereafter evaluate the result in order to choose the right treatment of the patient.
  • the system according to the present invention could also be implemented as an all in one digital stethoscope.
  • the stethoscope would therefore automatically perform the analysis described in FIG. 3 when a heart sound has been recorded.
  • An advantage of this embodiment is that most clinicians are familiar with a digital stethoscope and could therefore easily learn to use the stethoscope to diagnose whether or not the patient has a coronary stenosis.
  • FIG. 5 illustrates another embodiment of the method described in FIG. 3 .
  • relevant segments are selected for further analysis.
  • a part of the diastolic segment is selected for further analysis as the murmur due to stenosis is most likely to be audible in the diastolic segment.
  • the diastolic segment comprising respiration sounds is discarded ( 501 ). This is done by calculating the energy level of the diastolic segment in the frequency band 200-440 Hz and comparing this energy level with the median energy level of the entire diastolic segment. The diastolic segment would be discarded if the energy level of the 200-440 Hz frequency band is a factor 1.1 larger than the energy level in the entire diastolic segment.
  • the remaining diastolic segments are hereafter divided into sub-segments ( 502 ) with a duration of 37.5 ms or 300 samples. This is done because the blood flow in the coronary artery is not constant during a diastole, and the murmurs due to stenosis would therefore not be constant.
  • the variance of the signal in all sub-segments is then calculated and the sub-segments with a variance larger than 1.3 of the median variance of all sub-segments are then discarded ( 503 ). Hereby sub-segments comprising high noise spikes are removed.
  • none stationary sub-segments are removed. This is done by dividing the sub-segment into sub-sub-segments with a duration of 3.75 ms or 30 samples and then calculate the variance of each sub-sub-segments. Thereby an outline of the variance throughout the sub-segment is constructed. The variance of the outline is then calculated and the sub-segment is removed if the variance of the outline is larger than 1.
  • the remaining sub-segments are thereafter used in step ( 308 ) and ( 309 ) as described in FIG. 3 in order to extract parameters describing different properties of the cardiovascular signal. Thereafter the median of each parameter is calculated using the values of the parameter form each sub-segment ( 505 ). The median of each parameter is thereafter used in the multiparametric discriminant function as described in FIG. 3 . In this embodiment the following parameters are used: the mobility, the power-ratio and the pole-amplitude of a 3 pole in an AR model of order 6 .
  • FIG. 6 illustrates a flow diagram of the segmentation method ( 306 ) according to the present invention used to automatic divide a heart sound ( 601 ) into sub-segments.
  • the heart sound ( 601 ) has been recorded by a stethoscope and the signal has been digitized in order to digitally process the signal.
  • the graph shows the amplitude (A) of the sound intensity as a function of time (t).
  • the heart sounds reflect events in the cardiac cycle; the deceleration of blood, turbulence of the blood flow and the closing of valves.
  • the closing of the valves is typically represented by two different heart sounds, the first (S 1 ) and the second (S 2 ) heart sound.
  • the first and second heart sounds are illustrated in the figure, and (S 1 ) marks the beginning of systole, which is the part of cardiac cycle in which the heart muscle contracts, forcing the blood into the main blood vessels, and the end of the diastole which is the part of the heart cycle during which the heart muscle relaxes and expands. During diastole, blood fills the heart chambers.
  • the purpose of the segmentation method is to classify the recorded heart sound into systolic, diastolic and noise segments.
  • the illustrated method includes steps of noise reduction ( 602 ) followed by envelope creation ( 603 ).
  • the noise reduction could be implemented as a high-pass filer followed by removal of high amplitude friction noise spikes due to external noise like movement of the stethoscope during recording and thereafter a low pass filter.
  • the purpose of the envelope creation is to enhance the trend of the signal.
  • the envelope is in this embodiment created by calculating the Shannon energy of the signal:
  • Candidates S 1 's and S 2 's are then detected ( 605 ) using the time intervals extracted above and a threshold ( 614 ) on the envelope ( 613 ).
  • a threshold 614
  • a minimum requirement is applied to the candidate segments, which effectively removes some of the erroneously detected noise spikes.
  • the segmentation method performs a test for missing S 1 and S 2 sounds ( 606 ). If it can be determined that some segments are missing, the threshold procedure is rerun ( 607 ) using lower local thresholds.
  • interval parameters and frequency parameters for each segment are then extracted ( 608 ).
  • the parameters aid in the classification of the sounds into systolic segments and diastolic segments.
  • the interval parameters are four Boolean parameters extracted for each sound by comparing the time duration to the previous sound and to the next sound with the time intervals extracted using the autocorrelation.
  • the parameters are:
  • the frequency parameter divides the sounds into low frequency and high frequency sounds by calculating the median frequency of the sound. This is useful information as the first heart sound is expected to be a low frequency sound and the second heart sound is expected to be a high frequency sound.
  • the parameters are parsed into a Bayesian network where the probability of a segment being a S 1 , S 2 and noise sound is computed ( 609 ).
  • the figure illustrates a bar chart ( 615 ) of the probability calculated for each sound in the heart signal ( 601 ).
  • Each sound would typically have one dominating probability indicating the type (S 1 , S 2 or noise) of the sound. Thereby all sounds are classified into S 1 , S 2 and noise sounds.
  • the probability of the three types would in some cases be more or less equal and in such cases it is not possible to classify the sound into a S 1 , S 2 or noise sound using the Bayesian network.
  • the probabilities are used in the last step ( 610 ) to divide and verify the heart signal into systole and diastole segments. This is done by using the position of the identified S 1 and S 2 sounds to mark the beginning of a systolic and diastolic sound segment respectively
  • the final result of the method ( 611 ) is the beginnings and ends of all identified systoles and diastoles. Therefore a “train” ( 616 ) of alternating systoles ( 617 ) and diastoles ( 618 ) can be created. Once the systoles and diastoles have been identified they can be used in further data handling, e.g. to extract further parameters from these segments and thereafter use the parameters to classify the medical condition of the recorded heart sound.
  • FIG. 7 illustrates the relationship between the envelope autocorrelation and the cardiac cycle, and how the intervals between heart sounds S 1 and S 2 can be found from the autocorrelation.
  • FIG. 7 a illustrates the envelope autocorrelation with the normalized autocorrelation at the y-axis (NA) and the displacement (m) of the shifted envelope at the x-axis.
  • FIG. 7 b illustrates the displacement (m 1 ) when the shifted envelope ( 701 ) is displaced by the duration of the systole corresponding to the unshifted envelope ( 702 ).
  • the y-axis shows the amplitude (A) of the envelope and the x-axis the time (t).
  • the S 1 's in the displaced envelope are multiplied by the S 2 's in the unshifted envelopes resulting in the first peak ( 703 ) seen in the autocorrelation.
  • FIG. 7 c illustrates the displacement (m 2 ) when the shifted envelope ( 701 ) is displaced by the duration of the diastole corresponding to the unshifted envelope ( 702 ).
  • the displaced S 2 's are multiplied by the S 1 's in the unshifted envelope resulting in the second peak ( 704 ) seen in the autocorrelation.
  • FIG. 6 b illustrates the displacement (m 3 ) when the shifted envelope ( 701 ) is displaced by the duration of the cardiac cycle corresponding to the unshifted envelope ( 702 ).
  • the S 1 's in the displaced envelope are multiplied by the S 1 's in the unshifted envelope, and the S 2 's in the displaced envelope are multiplied by the S 2 's in the unshifted envelope.
  • the dominating peak ( 705 ) in the autocorrelation is produced.
  • the interval between the heart sounds could therefore be found by measuring the distance between the peaks in the autocorrelation as described above.
  • FIG. 8 illustrates the implementation of the Bayesian network used to calculate the probability of a sound of being an S 1 , S 2 and noise sound in step ( 809 ).
  • the basic concept in the Bayesian network is the conditional probability and the posterior probability.
  • the conditional probability describes the probability of the event a given the event b.
  • Equation [8.3] only describes the relation between one parent and one child, but since the event a can be the combination of several events ⁇ a 1 , a 2 , , , a n ⁇ the equation can be expanded to:
  • equation [8.5] can be generalized to:
  • Equation [8.6] is useful in determining the probability of the event b if the states of all a events are known and if all a events are conditionally independent.
  • a Bayesian network based on equation [8.6] is called a naive Bayesian network because it requires conditional independency of the children.
  • the task for the Bayesian network is to evaluate the type of each detected sound above the detection threshold. For each of these sounds, the posterior probability of being an S 1 sound, an S 2 sound or a noise component is calculated and the Bayesian network is constructed using one parent and five children.
  • the parent is a sound above the envelope threshold ( 801 ), and the children are the five parameters described above: Frequency ( 802 ), AfterSys ( 803 ), AfterDia ( 804 ), BeforeSys ( 805 ) and BeforeDia ( 806 ).
  • the prior probability for the different states of a sound type P(S) and the conditional probabilities must be known, i.e.
  • This posterior probability requires definition of P(S), P(AfterSys
  • the prior probability that a sound is an S 1 , S 2 or a noise component changes between recordings.
  • this optimal condition cannot be assumed for real signals, and noise sounds would be detected. This will increase the prior probability that a given sound is noise.
  • the S 1 sounds will normally be followed by S 2 sounds after an interval of duration equal to the systole. Deviations from this can also occur, e.g. when S 1 is the last sound in the recording, or if S 2 is missing because it is not detected by the threshold. It may also occur that a weak (below threshold) S 2 is detected because noise occurs in the tolerance window associated with those sounds.
  • the probability that “AfterSys” is false if the sound is an S 1 sound may thus be calculated as
  • EndSound is an event describing that the sound is the last sound in the recording.
  • SingleSound describes that S 1 is not followed by S 2 as the next S 2 sound is not detected due to sub-threshold amplitude.
  • NoiselnWin describes noise occurrence in the window, where the S 2 sound was expected. The conditional probability that “AfterSys” is true given that the examined sound is an S 1 sound is given by:
  • the examined sound is an S 2 sound it is not likely that any sound occurs after an interval corresponding to the systolic duration since the next S 1 sound will occur after the duration of the diastole.
  • An exception is if a noise sound occurs in the window P(NoiseInWin) or if the systole and diastole durations are equal. If the duration of the diastole is equal to the duration of the systole, the S 1 sound which follows the S 2 sound after the duration of a diastole occurs in both the systole tolerance window and in the diastole tolerance window. This will happen if the heart rate of the subject is high.
  • conditional probability that a sound does not occur after a systole duration, if the examined sound is an S 2 is the opposite of the conditional probability that it does occur:
  • conditional probability that a detected noise sound is followed by another sound after the systole duration is based on the probability that a sound of any kind is present in a segment with the length of the used tolerance window. This can be estimated from the ratio of the tolerance window length multiplied by the number of detected sounds minus one to recording length.
  • the conditional probability that a noise sound is followed by another sound after a systole duration is therefore:
  • N sound is the number of sounds within the recording
  • Sys tot is the duration of a systole
  • RecLength is the length of the recording.
  • S) are based on the same assumptions used to define P(AfterSys
  • equation [8.6] is used by the Bayesian network to calculate the posterior probabilities for all detected sounds. This way, three probabilities are calculated for each sound that reflect how likely it is that the current sound is a given type.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Acoustics & Sound (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Mathematical Physics (AREA)
  • Fuzzy Systems (AREA)
  • Vascular Medicine (AREA)
  • Cardiology (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
US12/308,755 2006-06-26 2006-06-26 Multi parametric classification of cardiovascular sounds Abandoned US20100145210A1 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/DK2006/000374 WO2008000254A1 (fr) 2006-06-26 2006-06-26 Classification multiparamétrique de bruits cardio-vasculaires

Publications (1)

Publication Number Publication Date
US20100145210A1 true US20100145210A1 (en) 2010-06-10

Family

ID=37680714

Family Applications (2)

Application Number Title Priority Date Filing Date
US12/308,755 Abandoned US20100145210A1 (en) 2006-06-26 2006-06-26 Multi parametric classification of cardiovascular sounds
US12/308,754 Active 2029-07-02 US8684943B2 (en) 2006-06-26 2007-06-26 Multi parametric classfication of cardiovascular sound

Family Applications After (1)

Application Number Title Priority Date Filing Date
US12/308,754 Active 2029-07-02 US8684943B2 (en) 2006-06-26 2007-06-26 Multi parametric classfication of cardiovascular sound

Country Status (11)

Country Link
US (2) US20100145210A1 (fr)
EP (1) EP2034900B1 (fr)
JP (1) JP5576116B2 (fr)
CN (1) CN101478921B (fr)
AU (1) AU2007264203B2 (fr)
BR (1) BRPI0713474A2 (fr)
CA (1) CA2654252C (fr)
DK (1) DK2034900T3 (fr)
HK (1) HK1129821A1 (fr)
RU (1) RU2449730C2 (fr)
WO (2) WO2008000254A1 (fr)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090177107A1 (en) * 2005-04-13 2009-07-09 Marie A. Guion-Johnson Detection of coronary artery disease using an electronic stethoscope
US20110137210A1 (en) * 2009-12-08 2011-06-09 Johnson Marie A Systems and methods for detecting cardiovascular disease
US8702628B2 (en) 2010-03-18 2014-04-22 Panasonic Corporation Physiological sound examination device
US8870791B2 (en) 2006-03-23 2014-10-28 Michael E. Sabatino Apparatus for acquiring, processing and transmitting physiological sounds
WO2013154655A3 (fr) * 2012-04-09 2015-01-22 Geissler Companies, Llc Système et procédé de diagnostic de maladies bovines faisant appel à une analyse d'auscultation
US20150157218A1 (en) * 2013-12-05 2015-06-11 Comsats Institute Of Information Technology Heart murmur extraction and heart impairments identification using fuzzy controller
US20150164466A1 (en) * 2010-08-25 2015-06-18 Diacoustic Medical Devices (Pty) Ltd System and method for classifying a heart sound
US20160364963A1 (en) * 2015-06-12 2016-12-15 Google Inc. Method and System for Detecting an Audio Event for Smart Home Devices
US10064580B2 (en) 2008-11-07 2018-09-04 Intervet Inc. System and method for determining antibiotic effectiveness in respiratory diseased animals using auscultation analysis
EP3242600A4 (fr) * 2015-01-07 2018-11-14 Children's National Medical Center Appareil et procédé de détection de souffles cardiaques
WO2019027693A1 (fr) * 2017-07-31 2019-02-07 University Of Pittsburgh-Of The Commonwealth System Of Higher Education Système de détection automatique et précoce de défaillance de volet libre
US10813586B2 (en) 2013-02-06 2020-10-27 Intervet Inc. System and method for determining antibiotic effectiveness in respiratory diseased animals using auscultation analysis
US11284827B2 (en) 2017-10-21 2022-03-29 Ausculsciences, Inc. Medical decision support system
US11389080B2 (en) 2007-11-28 2022-07-19 Intervet Inc. System and method for diagnosis of bovine diseases using auscultation analysis
US11412944B2 (en) * 2015-06-15 2022-08-16 Cvr Medical Corporation Non-invasive method for measuring sound frequencies created by vortices in a carotid artery

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL172754A0 (en) 2005-12-22 2006-04-10 Menashe Shahar Urethral blockage diagnosis
WO2008000254A1 (fr) * 2006-06-26 2008-01-03 Coloplast A/S Classification multiparamétrique de bruits cardio-vasculaires
US20110125061A1 (en) * 2006-12-20 2011-05-26 P Square Medical Ltd. Monitoring Conditions Of A Patient's Urinary System
ES2293860B1 (es) 2007-10-01 2009-02-01 Hierros Y Aplanaciones, S.A. (Hiasa) Sistema de anclaje fusible directo para poste de soporte de una barrera de seguridad o pretil en un terreno rigido.
WO2009098312A1 (fr) * 2008-02-06 2009-08-13 Capis Sprl Méthode et dispositif pour la détermination de la bande de fréquence d'un souffle au cœur
US8560492B2 (en) * 2008-10-20 2013-10-15 Siemens Corporation Robust filtering and prediction using switching models for machine condition monitoring
CN102271589A (zh) * 2008-12-30 2011-12-07 皇家飞利浦电子股份有限公司 用于处理心音信号的方法和系统
EP2462871A1 (fr) * 2010-12-13 2012-06-13 Acarix A/S Système, stéthoscope et procédé pour indiquer le risque de maladie coronarienne
CN103052355B (zh) * 2011-03-25 2015-02-04 松下电器产业株式会社 生物体声音处理装置及生物体声音处理方法
WO2013079073A1 (fr) 2011-12-01 2013-06-06 Aalborg Universitet Surveillance multifactorielle par télématique de santé pour une grossesse et une naissance
JP6031863B2 (ja) * 2012-07-19 2016-11-24 株式会社Ｊｖｃケンウッド 心音分析装置、心音分析方法および心音分析プログラム
JP6036174B2 (ja) * 2012-10-31 2016-11-30 株式会社Ｊｖｃケンウッド 心音情報処理装置、心音情報処理方法および心音情報処理プログラム
JP6036242B2 (ja) * 2012-12-06 2016-11-30 株式会社Ｊｖｃケンウッド 心音情報処理装置、心音情報処理方法およびプログラム
EP3019077B1 (fr) * 2013-07-11 2020-10-21 Coala-Life AB Classification des bruits du coeur
US10271737B2 (en) 2014-09-18 2019-04-30 National Central University Noninvasive arterial condition detecting method, system, and non-transitory computer readable storage medium
EP3078329A1 (fr) 2015-04-09 2016-10-12 Acarix A/S Indication du risque de coronaropathie
JP6604113B2 (ja) * 2015-09-24 2019-11-13 富士通株式会社 飲食行動検出装置、飲食行動検出方法及び飲食行動検出用コンピュータプログラム
JP6468986B2 (ja) * 2015-10-27 2019-02-13 日本電信電話株式会社 ノイズ判定装置、方法、およびプログラム
WO2017075601A1 (fr) * 2015-10-30 2017-05-04 The Johns Hopkins University Dispositifs, algorithmes, systèmes, et procédés de stéthoscope électronique programmable
CN106530927A (zh) * 2016-12-28 2017-03-22 侯江昊 一种展示心脏杂音产生原理的装置
US11045163B2 (en) 2017-09-19 2021-06-29 Ausculsciences, Inc. Method of detecting noise in auscultatory sound signals of a coronary-artery-disease detection system
CN108323158A (zh) * 2018-01-18 2018-07-24 深圳前海达闼云端智能科技有限公司 心音的识别方法及云系统
CN111914655B (zh) * 2020-07-06 2022-07-12 电子科技大学 一种基于显著信号子段提取的脉象分类方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4792145A (en) * 1985-11-05 1988-12-20 Sound Enhancement Systems, Inc. Electronic stethoscope system and method
US5010889A (en) * 1988-02-04 1991-04-30 Bloodline Technology Intelligent stethoscope
US5036857A (en) * 1989-10-26 1991-08-06 Rutgers, The State University Of New Jersey Noninvasive diagnostic system for coronary artery disease
US5218969A (en) * 1988-02-04 1993-06-15 Blood Line Technology, Inc. Intelligent stethoscope
US5550902A (en) * 1994-08-17 1996-08-27 American Telecare, Inc. Remote stethoscope signal processing system
US5595182A (en) * 1994-10-24 1997-01-21 Transonic Systems, Inc. Cardiovascular measurements by sound velocity dilution
US5841846A (en) * 1996-09-16 1998-11-24 American Telecare, Inc. Digital telephonic system for stethoscope signal processing
US6048319A (en) * 1998-10-01 2000-04-11 Integrated Medical Systems, Inc. Non-invasive acoustic screening device for coronary stenosis
US6053872A (en) * 1996-12-18 2000-04-25 Aurora Holdings, Llc Cardiac sonospectrographic analyzer
US20020099286A1 (en) * 2001-01-16 2002-07-25 Sandler Richard H. Acoustic detection of vascular conditions
US20040260188A1 (en) * 2003-06-17 2004-12-23 The General Hospital Corporation Automated auscultation system
US20050090755A1 (en) * 2003-10-22 2005-04-28 Guion Marie A. Analysis of auscultatory sounds using single value decomposition
US20080119749A1 (en) * 2006-11-20 2008-05-22 Cardiac Pacemakers, Inc. Respiration-synchronized heart sound trending
US20100160807A1 (en) * 2006-06-26 2010-06-24 Samuel Emil Schmidt Multi Parametric Classfication of Cardiovascular Sound

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1317708A1 (ru) * 1985-03-29 1988-03-15 Предприятие П/Я М-5629 Устройство дл анализа кардиоритма
DE3732122C2 (de) * 1987-09-24 1996-02-22 Dieter Prof Dr Ing Barschdorff Vorrichtung zur Bewertung einer periodischen Folge von Herzgeräuschsignalen
US5146414A (en) * 1990-04-18 1992-09-08 Interflo Medical, Inc. Method and apparatus for continuously measuring volumetric flow
JP3503962B2 (ja) 1993-07-16 2004-03-08 富士通テン株式会社 プログラム自動配置装置
US20030220584A1 (en) * 2002-04-19 2003-11-27 Southwest Research Institute Headset for measuring physiological parameters
JP4759727B2 (ja) * 2004-03-19 2011-08-31 国立大学法人山口大学 デジタル聴診解析システム

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4792145A (en) * 1985-11-05 1988-12-20 Sound Enhancement Systems, Inc. Electronic stethoscope system and method
US5010889A (en) * 1988-02-04 1991-04-30 Bloodline Technology Intelligent stethoscope
US5218969A (en) * 1988-02-04 1993-06-15 Blood Line Technology, Inc. Intelligent stethoscope
US5036857A (en) * 1989-10-26 1991-08-06 Rutgers, The State University Of New Jersey Noninvasive diagnostic system for coronary artery disease
US5550902A (en) * 1994-08-17 1996-08-27 American Telecare, Inc. Remote stethoscope signal processing system
US5595182A (en) * 1994-10-24 1997-01-21 Transonic Systems, Inc. Cardiovascular measurements by sound velocity dilution
US5841846A (en) * 1996-09-16 1998-11-24 American Telecare, Inc. Digital telephonic system for stethoscope signal processing
US6053872A (en) * 1996-12-18 2000-04-25 Aurora Holdings, Llc Cardiac sonospectrographic analyzer
US6048319A (en) * 1998-10-01 2000-04-11 Integrated Medical Systems, Inc. Non-invasive acoustic screening device for coronary stenosis
US20020099286A1 (en) * 2001-01-16 2002-07-25 Sandler Richard H. Acoustic detection of vascular conditions
US20040260188A1 (en) * 2003-06-17 2004-12-23 The General Hospital Corporation Automated auscultation system
US20050090755A1 (en) * 2003-10-22 2005-04-28 Guion Marie A. Analysis of auscultatory sounds using single value decomposition
US20100160807A1 (en) * 2006-06-26 2010-06-24 Samuel Emil Schmidt Multi Parametric Classfication of Cardiovascular Sound
US20080119749A1 (en) * 2006-11-20 2008-05-22 Cardiac Pacemakers, Inc. Respiration-synchronized heart sound trending

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090177107A1 (en) * 2005-04-13 2009-07-09 Marie A. Guion-Johnson Detection of coronary artery disease using an electronic stethoscope
US10039520B2 (en) 2005-04-13 2018-08-07 Aum Cardiovascular, Inc Detection of coronary artery disease using an electronic stethoscope
US8920343B2 (en) 2006-03-23 2014-12-30 Michael Edward Sabatino Apparatus for acquiring and processing of physiological auditory signals
US8870791B2 (en) 2006-03-23 2014-10-28 Michael E. Sabatino Apparatus for acquiring, processing and transmitting physiological sounds
US11357471B2 (en) 2006-03-23 2022-06-14 Michael E. Sabatino Acquiring and processing acoustic energy emitted by at least one organ in a biological system
US11389080B2 (en) 2007-11-28 2022-07-19 Intervet Inc. System and method for diagnosis of bovine diseases using auscultation analysis
US10064580B2 (en) 2008-11-07 2018-09-04 Intervet Inc. System and method for determining antibiotic effectiveness in respiratory diseased animals using auscultation analysis
US20110137210A1 (en) * 2009-12-08 2011-06-09 Johnson Marie A Systems and methods for detecting cardiovascular disease
US8702628B2 (en) 2010-03-18 2014-04-22 Panasonic Corporation Physiological sound examination device
US20150164466A1 (en) * 2010-08-25 2015-06-18 Diacoustic Medical Devices (Pty) Ltd System and method for classifying a heart sound
US9168018B2 (en) * 2010-08-25 2015-10-27 Diacoustic Medical Devices (Pty) Ltd System and method for classifying a heart sound
WO2013154655A3 (fr) * 2012-04-09 2015-01-22 Geissler Companies, Llc Système et procédé de diagnostic de maladies bovines faisant appel à une analyse d'auscultation
AU2013246522B2 (en) * 2012-04-09 2017-05-11 Intervet International B.V. System and method for diagnosis of bovine diseases using auscultation analysis
US10813586B2 (en) 2013-02-06 2020-10-27 Intervet Inc. System and method for determining antibiotic effectiveness in respiratory diseased animals using auscultation analysis
US20150157218A1 (en) * 2013-12-05 2015-06-11 Comsats Institute Of Information Technology Heart murmur extraction and heart impairments identification using fuzzy controller
US9167971B2 (en) * 2013-12-05 2015-10-27 Comsats Institute Of Information Technology Heart murmur extraction and heart impairments identification using fuzzy controller
EP3242600A4 (fr) * 2015-01-07 2018-11-14 Children's National Medical Center Appareil et procédé de détection de souffles cardiaques
US10251562B2 (en) 2015-01-07 2019-04-09 Children's National Medical Center Apparatus and method for detecting heart murmurs
US9965685B2 (en) * 2015-06-12 2018-05-08 Google Llc Method and system for detecting an audio event for smart home devices
US10621442B2 (en) 2015-06-12 2020-04-14 Google Llc Method and system for detecting an audio event for smart home devices
US20160364963A1 (en) * 2015-06-12 2016-12-15 Google Inc. Method and System for Detecting an Audio Event for Smart Home Devices
US11412944B2 (en) * 2015-06-15 2022-08-16 Cvr Medical Corporation Non-invasive method for measuring sound frequencies created by vortices in a carotid artery
US20200222026A1 (en) * 2017-07-31 2020-07-16 University Of Pittsburgh-Of The Commonwealth System Of Higher Education System for automatic and early detection flap failure
WO2019027693A1 (fr) * 2017-07-31 2019-02-07 University Of Pittsburgh-Of The Commonwealth System Of Higher Education Système de détection automatique et précoce de défaillance de volet libre
US11877885B2 (en) * 2017-07-31 2024-01-23 University of Pittsburgh—of the Commonwealth System of Higher Education System for automatic and early detection flap failure
US11284827B2 (en) 2017-10-21 2022-03-29 Ausculsciences, Inc. Medical decision support system
US11896380B2 (en) 2017-10-21 2024-02-13 Ausculsciences, Inc. Medical decision support system

Also Published As

Publication number Publication date
WO2008000259A3 (fr) 2008-04-10
DK2034900T3 (en) 2016-06-20
EP2034900B1 (fr) 2016-03-09
RU2449730C2 (ru) 2012-05-10
US8684943B2 (en) 2014-04-01
AU2007264203A1 (en) 2008-01-03
JP2009540971A (ja) 2009-11-26
HK1129821A1 (zh) 2009-12-11
WO2008000259A2 (fr) 2008-01-03
AU2007264203B2 (en) 2014-01-30
CA2654252C (fr) 2016-08-09
JP5576116B2 (ja) 2014-08-20
US20100160807A1 (en) 2010-06-24
RU2009102245A (ru) 2010-08-10
BRPI0713474A2 (pt) 2012-02-28
CA2654252A1 (fr) 2008-01-03
CN101478921A (zh) 2009-07-08
WO2008000254A1 (fr) 2008-01-03
CN101478921B (zh) 2011-05-25
EP2034900A2 (fr) 2009-03-18

Similar Documents

Publication Publication Date Title
US8684943B2 (en) Multi parametric classfication of cardiovascular sound
Kachuee et al. Cuffless blood pressure estimation algorithms for continuous health-care monitoring
Schmidt et al. Acoustic features for the identification of coronary artery disease
US9125574B2 (en) System and method for acoustic detection of coronary artery disease and automated editing of heart sound data
Papadaniil et al. Efficient heart sound segmentation and extraction using ensemble empirical mode decomposition and kurtosis features
EP2651304B1 (fr) Système, stéthoscope et procédé pour indiquer le risque d'une coronaropathie
US9198634B2 (en) Medical decision support system
US5109863A (en) Noninvasive diagnostic system for coronary artery disease
US5036857A (en) Noninvasive diagnostic system for coronary artery disease
US5638823A (en) System and method for noninvasive detection of arterial stenosis
US8469896B2 (en) Method for segmenting a cardiovascular signal
Haghighi-Mood et al. A sub-band energy tracking algorithm for heart sound segmentation
WO1997007733A9 (fr) Systeme et procede de detection non invasive de la stenose arterielle
JP6689285B2 (ja) 被験者のバイタルサインを監視する際に用いる加速度計信号を処理する処理装置、システム、及び方法
Schmidt et al. Coronary artery disease and low frequency heart sound signatures
Balogh et al. Application of phonocardiography on preterm infants with patent ductus arteriosus
Shokouhmand et al. Diagnosis of peripheral artery disease using backflow abnormalities in proximal recordings of accelerometer contact microphone (ACM)
Guedes et al. Heart sound analysis for blood pressure estimation
Kumar Automatic heart sound analysis for cardiovascular disease assessment
Mohan et al. Real-time signal processing of photoplethysmographic signals to estimate the on-demand and continuous heart rate by spectral analysis
Cherif et al. Comparison between analysing wavelets in continuous wavelet transform based on the fast Fourier transform: application to estimate pulmonary arterial hypertension by heart sound
Kumar et al. Fourier Analysis Based Respiration Rate Estimation Using Corrupted Photoplethysmogram Signal
Jain Processing of Heart Sound Signal to Monitor Cardiovascular Functions in Real-life Scenario
Turkoglu et al. A computer aided method for diagnose of the heart mitral valve diseases

Legal Events

Date Code Title Description
AS Assignment

Owner name: COLOPLAST A/S,DENMARK

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:GRAFF, CLAUS;SCHMIDT, SAMUEL EMIL;STRUIJK, JOHANNES;SIGNING DATES FROM 20080320 TO 20090327;REEL/FRAME:023344/0108

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION