METHOD AND SYSTEM FOR HEART SOUND IDENTIFICATION

• US 20090192401A1
• Filed: 03/07/2008
• Published: 07/30/2009
• Est. Priority Date: 01/25/2008
• Status: Abandoned Application

• Current Assignee: Texas Instruments, Inc.
• Original Assignee: Texas Instruments, Inc.
• Inventors: Ravindran, Sourabh
• Application Number:
  US12/044,807
• Publication Number:
  US 20090192401A1
• Time in Patent Office: Days
• US Class Current: 600/528
• CPC Class Codes: A61B 7/04 Electric stethoscopes microphones, acoustic transducers therefor H04R

Abstract

Methods, systems, and computer readable media are provided for identification of heart sound components in an audio signal of heart sounds. Time domain kurtosis and frequency domain kurtosis are used to distinguish peaks corresponding to the primary heart sounds, S₁ and S₂, from murmur peaks. Timing based error correction may also be used to verify that appropriate peaks corresponding to the primary heart sounds are identified.

Claims

1. A method for identification of heart sound components comprising:

receiving an audio signal comprising heart sounds;

identifying a first peak corresponding to S₁ within a first search window of the audio signal, wherein identifying the first peak comprises distinguishing the first peak from a murmur peak using time domain kurtosis;

identifying a second peak corresponding to S₁ within a second search window of the audio signal, wherein identifying the second peak comprises distinguishing the second peak from a murmur peak using time domain kurtosis;

identifying a third peak corresponding to S₂ between the first peak and the second peak, wherein identifying the third peak comprises using frequency domain kurtosis to determine whether another peak that may be the third peak is a murmur peak; and

storing a location of the first peak as a first S₁ location, storing a location of the second peak as a second S₁ location, and storing a location of the third peak as an S₂ location.

2. The method of claim 1, further comprising:

verifying that a first distance between the first peak and the third peak is smaller than a second distance between the third peak and the second peak.

3. The method of claim 1, further comprising:

locating a fourth peak and a fifth peak that may correspond to S₁ and S₂ between the first peak and the third peak; and

distinguishing the fourth peak and the fifth peak from murmurs using time domain kurtosis.

4. The method of claim 3, further comprising:

when the fourth peak and the fifth peak correspond to S₁ and S₂,

reducing a size of the first search window;

reducing a length between the first peak identified in the first search window and a beginning of the second search window; and

repeating the identifying a first peak, the identifying a second peak, and the identifying a third peak.

5. The method of claim 1, further comprising:

performing timing based error correction to verify that the second peak corresponds to S₁ and the third peak corresponds to S₂.

6. The method of claim 5, wherein performing timing based error correction comprises:

when a distance between the first peak and the second peak is not within a predetermined percentage of an average distance between two consecutive S₁ peaks, locating another peak corresponding to S₁, wherein a distance between the first peak and the another peak is within the predetermined percentage of the average distance, and identifying the another peak as the second peak.

7. The method of claim 5, wherein performing timing based error correction comprises:

when a distance between the third peak and a fourth peak corresponding to S₂ is not within a predetermined percentage of an average distance between two consecutive S₂ peaks, locating another peak corresponding to S₂, wherein a distance between the third peak and the another peak is within the predetermined percentage of the average distance, and identifying the another peak as the third peak.

8. The method of claim 5, wherein performing timing based error correction comprises:

when a distance between the first peak and the third peak is not within a predetermined percentage of an average distance between an S₁ peak and a subsequent S2 peak, locating another peak corresponding to S₂, wherein a distance between the first peak and the another peak is within the predetermined percentage of the average distance, and identifying the another peak as the third peak.

9. The method of claim 5, wherein performing timing based error correction comprises:

when a distance between the third peak and the second peak is not within a predetermined percentage of an average distance between an S₂ peak and a subsequent S₁ peak, locating another peak corresponding to S₁, wherein a distance between the third peak and the another peak is within the predetermined percentage of the average distance, and identifying the another peak as the second peak.

10. The method of claim 1, further comprising:

determining a heart rate using the stored S₁ locations.

11. The method of claim 1, wherein distinguishing the first peak from a murmur further comprises:

computing a first kurtosis of a segment of the audio signal that is a predetermined number of milliseconds on either side of the first peak;

computing a second kurtosis of a segment of the audio signal that is the predetermined number of milliseconds before the first peak; and

computing a third kurtosis of a segment of the audio signal that is the predetermined number of milliseconds after the second peak,

wherein when the first kurtosis is greater than a first predetermined value or an absolute difference between the second kurtosis and the third kurtosis is greater than a second predetermined value, the first peak is not a murmur.

12. The method of claim 1, wherein using frequency domain kurtosis further comprises:

computing a first kurtosis of a magnitude of a Fourier transform of a segment of the audio signal beginning at a location of the third peak; and

computing a second kurtosis of a magnitude of a Fourier transform of a segment of the audio signal beginning at a location of the another peak,

wherein a length of the segments is the nearest power of two to the length in samples that equals 50 ms, and

wherein when an absolute difference between a geometric mean of the first kurtosis and the second kurtosis and an arithmetic mean of the first kurtosis and the second kurtosis is greater than a predetermined value and the first kurtosis is greater than the second kurtosis, the another peak is determined to be a murmur peak.

13. The method of claim 1, wherein

distinguishing the first peak from a murmur further comprises:

identifying a murmur; and

storing a location of the murmur, and

wherein the location of the murmur and the stored S₁ locations and stored S₂ location are used to determine a type of the murmur.

14. The method of claim 5, wherein performing timing based error correction further comprises identifying a new peak as the second peak, wherein identifying the new peak as the second peak indicates a possible murmur interfering with the second peak.

15. The method of claim 5, wherein performing timing based error correction further comprises identifying a new peak as the third peak, wherein identifying the new peak as the third peak indicates a possible murmur interfering with the third peak.

16. A system comprising:

a processor;

a display operatively connected to the processor;

a memory operatively connected to the processor; and

instructions stored in the memory that are executable by the processor to identify heart sound components by;

receiving an audio signal comprising heart sounds;

identifying a first peak corresponding to S₁ within a first search window of the audio signal, wherein identifying the first peak comprises distinguishing the first peak from a murmur peak using time domain kurtosis;

identifying a second peak corresponding to S₁ within a second search window of the audio signal, wherein identifying the second peak comprises distinguishing the second peak from a murmur peak using time domain kurtosis;

identifying a third peak corresponding to S₂ between the first peak and the second peak, wherein identifying the third peak comprises using frequency domain kurtosis to determine whether another peak that may be the third peak is a murmur peak; and

storing a location of the first peak as a first S₁ location, storing a location of the second peak as a second S₁ location, and storing a location of the third peak as an S₂ location,

wherein the first S₁ location, the second S₁ location, and the S₂ location are shown in a phonocardiogram on the display.

17. The system of claim 16, wherein the instructions further identify heart sound components by:

performing timing based error correction to verify that the second peak corresponds to S₁ and the third peak corresponds to S₂.

18. The system of claim 16, wherein the system is one selected from a group consisting of a digital stethoscope, a personal computer, a laptop computer, a server, a mainframe, a personal digital assistant, a mobile phone, an iPod, and an MP3 player.

19. A computer readable medium storing instructions for identifying heart sound components, the instructions comprising functionality for:

receiving an audio signal comprising heart sounds;

identifying a first peak corresponding to S₁ within a first search window of the audio signal, wherein identifying the first peak comprises distinguishing the first peak from a murmur peak using time domain kurtosis;

identifying a second peak corresponding to S₁ within a second search window of the audio signal, wherein identifying the second peak comprises distinguishing the second peak from a murmur peak using time domain kurtosis;

identifying a third peak corresponding to S₂ between the first peak and the second peak, wherein identifying the third peak comprises using frequency domain kurtosis to determine whether another peak that may be the third peak is a murmur peak; and

storing a location of the first peak as a first S₁ location, storing a location of the second peak as a second S₁ location, and storing a location of the third peak as an S₂ location.

20. The computer readable medium of claim 19, wherein the instructions further comprise functionality for:

performing timing based error correction to verify that the second peak corresponds to S₁ and the third peak corresponds to S₂.