Bayesian denoising framework of phonocardiogram based on a new dynamical model

In this paper, we introduce a model-based Bayesian denoising framework for phonocardiogram (PCG) signals. The denoising framework is founded on a new dynamical model for PCG, which is capable of generating realistic synthetic PCG signals. The introduced dynamical model is based on PCG morphology and is inspired by electrocardiogram (ECG) dynamical model proposed by McSharry et al. and can represent various morphologies of normal PCG signals. The extended Kalman smoother (EKS) is the Bayesian filter that is used in this study. In order to facilitate the adaptation of the denoising framework to each input PCG signal, the parameters are selected automatically from the input signal itself. This approach is evaluated on several PCGs recorded on healthy subjects, while artificial white Gaussian noise is added to each signal, and the SNR and morphology of the outputs of the proposed denoising approach are compared with the outputs of the wavelet denoising (WD) method. The results of the EKS demonstrate better performance than WD over a wide range of PCG SNRs. The new PCG dynamical model can also be employed to develop other model-based processing frameworks such as heart sound segmentation and compression.     

A robust heart sounds segmentation module based on S-transform

This paper presents a new module for heart sounds segmentation based on S-transform. The heart sounds segmentation process segments the PhonoCardioGram (PCG) signal into four parts: S1 (first heart sound), systole, S2 (second heart sound) and diastole. It can be considered one of the most important phases in the auto-analysis of PCG signals. The proposed segmentation module can be divided into three main blocks: localization of heart sounds, boundaries detection of the localized heart sounds and classification block to distinguish between S1 and S2. An original localization method of heart sounds are proposed in this study. The method named SSE calculates the Shannon energy of the local spectrum calculated by the S-transform for each sample of the heart sound signal. The second block contains a novel approach for the boundaries detection of S1 and S2. The energy concentrations of the S-transform of localized sounds are optimized by using a window width optimization algorithm. Then the SSE envelope is recalculated and a local adaptive threshold is applied to refine the estimated boundaries. To distinguish between S1 and S2, a feature extraction method based on the singular value decomposition (SVD) of the S-matrix is applied in this study. The proposed segmentation module is evaluated at each block according to a database of 80 sounds, including 40 sounds with cardiac pathologies.     

Hemodynamics-driven mathematical model of first and second heart sound generation

We propose a novel, two-degree of freedom mathematical model of mechanical vibrations of the heart that generates heart sounds in CircAdapt, a complete real-time model of the cardiovascular system. Heart sounds during rest, exercise, biventricular (BiVHF), left ventricular (LVHF) and right ventricular heart failure (RVHF) were simulated to examine model functionality in various conditions. Simulated and experimental heart sound components showed both qualitative and quantitative agreements in terms of heart sound morphology, frequency, and timing. Rate of left ventricular pressure (LV dp/dt_max) and first heart sound (S1) amplitude were proportional with exercise level. The relation of the second heart sound (S2) amplitude with exercise level was less significant. BiVHF resulted in amplitude reduction of S1. LVHF resulted in reverse splitting of S2 and an amplitude reduction of only the left-sided heart sound components, whereas RVHF resulted in a prolonged splitting of S2 and only a mild amplitude reduction of the right-sided heart sound components. In conclusion, our hemodynamics-driven mathematical model provides fast and realistic simulations of heart sounds under various conditions and may be helpful to find new indicators for diagnosis and prognosis of cardiac diseases.New & noteworthyTo the best of our knowledge, this is the first hemodynamic-based heart sound generation model embedded in a complete real-time computational model of the cardiovascular system. Simulated heart sounds are similar to experimental and clinical measurements, both quantitatively and qualitatively. Our model can be used to investigate the relationships between heart sound acoustic features and hemodynamic factors/anatomical parameters.     

Segmentation of heart sounds based on dynamic clustering

The heart sound signal is first separated into cycles, where the cycle detection is based on an instantaneous cycle frequency. The heart sound data of one cardiac cycle can be decomposed into a number of atoms characterized by timing delay, frequency, amplitude, time width and phase. To segment heart sounds, we made a hypothesis that the atoms of a heart sound congregate as a cluster in time–frequency domains. We propose an atom density function to indicate clusters. To suppress clusters of murmurs and noise, weighted density function by atom energy is further proposed to improve the segmentation of heart sounds. Therefore, heart sounds are indicated by the hybrid analysis of clustering and medical knowledge. The segmentation scheme is automatic and no reference signal is needed. Twenty-six subjects, including 3 normal and 23 abnormal subjects, were tested for heart sound signals in various clinical cases. Our statistics show that the segmentation was successful for signals collected from normal subjects and patients with moderate murmurs.     

Heart energy signature spectrogram for cardiovascular diagnosis

A new method and application is proposed to characterize intensity and pitch of human heart sounds and murmurs. Using recorded heart sounds from the library of one of the authors, a visual map of heart sound energy was established. Both normal and abnormal heart sound recordings were studied. Representation is based on Wigner-Ville joint time-frequency transformations. The proposed methodology separates acoustic contributions of cardiac events simultaneously in pitch, time and energy. The resolution accuracy is superior to any other existing spectrogram method. The characteristic energy signature of the innocent heart murmur in a child with the S3 sound is presented. It allows clear detection of S1, S2 and S3 sounds, S2 split, systolic murmur, and intensity of these components. The original signal, heart sound power change with time, time-averaged frequency, energy density spectra and instantaneous variations of power and frequency/pitch with time, are presented. These data allow full quantitative characterization of heart sounds and murmurs. High accuracy in both time and pitch resolution is demonstrated. Resulting visual images have self-referencing quality, whereby individual features and their changes become immediately obvious.     

An algorithm for FHR estimation from foetal phonocardiographic signals

The long-term foetal surveillance is often to be recommended. Hence, the fully non-invasive acoustic recording, through maternal abdomen, represents a valuable alternative to the ultrasonic cardiotocography. Unfortunately, the recorded heart sound signal is heavily loaded by noise, thus the determination of the foetal heart rate raises serious signal processing issues. In this paper, we present a new algorithm for foetal heart rate estimation from foetal phonocardiographic recordings. A filtering is employed as a first step of the algorithm to reduce the background noise. A block for first heart sounds enhancing is then used to further reduce other components of foetal heart sound signals. A complex logic block, guided by a number of rules concerning foetal heart beat regularity, is proposed as a successive block, for the detection of most probable first heart sounds from several candidates. A final block is used for exact first heart sound timing and in turn foetal heart rate estimation. Filtering and enhancing blocks are actually implemented by means of different techniques, so that different processing paths are proposed. Furthermore, a reliability index is introduced to quantify the consistency of the estimated foetal heart rate and, based on statistic parameters; [,] a software quality index is designed to indicate the most reliable analysis procedure (that is, combining the best processing path and the most accurate time mark of the first heart sound, provides the lowest estimation errors). The algorithm performances have been tested on phonocardiographic signals recorded in a local gynaecology private practice from a sample group of about 50 pregnant women. Phonocardiographic signals have been recorded simultaneously to ultrasonic cardiotocographic signals in order to compare the two foetal heart rate series (the one estimated by our algorithm and the other provided by cardiotocographic device). Our results show that the proposed algorithm, in particular some analysis procedures, provides reliable foetal heart rate signals, very close to the reference cardiotocographic recordings.     

A method for computerized modification of certain natural animal sounds for communication study purposes

A new method for computerized modification of sound signals is presented. With digital signal processing in the time domain it is possible to alter the amplitude, the frequency and the time scale of natural sounds independently. The method can be applied to natural sounds with reasonably pure tonal quality.     

Comparison of first and second heart sounds after mechanical heart valve replacement

In this article, the spectral features of first heart sounds (S1) and second heart sounds (S2), which comprise the mechanical heart valve sounds obtained after aortic valve replacement (AVR) and mitral valve replacement (MVR), are compared to find out the effect of mechanical heart valve replacement and recording area on S1 and S2. For this aim, the Welch method and the autoregressive (AR) method are applied on the S1 and S2 taken from 66 recordings of 8 patients with AVR and 98 recordings from 11 patients with MVR, thereby yielding power spectrum of the heart sounds. Three features relating to frequency of heart sounds and three features relating to energy of heart sounds are obtained. Results show that in comparison to natural heart valves, mechanical heart valves contain higher frequency components and energy, and energy and frequency components do not show common behaviour for either AVR or MVR depending on the recording areas. Aside from the frequency content and energy of the sound generated by mechanical heart valves being affected by the structure of the lungs–thorax and the recording areas, the pressure across the valve incurred during AVR or MVR is a significant factor in determining the frequency and energy levels of the valve sound produced. Though studies on native heart sounds as a non-invasive diagnostic method has been done for many years, it is observed that studies on mechanical heart valves sounds are limited. The results of this paper will contribute to other studies on using a non-invasive method for assessing the mechanical heart valve sounds.     

Time-frequency analysis of phonocardiogram signals using wavelet transform: a comparative study

Analysis of phonocardiogram (PCG) signals provides a non-invasive means to determine the abnormalities caused by cardiovascular system pathology. In general, time-frequency representation (TFR) methods are used to study the PCG signal because it is one of the non-stationary bio-signals. The continuous wavelet transform (CWT) is especially suitable for the analysis of non-stationary signals and to obtain the TFR, due to its high resolution, both in time and in frequency and has recently become a favourite tool. It decomposes a signal in terms of elementary contributions called wavelets, which are shifted and dilated copies of a fixed mother wavelet function, and yields a joint TFR. Although the basic characteristics of the wavelets are similar, each type of the wavelets produces a different TFR. In this study, eight real types of the most known wavelets are examined on typical PCG signals indicating heart abnormalities in order to determine the best wavelet to obtain a reliable TFR. For this purpose, the wavelet energy and frequency spectrum estimations based on the CWT and the spectra of the chosen wavelets were compared with the energy distribution and the autoregressive frequency spectra in order to determine the most suitable wavelet. The results show that Morlet wavelet is the most reliable wavelet for the time-frequency analysis of PCG signals.     

Flowers respond to pollinator sound within minutes by increasing nectar sugar concentration

Can plants sense natural airborne sounds and respond to them rapidly? We show that Oenothera drummondii flowers, exposed to playback sound of a flying bee or to synthetic sound signals at similar frequencies, produce sweeter nectar within 3 min, potentially increasing the chances of cross pollination. We found that the flowers vibrated mechanically in response to these sounds, suggesting a plausible mechanism where the flower serves as an auditory sensory organ. Both the vibration and the nectar response were frequency‐specific: the flowers responded and vibrated to pollinator sounds, but not to higher frequency sound. Our results document for the first time that plants can rapidly respond to pollinator sounds in an ecologically relevant way. Potential implications include plant resource allocation, the evolution of flower shape and the evolution of pollinators sound. Finally, our results suggest that plants may be affected by other sounds as well, including anthropogenic ones.     

Novel PCG Analysis Method for Discriminating Between Abnormal and Normal Heart Sounds

A novel approach for separation among normal and heart murmurs sounds based on Phonocardiogram (PCG) analysis is introduced in this paper. The purpose of this work is to find the appropriate algorithm able to detect heart failures. Different features have been extracted from time and frequency domains. After the normalization step, the Principal Component Analysis algorithm is used for data reduction and compression. Support Vectors Machine (SVM), and k-Nearest Neighbors (kNN) classifiers were used with different kernels and number of neighbors in the classification step. Simulation results obtained from different databases are compared. The developed system gave good results when applied to different datasets. The accuracy of 96%, and 100% for the first, and the second dataset respectively were obtained. The algorithm shows its effectiveness in separation between normal and pathological cases.     

Probing real sensory worlds of receivers with unsupervised clustering

The task of an organism to extract information about the external environment from sensory signals is based entirely on the analysis of ongoing afferent spike activity provided by the sense organs. We investigate the processing of auditory stimuli by an acoustic interneuron of insects. In contrast to most previous work we do this by using stimuli and neurophysiological recordings directly in the nocturnal tropical rainforest, where the insect communicates. Different from typical recordings in sound proof laboratories, strong environmental noise from multiple sound sources interferes with the perception of acoustic signals in these realistic scenarios. We apply a recently developed unsupervised machine learning algorithm based on probabilistic inference to find frequently occurring firing patterns in the response of the acoustic interneuron. We can thus ask how much information the central nervous system of the receiver can extract from bursts without ever being told which type and which variants of bursts are characteristic for particular stimuli. Our results show that the reliability of burst coding in the time domain is so high that identical stimuli lead to extremely similar spike pattern responses, even for different preparations on different dates, and even if one of the preparations is recorded outdoors and the other one in the sound proof lab. Simultaneous recordings in two preparations exposed to the same acoustic environment reveal that characteristics of burst patterns are largely preserved among individuals of the same species. Our study shows that burst coding can provide a reliable mechanism for acoustic insects to classify and discriminate signals under very noisy real-world conditions. This gives new insights into the neural mechanisms potentially used by bushcrickets to discriminate conspecific songs from sounds of predators in similar carrier frequency bands.     

SOUND ANALYSIS AND SYNTHESIS USING AN AMIGA MICRO-COMPUTER

ABSTRACT

This note presents a program for analysis and synthesis of sounds designed to be used with an Amiga micro-computer. The user is in real-time control of the capture of the signal with reference to both time and amplitude parameters. Subsequently the chosen signal can be edited and analysed. Sonograms can be obtained at three resolutions, both in time and frequency. Measurements of time can be made on the signal, of frequency on the sonogram and of amplitude on sections. Three functions for the extraction of the fundamental are available. A synthetic sound may be obtained, either by modifying the amplitude and frequency modulation of an existing sound, or by creating a new sound by additive synthesis. All the commands can be selected by using a mouse, and the program is user-friendly.     

Segmentation of cardiac sound signals by removing murmurs using constrained tunable-Q wavelet transform

The automatic segmentation of cardiac sound signals into heart beat cycles is generally required for the diagnosis of heart valve disorders. In this paper, a new method for segmentation of the cardiac sound signals using tunable-Q wavelet transform (TQWT) has been presented. The murmurs from cardiac sound signals are removed by suitably constraining TQWT based decomposition and reconstruction. The Q-factor, redundancy parameter and number of stages of decomposition of the TQWT are adapted to the desired statistical properties of the murmur-free reconstructed cardiac sound signals. The envelope based on cardiac sound characteristic waveform (CSCW) is extracted after the removal of low energy components from the reconstructed cardiac sound signals. Then the heart beat cycles are derived from the original cardiac sound signals by mapping the required timing information of CSCW which is obtained using established methods. The experimental results are included in order to show the effectiveness of the proposed method for segmentation of cardiac sound signals in comparison with other existing methods for various clinical cases.     

基于MFCC和GMM的昆虫声音自动识别

昆虫的运动、 取食、 鸣叫都会发出声音, 这些声音存在种内相似性和种间差异性, 因此可用来识别昆虫的种类。基于昆虫声音的昆虫种类自动检测技术对协助农业和林业从业人员方便地识别昆虫种类非常有意义。本研究采用了语音识别领域里的声音参数化技术来实现昆虫的声音自动鉴别。声音样本经预处理后, 提取梅尔倒谱系数（Mel frequency cepstrum coefficient, MFCC）作为特征, 并用这些样本提取的MFCC特征集训练混合高斯模型（Gaussian mixture model, GMM）。最后用训练所得到的GMM对未知类别的昆虫声音样本进行分类。该方法在包含58种昆虫声音的样本库中进行了评估, 取得了较高的识别正确率（平均精度为98.95%）和较理想的时间性能。该测试结果证明了基于MFCC和GMM的语音参数化技术可以用来有效地识别昆虫种类。     

Non-invasive diagnosis methods of coronary disease based on wavelet denoising and sound analyzing

The heart sound is the characteristic signal of cardiovascular health status. The objective of this project is to explore the correlation between Wavelet Transform and noise performance of heart sound and the adaptability of classifying heart sound using bispectrum estimation. Since the wavelet has multi-scale and multi-resolution characteristics, in this paper, the heart sound signal with different frequency ranges is decomposed through wavelet and displayed on different scales of the resolving wavelet result. According to distribution features of frequency of heart sound signals, the interference components in heart sound signal can be eliminated by selecting reconstruction coefficients. Comparing de-noising effects of four wavelets which are haar, db6, sym8 and coif6, the db6 wavelet has achieved an optimal denoising effect to heart sound signals. The de-noising result of contrasting different layers in the db6 wavelet shows that decomposing with five layers in db6 provide the optimal performance. In practice, the db6 wavelet also shows commendable denoising effects when applying to 51 clinical heart signals. Furthermore, through the clinic analyses of 29 normal signals from healthy people and 22 abnormal heart signals from coronary heart disease patients, this method can fairly distinguish abnormal signals from normal signals by applying bispectrum estimation to denoised signals via ARMA coefficients model.     

Estimation of Respiratory Nasal Pressure and Flow Rate Signals Using Different Respiratory Sound Features

BackgroundRespiratory sounds are associated with the flow rate, nasal flow pressure, and physical characteristics of airways. In this study, we aimed to develop the flow rate and nasal flow pressure estimation models for the clinical application, and find out the optimal feature set for estimation to achieve the optimal model performance.MethodsRespiratory sounds and flow rate were acquired from nine healthy volunteers. Respiratory sounds and nasal flow pressure were acquired from twenty-three healthy volunteers. Four types of respiratory sound features were extracted for flow rate and nasal flow pressure estimation using different estimation models. Effects of estimations using these features were evaluated using Bland-Altman analysis, estimation error, and respiratory sound feature calculation time. Besides, expiratory and inspiratory phases divided estimation errors were compared with united estimation errors.ResultsThe personalized logarithm model was selected as the optimal flow rate estimation model. Respiratory nasal flow pressure estimation based on this model was also performed. For the four different respiratory sound features, there is no statistically significant difference in flow rate and pressure estimation errors. LogEnvelope was, therefore, chosen as the optimal feature because of the lowest computational cost. In addition, for any type of respiratory sound feature, no statistically significant difference was observed between divided and united estimation errors (flow rate and pressure).ConclusionRespiratory flow rate and nasal flow pressure can be estimated accurately using respiratory sound features. Expiratory and inspiratory phases united estimation using respiratory sounds is a more reasonable estimation method than divided estimation. LogEnvelope can be used for this united respiratory flow rate and nasal flow pressure estimation with minimum computational cost and acceptable estimation error.     

长江江豚声信号及其声行为的初步研究

长江江豚的声信号可分为两大类,即高频脉冲信号和低频连续信号。高频脉冲信号可能与回声定位有关,而低频连续信号可能与通信和情感表达有关。不管是高频脉冲信号还是低频连续信号,在豚处于自由状态时,夜间的发声次数要多于白天。       

A Humanoid Method for Extracting Abnormal Engine Sounds from Engine Acoustics Based on Adaptive Voiterra Filter

The improvement of SNR (Signal-to-Noise Ratio) of abnormal engine sounds is of great help in improving the accuracy of engine fault diagnosis.By imitating the way that human technicians use to distinguish abnormal engine sounds from engine acoustics,a humanoid abnormal sound extracting method is proposed.By implementing adaptive Volterra filter in the canonical Adaptive Noise Cancellation (ANC) system,the proposed method is capable of tracing the engine baseline sound which exhibits an intrinsic nonlinear dynamics.Besides,by introducing a template noise tailored from the records of engine baseline sound and taking it as virtual input of the adaptive Volterra filter,the priori knowledge of engine baseline sound,such as inherent correlation,periodicity or phase information,and stochastic factors,is taken into consideration.The hybrid simulations prove that the proposed method is functional.Since the method proposed is essentially a single-sensor based ANC,hopefully,it may become an effective way to extricate the dilemma that canonical dual-sensor based ANC encounters when it is used in extracting fault-featured signals from observed signals.