AUCTORES
Globalize your Research
Research Article | DOI: https://doi.org/10.31579/2641-0419/216
*Corresponding Author: S.M. DEBBAL, Genie -Biomedical Laboratory (GBM), Department of electronic Faculty of science engineering, university Aboubekr Belkaid, BP 119, Tlemcen 13000, Algeria.
Citation: S.M. DEBBAL. (2021) Heart Cardiac Sounds analysis using the Wigner distribution (WD) Method. J. Clinical Cardiology and Cardiovascular Interventions, 4(15); Doi:10.31579/2641-0419/216
Copyright: © 2021 S.M. DEBBAL, This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Received: 24 August 2021 | Accepted: 07 September 2021 | Published: 14 September 2021
Keywords: phonocardiogram; signal processing; sounds; heart cardiac’s sounds; time frequency analysis; wigner distribution
This paper is concerned a “The Wigner distribution (WD)” analysis of the Heart cardiac (or phonocardiogram signals: PCG). The Wigner distribution (WD) and the corresponding WVD (Wigner Ville Distribution) have shown good performances in the analysis of non-stationary and quantitative measurements of the time-frequency PCG signal characteristics. It is shown that these transforms provides enough features of the PCG signals that will help clinics to obtain diagnosis.
Heartbeat sound analysis by auscultation is still insufficient to diagnose some heart diseases. It does not enable the analyst to obtain both qualitative and quantitative characteristics of the phonocardiogram signals [1], [2]. Abnormal heartbeat sounds may contain, in addition to the first and second sounds, S1 and S2, murmurs and aberrations caused by different pathological conditions of the cardiovascular system [2]. Moreover, in studying the physical characteristics of heart sounds and human hearing, it is seen that the human ear is poorly suited for cardiac auscultation [3]. Therefore, clinic capabilities to diagnose heart sounds are limited.
The sound emitted by a human heart during a single cardiac cycle consist of two dominant events, known as the first heart sound S1 and the second heart sound S2 (Figure1), S1 relates to the closing of the mitral and tricuspid valves whilst S2 is generated by the halting of the aortic and pulmonary valves leaflets [1]. S1 corresponds in timing to the QRS complex in ECG (Electrocardiogram) and S2 follows the systolic pause in the normal cardiac cycle. Heartbeat sound analysis by auscultation only is still insufficient to diagnose some heart diseases. It does not enable the analyst to obtain both qualitative and quantitative characteristics of S1 and S2 of the phonocardiogram [1-2]. Moreover, in studying the physical characteristics of heart sounds and human hearing, it is seen that the human ear is poorly suited for cardiac auscultation [3]. Therefore, clinic capabilities to diagnose heart sounds are limited. In this paper we are interested in the study of the physical characteristics of the second heart sound S2 which consist of two major components in the spectrum of the signal. One of these components A2 is due to the closure of the aortic valve and the other P2 is due to the closure of the pulmonary valve.
The aortic component is loudest than the pulmonic component. It is discernible at all the auscultation sites.It is best heard at the right base, with the diaphragm of the chest piece firmly pressed, whereas the pulmonic on a may only be audible at the left base, with the diaphragm of the chest piece firmly pressed.
The aortic component has higher frequency contents and generally precedes the pulmonary component because in normal heart activity the aortic valve closes before the pulmonary valve. The difference of time occurrence between these valves activities is known in medical community [4 - 6] as split. However in many diseases [1] this order of time occurrence may be reverse and its delay varies.
The characteristics of the PCG signal and other features such as heart sounds S1 and S2 location; the number of components for each sound; their frequency content; their time interval; all can be measured more accurately by digital signal processing techniques.
The Wigner distribution (WD) and the corresponding WVD (Wigner Ville Distribution) have shown good performances in the analysis of non-stationary signals. This comes from the ability of the WD to separate signals in both time and frequency directions. One advantage of the WD over the STFT is that it does not suffer from the time-frequency trade-off problem. On the other hand, the WD has a disadvantage since it shows cross-terms in its response. These cross-terms are due to the nonlinear behaviour of the WD, and bear no physical meaning. One way to remove these cross-terms is by smoothing the time-frequency plane, but this will be at the expense of decreased resolution in both time and frequency [7].
The WD was applied to heart sound signal it shows no success in displaying or separating the signal components in both the time and frequency direction [6], although it provides high time-and frequency- resolution in simple monocomponent signal analysis [8]
The WD of two signals, x(t), y(t), is defined via,
The auto WD is obtained when x(t) = y(t)=s(t), it is a bilinear function of the signal s(t). The WD, and can also be expressed as
Where t and are respectively, the time and frequency variables, and * denotes the complex conjugate.
The WD had shown good applications in the analysis of non-stationary signal [9-11]. This comes from the ability of this method to separate signals in both time and frequency directions. The WD has a disadvantage that it is limited by the appearance of cross-terms. These cross-terms are due to the nonlinear behaviour of the WD, and bear no physical meaning. One way to remove these cross-terms is by smoothing the time-frequency plane [7], but this will be at the expense of decreased resolution in both time and frequency.
3. Results and discussion of the Wigner distribution application
Figure2 shows the WD application of the cardiac sound of the normal case (2a), the coarctation of the aorta (2b) and the innocent murmur (2c). One can notice here that the two principal components (A2 and P2) start to appear in the presence cross-terms. The WD results may be improved by increasing the sampling rate of original signal,, but it still suffers from the cross-terms problem because of the nonlinearity of the WD.
However the Wigner distribution (WD) has shown good performances in the analysis of non-stationary signals. This comes from the ability of the WD to separate signals in both time and frequency directions. One advantage of the WD over the STFT is that it does not suffer from the time-frequency trade-off problem. On the other hand, the WD has a disadvantage since it shows cross-terms in its response. These cross-terms are due to the nonlinear behaviour of the WD, and bear no physical meaning. One way to remove these cross-terms is by smoothing the time-frequency plane, but this will be at the expense of decreased resolution in both time and frequency
Figure2 shows the WD application of the cardiac sound of the normal case (2a), the coarctation of the aorta (2b) and the innocent murmur (2c). One can notice here that the two principal components (A2 and P2) start to appear in the presence cross-terms. The WD results may be improved by increasing the sampling rate of original signal,, but it still suffers from the cross-terms problem because of the nonlinearity of the WD.
However the Wigner distribution (WD) has shown good performances in the analysis of non-stationary signals. This comes from the ability of the WD to separate signals in both time and frequency directions. One advantage of the WD over the STFT is that it does not suffer from the time-frequency trade-off problem. On the other hand, the WD has a disadvantage since it shows cross-terms in its response. These cross-terms are due to the nonlinear behaviour of the WD, and bear no physical meaning. One way to remove these cross-terms is by smoothing the time-frequency plane, but this will be at the expense of decreased resolution in both time and frequency
The cardiac (heartbeat sound) cycle of phonocardiogram (PCG) is characterized by transients and fast changes in frequency as time progresses. It was shown that basic frequency content of PCG signal can be easily provided using FFT technique. However, time duration and transient variation cannot be resolved; the Wigner Distribution (WD) is a suitable technique to analyze such a signal. It was also shown that the graphic representation that provides a quantitative analysis simultaneously in time and frequency. It is therefore very helpful in extracting clinically useful information.
The measurement of the time difference between the A2 and P2 components in the sound S2, the number of major components of the sounds S1 and S2 and the frequency range and duration for all these components and sounds can be accurately achieved for Wigner Distribution (WD) simultaneously as was clearly illustrated.
It is shown that the WD Technic provides more information of the PCG signals with murmurs that will help physicians to obtain qualitative and quantitative measurements of the time and the time-frequency PCG signal characteristics and consequently aid to diagnosis.