1. Introduction
With the advancement of electrical and communication technologies, various wearable devices have been developed in a variety of areas including communications, sports, medical diagnosis, and safety. The number of wearable devices will reach up to 485 million in 2018, according to a Business Intelligence report in 2013 [
1]. This number is almost 28% of the smart phone market. After the launch of the smart watch, the public has become more interested in wearable devices. In the beginning, most wearable devices were developed for fitness. After smart watch, the application areas expanded from fitness to diverse areas. Healthcare integrating wearable devices with bio-sensing technology are expected to become one of the biggest markets of the Internet of Things [
2].
For healthcare, various types of bio-signals are measured such as electroencephalography (EEG), electrocardiography (ECG), electrooculography (EOG), and photoplethysmography (PPG) signals. EEG and ECG are monitoring methods to record electrical activity of the brain and the heart, respectively. EOG is a technique for measuring the corneo-retinal standing potential that exists between the front and the back of the human eye. EEG, ECG, and EOG are measured with a differential amplifier, which registers the difference between two electrodes attached to the skin. These three measurement methods are very uncomfortable for subjects and require knowledge about the human body.
Because of their convenience and portability of measurement, PPG-related signals such as PPG and the second derivative of photoplethysmography (SDPPG) signals are widely used for wearable devices [
3]. PPG and SDPPG can monitor heart rates and cardiac cycles by optically detecting blood volume changes in the microvascular bed of tissues. The change in volume caused by the pressure pulse is detected by illuminating the skin with the light from an LED (Light emitting diode), and then measuring the amount of light either transmitted or reflected to a photodiode, as shown in
Figure 1 (Near infrared light (600–1000 nm) is suitable for PPG measurement because it permeates the skin and is absorbed into hemoglobin [
4]. Therefore, a PPG uses an LED that emits near infrared light to the skin.). For example, more blood absorbs more light so that the receiver receives weaker signal strength. Since PPG and SDPPG only use light signals, they are easier and more comfortable to measure than other kinds of bio-signals. Hence various types of wearable devices with PPG have been developed [
5]. However, vulnerability to noise is a major drawback of the PPG-related signals.
In PPG measurement, a signal includes two types of noise: background noise and movement noise [
6]. The background noise is caused by electrical noise, such as thermal noise and electromagnetic interference in cables. Movement noise is caused by voluntary or involuntary movements of a subject and affects a wide frequency range. Such noises in measurement distort PPG signal shapes, which can result in wrong diagnoses. Although background noise can be filtered out using a low-pass filter, movement noise cannot be cancelled in this manner [
5].
In previous work, there have been efforts to reduce movement noise [
6,
7,
8,
9,
10,
11,
12]. Patterson et al. [
7] studied the best location to measure PPG signals, but movement effect was not considered. Motion artifacts were studied using wavelet transform [
8]. Two PPG signals were measured at two different locations with identical PPG sensors. One was used as a reference and the other included movement noise. By comparing the two signals, the properties of movement noises were studied but there was no solution for removal of the movement noise. Kim et al. [
6] analyzed the relationship between the accelerometer-measured signal and the PPG signal. Based on this analysis and the least mean squares (LMS), a movement noise cancellation algorithm was proposed. However, the proposed algorithm improves only the peak detection accuracy of PPG signals and needs additional hardware. The adaptive step-size LMS (AS-LMS) adaptive filter [
9] was proposed to reduce motion artifacts using a LMS adaptive filter without additional hardware, but the scheme considered only limited cases such as horizontal and vertical movement and therefore cannot be applied for random movement. The periodic moving average filter (PMAF) [
10] was proposed to remove movement noise by considering the periodicity of the PPG signals. PMAF segments the PPG signal into periods and resamples each period. The algorithm mitigates a noise at an instant period by averaging multiple signals, but the noise can be propagated to the next signals, which distorts the signal shape. An adaptive spectrum noise cancellation (ASNC) approach for motion artifacts removal in PPG signals is presented by Yang et al. [
11]. However, the proposed ASNC utilizes the onboard accelerometer and gyroscope sensors to detect and remove the artifacts adaptively. Peng et al. [
12] proposed a method to extract the clean PPG signals from the motion artifacts corrupted PPG signals. The method combined temporally constrained independent component analysis (cICA) and adaptive filters. The cICA extracts the clean PPG signal, and the adaptive filters recover the amplitude information of the PPG signal. But when motion artifacts and PPG signals have dependency, the performance of the proposed method is degraded, i.e., the PPG signal shapes are distorted.
All of the previous studies mainly considered movement noise reduction for measuring heart rates from only PPG signals [
6,
7,
8,
9,
10,
11,
12] even though SDPPG signals are more useful for diagnosing diseases [
13]. The information can be extracted from the signal shapes, so if the shape of SDPPG signals is distorted the applications are limited for diagnosis. Hence, movement noise reduction is important for SDPPG signals.
To remove movement noise generally, and to maintain signal fidelity, we apply signal processing and diversity combining schemes for multiple SDPPG signals. The heart rate is almost periodic with a subject-dependent shape; therefore, a wavelet transform that has multiple bases and changes the time-frequency resolution [
14] can detect movement noise effectively. Moreover, the human heart generates a heartbeat whereby PPG signals can be measured at different locations, such as fingers, ears, and toes, as shown in
Figure 2. Hence, diversity combining schemes, as with communications, can improve the received PPG signal quality. The quality of each channel is measured by the noise level using the wavelet transform, and the algorithm chooses the best signal. This work can also be applicable to other types of bio-signals that have periodicity and diversity properties for healthcare and disease diagnosis using wearable devices.
Compared to earlier work [
15], this paper provides more detail of analysis and combining schemes, improves the proposed algorithm, and verifies the algorithm through more simulations and comparisons. The noise detection using wavelet transform and the noise reduction using a diversity combining scheme are explained in detail. The proposed algorithm is verified by comparing with previous algorithms in various environments.
The rest of this paper is organized as follows.
Section 2 describes the properties of SDPPG signals and problems in measurement.
Section 3 overviews the wavelet transform and diversity combining for application in analyzing SDPPG signals. The noise reduction algorithm with wavelet transform and selection combining is proposed in
Section 4.
Section 5 presents experiment results. Finally, this paper is concluded in
Section 6.
4. Proposed Algorithm
In this paper, an algorithm is proposed to detect movement noise and mitigate the impact of mobility on SDPPG measurement by using multipath signals and the wavelet transform. First, SDPPG signals are measured at different locations on a body, such as fingers and ears, and background noise is removed using a low-pass filter. The algorithm determines a reference signal using the measured SDPPG signals. Time offsets between received signals are aligned in order to compare the measured signals. To determine whether the measured signals include movement noise, the wavelet transform is applied to each signal. Based on the composed signals in the frequency domain, the movement noise is detected and removed. The algorithm procedure is summarized in
Figure 7 and the details of each step will be explained in the following subsections.
4.1. Multi-channel SDPPG Measurement
SDPPG signals can be measured at multiple locations on body details such as fingers, toes, ears, and wrists even if the heart is the only source. Measuring locations that are independently influenced by movement are chosen since the more independent multipath signals bring the better signal combining performance in wireless communications [
18]. For example, two fingers on the same hand may be identically influenced by arm motion and have strongly correlated movement noise. Hence, a toe and a finger, or a finger and an ear would be a better selection to measure SDPPG signals.
4.2. Selection of Mother Wavelet
Because of ambient and electric noises, measurement signals inherently include a background noise. The electric noise is generally included in the high-frequency band, and interferes with detection of PPI. Such a noise can be removed by a low-pass filter [
6]. Since human heart rates are in the range of 40 to 200 per minute [
19] and the 99% power of SDPPG is below 10 Hz described in
Section 2, a low-pass filter with a cutoff frequency of 20 Hz can remove such a noise.
After removing background noise from the measured signal, the algorithm determines a reference signal of a single period to select the mother wavelet for wavelet transform. To that end, the SDPPG signal of a few periods is averaged, which is set to a reference signal.
Since an SDPPG signal is dependent on the subject, by cross-correlation between the reference signal and mother wavelets, the algorithm chooses the mother wavelet that is the most similar to the reference signal among mother wavelets such as Harr, Daubechies, or Symlet.
4.3. Normalization and Time-Synchronization
Even if the measured SDPPG signals originate from one source (heart), the received signals may have different amplitudes and traveling times owing to the different paths taken, as shown in
Figure 2. The differences in traveling time causes phase differences between received signals, as in
Figure 8a. Such amplitude attenuation and time deviation induce errors in detection and decision during signal processing. To adjust the amplitudes, the signals are scaled to have the same amplitude as the strongest signal. To compensate for time differences between signals, cross-correlation between channel signals is performed to obtain the time difference, as shown in
Figure 9. The figure has a maximum value at 0.23 s, which means that the two signals had a 0.23-s time difference due to the travelling paths. Based on time difference computation, the algorithm aligns the times of the two signals, as shown in
Figure 8b.
4.4. Wavelet Transform of SDPPG Signals
With the selected mother wavelet, the wavelet transform is performed at each channel. As described in
Section 3, signals are iteratively decomposed up to level
N, which depends on the sampling rate, as in
Figure 10. For example, if the sampling rate of an SDPPG signal is 256 Hz, then the level
N is set to 6 so that the frequency band of the decomposed approximation components becomes between 0 and 4 Hz [
14].
Table 1 shows the frequency band according to the level of the wavelet transform at a 256 Hz sampling rate.
A signal period is detected using approximation components,
AN in level
N, for which the frequency range is between 0 and 4 Hz, whereas movement noise is detected using detail components
DN-3 in a level N-3, for which the frequency range is between 32 and 64 Hz. Since measured signals are low-pass filtered out with a cutoff frequency of 20 Hz, as explained in
Section 2 and
Section 4.2, signals higher than the cutoff frequency are mainly related to movement noise at the sampling rate of 256 Hz. In addition, the heart rate is obtained from
AN by computing PPIs.
4.5. Detection and Cancellation of Movement Noise
The algorithm detects noise in SDPPG signals during periods of pulse waves. Instead of continuous time, we consider time slots of which length is the pulse wave period. The length of the time slots can vary according to subjects. Even if signals are measured for a single subject, the duration of the pulse wave can vary according to physical conditions, so a fixed slot time cannot be used for this work. Hence, the algorithm first defines time slots of a pulse wave from the approximation coefficients that contains frequencies between 0 and 4 Hz, and determines whether the movement noise is included in each time period using the detail coefficient that contains a frequency range between 32 and 64 Hz. Each time slot is defined as the time from a peak time to the next peak time of
AN. The procedure for movement noise detection is summarized in
Figure 11.
After determining a time slot, to decide whether movement noise exists the algorithm checks whether the period of the time slot is less than a threshold (
Thp). In this paper, the threshold is set to 0.3 (In [
17], the heart rate is 40–200 beats per minute, which means that the shortest period is 0.3 s.), since the typical heart rate is between 0 and 4 Hz. If the period is less than 0.3 s, i.e., the heart rate is greater than 200 beats per minute, it is decided that the subject is abnormal or the signal is contaminated. If the period of the measured signal is greater than or equal to 0.3 s, the algorithm computes the power of the detail components between 32 and 64 Hz to detect a movement noise.
To determine how much noise is included in a signal for any duration, for a given period the signal-plus-noise-to-signal (
SNS) is defined as
As a signal and a noise cannot be separately measured, the algorithm is not able to directly use the signal quality used in communications, which is defined as the signal-to-noise ratio. SNSiS and SNSiD(N-3) denote SNS of the total signal and SNS of the detail component of level N − 3 at channel i, respectively.
If the
SNS of detail component of level
N − 3, whose frequency range is between 32 and 64 Hz, is greater than a threshold
ThD(N-3), the algorithm decides that the SDPPG signal includes movement noise during the period. Otherwise, the SDPPG signal has no noise. If only one channel has no noise among the channels, the algorithm selects the channel signal at the period without further computation. Otherwise, the algorithm calculates
SNSs of the measured signals for all channels and chooses the best channel signal as the selection combining scheme [
18]. Since the
SNS closest to one means that the signal has the highest fidelity to the reference signal, the algorithm chooses the best channel such that
Instead of SNSiS, alternative measures for selecting the channel can be used, such as SNSiD(N−3) of the SDPPG signal instead of its detail components.
6. Conclusions
Among a variety of bio-signals, PPG-related signals such as PPG and SDPPG are widely used for health care applications because of their easy and comfortable measurement. However, PPG and SDPPG signals are vulnerable to movement noise in measurement. Compared to PPG, an SDPPG signal has more information for diagnosing diseases, but is more sensitive to movement noise. Deterioration of SDPPG can result in incorrect diagnosis.
To remove movement noise, an algorithm was developed using wavelet transform and a selection combining technique. Since SDPPG is almost periodic with subject-dependent pulse shape and a band-limited signal, wavelet transform is used for detecting movement noises. SDPPG signals originating from one source (heart) can be measured at different locations so that a diversity combining scheme can be adopted to reduce movement noise. Contrast to previous algorithms, the proposed algorithm improves not only PPI performance but also the shape of SDPPG, which results in better diagnosis. In experiments, movement noise could be reduced by up to 70.89%. The proposed approach is also applicable to PPG and other bio-signals that have periodicity and measurement diversity.