Design methodology of a new wavelet basis function for fetal phonocardiographic signals.
Continuous and long-term fetal monitoring has become an essential approach for better accuracy in diagnosis . Variations in fetal heart rate (FHR) provide up-to-date information about the fetal health status . Doppler ultrasound based fetal cardiotocography (fCTG) is currently being used as a technique for recording and analysis of the FHR. In this technique, a graph of fetal heart rate (cardio-) and uterine con tractions (-toco-) are recorded during pregnancy . This technique cannot be used for long-term monitoring of the fetus due to several reasons  such as the cost and complexity of the monitoring instrument, ultrasound radiation exposure, and gel application.
Fetal phonocardiography, recording of vibroacoustic (fPCG) signals from the maternal abdominal surface, may become an important alternative to fCTG . However, these fPCG signals are heavily contaminated by noise from various sources [6-11]. Hence to extract important diagnostic information such as FHR, frequency/pitch, intensity, timing, and energy of fetal heart sound, this technique requires matured signal processing strategies.
Researchers from biomedical signal processing community have applied several techniques for denoising of fPCG signals. Varady (2001) presented a wavelet-based denoising method for phonocardiographic signals using two-channel signal recording and an adaptive cross-channel coefficient thresholding technique . Messer et al. (2001) attempted to answer about which wavelet families, levels of decomposition, and thresholding techniques best removes the noise in a PCG . Jianfeng et al.  used a modified spectral subtraction algorithm to remove the unwanted stationary background noise from the noisy fetal heart sound. The AM/FM modulation technique is employed to make the fetal heart sound more audible so that both pregnant women and gynecologists can identify the rhythmic fetal heart beat sounds easily . Mittra et al.  compared and analyzed the performance of various digital noise cancellation techniques for fetal heart sound. Their experimental results show that adaptive filtering using Recursive Least Squares (RLS) algorithm is an appropriate methodology for the fetal phonocardiographic signals denoising implementation . Tosanguan et al. (2008) proposed a 2-channel signal processing technique, termed interference suppression via spectral comparison (ISSC) and aims at improving the quality of the recorded heart sound or the PCG data . Zhao et al. (2009) developed a denoising method for heart sound signal using improved thresholding function in wavelet domain . Chourasia et al. (2010) presented a technique for denoising of fPCG signals using wavelet transform . In another work, Chourasia et al. (2012) also developed a methodology for removal of unwanted noise from the fPCG signal using non negative matrix factorization .
The above-mentioned researchers have used conventional filtering techniques and wavelet transform for denoising the fPCG signals. The conventional techniques also impact on the useful signals and hence may result in loss of information of diagnostic importance. Additionally, the passed band may still contain noise. In view of these limitations of conventional methods of denoising, wavelet thresh old denoising has been applied to denoise heart sound recordings. The wavelet based denoising method try to pre serve the signal by operating only on those selected regions of the bandwidth that need filtering. This technique requires appropriate selection of wavelet family, level of decomposition, and the method to be used for calculation of threshold. The existing works are limited only to employ the existing wavelet family from the bank of wavelet transform.
This present work proposes a novel algorithm for denoising of fPCG signals using wavelet transform (WT). In this approach, a new wavelet family and its mother wavelet are developed. For this purpose, a quadratic mirror filter (QMF) bank  is designed based on the characteristics of the fPCG signals. This filter bank requires a low-pass and a high pass filter in decomposition (analysis) phase and reconstruction (synthesis) phase. Appropriate denoising algorithm and thresholding rule have been selected and used with the developed mother wavelet. The developed wavelet family provides better results when compared with the existing wave lets. The obtained denoised fPCG signals retain the vital diagnostics information contained in the original signals.
The rest of the paper is organized as follows. A brief introduction about fetal phonocardiographic signals and wavelet theory has been provided in Sections 2 and 3, respectively. Section 4 discusses the selection of appropriate denoising algorithm and thresholding rule for fPCG signals. Design methodology of filter bank for new wavelet has been explained in Section 5. Finally, the experimental results are described in Section 6, followed by conclusion and discussion in Section 7.
2. Fetal Phonocardiographic Signals
Fetal phonocardiography is the recording of natural vibroacoustic signals from the maternal abdominal surface. The fPCG signal carries valuable information concerning the physiological state of unborn . It is also capable of recognizing additional potential dysfunctional signals of the fetal heart such as those related to cardiac murmurs, split effect, and breathing movements, which can be detected from the analysis of abdominal sound (fPCG) signals. However these cardiac anomalies are impossible to detect with the widely accepted fCTG technique due to its principle of operation . Additionally, phonocardiography provides a permanent and lasting record of fetal heart sounds which by compari son at a later stage may prove to be of great prognostic importance. The fPCG technique has been introduced earlier but overlooked by biomedical scientists and medical experts mainly due to its low signal-to-noise ratio (SNR) . The fPCG signals are a linear summation of 
(i) fetal heart sound,
(ii) internal noise,
(iii) external noise.
The fetal heart sound is a signal produced by mechanical activity of the fetal heart. The fetal heart is basically divided into two pairs of chambers and has four valves: the mitral and tricuspid valves. In the fetal cardiac cycle, when the ven tricles begin to contract, the blood attempts to flow back into the lower pressure atrial chambers. This reverse flow of blood is arrested by the shutting of the mitral and tricuspid valves, which produces the first heart sound (S1). Whenever the pressure in the ventricular chambers becomes too high for the pulmonary valves to withstand, they open, and the pre ssurized blood is rapidly ejected into the arteries. While the ventricles are being evacuated, the pressure of the remaining blood decreases with respect to that in the arteries. This pressure gradient causes the arterial blood to flow back into the ventricles. The pulmonary valves arrest this reverse flow by shutting, which gives rise to the second heart sound (S2) . The frequency spectrum of fetal heart sound lies below 200 Hz. Figures 1(a) and 1(b) show a typical fetal phono-cardiographic (fPCG) signal and its frequency spectrum, respectively.
The internal noise is a random signal caused due to mat ernal respiratory sounds, acoustic noise produced by the fetal movement, maternal digestive sound, maternal respiratory sound, maternal heart sound, and placental blood turbulence. These noises are of low amplitude with main frequency com ponents from 0 to 25 Hz .
Similarly, the external noise is a combination of shear noise from movement of the sensor during recording and ambient noise originating from the environment such as sound produced by fan, air conditioner, and hue and cry of the nearby people. It is comparatively of high amplitude and frequency (100-20000 Hz) .
3. Theoretical Background
3.1. Wavelet Transform. The WT is a two-dimensional time scale processing method for nonstationary signals with adequate scale values and shifting in time. The major advantage of the WT is that it has a varying window size, being broad at low frequencies and narrow at high frequencies, thus leading to best possible time-frequency resolution in all frequency ranges [25,26]. The analysis of nonstationary signals requires proper location of transitions or discontinuities and identification of their long-term behavior. The WT represents a time function in terms of simple and fixed building blocks derived from mother wavelet by translation and dilation operations. The translation is the shifting of mother wavelet along the time axis while dilation or scaling starches or compresses it.
The WT can be categorized as continuous wavelet trans form (CWT) or discrete wavelet transform (DWT). The CWT is defined as the convolution between the original signal s(t) and a wavelet y(t) which can be calculated by
[CWT.sub.[psi]] (a,b) = [[integral].sup.+[infinity].sub.-[infinity]] s (t) [[psi].sup.*.sub.a,b] (t) dt = 1/[square root of (a)] [[integral].sup.+[infinity].sub.-[infinity]] s (t) [[psi].sup.*] (t - b/a)dt, (1)
where [CWT.sub.[psi]] (a, b) is a continuous wavelet transform, "s(t)" is a signal under study, "a" is a scale coefficient connected with stretching or compression of signal in time, "b" is a shift connected with time location, "[[psi].sup.*.sub.a,b](t)" is a wavelet function or mother wavelet representing a wavelet family, "*" denotes the complex conjugation, and the factor 1 [square root of a] is used for energy normalization purposes so that the transformed signal will have the same energy at every scale. In CWT, the scaling parameter a and translation parameter b change with time continuously. Hence wavelet coefficients are calculated for every possible scale, which requires huge processing power and results in large amount of data.
The Discrete Wavelet Transform (DWT) coefficients are usually sampled from the CWT on a dyadic grid, choosing parameters of translation b = [k.sup.*] [2.sup.-j] and scale a = [2.sup.-j]. Where j, k [member of] Z a set of positive integers and k = 0, 1, ..., n - 1, n represents the number of samples. These dilation and translation parameters are discretized leading to the DWT. After discretization, the wavelet function is defined as
[DWT.sup.[psi]] (j, k) = [[integral].sup.+[infinity].sub.-[infinity]] s(t) [[psi].sup.*.sub.j,k](t)dt. (2)
Here [[psi].sup.*.sub.j,k](t) is the dilated and translated version of the wavelet function and given as
[[psi].sub.j,k] (*) = [2.sup.j/2] [psi]([2.sup.j]t - k) , (3)
where [psi] is called as mother wavelet and [[psi].sub.j,k] is called as daughter wavelet. The level j determines how many wavelets are needed to cover the mother wavelet, and the number k determines the position of the wavelet and gives the indication of time. DWT analyzes the signal by decomposing it into its coarse and detail information, which is accomplished by using successive high-pass and low-pass filtering operations, on the basis of the following equations:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)
where [y.sub.high](k) and [y.sub.low](k) are the outputs of the high-pass and low-pass filters with impulse response h and g, respectively, after downsampling by 2 [27, 28]. The coefficients of low-pass filter are called "approximation" ([c.sub.j,k]) and coefficients from high-pass filter are called "detail" ([b.sub.j,k]) wavelet coefficients. The detail coefficients are defined by the following equation:
[b.sub.j, k] = [integral] s(t) [[psi].sup.*.sub.j, k](*)dt, (5)
where [[psi].sub.j,k] are wavelet functions given by
[[psi].sub.j,k](t) = [1/[square root ([2.sup.j])] [psi] (t - k[2.sup.j]/[2.sup.j]). (6)
Similarly, the approximate coefficients are
[c.sub.j,k] = [integral] S(t) [[PHI].sup.*.sub.j,k] (t)dt, (7)
where [[PHI].sub.j,k] are called scaling functions as follows:
[[PHI].sub.j,k](t) = 1/[square root of ([2.sup.j])] [PHI] (t - k[2.sup.j]/[2.sup.j])). (8)
The discrete inverse transform is found by adding the translated, dilated wavelets, weighted by the coefficients:
f(t) = [summation over (j,k)] [b.sub.j,k][[psi].sub.j,k](t). (9)
The DWT gives a multiresolution description of a signal which is very useful in analyzing real-time signals .
The general wavelet denoising procedure is described in the following steps:
(i) decomposition of the fPCG signal using DWT to obtain the approximation and detail coefficients,
(ii) thresholding of these decomposed coefficients using an appropriate denoising algorithm,
(iii) reconstruction of the fPCG signal from these thresholded coefficients using the inverse transform (IDWT) .
3.2. Denoising Algorithm. The denoising algorithm uses statistical regression of noisy coefficients over time to obtain a nonparametric estimation of the reconstructed signal with out noise. The thresholding algorithms commonly employed for denoising of the nonstationary signals are 
(i) universal threshold,
(ii) minimax threshold,
(iii) rigorous Stein's Unbiased Risk Estimate (SURE) threshold.
(i) Universal Threshold (Sqtwolog). The universal threshold denoising algorithm is a fixed threshold method which can be calculated by
[lambda] = [sigma] [square root of 2log(n)], (10)
where n denotes the length of the signal and [sigma] is the standard deviation.
(ii) Minimax Threshold (Minimaxi). Minimax threshold also uses fixed threshold and it yields minimax performance for Mean Square Error (MSE) against an ideal procedures. This threshold level depends on the noise and signal relationships in the input data and it is given by [lambda] = [sigma][[lambda].sub.n], where [[lambda].sub.n] is determined by a minimax rule such that the maximum risk of estimation error across all locations of the data is minimized.
(iii) Rigorous SURE Threshold (Rigrsure). The denoising algorithms described previously use global thresholds. That is, the computed threshold is applied to all wavelet coefficients. The rigorous SURE threshold algorithm describes a scheme that uses a threshold value [[lambda].sub.j] at each resolution level j of the wavelet coefficients. This algorithm is also known as SureShrink and uses the Stein's Unbiased Risk Estimate (SURE) criterion to get an unbiased estimate.
These denoising algorithms can be divided into linear and nonlinear methods. The linear method is independent of the size of empirical wavelet coefficients, and therefore the size of the coefficient by itself is not taken into account. It assumes that signal noise can be found mainly in fine scale coefficients and not in coarse scales. The nonlinear method is based on the idea that the signal noise can be found in every coefficient and is distributed over all scales.
3.3. Wavelet Thresholding Rules. Wavelet thresholding is a signal estimation technique that exploits the capabilities of wavelet transform for signal denoising . This method has been researched extensively due to its effectiveness and simplicity. Any signal s(t) can be represented by the summation of the original signal x(t) and the noise n(t) as follows:
s(t) = x(t) +n(t). (11)
After performing the wavelet transform
[S.sub.j,k] = [X.sub.j,k] + [N.sub.j,k], (12)
where [S.sub.jk] is the kth wavelet coefficient in the scale j. There are two ways of thresholding with threshold [lambda]; the shapes of these thresholding operators are illustrated in Figure 2.
Hard Thresholding (h). In hard thresholding, those wavelet coefficients with absolute values below or at the threshold level ([lambda]) are affected only and they are replaced by zero value whereas others are kept unchanged
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)
Soft Thresholding (s). In soft thresholding, coefficients above threshold level ([lambda]) are also modified; they are reduced by particular value of the threshold
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)
Hard thresholding maintains the scale of the signal but introduces ringing and artifacts after reconstruction due to a discontinuity in the wavelet coefficients. Soft thresholding eliminates this discontinuity resulting in smoother signals but slightly decreases the magnitude of the reconstructed signal.
4. Selection of Denoising Algorithm and Thresholding Rule for fPCG Signals
The presented work also contributes in the selection of suit able algorithm and thresholding rule for denoising of fPCG signals. Based on the discussion in Sections 3.2 and 3.3, there are three denoising algorithms, namely, universal threshold, minimax threshold, and rigorous SURE threshold and two thresholding rules which can be employed in denoising of fPCG signals. These denoising algorithms and thresholding rules are practically implemented and their results are com pared. The mean squared error (MSE) is used to evaluate the performance of all denoising algorithms and thresholding rules in denoising the fPCG signals. It can be obtained using following expression:
MSE = [[summation].sup.n.sub.i=1][(s - [s.sub.e]).sup.2.sub.i]/n, (15)
where n denotes the length of the signal, s represents the original signal and [s.sub.e] is the estimated signal obtained from the denoised wavelet coefficients.
Figure 3(a) shows the waveform of one simulated fPCG signal as an example. This signal is used as a reference signal in the selection process of appropriate denoising algorithms and wavelet threshold. Figure 3(b) is a test signal generated by adding simulated stationary random noise in the original fPCG signal which is shown in Figure 3(a). The noise is generated by recording it, first, in actual conditions. After recording, the noise signal is simulated with similar average characteristics. The simulated noise is then mixed with the fPCG signal shown in Figure 3(a) and used as an input to evaluate the performance of the selection of appropriate denoising algorithm and thresholding rule.
The test signal so obtained is analyzed with DWT based multiresolution analysis. The signal is decomposed to five levels using fourth order Coiflets wavelet. This selection of mother wavelet is based on the fact that it possesses all the properties needed for analysis of the fPCG signals .
All the three algorithms with soft or hard thresholding rule are applied for denoising of the fPCG signal. The resultant waveforms from this implementation are shown in Figure 4, and the best estimations obtained are depicted in Table 1. In Figures 3 and 4, the x-axis represents number of samples and the y-axis represents relative amplitude.
Table 1 shows a comparison of three denoising algorithms with soft or hard thresholding rule. The rigorous SURE threshold algorithm with soft thresholding rule yields the best estimation with considerably smaller MSE as compared to the other algorithms for denoising of fPCG signals.
5. Design of Filter Bank for New Wavelet
Since it has been discussed earlier in Section 1 that the effective implementations of wavelet transform in denoising the fPCG signals requires appropriate wavelet basis function, we propose to design a new wavelet function which will be based on the characteristics of the fPCG signals and hence will improve the performance of its denoising. The main properties to be verified in designing a family of wavelet are existence of scaling function, symmetricity and compactly support of wavelet and scaling functions, availability of filter bank, orthogonality of filter bank, and smoothness of wavelet function. With these properties, the wavelet family may be called as suitable wavelet for the analysis of intended signal. This wavelet will lead to maximization of wavelet coefficient values to produce the highest local maxima of the signal in wavelet domain. It also produces the possibility of best characterization of frequency content of that signal .
In view of these considerations, we developed a new wavelet function which is orthogonal and named as "fetal." In orthogonal wavelet analysis, the number of convolutions at each scale is proportional to the width of the wavelet basis at that scale. This produces a wavelet spectrum that contains discrete "blocks" of wavelet power and is useful for signal processing as it gives the most compact representation of the signal . Conversely, a nonorthogonal analysis is highly redundant at large scales, where the wavelet spectrum at adjacent times is highly correlated. The new wavelet also has a small number of coefficients in high-pass subbands and allows the signal singularities, transitions, and edges intact in the low-pass subband.
To synthesize new wavelet it requires a two-channel filter bank which has a low-pass and a high-pass filter in decomposition (analysis) phase and reconstruction (synthesis) phase. This two-channel multirate filter bank consists of filters which process the input signal at half of its original rate. A block diagram of two-channel filter bank is shown in Figure 5. The signal s(n) gets filtered by Go and Ho filters and down-sampled. These filters are called analysis or decomposition filters. The output of these filters contains the signal at halfrate. These are called subbands of the signal. Each subband can be further divided into smaller subbands using the same filter bank. After being processed, the signal is upsampled and filtered using [G.sub.1] and [H.sub.1] filters, which are called synthesis or reconstruction filters. For perfect reconstruction, in a two-channel filter bank, the downsampled signal that contains only the even samples is given to the first channel and the down-sampled signal that contains only the odd samples is given to the second channel. This separation of signal into even and odd components is called polyphase representation of the signal.
Four important characteristics are taken into account while designing the filter bank: (i) perfect reconstruction, (ii) orthogonality of the filter bank and the underlying wavelet based multiresolution structure, (iii) flatness of the filters and vanishing moments in the wavelets, and (iv) smoothness of the wavelets. The output of a perfect reconstruction two channel filter bank is 
Y(z) = [1/2] [[G.sub.o](z)[G.sub.1](z) + [H.sub.o](z)[H.sub.1](z)] X(z) + [1/2] [[G.sub.o](-z) [G.sub.1](z) + [H.sub.o](-z) [H.sub.1](z)] X (-z), (16)
where X(z), G(z), and H(z) are the z-transform of the input signal, analysis filters and synthesis filters respectively. For perfect reconstruction: firstly, the alias term X(-z) must be zero, hence
[G.sub.o](-z)[G.sub.1](z) + [H.sub.o](-z)[H.sub.1](z) = 0. (17)
This can be achieved by letting
[G.sub.1](z) = [H.sub.o](-z), [H.sub.1](z) = -[G.sub.o](-z). (18)
Second, the distortion term must be constant or a pure delay time; that is,
[G.sub.o](z)[G.sub.1](z) + [H.sub.o](z)[H.sub.1](z) = 2[z.sup.-1], (19)
where l denotes a time delay.
Equation (18) can also be written as
[H.sub.o](z) = [G.sub.1](-z), [H.sub.1](z) = -[G.sub.o](-z). (20)
Substituting this in (19)
[G.sub.o](z)[G.sub.1](z) - [G.sub.1](-z) [G.sub.o](-z) = 2[z.sup.-1], (21)
[P.sub.o](z) - [P.sub.o](-z) = 2[z.sup.-1], (22)
where [P.sub.o](z) denotes the product of two low-pass filters [G.sub.o](z) and [G.sub.1](z). Equation (21) indicates that all odd terms of product of two low-pass filters must be zero for order l, where I must be odd and even order terms are arbitrary . Hence it can be written as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)
Hence it can be concluded that the design process of two channel filter bank for a new wavelet is reduced into two steps as follows:
(1) design of filter [P.sub.o](z) which satisfies (23),
(2) factorize [P.sub.o](z) into [G.sub.o](z) and [G.sub.1](z).
In this work, the filter [P.sub.o](z) is designed based on the characteristics of fPCG signals. The common requirements of this design are linear phase, minimum phase, and orthogonality of the filter. A fourth-order low-pass Butterworth filter is chosen because the transition width requirement is not stringent for the given cut-off frequency. This will also help in reducing the computational complexity. The Butterworth filter satisfies the conditions for perfect reconstruction. It has linear frequency response in the pass band as compared to Chebyshev Type I/Type II and elliptic filters. The filter [P.sub.o](z) is factorized into Go (z) and G1(z), and then the coefficients of Ho (z) and H1(z) are derived using (20) [32, 36]. Figure 6 shows the impulse response of the four filters computed for construction of filter bank of the new wavelet "fetal".
The wavelet and scaling functions are then derived from the coefficients of these filters using (6) and 8), respectively. Figure 7 shows the wavelet and scaling functions of the "fetal" wavelet.
With these wavelet and scaling functions, the wavelet and scaling coefficients for multiresolution analysis are obtained. The developed wavelet "fetal" is now ready to use. All the discrete analysis functions, including dwt, idwt, and wave-dec. can operate on the new wavelet. Similarly, all the continuous analysis functions, including cwt, wscalogram, and the corresponding GUI tools, can also operate on the new wavelet. In this work, this new wavelet "fetal" is used for intended denoising of the fPCG signals.
6. Experimental Results
The newly developed wavelet "fetal" is used for denoising of real-time fetal heart sound signals which are of widely distinct nature and quality. A Simulink model is developed for computer based denoising of these signals . Figure 8 shows the model for wavelet denoising of the fPCG signal through which the filters of developed wavelet family and selected threshold criterion are implemented.
In this model, the fPCG signals are fetched from the work space and applied to the Dyadic analysis filter bank. These signals carry fetal heart sound and a damped version of simulated maternal organs' sounds along with the external noise. The analysis filter bank decomposes the fPCG signals into a collection of subbands with smaller bandwidths and slower sample rates. This bank uses a series of high-pass and low-pass FIR filters to repeatedly divide the input frequency range. The fPCG signals are decomposed to the 3 levels by newly developed "fetal" wavelet with the analysis filter bank. These decomposed wavelet coefficients consist of details in the input fPCG signals. The denoising of fPCG signals is carried out by using selected algorithm (Rigorous SURE) and thresholding rule (Soft threshold). The synthesis filter bank reconstructs the signal decomposed by the analysis filter bank block. This bank takes in subbands of this signal and uses them to reconstruct the signal by using a series of high-pass and low-pass FIR filters. The reconstructed signals have a wider bandwidth and faster sample rate than the input subbands. The waveforms of input signal and denoised outputs are displayed through "Time Scope" blocks.
For performance evaluation of the developed wavelet, five test fPCG signals (S1-S5) are generated using a signal simulation module (SSM) . These test signals are denoised using some existing competitive wavelets families and the developed wavelet. The competitive wavelets are the one that possess required properties of orthogonality and a scaling function [PHI] such as Daubechies, Symlets, and Coiflets. The results of this comparison are computed using MSE equation (15) and presented in Table 2. As explained earlier, all the three denoising algorithms with soft or hard threshold are used for denoising of the fPCG signals.
The results in Table 2 show that the developed wavelet "fetal" along with rigorous SURE algorithm and soft thresh old provides the best performance for denoising the fPCG signals.
The developed wavelet family was also implemented for denoising of a real fPCG signal. This signal was recorded using a specially developed wireless data acquisition system . The recording of fPCG signal was obtained in a quiet room with the help of a medical expert and a trained nurse. A pregnant woman visiting for antenatal care was requested to contribute in the real-time testing of the developed system. The subject was having 34 weeks of gestation and singleton pregnancy. She was asked to lie in the supine position, with the head resting on a pillow. The phonocardiographic sensor was positioned on the abdominal surface of the mother and adjusted to acquire the maximum intensity of the signal. An abdominal belt was used to fasten the sensor. The signal was sampled with sampling frequency of 8000 Hz, 16 bit resolution and saved for further processing.
The recorded fPCG signal is now fetched from the memory through the model, as shown in Figure 8, for denoising. Figure 9 shows the waveforms of original and denoised version of this signal. In these waveforms, x-axis represents the time in seconds, whereas y-axis represents amplitude of signal in volts.
7. Discussion and Conclusion
The fPCG signals are of very low amplitude and contain poor signal-to-noise ratio. The main sources of noise are maternal biological activities, external noises such as sound produced by electrical appliances, movement of transducer, and so forth. These noises show overlapping spectra with the actual fPCG signals. Hence conventional noise removal techniques are not suitable for denoising of these signals. During denoising, care has to be practiced to preserve the features contained in the original signal. These preserved features are relevant and necessary for an appropriate diagnosis about fetal health.
In this paper design of a new wavelet basis function for denoising of fPCG signals has been carried out. The key features of newly designed family are its speed of convergence at infinity to 0, regularity, and orthogonality. It also has a small number of coefficients in high-pass subbands and allows the signal singularities, transitions, and edges intact in the low pass subband. It has been found that the combination of optimal perfect reconstruction filter bank and appropriate denoising algorithm can improve the performance of denoising. The experimental results revealed suitability of the newly developed wavelet to be the most appropriate wavelet basis function for denoising of fPCG signals in terms of MSE. The resultant denoised fPCG signals will preserve physiological information of diagnostic importance regarding health status of the unborn.
The authors of this paper would like to thank the experts and technical staff of Government Women Hospital, Gondia (M.S.), India, for their kind support. They also thank the pregnant women who contributed in the clinical experimen tation.
 V. S. Chourasia and A. K. Tiwari, "A review and comparative analysis of recent advancements in fetal monitoring tech niques," Critical Reviews in Biomedical Engineering, vol. 36, no. 5-6, pp. 335-373, 2008.
 M. Murray, Antepartal and Intrapartal Fetal Monitoring, Spring er, 3rd edition, 2007.
 J. Karin, M. Hirsch, O. Segal, and S. Akselrod, "Non invasive fetal ECG monitoring," in Computers in Cardiology, pp. 365 368, IEEE Computer Society Press, Bethesda, Md, USA, 1994.
 F. Kovacs, M. Torok, and I. Habermajer, "A rule-based phonocardiographic method for long-term fetal heart rate monitoring, " IEEE Transactions on Biomedical Engineering, vol. 47, no. 1, pp. 124-130, 2000.
 C. Jianfeng, K. Phua, S. Ying, and L. Shue, "A portable phonocardiographic fetal heart rate monitor," in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 2141-2144, Island of Kos, Greece, May 2006.
 H. G. Goovaerts, O. Rompelman, and H. P. Van Geijn, "A transducer for detection of fetal breathing movements," IEEE Trans actions on Biomedical Engineering, vol. 36, no. 4, pp. 471-478, 1989.
 J. T. E. McDonnell, "Knowledge-based interpretation of foetal phonocardiographic signals," IEE Proceedings F, vol. 137, no. 5, pp. 311-318, 1990.
 M. Godinez, A. Jimenez, R. Ortiz, and M. Pena, "On-line fetal heart rate monitor by Phonocardiography," in Proceedings of the 25th Annual International Conference of the IEEE EMBS, pp. 3141-3144, Cancun, Mexico, September 2003.
 R. Acharyya, N. L. Scott, P. Teal, E. Deuss, and J. Flierl, "Signal separation for non-invasive monitoring of foetal heartbeat," in Proceedings of the 1st International Conference on BioMedical Engineering and Informatics, pp. 497-501, May 2008.
 A. Jimenez-Gonzelez and C. J. James, "Extracting sources from noisy abdominal phonograms: a single-channel blind source separation method," Medical & Biological Engineering& Computing, vol. 47, no. 6, pp. 655-664, 2009.
 F. Kovacs, C. Horveth, A. T. Balogh, and G. Hosszu, "Extended non-invasive fetal monitoring by detailed analysis of data measured with phonocardiography," IEEE Transactions on Biomedical Engineering, vol. 58, no. 1, pp. 64-70, 2011.
 P. Varady, "Wavelet-based adaptive denoising of phonocardiographic records, "in Proceedings of the 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society, vol. 2, pp. 1846-1849, 2001.
 S. R. Messer, J. Agzarian, and D. Abbott, "Optimal wavelet denoising for phonocardiograms," Microelectronics Journal, vol. 32, no. 12, pp. 931-941, 2001.
 A. K. Mittra, A. Shukla, and A. S. Zadgaonkar, "System simulation and comparative analysis of foetal heart sound de-noising techniques for advanced phonocardiography," International Journal of Biomedical Engineering and Technology, vol. 1, no. 1, pp. 73-85, 2007
 T. Tosanguan, R. J. Dickinson, and E. M. Drakakis, "Modified spectral subtraction for de-noising heart sounds: interference suppression via spectral comparison," in Proceedings of the IEEE-BIOCAS Biomedical Circuits and Systems Conference, pp. 29-32, Baltimore, Md, USA, November 2008.
 X. M. Zhao and G. T. Cao, "A novel de-noising method for heart sound signal using improved thresholding function in wavelet domain," in Proceedings of the International Conference on Future BioMedical Information Engineering, pp. 65-68, Sanya, China, December 2009.
 V S. Chourasia and A. K. Mittra, "Wavelet-based denoising of fetal phonocardiographic signals," International Journal of Medical Engineering and Informatics, vol. 2, no. 2, pp. 139-150, 2010.
 V S. Chourasia, A. K. Tiwari, R. Gangopadhyay, and K. A. Akant, "Foetal phonocardiographic signal denoising based on non-negative matrix factorization," Journal of Medical Engineering & Technology, vol. 36, pp. 57-66, 2012.
 D. Esteban and C. Galand, "Application of quadrature mirror filters to split band voice coding schemes," in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '77), pp. 191-195, May 1977.
 N. Colley, D. G. Talberg, and D. P. Southall, "Biophysical profile in the fetus from a phonocardiographic sensor," European Journal of Obstetrics & Gynecology and Reproductive Biology, vol. 23, no. 5-6, pp. 261-266, 1986.
 A. Di Lieto, M. De Falco, M. Campanile et al., "Regional and international prenatal telemedicine network for computerized antepartum cardiotocography," Telemedicine and e-Health, vol. 14, no. 1, pp. 49-54, 2008.
 B. H. Tan and M. Moghavvemi, "Real time analysis of fetal phonocardiography," in Proceedings of the TENCON Proceedings, vol. 2, pp. 135-140, Kuala Lumpur, Malaysia, September 2000.
 N. Colley, N. G. Abraham, P. Fayers, D. G. Talbert, W. L. Daviesa, and D. P. Southall, "The fetal phonogram: a measure of fetal activity," Lancet, vol. 327, no. 8487, pp. 931-935,1986.
 D. G. Talbert, W. L. Davies, and F. Johnson, "Wide bandwidth fetal phonography using a sensor matched to the compliance of the mother's abdominal wall," IEEE Transactions on Biomedical Engineering, vol. 33, no. 2, pp. 175-181,1986.
 X. C. Yin, P. Han, J. Zhang, F. Q. Zhang, and N. L. Wang, "Application of wavelet transform in signal denoising," in Proceedings of the International Conference on Machine Learning and Cybernetics, vol. 1, pp. 436-441, November 2003.
 I. Daubechies, Ten Lectures on Wavelets, Society for Industrial and Applied Mathematics (SIAM), New York, NY, USA, 1992.
 M. R. Raghuveer and A. S. Bopardikar, Wavelet Transforms: Introduction to Theory and Applications, Pearson Education, Indian Branch, Singapore, 3rd edition, 2002.
 A. Daamouche, L. Hamami, N. Alajlan, and F. Melgani, "A wave let optimization approach for ECG signal classification," Bio medical Signal Processing and Control, vol. 7, no. 4, pp. 342-349, 2012.
 Md. A. Kabir and C. Shahnaz, "Denoising of ECG signals based on noise reduction algorithms in EMD and wavelet domains," Biomedical Signal Processing and Control, vol. 7, no. 5, pp. 481-489, 2012.
 V. S. Chourasia and A. K. Mittra, "Selection of mother wavelet and denoising algorithm for analysis of foetal phonocardiographic signals," Journal of Medical Engineering and Technology, vol. 33, no. 6, pp. 442-448, 2009.
 D. L. Donoho and I. M. Johnstone, "Adapting to unknown smoothness via wavelet shrinkage," Journal of the American Statistical Association, vol. 90, no. 432, pp. 1200-1224,1995.
 G. Strang and T. Nguyen, Wavelets and Filter Banks, Wellesley, Cambridge, Mass, USA, 1997
 I. Daubechies, "Wavelet transform, time-frequency localization and signal analysis," IEEE Transactions on Information Theory, vol. 36, no. 5, pp. 961-1005, 1990.
 P. P. Vaidyanathan, "Multirate digital filters, filter banks, poly phase networks, and applications: a tutorial," Proceedings of the IEEE, vol. 78, no. 1, pp. 56-93,1990.
 T. Q. Nguyen and P P Vaidyanathan, "Two-channel perfect reconstruction FIR QMF structures which yield linear-phase analysis and synthesis filters," IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 5, pp. 676-690,1989.
 M. Vetterli, "Wavelets and filter banks: theory and design," IEEE Transactions on Signal Processing, vol. 40, pp. 2207-2232, 1992.
 Simulink, "Simulation and model based design user's guide," Mathwork Incorporation, 2010, http://www.mathworks.com/.
 V. S. Chourasia and A. K. Tiwari, "Development of a signal simulation module for testing of phonocardiography based pre natal monitoring systems," in Proceedings of the IEEE India Council Conference (INDICON '09), pp. 1-4, Ahmadabad, India, December 2009.
 V. S. Chourasia and A. K. Tiwari, "Wireless data acquisition for fetal phonographic signals using Bluetooth," International Journal of Computers in Healthcare, vol. 1, no. 3, pp. 240-253, 2012.
Vijay S. Chourasia (1) and Anil Kumar Tiwari (2)
(1) Manoharbhai Patel Institute of Engineering & Technology, Kudwa, Gondia, Maharashtra 441 614, India
(2) Indian Institute of Technology, Rajasthan MBM College Campus, Old Residency Road, Ratanada, Jodhpur, Rajasthan 342011, India
Correspondence should be addressed to Vijay S. Chourasia; email@example.com
Received 5 March 2013; Accepted 24 April 2013
Academic Editors: Z.-L. Sun and V. Valimaki
TABLE 1: Comparison of denoising algorithms. Denoising algorithm Code MSE Universal threshold algorithm with soft sqtwolog(s) 0.8331 thresholding rule Universal threshold algorithm with hard sqtwolog(h) 0.843 thresholding rule Minimax threshold algorithm with soft minimaxi(s) 0.6295 thresholding rule Minimax threshold algorithm with hard minimaxi(h) 0.7892 thresholding rule Rigorous SURE threshold algorithm with rigrsure(s) 0.575 soft thresholding rule Rigorous SURE threshold algorithm with rigrsure(h) 0.8165 hard thresholding rule TABLE 2: Comparison of different mother wavelets for denoising of fPCG signals in terms of MSE. Wavelet family db5 Algorithm [down arrow] S1 S2 S3 S4 Signal [right arrow] (1) Sqtwolog(s) 0.8331 0.9124 1.099 0.3274 (2) Sqtwolog(h) 0.843 0.9124 1.099 0.3274 (3) Minimaxi(s) 0.6295 0.9073 1.099 0.3274 (4) Minimaxi(h) 0.7892 0.8953 1.099 0.7691 (5) Rigrsure(s) 0.575 0.8875 1.1001 0.4331 (6) Rigrsure(h) 0.8165 0.8055 1.2945 0.7476 Wavelet family sym7 Algorithm [down arrow] S1 S2 S3 S4 Signal [right arrow] (1) Sqtwolog(s) 0.5608 0.6064 1.1245 0.2663 (2) Sqtwolog(h) 0.5629 0.5934 1.1245 0.2663 (3) Minimaxi(s) 0.5615 0.602 1.1245 0.2877 (4) Minimaxi(h) 0.836 0.8363 1.5005 0.8267 (5) Rigrsure(s) 0.5033 0.5979 1.1673 0.4122 (6) Rigrsure(h) 0.8431 0.8394 1.5005 0.8369 Wavelet family coif4 Algorithm [down arrow] S5 S1 S2 S3 Signal [right arrow] (1) Sqtwolog(s) 0.3328 0.5826 0.5744 1.1137 (2) Sqtwolog(h) 0.3328 0.5906 0.5666 1.1137 (3) Minimaxi(s) 0.3328 0.5371 0.5716 1.1137 (4) Minimaxi(h) 0.7779 0.8199 0.7972 1.1137 (5) Rigrsure(s) 0.4396 0.4661 0.5531 1.081 (6) Rigrsure(h) 0.7561 1.0091 0.8633 1.2254 Wavelet family Fetal (New) Algorithm [down arrow] S5 S1 S2 S3 Signal [right arrow] (1) Sqtwolog(s) 0.2769 0.5723 0.5766 1.1054 (2) Sqtwolog(h) 0.2769 0.6001 0.5613 1.1054 (3) Minimaxi(s) 0.2837 0.5298 0.5743 1.1054 (4) Minimaxi(h) 0.8224 0.8263 0.7889 1.1054 (5) Rigrsure(s) 0.4436 0.4212 0.539 0.9998 (6) Rigrsure(h) 0.8353 1.1776 0.8519 1.2173 Wavelet family Algorithm [down arrow] S4 S5 Signal [right arrow] (1) Sqtwolog(s) 0.2615 0.2817 (2) Sqtwolog(h) 0.2615 0.2817 (3) Minimaxi(s) 0.2712 0.2817 (4) Minimaxi(h) 0.8162 0.8092 (5) Rigrsure(s) 0.2539 0.2669 (6) Rigrsure(h) 1.0062 0.8851 Wavelet family Algorithm [down arrow] S4 S5 Signal [right arrow] (1) Sqtwolog(s) 0.2712 0.2865 (2) Sqtwolog(h) 0.2712 0.2865 (3) Minimaxi(s) 0.2712 0.2865 (4) Minimaxi(h) 0.8133 0.8324 (5) Rigrsure(s) 0.2469 0.2518 (6) Rigrsure(h) 1.1123 0.8732
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Chourasia, Vijay S.; Tiwari, Anil Kumar|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2013|
|Previous Article:||Static testing of a bridge using an interferometric radar: the case study of "Ponte degli Alpini," Belluno, Italy.|
|Next Article:||Are preoperative Kattan and Stephenson nomograms predicting biochemical recurrence after radical prostatectomy applicable in the Chinese population?|