Open access peer-reviewed chapter - ONLINE FIRST

Use of Transforms in Biomedical Signal Processing and Analysis

By Ette Harikrishna and Komalla Ashoka Reddy

Submitted: September 14th 2020Reviewed: May 4th 2021Published: May 29th 2021

DOI: 10.5772/intechopen.98239

Abstract

Biomedical signals like electrocardiogram (ECG), photoplethysmographic (PPG) and blood pressure were very low frequency signals and need to be processed for further diagnosis and clinical monitoring. Transforms like Fourier transform (FT) and Wavelet transform (WT) were extensively used in literature for processing and analysis. In my research work, Fourier and wavelet transforms were utilized to reduce motion artifacts from PPG signals so as to produce correct blood oxygen saturation (SpO2) values. In an important contribution we utilized FT for generation of reference signal for adaptive filter based motion artifact reduction eliminating additional sensor for acquisition of reference signal. Similarly we utilized the transforms for other biomedical signals.

Keywords

• Fourier transform
• biomedical signals
• electrocardiogram signal
• photoplethysmographic signal
• wavelet transform

1. Introduction

The essence of mathematical design cannot be ignored in the analysis of real world engineering applications i.e. the research in engineering and mathematics is a two way parallel track that interrelates and coordinates towards value added research. In specific, the use of transforms in the field of electrical, electronic and communication engineering is unimaginable. In the present scenario of Covid-19 pandemic, world is looking to sustainable development of biomedical devices for critical monitoring and efficient vaccination for human survival [1, 2, 3, 4]. In general, the Fourier transform (FT) is a mathematical tool which transforms the time domain signal into a frequency domain representation used in analysis of biomedical, wireless communication, signal and image processing applications. In literature, many researchers had used this tool in frequency domain analysis of all biomedical signals like electrocardiogram (ECG), photoplethysmographic (PPG) and blood pressure (BP).

In continuation to FT, different transforms were developed to analyze and design of various applications based on the requirement [5, 6]. In general, the FT is used in analysis of stationary signals; the wavelet transform (WT) is a mathematical tool used in analysis of both stationary and non-stationary signals. Discrete wavelet transform (DWT) used in enormous application in various engineering fields.

So, in this chapter we addressed some of the research challenges in ECG and PPG signal processing using Fourier and Wavelet transforms.

2. Transforms and its use in biomedical signal analysis

2.1 Transforms

Jean-Baptiste Joseph Fourier a French mathematician had developed theoretical and mathematical framework for Fourier analysis and harmonic analysis, which was laid down foundation to other transforms and important applications. Forward and inverse transform of a continuous time Fourier transform (CTFT) for a time domain signal x(t)is defined as

Xjω=xtejϖtdtE1
xt=12πXjωejωtdωE2

Here, ejϖtis called as a basis function for CTFT.

For ease of processing the continuous time signal is converted into discrete time, then the corresponding forward and inverse transform of a discrete Fourier transform (DFT) for a discrete time signal x(n)is defined as

Xk=n=0N1xnej2πnkN;fork=1,2,3NE3
xn=1Nk=0N1Xkej2πnkN;forn=1,2,3NE4

Here, ej2πnkNis called as a basis function for DFT.

Similarly, the continuous wavelet transform (CWT) is defined as

Fab=1aψtabftdtE5

where, f(t)is a time domain signal, ψtis a basis function also called as mother wavelet; F(a, b) is WT of a signal f(t); ‘a’ and ‘b’ are shifting and scaling parameters respectively;

Discrete WT (DWT) uses a series of low pass filters (LPF) and high pass filters (HPF) to decompose the signal of interest into different scales as approximate (Aj) and detailed coefficients (Dj). The output coefficients of the LPF are called approximations while the output coefficients of the HPF are called details.

A 3-stage DWT decomposition tree is illustrated in Figure 1.

After first level of decomposition A1 and D1 will be the outputs, D1 will be stored in Cmatrix as shown in Figure 1. A1 will be further decomposed in 2nd level of decomposition as D2, A2. D2 will be stored in Cmatrix, and A2 will be decomposed in 3rd level of decomposition as D3, A3. Then A3 and D3 will be stored in Cmatrix, i.e. Cmatrix consists of concatenated A3, D3, D2, D1 coefficients. Lmatrix stores the length of each corresponding approximate and detailed coefficient as in Cmatrix.

C matrix gives the concatenated approximate and detail coefficients. Lmatrix gives length of each coefficient. Cand Lmatrices along with suitable filters will be used to get back the original time domain signal f(t).

2.2 Biomedical signals

The Electrocardiogram (ECG) is an electrical manifestation of contractile activity of the heart, is a representation of instantaneous electrical activity of the heart during its contraction and relaxation over a period of time. A standard 12 leads ECG is recorded with the help of surface electrodes placed on the limbs and chest. In general, ECG is a widely accepted diagnosis tool in clinical validations of heart related diseases. Figure 2 shows a typical ECG signal, marked with all the characteristic waves and durations such as P wave, QRS complex wave, T wave and ST segment [7].

Monitoring of blood oxygen saturation (SpO2) is one of the important parameter to know the health status of Covid-19 affected patient. Pulse oximeter uses sensor probes to record photoplethysmographic (PPG) signals so as to estimate the SpO2 values. Typical recorded PPG signal is shown in Figure 3.

2.3 Biomedical signal processing using transforms

Fourier transform (FT) is used to analyze the behavior of biomedical signals in frequency domain. In Matlab FFT command can be used to get the frequency domain signal.

Following is the sample code to plot time and frequency domain signals.

ecg_signal_fft=fft(ecg_signal);

figure(1)

subplot(211)

plot(ecg_signal)

subplot(212)

plot(ecg_signal_fft)

Some additional modifications will be done in the program to get the following plots.

In the bottom trace of Figure 4, it can be seen that the frequency domain components of ECG signal are extending from 0 to 100 Hz, the main source of noise 60 Hz power line interference (PLI) the spike at 60 Hz can be observed. So, the complete frequency domain behavior of signal can be computed using fft command.

Similarly, the frequency domain components of PPG signal is shown in Figure 5. It can be seen from the bottom trace that frequency components present in PPG signal are pulse rate or heart rate component and MA noise component. Likewise it can be continued for many biomedical signals to see the frequency domain behavior of the signal.

So, use of Fourier transform in de-noising of signal can be described as shown in Figure 6.

In general adaptive filter provide a viable solution when signal and noise are in same frequency range. Adaptive filter requires a two input signals [8].

For example, in a power line interference cancelation from ECG signal, one is recorded noisy ECG signal and other is power line noise. So, here a synthetic noise reference signal will be generated using Fourier transform which will be used as another input to the adaptive filter will potentially eliminate additional sensor for acquisition shown in Figure 7.

A synthetic noise reference signal is generated for use in adaptive filtering without using any extra hardware. It is generated from the motion corrupted signal in the following way. The frequency spectrum of MA corrupted PPG signal consists of various frequency components, the pulsatile (0.5–4 Hz), respiratory activity (0.2–0.35 Hz) and MA noise component (0.1 Hz or more) information. By setting the co-efficients of cardiac and respiratory activity frequency components in the spectrum of MA corrupted PPG to zero, a modified spectrum corresponding to noise is obtained. By applying inverse Fourier transform to this modified spectrum, a synthetic noise reference signal is generated. The corresponding adaptive filter is shown in Figure 8.

With the help of LMS adaptive algorithm, MA noise is removed by estimating the synthetic noise reference signal and adapting the filter coefficients based on filter order. The necessary equations to implement the proposed method are given below:

Ŝn=Sn+NnN̂nE6
N̂n=i=0LwiNRniE7
win+1=win+2μSnNRniE8

where i: 0,1,2, …, L, L: filter order, S(n) + N(n): MA corrupted PPG signal, Ŝ(n) =e(n): MA reduced PPG signal, N̂(n): estimated synthetic noise reference signal, and NR(n): the synthetic noise reference signal.

The result of above methodology is presented in Figure 9, below. Figure 9(b1)represents the generated synthetic noise reference signal which potentially eliminates the additional sensor for data acquisition.

The biomedical signals such as ECG and PPG signals are quasi periodic signals i.e. the period of the signal continuously changes with time, but it is a periodic signal. In general, the pure periodic signals are stationary in nature means its period will be constant irrespective of time. So, Fourier transform is not sufficient to analyze the quasi-periodic signals. Wavelet transform will provide a viable solution to the same [9, 10].

The general de-noising procedure follows the steps described below and shown in Figure 10.

Decomposition:Choose a wavelet and choose a convenient level N for decomposition. Compute the wavelet decomposition of the signal s at level N.

Thresholding detail coefficients:For each level from 1 to N, select a threshold and apply soft or hard thresholding to the detail coefficients.

Reconstruction:Perform the wavelet reconstruction using the original approximation coefficients of level N and the modified detail coefficients of levels from 1 to N.

There are two important steps: how to choose the threshold, and how to perform the thresholding [10]. In hard thresholding process, the elements whose absolute values are lower than the threshold will be set to zero. Soft thresholding is an extension of hard thresholding, first setting to zero the elements whose absolute values are lower than the threshold, remaining coefficients are compressed.

The same results were presented in Figures 1113. Similar results were presented for PPG signal as shown in Figure 14.

3. Conclusions

Transforms like Fourier and wavelet transforms were used in biomedical signal analysis and processing. Fourier and wavelet transforms were utilized to reduce motion artifacts from PPG signals so as to produce correct blood oxygen saturation (SpO2) values. In an important contribution we utilized FT for generation of reference signal for adaptive filter based motion artifact reduction eliminating additional sensor for acquisition of reference signal.

chapter PDF

More

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Cite this chapter Copy to clipboard

Ette Harikrishna and Komalla Ashoka Reddy (May 29th 2021). Use of Transforms in Biomedical Signal Processing and Analysis [Online First], IntechOpen, DOI: 10.5772/intechopen.98239. Available from: