The present study suggests a localization method for the first (
S
1
) and the second (
S
2
) feature of heart sounds, based on an algorithm involving frequency filtering, energy detection, and interval regulation. Localization accuracy was evaluated by comparing the algorithm with the traditional Hilbert transform-based localization method. Results show that the sensitivity and the positive predictivity value of proposed method, respectively, were 97.27 % and 99.94 % in
S
1
detection and 94.99 % and 100 % in
S
2
detection.
1. Introduction
Auscultation, the action of listening to the internal sounds of the body using a stethoscope, is a simple and effective way to diagnosis of cardiac dysfunction and cardiovascular disease. Moreover, auscultation has been widely used for fundamental cardiac diagnoses because of its non-invasive and simple diagnostic characteristics. The major purpose of auscultation is listening to the heart sound caused by the closing of the heart valves. Heart sound contains many signals related to magnitude and frequency content, splitting, and the presence of additional systolic or diastolic murmur components. There are two dominant features of heart sounds: the first heart sound (
S
1
) and the second heart sound (
S
2
).
S
1
reflects the closure of mitral and the tricuspid valves, and
S
2
marks the closure of the aortic and pulmonary valves.
Heart sounds are usually investigated with
S
1
and
S
2
analysis. Because
S
1
and
S
2
occur at the onset of ventricular contraction and the end of ventricular systole, respectively, the left ventricular ejection time (LVET) and the preejection period (PEP) have been estimated by heart sound analysis
[1]
. Moreover, the diagnosis of valve dysfunction and hemodynamic research are widely known applications of heart sound analysis. For example, cardiac valve disorder has been classified by a neural network classifier
[2]
and blood pressure has been estimated with pattern analysis
[3]
using heart sounds. Furthermore, the relationship between
S
2
and aortic blood pressure has been shown
[4]
. Although heart sounds are a useful tool clinically, they are hard to analyze from the non-stationary characteristic. That is, making a cardiac disorder diagnosis based on heart sounds is still a considerable challenge
[5
-
12]
. Thus, the accurate localization of
S
1
and
S
2
is an important issue for heart sound analysis and clinical applications.
Previous localization methods for heart sounds can be summarized as filtering, transform-based de-noising, and feature classification. In particular, wavelet de-noising has been researched to enable accurate heart sound localization
[13
-
15]
. There have been several approaches to
S
1
and
S
2
localization, such as spectral tracking
[16]
, Shannon energy
[17
,
18]
, boundary modification method
[19]
, time-delayed neural network
[20]
, ECG-gated method
[7]
, entropy based detection
[21]
, homomorphic filtering
[8]
, complexity based segmentation
[22]
, hidden Markov model
[23]
, and wavelet transform
[24
,
25]
. Recently, wavelet de-noising and ECG-gated peak-peaking methods have been improved by an ensemble average method for more accurate localization
[26]
. Also, in 2012,
S
1
/
S
2
heart sounds detection based on music beat tracking has been researched
[27]
.
Most heart sound localization methods are based on the energy of the features. A problem with energy-based detection methods is the low amplitude of
S
1
and
S
2
components; they can easily be buried in background noise, murmurs, or lung sounds, and they are not reflected in the energy envelope of the heart sound
[22]
. Moreover, it has been shown that the amplitude of the heart sound waveform is changed by respiration
[28]
; the amplitude of the heart sound waveform is attenuated by inspiration. The ECG-gated localization method, a commonly used technique for the localization, is based on ECG features, such as the R-wave. ECG-gated localization provides relatively stable localization
[26]
. However, ECG should also be measured as a reference in the ECG-gated localization method. Moreover, for a variety of pathological conditions, there may be detection errors because of inconsistencies between the electrical and mechanical activities of the heart
[29]
.
The aim of the present study is the localization of
S
1
and
S
2
in heart sounds. In this study, features were discriminated from a phonocardiogram only. To improve localization performance, cascade integrator comb (CIC) filtering and an energy- and temporal-based localization method were developed. An experiment was performed in spontaneous and a 0.1 Hz or 0.25 Hz respiration-controlled environment to produce practical respiration effects and distortion. Statistical methods were used for a quantitative evaluation of the proposed algorithm.
2. Experimental Setup
- 2.1 Dataset
Eighteen young, healthy subjects (11 males, 7 females, mean age is 24.1, range 17-30 years old) participated. Prior to the experiment, the subjects were requested to provide information on their physical condition. None reported any taking medicine or signs of cardiovascular or respiratory disease. Each experiment was performed in a typical sports medicine laboratory at ambient room temperature, from 11 am to 6 pm. Drinking and smoking were prohibited during 24 hours and 2 hours before experiments, respectively. Anthropometric data, such as height and weight, were measured and are summarized in
Table 1
.
Subject demographic data (aSystolic Blood Pressure,bDiastolic blood pressure,cHeart rate (beats per minute)
Subject demographic data (a Systolic Blood Pressure, b Diastolic blood pressure, c Heart rate (beats per minute)
- 2.2 Devices and software
Heart sounds were recorded with a MP150 (Biopac Systems, inc., CA, USA), DA 100C acquisition module, and SS17L contact microphone. Heart sounds were measured on the chest surface between the 3
rd
and 4
th
ribs, and the ECG was measured on the chest surface with a Lead II configuration (
Fig. 1
). Acknowledge 3.9.1 (BIOPAC, USA) was used for real-time monitoring and data storage. For signal conditioning, localization, and statistical analyses, MATLAB 2008b (Mathworks, Inc., Natick, MA, USA) was used and an HEM-907(Omron, Kyoto, Japan) was used for blood pressure measurements.
Location of heart auscultation and ECG measurements
- 2.3 Experimental protocol
Heart sounds were measured, including spontaneous and controlled respiration effects. Respiration was controlled at 0.1 Hz (low rate) and 0.25 Hz (high rate) in a supine position using a metronome. In the sitting position, there was no respiration controlling. Experiments were carried out sequentially in the following order: spontaneous, 0.1 Hz controlled, and 0.25 Hz controlled respiration in the supine position, and spontaneous respiration in the sitting position. Subjects had a 2-min rest period between experiments. Data were sampled at a 1 kHz sampling rate for a 5-min period.
3. Signal Conditioning and Localization
- 3.1 Cascade integrator comb filter
CIC filtering has been used for biomedical signal conditioning and feature enhancement
[29]
. A CIC filter, which is a linear-phase filter, was implemented for signal conditioning of the heart sounds. Linear phase-filtering, which prevents phase and timing distortion, is important because the heart sound has non-stationary characteristic. A CIC filter is composed of a high-pass and a low-pass filter and its transfer function is described in (1) and (2). N
hpf
, N
lpf
, and M indicate the high-pass filter order, low-pass filter order, and number of cascaded stages, respectively; their values were 16, 8, and 2, respectively, in the filter design. N
hpf
, N
lpf
, and M had a 43.2 Hz center frequency and pass-band between 26.6 Hz and 63 Hz (3 dB point). The frequency characteristics of the CIC filter are presented in
Fig. 2(a)
.
- 3.2 Moving average filter
A general heart sound is described with a high-frequency signal. However, specific features, such as
S
1
and
S
2
, are decided by the value of the envelope. From the natural characteristics of heart sounds, the envelope detected by the Hilbert transform contains many tiny fluctuations that could have an effect on specific feature extraction. These fluctuations are also regarded as high frequency noise in feature extraction; thus, a moving-average filter (MAF) was used for envelope smoothing as a kind of low-pass filter. The frequency response of MAF is shown in (3). L is the MAF length and it was of 0.15 s duration (cut-off frequency is ~3 Hz).
Fig. 2(b)
shows the frequency characteristics of MAF.
Frequency characteristic of filters: (a) Cascade integrator comb filter. The center frequency is 43.2 Hz and the pass-band is 26.6 to 63 Hz; (b) Moving average filter (cut-off frequency is ~3 Hz)
- 3.3 Hilbert transform
The Hilbert transform is proper for time-varying signal analysis, used frequently for envelope detection with heart sounds
[31]
. The analytic function of
z
[
n
] contains a real-value function
x
[
n
] and a complex value function
, which is the Hilbert transform pair of
x
[
n
] . The analytical signal and envelope detection equations are described in (4) and (5), respectively.
- 3.4S1andS2Localization
The heart sound signal is processed sequentially, with CIC filtering, the Hilbert transformation, and movingaverage filtering. Then, the
S
1
and
S
2
peaks are detected from the filtered envelope of the heart sound.
S
1
peaks were detected by a zero-crossing method, based on signal derivation. For
S
2
localization, an energy based detection method was used within the
S
1
intervals. The
S
1
-gated window was determined from the (n-1)
th
position of
S
1
(PS1
n-1
) to the n
th
position of
S
1
(PS1
n
). A 0.2-s refractory period was also defined to reduce misdetection for the former and latter
S
1
location. Thus, the scanning window had a range between PS1
n-1
+ 0.2 and PS1
n
- 0.2. Within the scanning window,
S
2
was localized by taking the maximum value of the energy waveform of the heart sound.
The reference peak position was detected by an ECG-gated detection method and the Hilbert transformation. The
S
1
location was decided within a 250 ms detection window (P1) following the R peak of the ECG waveform. The
S
2
location was detected by the maximum energy within the specified range (P
2
). P
2
starts from the
S
1
location to the next R peak of the ECG. Each detected
S
1
and
S
2
was corrected manually (
Fig. 3
).
Localization method of the reference S1 and S2 positions
- 3.5 Peak correction
The envelope of the heart sound has two distinguishable peaks, at
S
1
and
S
2
. In using time difference-based classification method, the detection of both
S
1
and
S
2
should precede
S
1
and
S
2
classification. The detected peaks are typically located at the
S
1
and
S
2
positions. However, redundant peaks or missed peaks may still remain. In the pre-processing stage of
S
1
localization, mis-detected peak correction is directly related to detection accuracy. The proposed localization method is based on time difference characteristics between
S
1
and
S
2
. Because the
S
2
-
S
1
interval is generally longer than the
S
1
-
S
2
interval, the difference in peak time interval (PTI) can be described with an up-and-down curve. If we assume that every peak was detected correctly, then the PTI curve should be repeated with a regular up-and-down pattern. However, a PTI curve could show an irregular trend from missed or redundant peaks. In the missed peak case, PTI becomes higher than the normal upper limit (6) and redundant peaks causes a lower PTI value than the other lower peaks (7). To remove these irregularities, a peak correction method was used. Peak correction was carried out with two processes: False Positive Peak correction (FPPC) and False Negative Peak correction (FNPC). FPP correction modifies PTI values using (8); a redundant peak time interval at the nth PTI, PTIn, is removed by comparison between the former and the latter PTI value. In the FNP correction, in contrast, a new PTI is added using (9). The estimated
S
1
and
S
2
intervals are decided by average values of 50% higher or lower.
Moreover, both
S
1
and
S
2
should be detected after the signal conditioning procedures and classified with PTI analysis for
S
1
localization.
Fig. 4
shows the PTI analysis procedure. PTI can be represented with an up-and-down plot in the general case, and the up-to-down breakpoint indicates the
S
1
location.
Peak time interval (PTI) analysis and S1 localization: (a) Detected peaks and intervals between S1 and S2; (b) PTI and S1 localization. S1 is located at the up-to-down breakpoint.
- 3.6 Replace, add, and delete algorithm
As noted previously, a representative detection problem can be summarized with FPP and FNP. In
S
1
classification, a similar problem occurs. Typically, FNP means missed peaks or confusing
S
1
and
S
2
. Moreover, FPP means that there are redundant peaks around the actual
S
1
location. To improve this detection problem, a replace, add, and delete (RAD) algorithm is proposed.
Replace, add, and delete procedures are defined by PTI range. In the replace procedure, when
S
2
is detected as a substitute of
S
1
, PTI is between
Limit
h1
and
Limit
h2
(
Limit
h1
< PTI <
Limit
h2
). If PTI is longer than the replace interval (
Limith
2
< PTI), it is regarded as an add procedure. The delete procedure is invoked when PTI is shorter than Limitl (
Limit
l
> PTI). Classification criteria,
Limit
h1
,
Limit
h2
, and
Limit
l
, are determined using (10), (11) and (12), respectively, with empirical constant, and they are adaptively changed by the previous
S
1
interval.
Fig. 5
shows three detection problems and correction methods such as ‘Replace’, ‘Add’ and ‘Delete’ algorithm. In
Fig. 5
, dotted box shows the detection problem with
S
1
detection.
Fig. 5(a)
indicates a false detection of
S
1
when PTI is between
Limit
h1
and
Limit
h2
.
Fig. 5(b)
shows the corrected result of the ‘replace’ procedure for
(a)
. A new
S
1
is detected with the mean beat rate (MBR) and detection range (DR). MBR is decided with the average value of the previous
S
1
time difference and DR is 0.2 s, empirically. The second detection problem is shown in
Fig. 5(c)
. A
S
1
is missed when PTI is longer than
Limit
h2
.
Fig. 5(d)
is the result of the ‘add’ procedure. In add procedure,
S
1
is added by the same rule of replace method. However, in add procedure, the next
S
1
would not be deleted because a
S
1
is missed. The last detection problem is represented in
Fig. 5(e)
.
S
1
is detected again when PTI is shorter than
Limit
l
.
Fig. 5(f)
shows the result of the ‘delete’ procedure. In this procedure, the latter
S
1
is deleted between two adjacent
S
1
, and the former
S
1
is retained.
S1 detection problems and results of the Replace, Add, and Delete (RAD) algorithm. A solid circle indicates a newly added S1 and a dotted circle indicates a deleted S1: (a) False detection of S1; (b) The correction result by the ‘Replace’ procedure of (a); (c) Detection failure of S1; (d) The corrected result by the ‘Add’ procedure of (c); (e) Duplicated S1 detection; (f) The result of the ‘Delete’ procedure of (e).
Fig. 6
shows a flowchart of the RAD algorithm and
S
2
detection. The RAD algorithm is repeated until every
S
1
is corrected. The
S
2
localization depends on
S
1
detection and it is represented in
Fig. 5(b)
.
S
2
is determined by the maximum value position between two adjacent
S
1
.
Flowchart of the RAD algorithm and S2 detection: (a) RAD algorithm for S1 correction; (b) S2 localization procedure. MBR = mean beat rate, PI = peak interval, PPn = n-th peak position.
4. Evaluation
For localization performance evaluation, the true positive (TP), false positive (FP), false negative (FN) and failed detection (FD) were calculated. Then, we derived the sensitivity (SE), positive predictivity (+P), and failed detection rate (FDR) to evaluate the accuracy of the detection algorithm. Detected peaks were verified with a 200-ms approval range from the reference peak positions that were obtained by ECG-gated detection. Sensitivity was calculated using (13), +P using (14), and FDR using (15). Evaluations were carried out for both
S
1
and
S
2
localizations.
5. Results
- 5.1 Signal conditioning
The signal conditioning results are shown in
Fig. 7
.
Fig. 7(a)
represents the original heart sound waveform, and
(b)
,
(c)
, and
(d)
are the CIC-filtered signal, Hilbert-transformed signal, and moving-average-filtered signal, respectively. In
Fig. 7(d)
, the detected
S
1
and
S
2
peaks are described.
Signal conditioning procedure and digital filtering: (a) Original heart sound; (b) CIC-filtered result (M = 2, Nlpf = 16, Nhpf = 32); (c) Hilbert-transformation result; (d) Moving average integrator result (window size=0.15 s). ‘▼’ indicates a detected peak.
- 5.2 Peak correction and classification
The initially detected peaks after signal conditioning are processed with the FPPC and FNPC methods. The FPPC and FNPC results are represented in
Fig. 8
. Initially detected peaks are shown in
Fig. 8(a)
, and
Fig. 8(b)
shows the PTI of
(a)
. FPP and FNP could be detected with detection criteria (1) and (2), respectively, described as relatively low and high PTI values in
Fig. 8(b)
. The corrected peak and its PTI are presented in
Figs. 8(c)
and
8(d)
, respectively. After peak correction, the
S
1
and
S
2
locations are decided at the up-or-down peaks of the PTI curve.
Correction and classification of the detected peaks after signal conditioning: (a) Detected peaks by differential method; (b) Peak time intervals of (a). False positive peaks (FPP) and false negative peaks (FNP) were detected in this stage; (c) Corrected peaks. (d) Peak time intervals of (c). A longer time interval (UP peak) means a S1 and shorter time interval (Down peak) means a S2 peak.
- 5.3 Replace, Add, and delete correction
Classified
S
1
peaks were also corrected with the RAD algorithm.
Fig. 9
shows an example of
S
1
correction with the RAD algorithm.
Fig. 9(a)
and
(b)
show the initially detected
S
1
and its PTI. RAD correction was carried out using (8), (9), and (10).
Figs. 9(c)
and
(d)
,
Figs. 9(e)
and
(f)
, and
Figs. 9(g)
and
(h)
show the replace, add, and delete procedures and the PTI, respectively.
An example of the S1 correction method: (a) S1 detection result before RAD algorithm. There is a misdetected peak between the 3rd and 5th peaks; (b) Peak time interval (PTI) of (a). It was found to be a PTI>LimitHIGH1 peak; (c) ‘Replace’ procedure; (d) PTI of (c); (e) ‘Add’ procedure; (f) PTI of (e); (g) ‘Delete’ procedure; (h) PTI of (g).
- 5.4 Statistical evaluation
In the statistical analysis, SE, +P, and FDR were calculated to evaluate the classification performance at each stage. The results of
S
1
and
S
2
localization are presented in
Tables 2
and
3
. In
S
1
localization, SE was increased, 87.57% to 97.27%, and +P was increased, 99.72% to 99.94%. FDR decreased, 14.48% to 2.86%. Except for slight fluctuation in +P, each performance evaluation index was increased in every correction stage.
Results evaluating theS1detection algorithm at each stage. The method consists of CIC filtering, peak correction (PC), and the replace, add, and delete (RAD) algorithm
a Total number of reference S1 b Total number of detected peaks
Results evaluating theS2detection algorithm at each stage. The method consists of CIC filtering, peak correction (PC), and the replace, add, and delete (RAD) algorithm
a Total number of reference S2 b Total number of detected peaks
Performance for
S
2
localization was evaluated in the same way as
S
1
. In
S
2
localization, SE increased, 83.94% to 94.9%, and +P was 100% in each case. FDR decreased, 19.13% to 5.27%. Each performance evaluation index of
S
2
also increased in every correction stage. Overall localization performance for
S
1
was better than for
S
2
except in terms of +P.
6. Conclusion
- 6.1 Performance evaluation
The proposed method was validated by a quantitative study with 27,082 pulses, posture changes, and respiratory variation. Localization performance was most improved by CIC filtering; SE and FDR were improved by 6.21% and 7.8%, respectively, in the
S
1
case. In the
S
2
case, they were also improved, by 6.1% and 8.07%. +P was not changed significantly. Remaining misdetected peaks after frequency domain-filtering were corrected by peak correction and the RAD algorithm. These algorithms improved FDR by 1.05% and 1.95% in the
S
1
case. These results could be considered as slight improvement; however, it is meaningful to recover from missing or extremely distorted peaks. Because
S
2
detection depends largely on
S
1
results,
S
2
detection performance had a lower value than
S
1
results; however, +P was always 100%.
- 6.2 Limitation
A time-dependent classification algorithm is hard to adapt to an irregular heartbeat. For example, the heart sound of an arrhythmia patient could generate more errors because a correction criterion is affected by the previous PTI and specific constant. An error could even be propagated in certain cases.
In
S
2
localization, the
S
2
detection rate depends on the
S
1
detection result. Thus, failed detection of
S
1
could be a reason of the failure of
S
2
localization. The proposed algorithm focuses on a post-peak correction method of heart sounds; however, it does not provide a noise reduction method for the inherent problems of heart sounds. Major heart sound problems, such as lung sounds, murmurs, respiration effects, and signal distortion in female subjects, were not addressed significantly in this paper.
- 6.3 Concluding remarks
In this study, only heart sound data were used to localize
S
1
and
S
2
features, with no other reference signal. This is a novel study in heart sound localization because it may provide an efficient and accurate method for heart sound applications. From the research presented, it appears that the proposed first and second heart sound localization method can improve localization performance. It may be a useful method for robust peak detection in various applications, such as LVET/PEP detection and cardiac disorder diagnoses. Although there are still unsolved problem related to abnormal heart activity and optimization for robust detection, we believe that the proposed method can provide simple and reliable heart sound feature detection.
Acknowledgements
This work was supported by the Soonchunhyang University Research Fund (No. 20120676)This research was supported by the MSIP(Ministry of Science, ICT and Future Planning), Korea, under the CITRC(Convergence Information Technology Research Center) (IITP-2015-IITP-2015-H8601-15-1009) supervised by the IITP(Institute for Information & communications Technology Promotion)
BIO
Se Dong Min He received the M.S. and Ph.D. degrees in electrical and electronic engineering from the Department of Electrical and Electronics Engineering, Yonsei University, Seoul, in 2004 and 2010, respectively. He is currently an Assistant Professor at the Department of Medical IT Engineering, Soonchunhyang University, Asan, Korea. His research area includes biomedical signal processing, healthcare sensor application, and mobile healthcare technologies.
Hangsik Shin is the Assistant Professor at the Chonnam National University. He received the Ph.D. degree in electrical and electronic engineering from the departent of Electrical and Electronics Engineering, Yonsei University, in 2010. His research area include biomedical signal processing, physiological modeling, u-Healthcare and mobile Healthcare Technologies.
Lewis RP
,
Rittogers SE
,
Froester WF
,
Boudoulas H
1977
A critical review of the systolic time intervals
Circulation
56
146 -
58
DOI : 10.1161/01.CIR.56.2.146
Babaei S
,
Geranmayeh A
2009
Heart sound reproduction based on neural network classification of cardiac valve disorders using wavelet transforms of PCG signals
Comput Biol Med
39
8 -
15
DOI : 10.1016/j.compbiomed.2008.10.004
Bartels A
,
Harder D
1992
Non-invasive determination of systolic blood pressure by heart sound pattern analysis
Clin Phys Physiol Meas
13
249 -
56
DOI : 10.1088/0143-0815/13/3/004
Zhang XY
,
MacPherson E
,
Zhang YT
2008
Relations between the timing of the second heart sound and aortic blood pressure
IEEE Trans Biomed Eng
55
1291 -
7
DOI : 10.1109/TBME.2007.912422
Debbal SM
,
Bereksi-Reguig F
2007
Automatic measure of the split in the second cardiac sound by using the wavelet transform technique
Comput Biol Med
37
269 -
76
DOI : 10.1016/j.compbiomed.2006.01.005
El-Segaier M
,
Lilja O
,
Lukkarinen S
,
Sornmo L
,
Sepponen R
,
Pesonen E
2005
Computer-based detection and analysis of heart sound and murmur
Ann Biomed Eng
33
937 -
42
DOI : 10.1007/s10439-005-4053-3
Gupta C
,
Palaniappan R
,
Swaminathan S
,
Krishnan S
2007
Neural network classification of homomorphic segmented heart sounds
Applied Soft Computing Journal
7
286 -
297
DOI : 10.1016/j.asoc.2005.06.006
Hult P
,
Fjallbrant T
,
Hilden K
,
Dahlstrom U
,
Wranne B
,
Ask P
2005
Detection of the third heart sound using a tailored wavelet approach: method verification
Med Biol Eng Comput
43
212 -
7
DOI : 10.1007/BF02345957
Torry J
2003
Heart sound analysis comparing wavelet and autoregressive techniques
Computers in Cardiology
657 -
660
Turkoglu I
,
Arslan A
,
Ilkay E
2003
An intelligent system for diagnosis of the heart valve diseases with wavelet packet neural networks
Comput Biol Med
33
319 -
31
DOI : 10.1016/S0010-4825(03)00002-7
Voss A
,
Mix A
,
Hubner T
2005
Diagnosing aortic valve stenosis by parameter extraction of heart sound signals
Ann Biomed Eng
33
1167 -
74
DOI : 10.1007/s10439-005-5347-x
Ahlstrom C
,
Johansson A
,
Uhlin F
,
Lanne T
,
Ask P
2005
Noninvasive investigation of blood pressure changes using the pulse wave transit time: a novel approach in the monitoring of hemodialysis patients
J Artif Organs
8
192 -
7
DOI : 10.1007/s10047-005-0301-4
Johansson A
,
Ahlstrom C
,
Lanne T
,
Ask P
2006
Pulse wave transit time for monitoring respiration rate
Med Biol Eng Comput
44
471 -
8
DOI : 10.1007/s11517-006-0064-y
Iwata A
,
Ishii N
,
Suzumura N
,
Ikegaya K
1980
Algorithm for detecting the first and the second heart sounds by spectral tracking
Med Biol Eng Comput
18
19 -
26
DOI : 10.1007/BF02442475
Liang H
,
Lukkarinen S
,
Hartimo I
1997
Heart sound segmentation algorithm based on heart sound envelogram
Computers in Cardiology
105 -
108
Moukadem A.
,
Dieterlen A.
,
Hueber N.
,
Brandt C.
2011
Comparative study of heart sounds localization
Proc. SPIE 8068, Bioelectronics, Biomedical, and Bioinspired Systems V; and Nanotechnology V, 80680P
99 -
Liang H
,
Lukkarinen S
,
Hartimo I
1998
A boundary modification method for heart sound segmentation algorithm
Computers in Cardiology
593 -
595
Oskiper T
,
Watrous R
2002
Detection of the first heart sound using a time-delay neural network
Computers in Cardiology
537 -
540
Yadollahi A
,
Moussavi Z
2006
A robust method for heart sounds localization using lung sounds entropy
IEEE Transactions on Biomedical Engineering
53
497 -
502
DOI : 10.1109/TBME.2005.869789
Ricke A
,
Povinelli R
,
Johnson M
2005
Automatic Segmenatation of Heart Sound Signals Using Hidden Markov Models
Computers in Cardiology
32
953 -
Wang P
,
Kim Y
,
Ling L
,
Soh C
2005
First heart sound detection for phonocardiogram segmentation
Conf Proc IEEE Eng Med Biol Soc
5
5519 -
22
Ahlstrom C
,
Lanne T
,
Ask P
,
Johansson A
2008
A method for accurate localization of the first heart sound and possible applications
Physiol Meas
29
417 -
28
DOI : 10.1088/0967-3334/29/3/011
Barabasa C
,
Jafari M
2012
A Robust Method for S1/S2 Heart Sounds Detection Without ECG Reference Based on Music Beat Tracking, Electronics and Telecommunications (ISETC)
2012 10th International Symposium
307 -
310
Amit G
,
Shukha K
,
Gavriely N
,
Intrator N
2009
Respiratory modulation of heart sound morphology
Am J Physiol Heart Circ Physiol
296
H796 -
805
DOI : 10.1152/ajpheart.00806.2008
Haghighi-Mood A
,
Torry J
1995
A sub-band energy tracking algorithm for heart sound segmentation
Computers in Cardiology
501 -
504
Kohama T
,
Nakamura S
,
Hoshino H
1999
An efficient RR interval detection for ECG monitoring system
IEICE Transactions on Information and Systems E Series D
82
1425 -
1432
Baranek H
,
Lee H
,
Cloutier G
,
Durand L
1989
Automatic detection of sounds and murmurs in patients with lonescu-Shiley aortic bioprostheses
Medical and Biological Engineering and Computing
27
449 -
455
DOI : 10.1007/BF02441460