Next Article in Journal
Hierarchical Area-Based and Path-Based Heuristic Approaches for Multirobot Coverage Path Planning with Performance Analysis in Surveillance Systems
Previous Article in Journal
Comparing Three Methods of Selecting Training Samples in Supervised Classification of Multispectral Remote Sensing Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Signal-to-Noise Ratio Indicator for Wearable Bioimpedance Monitoring

Institute of Electronics and Computer Science (EDI), Dzerbenes 14, LV-1006 Riga, Latvia
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(20), 8532; https://doi.org/10.3390/s23208532
Submission received: 13 September 2023 / Revised: 3 October 2023 / Accepted: 13 October 2023 / Published: 17 October 2023
(This article belongs to the Section Electronic Sensors)

Abstract

:
Bioimpedance monitoring is an increasingly important non-invasive technique for assessing physiological parameters such as body composition, hydration levels, heart rate, and breathing. However, sensor signals obtained from real-world experimental conditions invariably contain noise, which can significantly degrade the reliability of the derived quantities. Therefore, it is crucial to evaluate the quality of measured signals to ensure accurate physiological parameter values. In this study, we present a novel wrist-worn wearable device for bioimpedance monitoring, and propose a method for estimating signal quality for sensor signals obtained on the device. The method is based on the continuous wavelet transform of the measured signal, identification of wavelet ridges, and assessment of their energy weighted by the ridge duration. We validate the algorithm using a small-scale experimental study with the wearable device, and explore the effects of variables such as window size and different skin/electrode coupling agents on signal quality and repeatability. In comparison with traditional wavelet-based signal denoising, the proposed method is more adaptive and achieves a comparable signal-to-noise ratio.

1. Introduction

Bioimpedance monitoring is a valuable technique to noninvasively assess various physiological parameters such as body composition, hydration levels, heart rate and breathing. However, signals acquired in real-world experimental conditions are often corrupted by noise, posing challenges to the reliability of derived quantities. Therefore, it is essential to evaluate the quality of measured signals to ensure the accuracy of the derived physiological parameter values. This includes applications of bioimpedance signals such as heart-rate monitoring [1], which is the focus of this article. Although heartbeat signals obtained from different individuals tend to be similar, the shapes of individual heartbeat signals are different and variable, and it is important to accurately distinguish between the authentic signal from the interfering noise in terms of heart rate.
Existing approaches for signal quality assessment in bioimpedance monitoring mainly focus on traditional signal-to-noise ratio (SNR) estimation techniques such as autocorrelation [2]. However, these universal methods fail to take into account the specific shape and features of bioimpedance signals; in short, the issue lies in the difficulty of identifying the reference signal that should be compared with the noise. With regard to our application—heart rate detection from the bioimpedance measurements—a problem is that heart rate signals are non-stationary and not normally distributed, as shown further in this paper. Therefore, they are challenging to the existing SNR estimation techniques.
The primary objective of this study is to propose an adaptive signal-to-noise ratio indicator for bioimpedance signals recorded on a wearable device, with the target application of heart-rate monitoring. We focus on signals from the wrist, as it is a typical location for wearable devices. Our contributions are:
  • Present a custom wearable device for bioimpedance monitoring.
  • Develop a SNR estimation method based on continuous wavelet transform (CWT) to identify wavelet ridges and measure their energy.
  • Validate the proposed method using bioimpedance signals obtained in a small-scale experimental trial, with photoplethysmogram (PPG) signals as the reference.
  • Compare the performance of the proposed adaptive method with a traditional wavelet-based signal denoising approach.
  • Compare the SNR depending on different electrode-to-skin contact materials and other parameters.
The study shows that our adaptive method accurately captures the SNR information. Unlike classical wavelet-based signal denoising and signal-to-noise ratio estimation methods, our adaptive method does not require setting fixed threshold levels or selecting proper threshold types. The PPG signals serve as a reference as they have relatively lower amounts of noise, and as such can validate the separation of bioimpedance signals in “signal” and “noise” components. Developing methods for accurate heart rate estimation based on this decomposition in the signal and noise is in our future plans.
This paper is structured as follows. Section 2 describes the related work. Section 3 describes the experimental platform used in this study. Section 4 explains the novel SNR estimation method and the experimental setup, Section 5 shows the results of this approach, and Section 6 concludes the paper.

2. Related Work

Various techniques are commonly employed to denoise biological signals. These include the application of highpass filters [3], lowpass filters [4,5], notch filters [6,7], adaptive filters [8,9], empirical mode decomposition [10,11], blind source separation methods [12] based on the Independent Component Analysis algorithm, and the Fourier decomposition method [13].
Adaptive filters eliminate noise by utilizing a correlated reference signal alongside the original signal. However, their suitability for real-time applications is limited due to the requirement of a reference signal [14]. Nevertheless, the combination of blind source separation methods and adaptive filtering can yield favorable outcomes by generating a reference signal without the need for an implantable sensor [14]. It is important to note, though, that blind source separation methods are computationally intensive, resulting in higher power consumption and processing costs. Consequently, they are not ideal for real-time monitoring.
To address the limitations of the aforementioned techniques and effectively denoise non-stationary signals while preserving the essential signal characteristics, the wavelet transform (WT) methods, namely, discrete wavelet transform [15,16], stationary wavelet transform [14,17] and synchrosqueezed wavelet transform [18] had been used. These approaches operate by decomposing the signals using wavelets and subsequently applying coefficient thresholding to remove noise. These methods have proven to be successful in eliminating the specified types of noise while maintaining the integrity of important signal attributes [14].
Thresholding plays a crucial role in wavelet transform-based signal denoising, as it determines the removal of WT coefficients below the threshold level. Therefore, careful selection of the threshold is essential. The use of the universal threshold rule, proposed by Donoho and Johnstone, has been widely adopted for signal denoising [19,20,21]. However, it should be noted that the universal threshold method is not suitable when the noise variance of the signal is unknown [18].
Attempts to overcome the limitations of conventional thresholding methods were carried out using metaheuristic algorithms. Li et al. [18] use the differential evolution algorithm to dynamically determine the optimal threshold, followed by the application of the soft threshold method to threshold the coefficients obtained from the wavelet transform, similar to the approach adopted in [22].

3. Experimental Platform

As part of the work, we have developed a prototype of a novel wearable device (Figure 1). The main novelty of the device lies in the combination of bioimpedance sensing with a nine-axial inertial measurement unit (IMU) sensor. Coupled with a next-generation, powerful and energy-efficient ARM Cortex-M33 dual-core microcontroller (MCU), this makes it possible to run on-board sensor fusion and activity recognition software [23], making the wearable device suitable for novel applications in health and behavior monitoring.
The main features of the wearable device are
  • nRF5340 System-on-Chip (SoC) with dual-core ARM Cortex-M33;
  • MAX30001 bioimpedance chip;
  • ICM-20994 nine-axial IMU chip, with accelerometer, gyroscope, and magnetometer sensors;
  • Li-Ion battery;
  • Power sourcing from the USB or from the battery;
  • LEDs and a vibration motor for feedback to the user;
  • External flash for additional on-board storage;
  • User button.
In the short-term future, we aim to develop a new version of the prototype with a smaller form factor. The current PCB size of 53 × 55 mm is meant to enable easy testing and debugging, rather than to be convenient for wearing.
For the experiments described in this paper, we connect two stainless steel electrodes to this wearable device.

3.1. Component Selection

3.1.1. Microcontroller and Radio

One of the most important design choices is to find a suitable microcontroller. The latest generation microcontrollers have 32-bit cores and offer much better performance/energy consumption tradeoff than the older 8-bit and 16-bit MCU. Microcontrollers several families were considered, all based on 32-bit ARM Cortex-M cores.
Table 1 and Table 2 show the comparison. EFR32 and nRF5340 have the best overall results across the variety of metrics that are important for our goals: low energy consumption, sufficient RAM and flash memory size, and high performance. From these two, nRF5340 is selected as the one with the better software ecosystem support. While the EFR is the absolute leader in energy efficiency from the MCU considered, it does not have as many existing hardware drivers and experienced programmers as the nRF MCU.

3.1.2. Bioimpedance Measurement Unit

The MAX30001 is a full analog front-end (AFE) solution for wearable applications that includes biopotential and bioimpedance measurements. For bioimpedance input, ECG input, and signal generator output, it features three distinct isolated channels. The biopotential channel detects the electrocardiogram (ECG) waveform, heart rate, and pacemaker edges, whereas the bioimpedance channel measures tissue impedance, such as respiration. ECG and bioimpedance AFE of clinical quality with high-resolution data converter. The 32-word ECG and 8-word bioimpedance first-in first-out buffers (FIFOs) allow the MCU to be turned off for 256 ms while still collecting data. Bioimpedance measurement parameters on the MAX30001 frontend are set as follows: the measurement frequency is 80 kHz, the injected current does not exceed 32  μ A, the gain is 80 v/v, low-pass filter 16 Hz, high-pass filter 7200 Hz, low noise mode selected. High input impedance, low noise, adjustable gain, multiple low-pass and high-pass filter choices, and a high resolution analog-to-digital converter are other features worth mentioning. Critcher et al. [24] provides a far more extensive comparison and examination of this IC, detailing the capabilities supplied by this component.
The existing measurement system uses a bipolar electrode measurement connection, which gives acceptably good results for our requirements. Although the tetrapolar configuration is recognized as more resistant to noise and has a better signal quality [25], the bipolar system is still used for convenience and simplicity [26]. Also, the SNR ratio in the tetrapolar configuration is not always unequivocally better [26].

3.2. Wearable Device Design

3.2.1. Design Overview

The wearable device includes the following essential components: nRF5340 central microcontroller, MAX30001 Bioimpedance measurement integrated circuit (IC), and ICM-20948 IMU sensor module.
nRF5340 is a wireless SoC with two ARM Cortex-M33 processors. It supports Bluetooth 5.3, high-speed SPI, QSPI, USB, and an operational temperature of up to 105 °C. The application processor is performance-optimized and may run at 128 or 64 MHz thanks to voltage-frequency scaling. It has 1 MB of flash memory, 512 kb of RAM, a floating-point unit (FPU), an 8 kb 2-way associative cache, and digital signal processing (DSP) instructions. The network processor runs at 64 MHz and is designed to be energy efficient (101 CoreMark/mA). It has a flash memory of 256 kb and a RAM of 64 kb. The SoC supports Bluetooth 5.3 and other protocols including NFC, ANT, IEEE 802.15.4 at the 2.4 GHz frequency band.
The sensor part consists of two data collection channels. The first is the IMU (ICM-209948) Motion Tracking system, which can transmit both the raw data and the results from signal processing on the IMU chip itself. The MAX30001 IC is the second sensor channel. Depending on where it is positioned and what signals are examined, the capacity to measure bioimpedance signals permits the detection of human breathing, heart rate, and other parameters. Both sensors support SPI data transmission, and system interrupts can be used to wake up the main microcontroller from sleep mode.
The device has two power options: micro USB input, and a lithium ion battery. An on-board switch selects the power supply mode. Additionally, 1.8 V voltage regulation IC is included, which then provides the main voltage to all other nodes on the board.
For reading the bioimpedance signals, two stainless steel electrodes are connected to the MAX30001 chip. The considerations behind the electrode design are described in [2].

3.2.2. Energy Consumption

The component selection allows to achieve a long battery life. Using high-accuracy multimeter (Keysight Technologies 34460A, Santa Clara, U.S.) measurements and values from datasheets, we found out that the wearable device with the measurement functionality enabled and Bluetooth Low Energy (BLE) active consumes 3.16 mA average current. When performing bioimpedance measurements without BLE, the current consumption is 2.7 mA. The nRF5340 module consumes approximately 1.5  μ A in its sleep mode, while the MAX30001 module consumes 0.6  μ A during sleep. During its active operational mode, the MAX30001 module typically consumes 163  μ A.

3.3. Software

In terms of an operating system for the device, we selected the Zephyr OS, a small-footprint OS designed for use on resource-constrained and embedded systems (https://docs.zephyrproject.org/latest/introduction/index.html (accessed on 14 October 2023)). The main benefits of Zephyr include an extensive ecosystem, inter alia direct support from the Nordic Semiconductors company itself, as well as a complete BLE network stack and support for a large number of hardware platforms, including nRF5340-based platforms. Zephyr also features primitives required to work with hardware in a platform-independent way, for example, abstractions for the SPI and I 2 C protocols, for general purpose input/output (GPIO) and interrupt management, and similar.

4. Methods

Methodology of the signal quality estimation applicable in human heart beat monitoring is schematically illustrated in Figure 2 and described below.

4.1. Band-Pass Filtering

To eliminate undesirable signals such as low-frequency and high-frequency hiss from the data, a standard bandpass filter is employed. The filter’s settings were established based on a frequency response range that was determined using the upper and lower limits of the heart rate, between 30 beats per minute and 300 beats per minute, assuming that these limits were unattainable. The filter is designed using an infinite impulse response (IIR) filter type, as specified by parameter bandpassiir in Matlab filter designer. The frequency range is defined by the following parameters: The lower frequency edge of the first stop-band set to 0.1 Hz, but the lower frequency edge of the pass-band set to 0.5 Hz. The upper frequency edge of the pass-band, set to 5 Hz, but the upper frequency edge of the second stop-band, set to 6 Hz.
To meet the desired frequency response, the filter has specific magnitude constraints, defined by stop-band attenuation. The amount of attenuation in the first stop-band is set to 40 dB and the amount of ripple allowed in the pass-band is set to 1 dB. The amount of attenuation in the second stop-band is set to 50 dB. The filter is designed using the elliptic filter design method. The MatchExactly parameter indicates that the filter design should match the pass-band frequencies exactly. Finally, the sample rate of the signal to be filtered is given as a constant value of 64 samples per second.

4.2. Wavelet-Based Signal Quality Estimation

Continuous wavelet transform (CWT) is a time-frequency analysis technique allowing us to observe changes of the signal’s frequency components with time. CWT is performed on the measured signal y ( t ) to extract the wavelet transform coefficients W y ( s , b ) at each scale parameter s and time instant b according to the formula
W y ( s , b ) = y ( t ) × 1 | s | × ψ * ( t b s ) d t
= y ( t ) × ψ s , b * d t
where ψ are wavelet functions—special functions satisfying conditions of zero mean and finite energy. Wavelet scale parameter s is related to signal frequency f according to this scale-to-frequency conversion
f = ω 0 2 π × s
where ω 0 is a central frequency of the wavelet function. CWT decomposition of the measured signal in time and scale (frequency) domains is conveniently viewed through wavelet scalogram (WS), which shows proportion of signal energy in time-frequency plane. The highest signal energy concentrates in so-called wavelet ridges, where high energy of a signal corresponds to a specific scale parameter (related to signal frequency) and is traced over the time instants of the signal. It is an instantaneous frequency of the signal. The WS is calculated as
W S i j = | W y ( s , b ) | 2 i = 1 M j = 1 N | W y ( s i , b j ) | 2
where i = 1 , , M is a number of scales and j = 1 , , N is a number of time instances. The wavelet ridges are identified from two conditions:
1.
Derivative of absolute values of the wavelet coefficient with respect to the scale parameter s as calculated for all time instants b of the signal is zero
d d s | W y ( s , b ) | = 0
2.
The second derivative of absolute values of the wavelet coefficient with respect to the scale parameter s as calculated for all time instants b of the signal is negative
d 2 d s 2 | W y ( s , b ) | < 0
Derivatives of W y ( s , b ) of the first and second order with respect to s are calculated. For clarity, these quantities are denoted with D 1 and D 2 , respectively. Time instants and respective scale parameters s corresponding to negative values of D 2 are found and are denoted with b * and s * , respectively. The identified wavelet ridges are W y ( s , b ) = W y ( s * , b * ) .
After the identification of ridges, their duration in number of samples is normalized to the duration of an entire signal.
b n = j = 1 N b j b L
where b L is the duration of signal. The rationale is that the wavelet ridge corresponding to the heartbeat should persist through all of the duration of the signal, while noisy ridges may not have such a stability in time. Then, the energy of each ridge is computed by summing up the wavelet scalogram coefficients corresponding to the ridge scale or along the time axis as j = 1 N | W y ( s k * , b j * ) | .
The energy of the wavelet ridge characterizes the intensity of a particular signal component. Next, a normalized energy is calculated for each ridge by dividing the ridge energy by the total energy of the signal, yielding a percentage of energy associated with each wavelet ridge.
e n = 100 × j = 1 N | W y ( s k * , b j * ) | i = 1 M j = 1 N | W y ( s i , b j ) |
Afterwards, this normalized energy is weighted by the duration of the wavelet ridge to obtain the weighted ridge energy as
e w = e n × b n
The weighted ridge energy takes into account both ridge stability in time and also intensity of the ridge. Hence, the signal component corresponding to heartbeat should have a very high weighted energy compared to other ridges.
A cornerstone of the proposed signal quality estimation method is estimation of noise content. Here, the median filter is applied to the weighted ridge energies. The choice of median filter instead of other filters is based on the distribution of e w values. As this distribution is not normal, since relatively few samples have high values, median filter is more robust to such outliers compared to, for example, mean filter. Median filter, however, requires us to set a window length in which the median operation is to be performed. The choice of window length is optimized using k-fold cross-validation technique in the following manner: The weighted energy data e w have been split into 10 folds, where 9 folds were used for applying the median filter to the data and obtaining filtered data e w M . f i l t . at pre-defined window length values. It constitutes a training data set. The remaining one fold was used for testing where a mean-squared error metric between the filtered data and data in this one fold (testing data set) was computed. This operation was repeated by cycling through all 9 folds and the one fold for testing was fixed. The following 9 iterations, however, cycled through the fold for testing so that eventually all 10 folds were used for testing once. This scheme is shown in Figure 3.
When the optimum value of window length has been chosen. The noise component of a signal is determined as a difference between the weighted ridge energy e w and a filtered e w with a median filter as n o i s e = i = 1 M | e w , i e w , i M . f i l t . | 2 .
The proposed signal quality indicator, also called signal-to-noise-ratio (SNR) is defined as
S N R = 10 × log 10 i = 1 M e w , i 2 n o i s e
The numerator expresses the total energy of the signal y ( t ) across all frequencies and time instants, whereas the denominator reflects the part of the signal related to noise—information-carrying signal components subtracted from the total signal. Following other definitions of the SNR, its proposed version is also expressed in decibels.

4.3. Experimental Measurements

The experiment was carried out with the wearable device shown in Figure 1 and described in Section 3, Experimental Platform. The complete measurement system also includes a BLE receiver module, a computer as a data receiving unit, as well as an impedance model of organic human tissue as Z tissue, shown in Figure 4. The prototype includes a pair of stainless steel electrodes that are supported by a layer of foam to improve flexibility and reduce pressure on the skin. The dimensions of prototype PCB are 55 mm × 52 mm × 30 mm, including the battery and electrodes. The initial electrode size was approximately 8 mm × 18.5 mm, but this was later changed to larger electrodes with dimensions of 18.5 mm × 44 mm to improve the contact area with the skin. The decision to use stainless steel electrodes was based on previous research on electrode materials and their properties, with particular emphasis on their effectiveness in heart rate detection [2].
The MAX30001 integrated circuit was configured to sample bioimpedance data at a rate of 32 or 64 samples per second. Multiple readings were collected and transmitted to a PC via the nRF5340 using Bluetooth Low Energy (BLE) protocol. A Python script was used to collect data at the user end. The collected data were then transferred to Matlab for analysis, where heart rate was extracted and analyzed over time.
Bioimpedance signatures were measured on nine individuals (subjects I to IX). Informed informed consent was obtained from all experimental subjects. Biological data of the test subjects are shown in Table 3. Each measurement session consisted of simultaneous bioimpedance and PPG measurements, where a PPG sensor was held between the thumb and an index finger. The wearable device was put on the wrist parallel to ulnar and radial arteries. The test subjects were sitting still during the measurement session. The effect of two different couplings between skin and electrodes was explored. A total of 10 measurements were performed using a gel “E.C.G. AND TENS GEL” [27] and 10 measurements were performed using hydrogel sheets “MODEL MA1017” [28] with size 4 × 6 cm, cut to precisely match the electrode size. Ambient temperature during the measurement sessions was approximately 19 °C. For clarity, only the results for subject I are shown in the following plots.
TENS GEL that we use in experiments [27] is specifically designed for this purpose, does not contain formaldehyde and is free from salts, completely hypoallergenic, water-soluble, resistant to drying, odourless, and easily removable. It does not stain, does not grease, does not ionize and does not oxidize over time. It does not contain active ingredients that could damage the electrodes in any way. The manufacturer did not provide much information about the hydrogel used in our experiments; however, similar hydrogel materials and applications are discussed in [29].

5. Results

5.1. Measured Signals

The recorded raw bioimpedance measurements are presented in Figure 5. Out of the 10 measurements, only three are presented—the first, the fifth and the last (10th) signal with the highest SNR, signal with the lowest SNR and signal obtained as an average of all 10 signals. These three signals are shown for both gel and hydrogel couplant cases. Bioimpedance values have a settling time—a trend to decrease over time and asymptotically tend to an equilibrium value. Bioimpedance starting values differ for different coupling agents. In the case of hydrogel, the starting bioimpedance value is about two to three times higher compared to measurements with gel. Moreover, the slopes of bioimpedance that decay over time are also different. In the case of gel coupling, the rate of decay is lower and small-amplitude oscillations can be seen for the fifth and tenth measurements. The signal obtained after placing the electrodes on the skin can be affected by the resulting capacitance of the skin layers. This can cause the signal to appear as a falling exponent, which can obscure the true nature of the signal and affect the accuracy of heart rate measurements, as mentioned in [30].
Bioimpedance signals processed with band-pass filtering are shown in Figure 6. This step has removed the slope of the raw signals as well as some noise outside the pass bands of the filter. For ease of comparison, these processed signals are also normalized to the range from 1 to 1. Periodic oscillations corresponding to heartbeat can be clearly seen. Zoom-in on the first 10 s of measurement is shown in Figure 7. A clear pattern of hemodynamic activity with the fiducial points can be seen.
Results of raw PPG signals measured simultaneously with a bioimpedance are shown in Figure 8. These signals correspond to the ones presented for the bioimpedance. Several observations can be made. Firstly, compared to bioimpedance, the values of PPG do not change significantly with time. Secondly, there is a pronounced baseline wander. Thirdly, apart from higher frequency oscillations associated with a heart beat, smaller frequency components are presumably related to respiration. The respiration frequency is in the range 0.15–0.4 Hz [31]. Band-pass filtering of PPG signals, as shown in Figure 9 has removed a wander of baseline level and the zoom-in on the first 10 s is depicted in Figure 10.

5.2. Signal Stationarity and Normality

A more straightforward method of signal quality assessment in periodic signals is based on the autocorrelation function. This method involves computing the autocorrelation function of the signal. The SNR indicator is computed as S N R = R ( τ = 0 ) 1 N i = 1 N R ( τ i ) , where R is the correlation coefficient and τ are time lags of the signal with itself. The idea is that by combining the fact that, theoretically, the autocorrelation of noise is zero [32] and diminishing of autocorrelation values with increasing time lags, the mean over all time lags except at τ = 0 corresponds to noise. Compared to the classical wavelet-based denoising method, this method requires less computational resources, since there is no need to decompose the signal into frequencies. However, the restriction of the autocorrelation method is that the signal has to be periodic, stationary and it should not be correlated with the noise.
The measured bioimpedance signals were tested regarding stationarity and normality in order to asses the appropriateness of the autocorrelation method for estimation of the SNR of bioimpedance signals. For this task, Kolmogorov–Smirnov, Anderson–Darling and Jarque–Bera statistical tests were conducted on the normalized bioimpedance signals to assess their normality. On the other hand, Augmented Dickey–Fuler test and Variance ratio test were carried out to test the stationarity of the recorded signals. Multiple tests were performed to increase the assurance of the judgment on the nature of the signals. The results on the signals measured on subject I are shown in Table 4 and Table 5. The tests were performed at a five percent significance level. All of the included signals are non-normal, since the associated p-value is less than α and statistic value is larger than the critical value. The same is true for non-stationarity tests. Hence, the measured bioimpedance signals are non-normal and non-stationary at α = 0.05 and it is not correct to apply simpler signal quality estimation methods, such as the autocorrelation method.

5.3. Noise Detection

Further, both band-pass-filtered bioimpedance and PPG measurements are processed with CWT. Scalograms of the bioimpedance signals are shown in Figure 11. A dominating wavelet ridge can be seen around 1 Hz on the frequency axis and can be traced for all time instants in all measurements conducted on subject I. It is clearly related to heart beat as the frequency if 1 Hz approximately matches the human heart beat at rest and the amplitude of modulus of CWT coefficients related to this signal component is much higher compared to other signal components. Relatively low amplitude CWT coefficients are observed at higher frequencies. Analogous results for PPG signals are illustrated in Figure 12. In this case, scalograms are dominated by a single signal component around 1Hz. These scalograms contain less noise compared to the ones of bioimpedance signals as is also evidenced by higher SNR indicator values.
Results of window length optimization for median filtering through 10-fold cross-validation are depicted in Figure 13. Here, window length values in a range from 1 to 15 were tested in a total of 1000 independent runs. A high scatter of values as indicated by the coefficient of variation indicates that no particular value is favorable over others. By increasing the number of runs from 100 to 1000, the coefficient of variation decreased by only approximately 1 percent. As there is no significance to a particular value of window length, it was decided to select the one that yields the maximum SNR.
Median filtering results are illustrated in Figure 14. For clarity, these results are shown only for a trial I (gel) of subject I as intermediate results. The blue bars denote the scale parameters satisfying the ridge conditions. However, the duration of the associated ridges are not necessarily high. The red bars correspond to the relative energy expressed as percentage from the total energy of a signal. An overlay of both plots gives a visual indication of a mutual relationship between the ridge duration and its energy. The bottom plot shows weighted ridge energy (multiplication of both quantities from the upper plot). Hence, a ridge energy that is weighted by duration (red bars) is obtained. The maximum weighted energy is at a scale parameter of 51, which, in turn, corresponds to 0.868 Hz through scale-to-frequency conversion. Blue bars, on the other hand, are obtained from applying a median filter to the weighted energy data. A logarithmic scale is used so that the small and large amplitudes are clearly visualized.

5.4. Signal Quality

The effect of window length of the median filter on estimation of the SNR was carried out. The SNR values versus window lengths ranging from 2 to 40 were calculated on measurements performed on subject I and shown in Figure 15. The first, fifth and tenth measurements are displayed both for gel and hydrogel couplants. It can be seen that the maximum SNR is consistently achieved at a window length of three. As the window length increases, the SNR values decrease and eventually stabilize where a further increase in window length does not have an impact on the SNR. In those cases, the weighted energy data are over-smoothed.
The final results of the obtained SNR values are shown in Figure 16. Here, the mean values of the SNR are illustrated for all test subjects for both couplants along with vertical errorbars. The errorbars were calculated as confidence intervals C I = t 1 α , ν × σ N and the final result can be presented as S N R = S N R ¯ ± C I .
The obtained mean SNR values range between 19 dB and about 24 dB. There are no significant differences in SNR results between gel and hydrogel couplants for all subjects—all differences are within the error bounds. Moreover, no significant differences are observed between test subjects. Hence, the choice of couplant has an impact mostly on practical aspects of measurement, such as availability, ease of application and personal comfort of the test subject rather than the actual quality of the measured signal. These bioparameters might show a correlation with specific features of the bioimpedance signals rather than the SNR values.

5.5. Classical Wavelet Transform Approach

The procedure of signal denoising using the classical wavelet transform approach is as follows. Signal decomposition into approximation and detail coefficients at different levels is performed using a discreet wavelet transform (DFT) with a selected mother wavelet function. The maximum level of decomposition depends on the length of the signal. A rule of thumb is to set it as equal or less than log 2 ( N ) , where N is length of the signal. The threshold level is calculated as
T = T M   ×   threshold_type  
where T M is the threshold multiplier, and “threshold_type” is the type of the threshold. Common choices include soft, hard or universal. The universal threshold is defined by M A D ( W x ) = m e d i a n ( | W x m e d i a n ( W x ) | ) 0.6745 . and M A D is median absolute deviation of wavelet decomposition coefficients W x of signal x representing a noise level in the signal.
A drawback of the classical wavelet transform approach for SNR estimation is that many parameters need to be provided by the user or found via an optimization process. These parameters include threshold type (soft, hard, universal or other), threshold multiplier and mother wavelet function. Even though the threshold type can be selected through trial and error, no particular guidelines exist for wavelet mother function and an appropriate value of threshold multiplier.
In the current study, we attempt to present a solution to the selection of the threshold multiplier value. At each level of wavelet decomposition, a threshold T calculated for each value of T M is applied to wavelet transform coefficients to obtain a reconstructed denoised signal x d . Then, a mean-squared error between the original signal and a reconstructed denoised signal is computed at each value of T M . Finally, the SNR is computed as
S N R ( T M ) = 10 × log 10 | | x d ( T M ) | | 2 | | x ( T M ) x d ( T M ) | | 2
where | | . | | is the norm. In this study, 100 linearly spaced threshold multiplier values ranging from 0.5 to 10 were tested to estimate an appropriate threshold multiplier. As a threshold multiplier increases, the threshold level increases as well, meaning more aggressive denoising as more wavelet decomposition coefficients are set to zero. As a consequence, the mean squared error (MSE) increases and the SNR decreases. This effect can be seen in Figure 17, where both of these quantities change monotonically. A breakpoint of T M is chosen where the MSE and SNR curves intersect. With this approach, SNR is equal to 33.4 dB and T M to 2.13. The disadvantage of this approach is that it is empirical as no strict guidelines exist on how to adequately choose T M .

6. Conclusions

The main goal of the study is the development of signal quality estimation technique based on continuous wavelet transform (CWT). This method requires identification of wavelet ridges, weighing their energy by a duration normalized to the duration of the whole signal and, lastly, performing median filtering of the weighted energies. Optimal median filtering is assured by selecting appropriate window length. The method is validated on nine test subjects, with different body mass indices, ages and levels of physical activity. The bioimpedance signals are measured with a wearable device fixed on a wrist. In addition, plethysmography (PPG) signals are measured from an index finger. We assess the effect of two different electrode couplants: gel and hydrogel.
In addition, this work also showcases a custom wearable device for human bio-impedance measurements. It has proven to be reliable according to the inspection of the measured signals. Nevertheless, it was necessary to carry out band-pass filtering to remove a negative slope and reveal the details of signal morphology.
The following findings must be highlighted:
  • Visualization of signal components in time-frequency plane aids in signal quality estimation. The most intensive signal components corresponding to heartbeat are traced through the whole duration of the signal. Other components, including noise, can also be traced in terms of their amplitude and duration as indicated by wavelet scalograms.
  • PPG signals are helpful to be correlated with bioimpedance signals to match the heartbeat component. PPG scalograms show that apart from the heartbeat component, signals are relatively clean. This means that the other components as seen in bioimpedance scalograms likely correspond to noise.
  • The measured bioimpedance signals are non-stationary and are not normally distributed, meaning that many existing methods for signal quality estimation such as autocorrelation are not applicable. The proposed method, on the other hand, is based on wavelet transform of the measured signal and, hence, can handle signal non-stationarity and non-normality.
  • The developed signal quality estimation method is fully adaptive. Compared to the classical wavelet-based signal denoising and signal-to-noise ratio estimation, there is no need to set fixed threshold levels or select proper threshold types. Instead, median filtering with an optimized window length is performed and noise is estimated from the available ridge information. SNR estimates obtained with the classical and the proposed methods are within error bounds.
  • The effect of coupling agent does not have a significant impact on signal quality as assessed by the SNR. All the differences are within the error bounds. Also, there is no evidence that physical activity levels, age and BMI of test subjects have a role in signal quality of the bioimpedance signals. Again, all differences are within error bounds.
The limitations of the applicability of the proposed SNR indicator include relatively high computational requirements. First, the signal needs to be pre-processed, as the low frequency slope had to be removed, since otherwise the energy of the slope’s frequency component might have been taken into account in calculation of the SNR. Secondly, computation of continuous wavelet transform with a desired resolution itself requires significant computational resources, especially for the median filtering via a cross-validation routine.

Author Contributions

Conceptualization, D.L., R.J. and A.E.; methodology, D.L., R.J. and A.E.; software, D.L. and R.J.; validation, D.L. and R.J.; investigation, D.L. and R.J.; writing—original draft preparation, D.L., R.J. and A.E.; writing—review and editing, D.L., R.J. and A.E.; visualization, D.L. and R.J.; supervision, A.E.; project administration, A.E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Latvian Council of Science project: “Body-Coupled Communication for Body Area Networks”, No: lzp-2020/1-0358. The wearable device was originally designed in the ERDF Activity 1.1.1.2 “Post-doctoral Research Aid” (No. 1.1.1.2/VIAA/2/18/282).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available upon request.

Acknowledgments

We thank Stefans Romanoss for help with the embedded software, and all the volunteers of the study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Mansouri, S.; Chabchoub, S.; Alharbi, Y.; Alshrouf, A.; Nebhen, J. A Real-Time Heart Rate Detection Algorithm Based on Peripheral Electrical Bioimpedance. IEEJ Trans. Electr. Electron. Eng. 2022, 17, 1054–1060. [Google Scholar] [CrossRef]
  2. Lapsa, D.; Janeliukštis, R.; Elsts, A. Electrode Comparison for Heart Rate Detection via Bioimpedance Measurements. In Proceedings of the 2022 Workshop on Microwave Theory and Techniques in Wireless Communications (MTTW), Riga, Latvia, 5–7 October 2022; IEEE: Piscataway, NJ, USA, 2022; pp. 41–46. [Google Scholar]
  3. Lenis, G.; Pilia, N.; Loewe, A.; Schulze, W.H.; Dössel, O. Comparison of baseline wander removal techniques considering the preservation of ST changes in the ischemic ECG: A simulation study. Comput. Math. Methods Med. 2017, 2017, 9295029. [Google Scholar] [CrossRef] [PubMed]
  4. Asgari, S.; Mehrnia, A. A novel low-complexity digital filter design for wearable ECG devices. PLoS ONE 2017, 12, e0175139. [Google Scholar] [CrossRef] [PubMed]
  5. Jin, Z.; Dong, A.; Shu, M.; Wang, Y. Sparse ECG denoising with generalized minimax concave penalty. Sensors 2019, 19, 1718. [Google Scholar] [CrossRef]
  6. Chen, B.; Li, Y.; Cao, X.; Sun, W.; He, W. Removal of power line interference from ECG signals using adaptive notch filters of sharp resolution. IEEE Access 2019, 7, 150667–150676. [Google Scholar] [CrossRef]
  7. Appathurai, A.; Carol, J.J.; Raja, C.; Kumar, S.; Daniel, A.V.; Malar, A.J.G.; Fred, A.L.; Krishnamoorthy, S. A study on ECG signal characterization and practical implementation of some ECG characterization techniques. Measurement 2019, 147, 106384. [Google Scholar] [CrossRef]
  8. Poungponsri, S.; Yu, X. An adaptive filtering approach for electrocardiogram (ECG) signal noise reduction using neural networks. Neurocomputing 2013, 117, 206–213. [Google Scholar] [CrossRef]
  9. Mohebbian, M.R.; Alam, M.W.; Wahid, K.A.; Dinh, A. Single channel high noise level ECG deconvolution using optimized blind adaptive filtering and fixed-point convolution kernel compensation. Biomed. Signal Process. Control 2020, 57, 101673. [Google Scholar] [CrossRef]
  10. Rakshit, M.; Das, S. An efficient ECG denoising methodology using empirical mode decomposition and adaptive switching mean filter. Biomed. Signal Process. Control 2018, 40, 140–148. [Google Scholar] [CrossRef]
  11. Kumar, S.; Panigrahy, D.; Sahu, P. Denoising of electrocardiogram (ECG) signal by using empirical mode decomposition (EMD) with non-local mean (NLM) technique. Biocybern. Biomed. Eng. 2018, 38, 297–312. [Google Scholar] [CrossRef]
  12. Vázquez, R.; Velez-Perez, H.; Ranta, R.; Dorr, V.; Maquin, D.; Maillard, L. Blind source separation wavelet denoising and discriminant analysis for EEG artefacts and noise cancelling. Biomed. Signal Process. Control 2012, 7, 389–400. [Google Scholar] [CrossRef]
  13. Singhal, A.; Singh, P.; Fatimah, B.; Pachori, R. An efficient removal of power-line interference and baseline wander from ECG signals by employing Fourier decomposition technique. Biomed. Signal Process. Control 2020, 57, 101741. [Google Scholar] [CrossRef]
  14. Kumar, A.; Tomar, H.; Mehla, V.K.; Komaragiri, R.; Kumar, M. Stationary wavelet transform based ECG signal denoising method. ISA Trans. 2021, 114, 251–262. [Google Scholar] [CrossRef] [PubMed]
  15. Kumar, A.; Komaragiri, R.; Kumar, M. Heart rate monitoring and therapeutic devices: A wavelet transform based approach for the modeling and classification of congestive heart failure. ISA Trans. 2018, 79, 239–250. [Google Scholar] [CrossRef]
  16. Kumar, A.; Kumar, M.; Komaragiri, R. Design of a biorthogonal wavelet transform based R-peak detection and data compression scheme for implantable cardiac pacemaker systems. J. Med. Syst. 2018, 42, 102. [Google Scholar] [CrossRef]
  17. Berwal, D.; Vandana, C.; Dewan, S.; Jiji, C.; Baghini, M.S. Motion artifact removal in ambulatory ECG signal for heart rate variability analysis. IEEE Sens. J. 2019, 19, 12432–12442. [Google Scholar] [CrossRef]
  18. Li, C.; Deng, H.; Yin, S.; Wang, C.; Zhu, Y. sEMG signal filtering study using synchrosqueezing wavelet transform with differential evolution optimized threshold. Results Eng. 2023, 18, 101150. [Google Scholar] [CrossRef]
  19. Donoho, D.; Johnstone, I. Adapting to unknown smoothness via wavelet shrinkage. J. Am. Stat. Assoc. 1995, 90, 1200–1224. [Google Scholar] [CrossRef]
  20. Sai, C.Y.; Mokhtar, N.; Arof, H.; Cumming, P.; Iwahashi, M. Automated classification and removal of EEG artifacts with SVM and wavelet-ICA. IEEE J. Biomed. Health Inform. 2018, 22, 664–670. [Google Scholar] [CrossRef]
  21. Singh, P.; Pradhan, G.; Shahnawazuddin, S. Denoising of ECG signal by non-local estimation of approximation coefficients in DWT. Biocybern. Biomed. Eng. 2017, 37, 599–610. [Google Scholar] [CrossRef]
  22. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  23. Elsts, A.; McConville, R. Are microcontrollers ready for deep learning-based human activity recognition? Electronics 2021, 10, 2640. [Google Scholar] [CrossRef]
  24. Critcher, S.; Freeborn, T.J. Localized bioimpedance measurements with the max3000x integrated circuit: Characterization and demonstration. Sensors 2021, 21, 3013. [Google Scholar] [CrossRef] [PubMed]
  25. Grimnes, S.; Martinsen, Ø.G. Sources of error in tetrapolar impedance measurements on biomaterials and other ionic conductors. J. Phys. Appl. Phys. 2006, 40, 9. [Google Scholar] [CrossRef]
  26. Yılmaz, K.; Adeyemi, A.; Antink, C.H.; Vehkaoja, A. Comparison of Electrode Configurations for Impedance Plethysmography Based Heart Rate Estimation at the Forearm. In Proceedings of the 2022 IEEE Sensors, Dallas, TX, USA, 30 October–2 November 2022; IEEE: Piscataway, NJ, USA, 2022; pp. 1–4. [Google Scholar]
  27. TENS GEL. Available online: https://en.morettispa.com/prodotto/gel-per-e-c-g-e-tens-260-g/ (accessed on 1 October 2023).
  28. Hydrogel Pads. Available online: https://www.amazon.com/UYGHHK-Trainer-Replacement-Stimulator-Training/dp/B07MK6K8HK (accessed on 1 October 2023).
  29. Gajewski, P.; Żyła, W.; Kazimierczak, K.; Marcinkowska, A. Hydrogel Polymer Electrolytes: Synthesis, Physicochemical Characterization and Application in Electrochemical Capacitors. Gels 2023, 9, 527. [Google Scholar] [CrossRef]
  30. Ishai, P.B.; Talary, M.S.; Caduff, A.; Levy, E.; Feldman, Y. Electrode polarization in dielectric measurements: A review. Meas. Sci. Technol. 2013, 24, 102001. [Google Scholar] [CrossRef]
  31. Shin, H.S.; Lee, C.; Lee, M. Adaptive threshold method for the peak detection of photoplethysmographic waveform. Comput. Biol. Med. 2009, 39, 1145–1152. [Google Scholar] [CrossRef]
  32. Jia, Y.; Pengfei, X. Noise cancellation in vibration signals using an oversampling and two-stage autocorrelation model. Results Eng. 2020, 6, 100136. [Google Scholar] [CrossRef]
Figure 1. The wearable device used in the experiments, fully assembled view.
Figure 1. The wearable device used in the experiments, fully assembled view.
Sensors 23 08532 g001
Figure 2. Signal quality estimation in bioimpedance-based human heart rate monitoring applications.
Figure 2. Signal quality estimation in bioimpedance-based human heart rate monitoring applications.
Sensors 23 08532 g002
Figure 3. k-fold cross-validation scheme with 10 folds to optimize window length selection for median filtering.
Figure 3. k-fold cross-validation scheme with 10 folds to optimize window length selection for median filtering.
Sensors 23 08532 g003
Figure 4. Bioimpedance wearable device equivalent circuit connected to a computer via Bluetooth.
Figure 4. Bioimpedance wearable device equivalent circuit connected to a computer via Bluetooth.
Sensors 23 08532 g004
Figure 5. Raw bioimpedance measurements for subject I. Measurements no. 1, 5 and 10 are shown.
Figure 5. Raw bioimpedance measurements for subject I. Measurements no. 1, 5 and 10 are shown.
Sensors 23 08532 g005
Figure 6. Band-pass filtered bioimpedance signals for subject I.
Figure 6. Band-pass filtered bioimpedance signals for subject I.
Sensors 23 08532 g006
Figure 7. Zoomed-in on the first 10 s of band-pass filtered bioimpedance signals for subject I.
Figure 7. Zoomed-in on the first 10 s of band-pass filtered bioimpedance signals for subject I.
Sensors 23 08532 g007
Figure 8. Raw PPG measurements for subject I. Measurements no. 1, 5 and 10 are shown.
Figure 8. Raw PPG measurements for subject I. Measurements no. 1, 5 and 10 are shown.
Sensors 23 08532 g008
Figure 9. Band-pass filtered PPG signals for subject I. These signals correspond to the presented filtered bioimpedance signals.
Figure 9. Band-pass filtered PPG signals for subject I. These signals correspond to the presented filtered bioimpedance signals.
Sensors 23 08532 g009
Figure 10. Zoom-in on the first 10 s of band-pass filtered PPG signals for subject I.
Figure 10. Zoom-in on the first 10 s of band-pass filtered PPG signals for subject I.
Sensors 23 08532 g010
Figure 11. Bioimpedance scalograms for subject I. High amplitude can be traced around 1 Hz. Lower amplitude and lower duration signal components are present.
Figure 11. Bioimpedance scalograms for subject I. High amplitude can be traced around 1 Hz. Lower amplitude and lower duration signal components are present.
Sensors 23 08532 g011
Figure 12. PPG scalograms for subject I. Dominating signal component is around 1 Hz.
Figure 12. PPG scalograms for subject I. Dominating signal component is around 1 Hz.
Sensors 23 08532 g012
Figure 13. Values of the optimized window length values as obtained in 1000 runs of 10-fold cross-validation for subject I, trial 1. Each * denotes the result of an individual run. High coefficient of variation suggests that selection of a particular value is not critical.
Figure 13. Values of the optimized window length values as obtained in 1000 runs of 10-fold cross-validation for subject I, trial 1. Each * denotes the result of an individual run. High coefficient of variation suggests that selection of a particular value is not critical.
Sensors 23 08532 g013
Figure 14. Normalized durations and energies of the ridges identified from the post-processed bioimpedance signals for subject I. Ridge durations normalized to a measurement length are shown as blue bars corresponding to the left y axis. Ridge energies normalized to the total signal energy are shown as red bars.
Figure 14. Normalized durations and energies of the ridges identified from the post-processed bioimpedance signals for subject I. Ridge durations normalized to a measurement length are shown as blue bars corresponding to the left y axis. Ridge energies normalized to the total signal energy are shown as red bars.
Sensors 23 08532 g014
Figure 15. SNR values versus window length for the 1st, 5th and 10th trial measurements performed on subject I. The maximum SNR is consistently achieved at window length 3 both, for gel and hydrogel couplants.
Figure 15. SNR values versus window length for the 1st, 5th and 10th trial measurements performed on subject I. The maximum SNR is consistently achieved at window length 3 both, for gel and hydrogel couplants.
Sensors 23 08532 g015
Figure 16. Mean SNR values for gel and hydrogel couplants for all test subjects. Error bars show confidence intervals.
Figure 16. Mean SNR values for gel and hydrogel couplants for all test subjects. Error bars show confidence intervals.
Sensors 23 08532 g016
Figure 17. Estimation of threshold multiplier using classical wavelet transform method (subject I, trial 8, gel). Threshold multiplier is selected at the crossing of MSE and SNR curves. With this approach the SNR is equal to 33.4 dB.
Figure 17. Estimation of threshold multiplier using classical wavelet transform method (subject I, trial 8, gel). Threshold multiplier is selected at the crossing of MSE and SNR curves. With this approach the SNR is equal to 33.4 dB.
Sensors 23 08532 g017
Table 1. Qualitative comparison of different MCU and SoC. nRF5340 (bolded in the table) is selected for the wearable device.
Table 1. Qualitative comparison of different MCU and SoC. nRF5340 (bolded in the table) is selected for the wearable device.
FamilyNameMax Clock FrequencyArchitectureRadio
nRF52nRF5284064-MHzCortex M4FBLE 5, IEEE 802.15.4, Propr. 2.4 GHz
nRF52nRF5283264-MHzCortex M4FBLE 5.1, Propr. 2.4 GHz
nRF53nRF5340128-MHzCortex M33BLE 5.2, IEEE 802.15.4, Propr. 2.4 GHz
TI SimpleLinkCC2652R48-MHzCortex M4FBLE 5.1, IEEE 802.15.4
MSP432msp432p401r48-MHzCortex M4FNo Built-in
STM32LSTM32L47548-MHzCortex M4FNo Built-in
ESP32ESP32240-MHzXtensa LX6BLE 5.0, WIFI, IEEE 802.11 b/g/n,
ESP32-C3ESP32-C3160-MHzRISC-VBLE 5.0, WIFI, IEEE 802.11 b/g/n,
EFR32EFR32BG2276.8-MHzCortex M33BLE 5.2, 802.15.4 g
SAM3XATSAM3X8E84-MHzCortex M3No Built-in
Table 2. Performance comparison of different MCU and SoC. nRF5340 (bolded in the table) is selected for the wearable device.
Table 2. Performance comparison of different MCU and SoC. nRF5340 (bolded in the table) is selected for the wearable device.
NameActive Current/Radio Current @0 dBAvg. mARAMFLASH
mAµA/MHzRx (mA)Tx (mA)RAM ret. (µA)@ 48 MHzkBkB
nRF528403.328524.64.83.1653.12561024
nRF528323.7585.45.31.957.664/32512/256
nRF53407.3573.84.22.357512 + 641024 + 256
CC2652R3.470.86.97.30.946980352
msp432p401r3.8480--0.35-64256
STM32L4758166.7--0.236-1281000
ESP3268283.310013010282.9400384
ESP32-C320125--5-400384
EFR32BG222.07273.64.11.427.332512
ATSAM3X8E70.89923--2.5-96512
Table 3. Biological data of test subjects.
Table 3. Biological data of test subjects.
SubjectGenderAge
(Years)
Mass
(kg)
Height
(m)
BMIActivity Level
Imale34761.8023.46high
IImale54941.8427.76low
IIImale22571.7817.99medium
IVmale30751.9519.72medium
Vmale39701.8320.90medium
VImale26831.8324.78medium
VIIfemale23701.7522.86low
VIIImale30881.7528.73high
IXmale241052.0325.48low
Table 4. Normality test of bioimpedance data for subject I at α = 0.05 .
Table 4. Normality test of bioimpedance data for subject I at α = 0.05 .
Kolmogorov–Smirnov TestAnderson–Darling TestJarque–Bera Test
p-ValueKSKS Crit.p-ValueADAD Crit.p-ValueJBJB Crit.
Gel #100.270.01970.00053.590.750.00128.265.98
Gel #200.320.01950.00056.370.750.00131.095.98
Gel #300.280.01970.00055.970.750.00137.295.98
Gel #400.250.01960.00056.100.750.00155.095.98
Gel #500.290.01970.00052.310.750.00120.115.98
Gel #600.260.01960.00053.900.750.00120.005.98
Gel #700.300.01930.00051.800.750.00115.055.98
Gel #800.300.02030.00055.210.750.00168.055.98
Gel #900.270.01980.00052.700.750.00121.725.98
Gel #1000.260.01980.00441.180.750.00610.715.98
H.gel #100.360.01680.00052.290.750.001240.365.98
H.gel #200.310.01910.00051.870.750.00610.415.98
H.gel #300.270.01930.00055.520.750.00312.575.98
H.gel #400.310.01960.00056.430.750.00128.125.98
H.gel #500.260.01950.00054.170.750.00137.505.98
H.gel #600.250.01960.00057.050.750.00155.155.98
H.gel #700.260.01950.000515.190.750.00189.205.98
H.gel #800.270.01940.00054.720.750.00131.165.98
H.gel #900.260.01990.00057.070.750.00141.095.98
H.gel #1000.270.01950.000512.840.750.00181.955.98
Table 5. Stationarity test of bioimpedance data for subject I at α = 0.05 .
Table 5. Stationarity test of bioimpedance data for subject I at α = 0.05 .
Augmented Dickey–Fuller TestVariance Ratio Test
p-ValueADFADF Crit.p-ValueVRVR Crit.
Gel #10.001−6.79−1.94035.141.96
Gel #20.001−7.51−1.94034.971.96
Gel #30.001−6.88−1.94035.261.96
Gel #40.001−7.43−1.94036.391.96
Gel #50.001−7.06−1.94033.301.96
Gel #60.001−7.26−1.94034.651.96
Gel #70.001−7.53−1.94036.121.96
Gel #80.001−7.43−1.94034.971.96
Gel #90.001−7.23−1.94033.791.96
Gel #100.001−7.68−1.94034.221.96
Hydrogel #10.001−7.76−1.94040.811.96
Hydrogel #20.001−6.55−1.94035.911.96
Hydrogel #30.001−6.44−1.94033.951.96
Hydrogel #40.001−6.11−1.94034.351.96
Hydrogel #50.001−6.92−1.94034.741.96
Hydrogel #60.001−6.18−1.94035.131.96
Hydrogel #70.001−6.86−1.94036.191.96
Hydrogel #80.001−6.76−1.94033.421.96
Hydrogel #90.001−5.97−1.94034.171.96
Hydrogel #100.001−6.34−1.94034.921.96
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lapsa, D.; Janeliukstis, R.; Elsts, A. Adaptive Signal-to-Noise Ratio Indicator for Wearable Bioimpedance Monitoring. Sensors 2023, 23, 8532. https://doi.org/10.3390/s23208532

AMA Style

Lapsa D, Janeliukstis R, Elsts A. Adaptive Signal-to-Noise Ratio Indicator for Wearable Bioimpedance Monitoring. Sensors. 2023; 23(20):8532. https://doi.org/10.3390/s23208532

Chicago/Turabian Style

Lapsa, Didzis, Rims Janeliukstis, and Atis Elsts. 2023. "Adaptive Signal-to-Noise Ratio Indicator for Wearable Bioimpedance Monitoring" Sensors 23, no. 20: 8532. https://doi.org/10.3390/s23208532

APA Style

Lapsa, D., Janeliukstis, R., & Elsts, A. (2023). Adaptive Signal-to-Noise Ratio Indicator for Wearable Bioimpedance Monitoring. Sensors, 23(20), 8532. https://doi.org/10.3390/s23208532

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop