Advanced
A Localization Method for First and Second Heart Sounds Based on Energy Detection and Interval Regulation
A Localization Method for First and Second Heart Sounds Based on Energy Detection and Interval Regulation
Journal of Electrical Engineering and Technology. 2015. Sep, 10(5): 2126-2134
Copyright © 2015, The Korean Institute of Electrical Engineers
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : February 14, 2014
  • Accepted : May 12, 2015
  • Published : September 01, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Se Dong Min
Dept. of Medical IT Engineering, Soonchunhyang University, Korea. (sedongmin@sch.ac.kr)
Hangsik Shin
Corresponding Author: Dept. of Biomedical Engineering, Chonnam National University, Korea. (hangsik.shin@jnu.ac.kr)

Abstract
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.
Keywords
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)
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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) .
PPT Slide
Lager Image
PPT Slide
Lager Image
- 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.
PPT Slide
Lager Image
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)
PPT Slide
Lager Image
- 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
PPT Slide
Lager Image
, which is the Hilbert transform pair of x [ n ] . The analytical signal and envelope detection equations are described in (4) and (5), respectively.
PPT Slide
Lager Image
PPT Slide
Lager Image
- 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 ).
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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 .
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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
PPT Slide
Lager Image
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
PPT Slide
Lager Image
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.
References
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
Debbal SM , Bereksi-Reguig F 2008 Computerized heart sounds analysis Comput Biol Med 38 263 - 80    DOI : 10.1016/j.compbiomed.2007.09.006
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
Shah SJ , Michaels AD 2006 Hemodynamic correlates of the third heart sound and systolic time intervals Congest Heart Fail 12 (Suppl 1) 8 - 13    DOI : 10.1111/j.0889-7204.2006.05767.x
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
Nigam V , Priemer R 2005 Accessing heart dynamics to estimate durations of heart sounds Physiol Meas 26 1005 - 1018    DOI : 10.1088/0967-3334/26/6/010
Ricke A , Povinelli R , Johnson M 2005 Automatic Segmenatation of Heart Sound Signals Using Hidden Markov Models Computers in Cardiology 32 953 -
Messer S , Agzarian J , Abbott D 2001 Optimal wavelet denoising for phonocardiograms Microelectronics Journal 32 931 - 941    DOI : 10.1016/S0026-2692(01)00095-7
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