1. Introduction
Physiological signal capturing devices have evolved rapidly in recent years. From being a curiosity a decade ago, today, almost any electronic devices (mainly cell phones, smart watches or wearable devices in general) have several electronic sensors capable of recording the individual’s physical activity and registering physiological signals such as body temperature or heart rate. The capture of physiological signals itself has been the subject of extensive study in recent years [
1], and the variety and importance of its applications has grown steadily throughout this time, especially after the COVID-19 pandemic [
2].
The advent of all kind of non-invasive sensors has enabled the development of many platforms for acquiring an increasingly wide range of physiological signals. Examples of commonly measured physiological signals are electrocardiogram (ECG) [
3,
4,
5,
6,
7], electromyogram (EMG) [
8,
9,
10,
11], photoplethysmography (PPG) [
12,
13], including heart rate (HR) and its variability (HRV) [
14,
15], body or skin temperature (ST) [
16,
17], skin conductance and impedance [
18,
19,
20], electrodermal activity (EDA) [
21,
22,
23], oxygen saturation (SpO2) [
24,
25,
26], blood pressure (BP) [
27], and respiration rate (RR) [
28], among others.
1.1. Signal Capturing and Wearables
Wearable devices are intelligent electronic platforms, incorporated into clothing or worn on the body, generally in the form of accessories (forehead bands, bracelets, patches, watches), which allow the collection of physiological data from the user and, in general, a more or less detailed analysis of these data, with fewer restrictions compared to traditional laboratory devices [
29]. Today, the use of wearable devices has become widespread, as have the signals that these devices are capable of recording and analyzing.
The first wearable device could be the pocket watch (16th century), which evolved into the wristwatch (1904) and the first digital watches (1969). Today, in the space of a wristwatch, almost all the capabilities of a cell phone can be integrated, including the ability to send and receive calls and messages and the use of apps, as well as different types of sensors that record everything from the subject’s physical activity to some of their basic physiological signals. In addition to the wristwatch [
30], other wearable devices have taken the form of rings [
17], earrings [
31], forehead bands [
32] and different clothing such as gloves, shoes and others [
33].
Some examples of successful wearable devices that are capable of storing, analyzing and sending physiological signals are the Samsung Gear S3 Frontier watch (
https://www.samsung.com/es/wearables/gear-s3/highlights/ Accessed on 1 August 2022) (capable of measuring position and motion and HR, and featuring barometer, GPS and light sensor), the MUSE and MUSE 2 forehead bands (
https://choosemuse.com/ Accessed on 1 August 2022) (which record and analyze EEG signal, as well as being able to measure HR, SpO2, motion and position) and the Oura rings (
https://ouraring.com/ Accessed on 1 August 2022) (with heart rate monitoring, activity detection and seven temperature sensors). For reasons that will become evident later, it is necessary to mention one of the most successful wearable devices on the market today, the Empatica E4 wristband (
https://www.empatica.com/en-eu/research/e4/ Accessed on 1 August 2022), which has four sensors: PPG sensor, EDA sensor, a three-axis accelerometer and an optical thermometer for ST measurement. This device has several working modes including data recording and transfer, streaming mode for real-time data visualization and E4 connect mode, which visualizes and manages physiological data (encrypted) in a cell phone app.
The development of platforms for recording, monitoring and analyzing physiological signals has continued growing for more than a decade [
1,
34] and is now a very active branch of sensor electronics [
35,
36,
37,
38]. The range of uses of captured physiological signals has expanded proportionally to the spread of this technology, from the extraction of health-related characteristics to the classification of emotions or moods. More specific applications of physiological signals may include, among others, disease diagnosis [
25,
39], remote healthcare [
40,
41,
42,
43,
44,
45,
46], human–machine interaction (e.g., to increase the autonomy of people with disabilities) [
47,
48], rehabilitation processes [
49], stress detection [
34,
50,
51], monitoring of daily activity [
52,
53] and emotion detection [
32,
38,
54,
55,
56,
57].
This paper will detail the development of a platform for the acquisition of physiological signals, which includes a total of five different biomedical signals: skin temperature (ST), blood oxygen saturation (SpO2), heart rate monitoring (HR), movement/position of user’s head (AHRS) and electrodermal activity (EDA). The sensors used for the different measurements have been chosen among the most advanced ones available on the market, trying to prioritize the accuracy of the measurements as well as their performance. For example, the EDA sensor can also measure bioimpedance and allows different types of measurement configuration, including AC and DC (unipolar and switched bipolar), which improves the signal-to-noise ratio. As will be explained later, the accuracy of the system has been improved by a measurement based on two-channel multiplexing. A complete AHRS system has been implemented, including a three-axis accelerometer, gyroscope and magnetometer. The sensor for measuring SpO2 and HR is also state-of-the-art, allowing high-precision measurements. In addition, a fisheye camera is incorporated for the study of facial expression. Although the present work focuses on the electronic part of the design, the main objective is to use the Multi-sensory Wearable Headband (MsWH) platform for emotion recognition. Indeed, in emotion recognition, valence and arousal (among other emotional parameters) can be obtained through the selected physiological signals [
58,
59] and facial expression [
55,
60,
61,
62]. However, this device also has the potential to be used to detect and discriminate acute psychological stress and the performance of different physical activities occurring at different times or concurrently [
50,
63].
1.2. Aims of This Paper
As mentioned above, this paper will detail the development of the Multi-sensory Wearable Headband (MsWH) platform, which is capable of recording and analyzing five different physiological signals (ST, SpO2, HR, IMU and EDA), in addition to capturing the user’s facial expression using a fisheye camera. The MsHW platform integrates sensors already calibrated by the manufacturers (ST, HR, SpO2 and IMU). For this reason, aspects such as device consumption, EDA signal reading and response time of the temperature sensors are dealt with in greater depth in this work. In these measurements, the tests have been performed with accurate measurement systems on test bench. On the other hand, as a validation of results, the available compatible signals such as EDA, ST and HR have been compared against the Empatica E4 wristband data.
In the following section (
Section 1.1), the technological differences of the developed platform compared to similar devices will be detailed. The present work focuses basically on the electronics required to capture the different signals, first treated in a descriptive way (
Section 1.1). The necessary power supplies, the conditioning stages and the connection to and between the different integrated circuits capable of recording and analyzing (totally or partially) the raw signals are described and detailed in
Section 3. In
Section 4, there is a brief description of the software required to receive and interpret these signals. The results obtained in the capture of these signals are shown in
Section 5, in which the data obtained by the device presented are compared with those obtained through the Empatica E4 wristband when possible. Finally, some conclusions and the future work are presented in
Section 6.
1.3. Related Works
This section will briefly detail some works on devices similar in one or more aspects to the one detailed here in order to highlight the differences of the proposed platform with respect to other pre-existing ones.
There is a wide variability of integrated circuits capable of performing accurate ST measurement. In [
32], the authors use the LMT70 integrated contact temperature sensor, which is suitable for portable devices and has an error as low as 0.05 °C. In [
23,
37], the MLX90614 chip is used. This integrated sensor uses infrared technology, presents an error in the ST measurement of 0.5 °C, communicates by SMBus and has a library for use with Arduino. As will be seen in
Section 3.1, the MAX30205 has been used in the MsHW platform, with an accuracy of 0.1 °C in the range of 37 to 39 °C and sends data via an I2C-compatible interface.
As will be seen in
Section 3.2, for HR and SpO2 measurement, we have used the MAX30101 integrated unit, which is a similar but more modern version than, for example, the MAX30100 integrated unit used in [
23].
The most significant advantage in terms of user motion measurement in the present work (
Section 3.3) is that MsHW is provided with a complete AHRS system. Indeed, the BNO055 has a three-axis accelerometer, gyroscope and magnetometer, unlike other [
55] platforms, where the latter does not exist. Thus, a complete inertial system is available that is capable of recording relative position and movements with high accuracy.
The EDA measure is currently one of the most widespread and active in research. Most devices present a solution for EDA measurement based on conventional analog circuitry (in the sense of circuitry built basically from amplifiers). For example, in [
31] a three-level amplifier stage (3CA) has been designed: a first unity gain stage with an AD8244, a preamplifier stage with an AD8222 and a 24-bit low-noise ADC (ADS1299), plus a final adaptive amplifier stage with an MSP432. In [
32], a direct impedance measurement is performed with a circuit where 2 low-power MCP6422 amplifiers are used. In [
22], a single AD8603 is used, while a solution more similar to the one presented here is given in [
37], with a double impedance measurement system through two MPC6004, plus a second stage that adjusts the output voltage to that of the ADC range used for signal processing.
In this work, a more compact solution has been chosen. The AD5941 is a measurement system that can be used, among other applications, for electrochemical measurements, battery impedance, glucose monitoring and bioimpedance measurements. It includes the analog inputs for the corresponding sensors, the conditioning of the input signals (amplifying, filtering), the generator of the excitation signals, and the converters and the microcontroller required for the calculation of the variable to be measured.
In conclusion, in the MsHW platform, we have tried to use state-of-the-art technologies in the precision measurement of the corresponding signals, with the idea of having highly accurate data that will be used later for emotion recognition. In this sense, it is known [
64] that multimodal emotional detection systems tend to be more accurate than unimodal ones. As concluded in [
65], the combination of methods for recording and analyzing physiological variables (especially if it is a combination of several signals) with deep learning algorithms can be an extremely powerful tool in the detection of emotional states in any of its applications.
2. Platform Overview
As mentioned in the previous section, the present work will detail the development of a platform that includes a total of five different physiological signals: skin temperature (ST), blood oxygen saturation (SpO2), heart rate (HR) monitoring, user movement/position (IMU) and electrodermal activity/bioimpedance (EDA). In addition, a fisheye camera is incorporated for facial recognition and facial expression study. The device is attached to the subject’s head by means of an elastic band and to the earlobe by means of a clip.In this sense, the developed platform is “wearable”, but it has been developed for data collection in the laboratory or under controlled conditions, and it is not for everyday use in free living or physical activity. On the forehead, on the same bracket that holds the camera, are located the power and communications control system (Wi-Fi, Bluetooth, and a serial-USB port), the main microcontroller, the chip for EDA, the gyroscope and accelerometers (IMU), and the LIPO battery that powers all the electronics. The platform is completed with a second board (oximeter) that is placed in the earlobe. The information from all these sensors can be preprocessed inside the device itself (ST, EDA, IMU, HR, SpO2) or the raw data can be sent to a microcontroller capable of analyzing it and extracting high-level characteristics.
The complete system can be seen in
Figure 1. The developed platform is divided into two main subsystems or PCBs.
Figure 1a shows a side view of the platform, which is attached to the subject’s forehead by means of an elastic strap. The image also reflects the approximate view area of the camera.
Figure 1b shows the electrodes of the EDA/bioimpedance subsystem. These electrodes (the silver brackets attached to the forehead and perfectly static by means of a strip of ethyl–vinyl acetate or EVA foam). In green, the temperature sensor can be seen.
Figure 1c shows an identical view to the previous one with the electrodes placed in their final position and fixed by means of the EVA foam. This way, they will be in contact with the subject’s forehead and always separated at the same distance. The second PCB, attached to the ear by means of a clamp, is shown in
Figure 1d. This PCB contains the SpO2 and HR measure systems, as well as the fastening system (in pink). In these images, we can be mainly see a (1) microcontroller, (2) camera, (3) LIPO battery, (4) EDA electrode, (5) pulse oximeter, (6) temperature sensor.
Each of these subsystems has its own power supply as well as the signal conditioning stages necessary for its correct operation. In addition, both modules have a microprocessor (ESP32) capable of analyzing in situ the signals from the different sensors and extracting the most notable characteristics from them, although if necessary, both devices can, in turn, dump the data in raw format to an external computer for further processing via Bluetooth and Wi-Fi connection.
The total weight of the platform is 168 grams, 96 grams of which correspond to the head attachment system. The PCBs of both subsystems can be seen in
Figure 2. On the left side,
Figure 2a,b show the top and bottom views of the PCB containing the RS232 system, the camera connection port (bottom view), the main microcontroller (wich includes Wi-Fi), the communication and battery control, the switched-mode power supplies that generate the supply voltages for the different sensors and their conditioning, and a USB port (on the right side of the board, top view). The size of this PCB is 38 mm × 124 mm.
Figure 2b,d, on the right, show the top and bottom faces of the PCB corresponding to the pulse oximeter, which is attached to the subject’s earlobe by means of a clamp, as shown in
Figure 1d. The size of this PCB is 18 mm × 20 mm. There is a third very simple PCB, which contains only the ST sensor and must be separated from the others since the sensor used needs to be in thermal contact with the subject’s forehead. This small PCB can be partially seen in
Figure 1b,c, number (6).
The following section will detail the measurement devices for each of the physiological signals implemented in the platform that has been developed.
3. Hardware Description
In this section, we detail the proposed platform design in terms of hardware and signal acquisition. The block diagram of the platform is shown in
Figure 3.
As can be seen in this figure, the platform is powered by lithium–polymer batteries with output voltages of 5.0 V and 3.3 V for the different devices. The EDA sensor, the IMU and digital camera are located inside the forehead board. The temperature sensor and pulse oximeter are located in the earlobe subcircuit. These two devices and the IMU communicate with the main microcontroller via I2C, while the EDA sensor sends its measurements via SPI and the camera uses the parallel port. The main microcontroller can, in turn, communicate with external devices via Wi-Fi or Bluetooth. The approximate battery life, with all sensors and camera at full power, is about 2 h. In this time, the amount of data collected, analyzed and sent can be massive.
Some characteristics of the different sensors and devices included in the platform are shown in
Table 1. Each subsystem will be detailed in the following sections.
3.1. Skin Temperature (ST)
Attached to the earlobe, the device contains the MAX30205 (
https://datasheets.maximintegrated.com/en/ds/MAX30205.pdf, accessed on 14 June 2022) skin temperature sensor, which, as advanced before, has a 16-bit (0.00390625 °C) temperature resolution, and a final clinical grade accuracy of 0.1 °C from 37 to 39 °C. This device converts analog temperature measurements to the corresponding digital data via a high-resolution sigma-delta analog-to-digital converter (ADC). The measurements obtained can be sent via a lockup-protected I2C-compatible two-wire serial interface. The sensor has a supply voltage range of 2.7 to 3.3 V and a low supply current of 600 µA.
The main difficulty in measuring temperature using a wearable device is making sure that the sensor is measuring the body temperature and not the temperature of the device itself, or whether the measurement is excessively influenced by the ambient temperature. These objectives can be achieved by maximizing skin contact, minimizing the thermal mass of the circuit (including the PCB), and minimizing the thermal exchange between the device and its enclosure with the environment.
Figure 4 depicts the detailed circuit diagram of the MAX30205.
3.2. Blood Oxygen Saturation (SpO2) and Heart Rate (HR) Monitoring
The block diagram of the pulse oximeter and heart rate monitoring are shown in
Figure 5. The monitored data are sent to a MAX32664 (
https://datasheets.maximintegrated.com/en/ds/MAX32664.pdf, accessed on 14 June 2022) micro-controller. It is a sensor hub family with integrated firmware and world-class algorithms for wearables. The MAX32664 seamlessly enables any desired sensor functionality and includes communication with Maxim’s optical sensors. It can deliver raw or already processed data to the outside world through a fast-mode slave I2C interface. A firmware bootloader is also provided.
The MAX32664 firmware version for the developed platform is the so-called Version C. It supports the MAX30101 high-sensitivity pulse oximeter and heart-rate sensor for wearable health. It also supports estimated blood pressure monitoring. The MAX30101 is an integrated module that includes internal LEDs (red and IR for sensing and green for ambient light cancellation), photodetectors, and other optical elements, with low-noise electronics. The MAX30101 provides a complete system solution to ease the design-in process for mobile and portable devices. As shown in
Figure 5, it requires a 1.8 V power supply (provided by a SPX3809), plus a separate 5.0 V power supply for the internal LEDs (obtained using a Buck-Boost REG71050 DC-DC regulator). As previously advanced, communication with the MAX32664 is performed through a standard I2C interface. The module can be shut down by software with zero standby current, which allows energy saving.
The complete circuit diagram of the pulse oximeter is shown in
Figure 6. The connections between the different integrated circuits, as well as other external elements necessary for their correct operation, can be observed.
3.3. Attitude and Heading Reference System (AHRS)
The subject’s orientation, position and movement (Attitude and Heading Reference System, AHRS) are monitored by the Inertial Measuring Unit (IMU) BNO055 (
https://www.bosch-sensortec.com/media/boschsensortec/downloads/datasheets/bst-bno055-ds000.pdf, accessed on 14 June 2022), which is an absolute orientation sensor. It is a System in Package (SiP), integrating a 14-bit three-axis accelerometer, a 16-bit three-axis gyroscope (with a range of ±2000 degrees per second), and a three-axis magnetometer. These three sensors communicate via SPI with a 32-bit Cortex M0+ microcontroller running Bosch Sensortec sensor fusion software. The BNO055 includes internal algorithms to constantly calibrate the gyroscope, accelerometer and magnetometer inside the device. The BNO055 is equipped with bidirectional I2C and UART interfaces. Once programmed to run with the HID-I2C protocol, it turns into a plug-and-play sensor hub solution. The corresponding circuit diagram is shown in
Figure 7.
The final image size is 640 × 480 pixels. In addition, the OV2640 module performs on-chip JPEG encoding, converting the pixels into a machine-readable format without additional computing from the MCU. The camera works at 12.5 frames per second.
3.4. Face and Gesture Recognition
To perform gesture and face recognition, which is a very important characteristic in emotion detection and classification, the OV2640 (
https://www.uctronics.com/download/cam_module/OV2640DS.pdf, accessed on 14 June 2022) camera module was chosen for its high sampling rate and variable pixel resolution. As mentioned above, the camera is fixed on a very light plastic support, which protrudes from the subject’s forehead (see
Figure 1). Thus, the viewing area remains approximately constant, regardless of the user’s movements, which is expected to increase the accuracy of the results of the facial expression recognition algorithms.
3.5. Electro Dermal Activity (EDA) and Bioimpedance
Human skin puts up a certain resistance to the flow of electric current, and therefore, it has an associated electric impedance that can be experimentally measured. The bioimpedance depends on multiple variables both prior to the measurement (such as the separation between the electrodes) and dependent on external stimuli (such as the dermal response to variations in body temperature or perspiration).
3.5.1. Electrodes Selection and Configuration
The bioimpedance measurement process is carried out by means of electrodes (from two to four) that can be placed on different parts of the body, separated by a certain distance. In turn, electrodes can be divided into two broad categories: dry and wet. The measurement procedure is different in each case. To make it user-friendly, dry electrodes have been chosen for the developed platform, while the electrical excitation of the skin is carried out by means of a low-voltage sinusoidal signal at 100 Hz frequency, which the integrated circuit itself is capable of generating. As can be seen in
Figure 1, the electrodes are placed on the forehead, close to the measuring device itself, located in the main PCB. This location of the electrodes has two advantages: firstly, the forehead is one of the areas of greatest electrodermal activity [
66], making measurements in this area particularly sensitive and/or accurate. Secondly, the electrodes remain in a more stable position with respect to, for example, the wrist, thus reducing the risk of electrical noise and unwanted signal variations, which cloud the final quality of the captured data.
In any case, the bioimpedance measurement must meet the IEC 60601 standard, which limits the current that can flow through the human body to a maximum of 10 µA for frequencies below 1 kHz. According to this standard, the sensors must be isolated by a decoupling capacitor
(where
x is the wire number) which guarantees that no DC current penetrates the body. The input current is limited by a resistor
. The connection and the final result of the electrode impedance are shown in
Figure 8.
As can be seen in
Figure 8, there are three main different ways to measure skin bioimpedance, depending on the type of measurement required in each case: through a two-wire, three-wire or four-wire connection. The values of the different decoupling capacitors (
to
) and limiting resistors (
to
) for both connections can be seen in the figure. In all measurement cases, the bioimpedance is essentially located between the electrodes connected to the
and
terminals of the measuring device (see
Figure 9, below), which are decoupled in
Figure 8 by capacitors
and
, respectively. In the four-wire measurement, the remaining electrodes are connected in the proximity of the previous ones, so that the detected impedance is essentially the same.
Although the two-wire and four-wire measurements mentioned above are the most commonly used, in our case, a three-wire connection has been chosen. In this case, a reference electrode is taken, while the skin impedance measurement is obtained between the reference electrode and one of the two additional ones, all of them placed on the individual’s forehead (see
Figure 1). Given the short distance between the electrodes, two similar but slightly different impedances (
and
) will be measured, whose value is multiplexed and analyzed in the microprocessor, seeking to take full advantage of the device’s capabilities (e.g., improving the data accuracy and its signal-to-noise ratio). The measurement process will be detailed in the next section.
3.5.2. Bioimpedance Measurement Process
As mentioned above, the bioimpedance measurement is carried out by the AD5941 integrated circuit. The AD5941 is a high-precision, low-power analog front end, which is specifically designed for portable applications that require high-precision, electrochemical-based measurement techniques, such as current, voltage or impedance measurements. It has been designed for skin and body impedance measurements, and it works as a complete bioelectric or biopotential measurement system.
EDA is measured by voltammetry. To measure an unknown impedance, an excitation signal is applied across this impedance. First, the voltage across the terminals of the unknown impedance is measured. In a second stage, the current flow through the unknown impedance is measured. This current is converted to a voltage via a Trans-Impedance Amplifier (TIA), and this voltage is measured by an Analog to Digital Converter (ADC). In a first stage of analysis, a Discrete Fourier Transform (DFT) is performed on the data measured by the ADC for the current and voltage values, thus obtaining the real and imaginary parts of and . Although there are other methods to perform this calculus, DFT has been chosen since it is implemented in the AD5941 itself, and therefore, no additional operation or programming is required.
As advanced in the previous section, in our case, the bioimpedance measurement is carried out by means of a three-wire connection detailed in
Figure 9. In this case, the electrode connected to input CEO (
) is used as the reference electrode, while a selector (taken from the matrix of selectors of the integrated unit) selects the other point of the skin among the two remaining electrodes (
and
in the figure, which are connected to inputs AIN1 and AIN3, respectively) to obtain two different impedance measurements using the procedure described above.
Under these circumstances, the module of both measured impedances can be obtained through the expression:
The excitation signal,
, is a sinusoidal wave of 1.1 V amplitude, oscillating at a frequency of 100 Hz. The use of an AC signal improves the overall characteristics of the recorded signal compared to that obtained with a DC or a switched DC measurement. Knowing the characteristics of the excitation signal, it is possible to calculate the isolation capacitors necessary for the correct connection of the electrodes (
Figure 8). For this purpose, a limiting resistor per channel is taken as follows:
In order to match the body impedance analysis measurement, the return path capacitor
for EDA must have a value of:
With these values for
and
, the isolation capacitor
is calculated by means of the expression:
where:
is the module of the impedance of capacitors and .
is the RMS value of a sinusoidal signal of 1.1 V amplitude.
is the maximum allowed AC current in the human body (IEC 60601).
is the frequency of the excitation signal.
The complete EDA measurement process can be seen in
Figure 10. As advanced, the excitation signal is a sinusoid of 1.1 V amplitude and 100 Hz frequency. As can be seen in the figure, the complete process takes 250 ms. During the first approximately 100 ms, with the system clock working at 16 MHz, the multiplexed reading of the two channels takes place. The rest of the time the system remains hibernated. Specifically, during the first 50 ms (five full signal periods), the excitation and impedance reading is carried out through the left channel (blue in the figure). Then, by changing the switch mode in the switch matrix of the device (see
Figure 9), the excitation and reading of the right channel (in orange in
Figure 10) is performed. Then, the device hibernates until the next lecture cycle, about 150 ms later. During the hibernation period, the clock works at 32 MHz. The collected data (which include the FFT of the multiplexed read signals) are stored in a FIFO memory, which has a certain size (a user-defined threshold). When the FIFO buffer is full, the system sends a warning signal to the microprocessor (FIFO interrupt). Before the next read cycle, the MCU must find the time to receive the data, which is an event it will mark with a new interrupt. In this way, the FIFO buffer remains empty for a new storage cycle.
Finally, the detailed circuit of connections and external elements connected to the AD5941 can be seen in
Figure 11.
One of the main advantages of this subsystem is the possibility of measuring both EDA and bioimpedance. On the other hand, the multiplexing of the measured impedances, and , leads to an increase in the performance of the device.
4. Software Description
The design of the application can be seen in
Figure 12. It works on a producer–consumer basis, connected by queues that act as temporary buffers for discontinuities, caused in the Wi-Fi interface due to connection quality or performance changes in the data transfer rate.
As can be seen in
Figure 12, the thread producers are camera capture and signal acquisition. Both work asynchronously but with the same time reference (timestamp). The camera capture thread uses the ESP32-camera driver (
https://github.com/espressif/esp32-camera, accessed on 14 June 2022). As advanced in
Section 3.4, the capture rate is 12.5 frames per second in continuous mode for VGA resolution. For each frame, dynamic memory is reserved, and the pointer is queued. The thread consumer (video streaming in
Figure 12) pulls the data out of the queue and sends them over the UDP socket to the target host. Finally, the memory is freed.
The second producer thread is the signal acquisition. This thread synchronizes and acquires signals from all the sensors of the platform. For this purpose, a synchronous acquisition start trigger is performed, and then, the transfers are attended with their different intervals based on the configuration of their internal buffers. In the case of the AD5941, the transfer is controlled by interruption every second with a buffer containing four samples. For the MAX30205 temperature, an acquisition is performed by polling every 250 ms. These two sensors are queued together in one structure, since they are acquired at the same sampling rate. In the same task, also by polling, the internal buffer of the SpO2 implemented in the MAX32664 processor is checked and the data are placed in a specific queue. The same operation is performed with the data obtained from the AHRS, BMO055. For each of the queues, a consumer thread is associated to send the data over UDP sockets on consecutive ports with respect to the base port to the IP address of the destination host.
The processor load is balanced between two cores. On Core 0, the video tasks are executed, and on Core 1, the acquisition threads are executed. The network threads are also configured in Core 1. This load balancing and the use of UDP sockets ensure a correct latency of the application. It should be noted that the application performance degrades if the TCP socket is used, not reaching 8 FPS compared to the stable performance of 12.5 FPS when using the UDP socket. In our UDP protocol, the timestamp and frame counter are included in the same way as in RTP (Real-time Transport Protocol), and the receiver must use them to detect packet loss and to accommodate out-of-service delivery.
In general, the platform is in the same local network as the server and the use of UDP does not require specific configuration in the routers.
Finally, there is a control command shell accessible via the TCP socket. These commands allow to start or stop the different data streams and to parameterize some additional configurations.
5. Experimental Results
In this section, the raw data from the different sensors and devices included in the developed platform are shown. Where possible, the equivalent results sampled by means of the Empatica E4 wristband are additionally shown. This comparison is of particular relevance for EDA, heart rate and body temperature. Although the Empatica has a three-axis accelerometer, it does not have a complete AHRS system. It is also unable to record facial expressions, so a comparison of results in the corresponding MsWH subsystems is not possible or relevant.
5.1. Power Consumption
This section provides the power consumption data for the different devices of the platform. The results are summarized in
Table 2, in which the experimental measurements of power consumption are detailed in addition to their theoretical values, which are extracted from the data sheets.
5.2. Sensors and Calibration
As mentioned above, the ST (MAX30205), HR and SpO2 (MAX30101-MAX32664) and IMU (BNO055) sensors are factory calibrated, with measurement accuracy guaranteed by the manufacturer. The same is true for the OV2640 that controls the fisheye camera.
The facial recognition camera can work with 66° and 160° (fisheye) lenses. The distance between the camera lens and the individual’s forehead is 170 mm. Optical distortion can be corrected easily, since the distance to the plane remains constant throughout the recording. However, if necessary, a calibration standard can be used to more accurately correct for lens distortion. The frames acquired by the camera have been analyzed with the Amazon Rekognition Emotion API (
https://docs.aws.amazon.com/rekognition/index.html, accessed on 14 June 2022), interpreting emotional expressions such as happiness, sadness and surprise with satisfactory results.
In the case of the ST and AHRS (IMU) measurements, the calibration process is included in the documentation detailed respectively in
Section 3.1 and
Section 3.3.
In the case of ST, HR and EDA, the signals obtained by the MsHW platform can be compared with those obtained by the Empatica E4 wristband. The results are shown in the following sections.
5.3. The Empatica E4 Wristband
The Empatica E4 wristband is a wearable wireless device designed for comfortable, continuous, real-time data acquisition in daily life. It contains four sensors:
Photoplethysmography (PPG) to provide blood volume pulse, from which heart rate, heart rate variability, and other cardiovascular features may be derived.
Electrodermal activity (EDA), used to measure sympathetic nervous system arousal and to derive features related to stress, engagement, and excitement.
Three-axis accelerometer, to capture motion-based activity (IMU).
Infrared thermometer, reading skin temperature (ST).
5.4. Data Capture Comparison
The comparison of the experimental data will be carried out in three different ways. Firstly, a comparison will be made in the measurement of skin temperature both from the point of view of the experimental results obtained and the response time (a relevant data, given the difference in the technologies used in both devices). Then, a comparison is made by obtaining the HR (or interval between heartbeats) in both devices. Finally, the results are obtained in EDA by measuring a set of resistances (or impedances) of known value, and the (more appropriate) direct measurement on the skin of an individual will be compared.
It should be noted that although both the developed platform and the Empatica E4 wristband have accelerometers, the results of the measurements with these sensors cannot be compared, since they are placed in different places on the body, and it is extremely difficult to design an experiment in which the movements of both devices are compatible and comparable.
5.4.1. Skin Temperature
For skin temperature measurement, the resolution of the E4 is of 0.02 °C, while the MsWH resolution is of only 0.0039 °C. However, the accuracy of both systems turns out to be similar (±0.2 °C for the E4 and ±0.1 °C for the MsWH, both within 36–39 °C). The final difference in skin temperatures measured through E4 and MsWH differs by only 0.1 °C. This result was obtained by fixing the sensors of both devices against each other and leaving them at rest in the laboratory, thus measuring the ambient temperature.
Another comparative variable of interest is the Step Response Time (SRT), i.e., the time required to obtain an accurate measurement when faced with an abrupt variation of the variable to be measured. Using the Matlab System Identification Toolbox, and through the data obtained by both temperature sensors, the step response of each device has been generated. These responses are shown in
Figure 13a,b.
As can be seen by comparing both subfigures, the SRTs of the compared devices are very different. This is due to the different technologies of the temperature sensors used in each case. In the E4 wristband, an optical thermometer is used. The absence of skin contact reduces the thermal inertia in the measurement, and therefore, the SRT is very low—less than 2 s, as can be seen in
Figure 13a. On the other hand, as mentioned in
Section 3.1, the MsWH uses a MAX30205, which requires direct contact with the skin and therefore has a higher thermal inertia. In this case, the SRT is almost 1 min, as can be seen in
Figure 13b.
Although very significant, this difference in response times is not important because it is the measurement of skin temperature, and variations in skin temperature are not usually significant over such short time intervals for system thermal inertia to be an issue.
Figure 14 shows a more detailed comparison of the time response between the Empatica E4 wristband and our MsWH system in which the faster response of the Empatica E4 wristband can be clearly seen.
5.4.2. Heart Rate
Previous studies [
67] suggest that the Empatica E4 is potentially usable as a tool for investigating HR and its variability (HRV), at least while the user/patient is at rest. Nevertheless, the HR capture algorithms used in the case of the MAX32664 sensor hub and the Empatica wristband are very different. In the first case, it is an algorithm implemented within the device itself with a refresh rate of 100 Hz. In the case of the E4 wristband, the output data are provided by the device itself, which are calculated through the E4 streaming server software and sent through the Bluetooth connection. The Empatica captures the Inter-Beat Interval (IBI), eliminating unclear minimums due to the presence of noise. The remaining data are used to calculate the instantaneous heart rate and are sent via Bluetooth.
This makes signal comparison between these devices a complex process. The data supplied by Empatica must be filtered and subsampled, and the phase difference between the signals must be corrected before proceeding to a comparison of results.
The HR raw results are shown in
Figure 15. The instantaneous heart rate of E4 (in red in the figure) has a higher variability than the data provided by MsHW (blue). This is due to the different calculation algorithms explained in the previous paragraph. Nevertheless, both measurements are comparable in the module, and variations are similar.
After processing the E4 wristband data (as explained before, subsampling, filtering, and synchronizing the E4 raw data), the HR measurement from both devices can be seen in
Figure 16, which represents a longer time sequence than that in
Figure 15. As can be seen in the figure, the measures of both systems are very similar, and the time delay is negligible. In this figure, the HR measurement of the MsWH is extracted directly from the raw data.
The difference between the MsWH and the E4 HR measurements can be calculated in percentage and plotted as a histogram, as has been shown in
Figure 17.
As can be seen in this figure, the Gaussian fit shows that the difference in the HR measurement between these two devices is approximately of −0.14%, which can be used to validate the HR measurement of our system, assuming that the measurement of the E4 wristband is also accurate.
5.4.3. Electro-Dermal Activity and Bioimpedance
The performance of the EDA and bioimpedance measurement systems has been verified by means of two simple experiments: in the first one, the correct calibration of the equipment is checked by measuring a set of seven precision resistors’ (, , , , , and ), with a tolerance, which is equivalent to conductance values between S and S. In the second experiment, the measurements are carried out over the skin of the same subject.
The results of the calibration check of both devices are shown in
Table 3. As mentioned above, we have a set of seven precision resistance values (or in other words, conductance values). These values are shown in the table (columns
and
) and will be taken as merely informative. The complete calibration results are represented in
Figure 18.
First, the resistors were measured with a high-performance multimeter (Fluke) to check their actual resistance and conductance values (columns
and
in
Table 3). Next, the conductance values were measured using first the E4 wristband and then with our platform, MsWH. The results of these measurements are shown in the
and
columns. Based on the Fluke conductance results, we can easily calculate the absolute and relative errors of both devices (columns
and
for absolute errors and
and
for relative errors).
As can be seen both in
Table 3 and
Figure 18, the smallest error in the accuracy of the measurement taken with the Empatica wristband is as large as the biggest error made with the MsWH, which is better calibrated.
Once this point is verified, we proceed to perform a set of simultaneous measurements with both devices under equivalent measurement conditions and on the same individual. For this purpose, the alternative electrodes of the E4 will be used to measure EDA on the fingers. In a first experiment, the E4 is placed on the left hand and the MsWH is placed on the right hand, and then, the experiment is repeated by reversing the positions. By proceeding in this way, an attempt is made to minimize the effect of possible local differences in conductance between the two hands. Changes in skin conductance turn out to be very slow, so the duration of the experiment should be relatively long (just over 18 min).
The measurement processes of the two devices are different. As explained in
Section 3.5.2, the MsWH uses an AC excitation signal at a frequency of 100 Hz. On the other hand, the E4 wristband uses a DC signal whose polarity is reversed between successive measurements. This is equivalent to saying that both devices measure skin conductance in different layers and may partly explain the experimental differences obtained. Similarly, the capacitance of the human body differs in measurements taken in DC and AC.
In
Figure 19, the measurement results taken with the MsWH are shown. Being an AC signal, the results can be shown both in modulus (clear blue) and in phase (orange). This figure should be considered only as a sample of the raw results collected by the platform since, in its state, it cannot be compared with the results of the E4 wristband.
The MsWH module and phase data are combined to obtain a single reading, which is the real part of the signal obtained. This signal can be compared with that obtained with the E4 wristband. The results are shown in
Figure 20.
As can be seen in
Figure 20, although the numerical results differ notably (in the figure, the right and left axes differ in scale), their relative behavior (in terms of general trend and more or less abrupt variations) is quite similar.
6. Conclusions and Future Work
In the present work, a new multi-sensory platform for the collection of physiological signals, that will be used in emotion recognition, has been presented. The MsWH platform has been developed using both sensors and complete state-of-the-art measurement systems, seeking the highest precision in the measurement of physiological signals and the greatest versatility in the measurement and analysis processes available on the market. It is a device designed to record experimental data, so it is not intended as a commercial or wearable device. Although bulky, the complete platform is very light (168 g), and the battery has an autonomy of approximately 2 h with all the devices at full performance, which allows long-term experiments to be carried out, in which the amount of data recorded and analyzed are massive. Although the information can be partially analyzed within the system itself, the device has a high transfer rate Wi-Fi connection (including raw data as well as image and/or video), which allows the collected results to be sent to an external device for more detailed analysis.
The platform is divided into two subsystems. The first of them is fixed to the subject’s forehead with an elastic band. In it, the temperature of the skin, the EDA/bioimpedance and the inertial system (in addition to the camera) are recorded. The second subsystem is placed in the earlobe and measures blood oxygen saturation and heart rate. The situation of all of these sensors is optimal for several reasons. The forehead is a sufficiently smooth and flat surface to be able to fix the thermal sensor in almost permanent contact. In addition, it is one of the areas of the body where the EDA appears with greater amplitude, thus facilitating its measurement. In addition, the platform that supports the camera barely moves, so the facial reactions of the subject are accurately recorded. With regard to the earlobe, it is one of the most common areas for measuring oxygen saturation in the blood, after the fingers. Although many of the registered variables can be obtained in a wristband wearable system (more comfortable to wear and allowing the recording of, for example, arm movements, perhaps more efficient for expressing emotions), the truth is that it is usually more of a drawback than an advantage in the recording of other variables such as temperature or EDA, where vibrations and the movement of sensors and electrodes can be a source of significant noise. By placing these measure systems on the user’s forehead, the electrical noise associated with movement is reduced, without losing the ability to record movements and emotional reactions with our Attitude and Heading Reference System. Extending the use of this platform to the measurement of physical parameters with changes in intensity or physical activity will require appropriate sensor calibration and sensor validation.
Leaving aside the work of the porthole camera (whose vision area remains practically fixed, with the advantages that this entails), there are five signals recorded in total. With regard to skin temperature, a sensor has been used that requires thermal contact with the subject. Although the precision achieved is practically identical to that of Empatica, our measurement shows a high thermal inertia. The possibility of replacing the sensor presented here with another one based on contactless IR technology, with similar precision but faster response, has been assessed. Regarding the collection of EDA data, as has been said, the position of the sensors offers advantages with respect to, for example, the wrist (Empatica), since they are less sensitive to user movements. Multiplexing a dual analog signal increases the accuracy and ability to calibrate the device, and the use of an AC signal offers inherent advantages over DC or switched DC measurements. Furthermore, by placing the inertial measuring unit on the user’s head, the aim is to prioritize all kinds of gesticulations and movements that facilitate emotional classification in the future based on the set of physiological signals that this device is capable of capturing.
Author Contributions
Conceptualization, D.A., J.P.d.L. and J.R.B.; methodology, D.A., J.P.d.L. and J.R.B.; software, D.A.; validation, D.A., J.P.d.L. and J.R.B.; formal analysis, D.A., J.P.d.L. and J.R.B.; investigation, D.A., J.P.d.L. and J.R.B.; resources, D.A., J.P.d.L. and J.R.B.; data curation, D.A. and J.P.d.L.; writing—original draft preparation, J.P.d.L.; writing—review and editing, D.A., J.P.d.L. and J.R.B.; visualization, D.A. and J.P.d.L.; supervision, J.R.B.; project administration, J.R.B.; funding acquisition, J.R.B. All authors have read and agreed to the published version of the manuscript.
Funding
This research has been partially supported by the Spanish Science, Innovation and University Ministry by the RTI2018-096986-B-C31 contract and the Aragonese Government by the AffectiveLab-T60-20R project.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Acknowledgments
Our thanks to Alvaro Rina Corrales, student at the La Almunia Polytechnical University School (EUPLA), who has developed part of the hardware described in this work.
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
The following abbreviations are used in this manuscript:
AHRS | Attitude and Heading Reference System |
API | Application Programming Interface |
BP | Blood Pressure |
DFT | Discrete Fourier Transform |
ECG | Electrocardiogram |
EDA | Electro-Dermal Activity |
EMG | Electromiogram |
ESP-IDF | Espressif IoT Development Framework |
EVA | Ethyl–Vinyl Acetate |
FIFO | First In First Out Memory |
GPS | Global Positioning System |
HID-I2C | Human Interface Device over I2C protocol |
HR | Heart Rate |
HRV | Heart Rate Variability |
IBI | Inter-Beat Interval |
IMU | Inertial Measurement Unit |
LIPO | Lithium Polymer |
MCU | Micro-Controller Unit |
MsWH | Multi-sensor Wearable Headband |
PCB | Printed Circuit Board |
PPG | Photoplethysmography |
RR | Respiration Rate |
SiP | System in Package |
SPI | Serial Peripheral Interface |
SpO2 | Oxygen Saturation |
SRT | Step Response Time |
ST | Skin Temperature |
TIA | Trans-Impedance Amplifier |
TCP | Transmission Control Protocol |
UART | Universal Asynchronous Receiver/Transmitter |
UDP | User Datagram Protocol |
References
- Bonato, P. Wearable Sensors/Systemsand Their Impact onBiomedical Engineering. IEEE Eng. Med. Biol. Mag. 2003, 22, 18–20. [Google Scholar] [CrossRef] [PubMed]
- Channa, A.; Popescu, N.; Skibinska, J.; Burget, R. The rise of wearable devices during the COVID-19 pandemic: A systematic review. Sensors 2021, 21, 5787. [Google Scholar] [CrossRef]
- Xu, Y.; Liu, G.Y. A method of emotion recognition based on ECG signal. In Proceedings of the 2009 International Conference on Computational Intelligence and Natural Computing, CINC 2009, Wuhan, China, 6–7 June 2009; pp. 202–205. [Google Scholar] [CrossRef]
- Yan, H.J.; Yoo, L. A Low-Power Portable ECG Touch Sensor with Two Dry Metal Contact Electrodes. J. Semicond. Technol. Sci. 2010, 10, 300–307. [Google Scholar] [CrossRef] [Green Version]
- Agrafioti, F.; Hatzinakos, D.; Anderson, A.K. ECG pattern analysis for emotion detection. IEEE Trans. Affect. Comput. 2012, 3, 102–115. [Google Scholar] [CrossRef]
- Kumar, V.; Upadhyaya, D. Improved Performance of Human Emotion Detection Using ECG Signal Processing. Int. J. Sci. Eng. Res. (IJSER) 2017, 5, 58–63. [Google Scholar]
- Mohammadi, Z.; Frounchi, J.; Amiri, M. Wavelet-based emotion recognition system using EEG signal. Neural Comput. Appl. 2017, 28, 1985–1990. [Google Scholar] [CrossRef]
- Cheng, B.; Liu, G. Emotion Recognition from Surface EMG Signal Using Wavelet Transform and Neural Network. In Proceedings of the 2nd International Conference on Bioinformatics and Biomedical Engineering (ICBBE), Shanghai, China, 16–18 May 2008; pp. 1363–1366. [Google Scholar]
- Yang, G.; Yang, S. Emotion Recognition of Electromyography Based on Support Vector Machine. In Proceedings of the 2010 Third International Symposium on Intelligent Information Technology and Security Informatics, Washington, DC, USA, 2–4 April 2010; pp. 298–301. [Google Scholar] [CrossRef]
- Perdiz, J.; Pires, G.; Nunes, U.J. Emotional state detection based on EMG and EOG biosignals: A short survey. In Proceedings of the ENBENG 2017—5th Portuguese Meeting on Bioengineering, Coimbra, Portugal, 16–18 February 2017; pp. 1–4. [Google Scholar] [CrossRef]
- Tam, S.; Bilodeau, G.; Brown, J.; G., G.T.; Campeau-Lecours, A.; Gosselin, B. A Wearable Wireless Armband Sensor for High-Density Surface Electromyography Recording. In Proceedings of the 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Berlin, Germany, 23–27 July 2019; pp. 6040–6044. [Google Scholar]
- Biswas, D.; Simões-Capela, N.; Van Hoof, C.; Van Helleputte, N. Heart Rate Estimation from Wrist-Worn Photoplethysmography: A Review. IEEE Sens. J. 2019, 19, 6560–6570. [Google Scholar] [CrossRef]
- Han, D.; Bashar, S.K.; Mohagheghian, F.; Ding, E.; Whitcomb, C.; McManus, D.D.; Chon, K.H. Premature atrial and ventricular contraction detection using photoplethysmographic data from a smartwatch. Sensors 2020, 20, 5683. [Google Scholar] [CrossRef] [PubMed]
- Appelhans, B.M.; Luecken, L.J. Heart rate variability as an index of regulated emotional responding. Rev. Gen. Psychol. 2006, 10, 229–240. [Google Scholar] [CrossRef] [Green Version]
- Xiao, N.; Yu, W.; Han, X. Wearable heart rate monitoring intelligent sports bracelet based on Internet of things. Measurement 2020, 164, 108102. [Google Scholar] [CrossRef]
- Song, C.; Zeng, P.; Wang, Z.; Zhao, H.; Yu, H. Wearable continuous body temperature measurement using multiple artificial neural networks. IEEE Trans. Ind. Inform. 2018, 14, 4395–4406. [Google Scholar] [CrossRef]
- Maijala, A.; Kinnunen, H.; Koskimäki, H.; Jämsä, T.; Kangas, M. Nocturnal finger skin temperature in menstrual cycle tracking: Ambulatory pilot study using a wearable Oura ring. BMC Women’s Health 2019, 19, 1–10. [Google Scholar] [CrossRef]
- Peres, P.J.A. Evaluating a Sensor of Skin Conductance to Assess Dental Anxiety. Ph.D. Thesis, University of Lisboa, Lisbon, Portugal, 2015; p. 93. [Google Scholar]
- Howell, N.; Devendorf, L.; Tian, R.; Galvez, T.V.; Gong, N.W.; Poupyrev, I.; Paulos, E.; Ryokai, K. Biosignals as social cues: Ambiguity and emotional interpretation in social displays of skin conductance. In Proceedings of the DIS 2016—2016 ACM Conference on Designing Interactive Systems: Fuse, Brisbane, QLD, Australia, 4–8 June 2016; pp. 865–870. [Google Scholar] [CrossRef] [Green Version]
- Gatti, E.; Calzolari, E.; Maggioni, E.; Obrist, M. Data Descriptor: Emotional ratings and skin conductance response to visual, auditory and haptic stimuli. Sci. Data 2018, 5, 180120. [Google Scholar] [CrossRef] [PubMed]
- Setz, C.; Arnrich, B.; Schumm, J.; Marca, R.L.; Tröster, G.; Ehlert, U. Discriminating stress from cognitive load using a wearable eda device. IEEE Trans. Inf. Technol. Biomed. 2010, 14, 410–417. [Google Scholar] [CrossRef] [PubMed]
- Pope, G.C.; Halter, R.J. Design and Implementation of an Ultra-Low Resource Electrodermal Activity Sensor for Wearable Applications. Sensors 2019, 19, 2450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Borrego, A.; Latorre, J.; Alcaniz, M.; Llorens, R. Reliability of the Empatica E4 wristband to measure electrodermal activity to emotional stimuli. In Proceedings of the International Conference on Virtual Rehabilitation, ICVR 2019, Tel Aviv, Israel, 21–24 July 2019. [Google Scholar] [CrossRef]
- Longmore, S.K.; Lui, G.Y.; Naik, G.; Breen, P.P.; Jalaludin, B.; Gargiulo, G.D. A Comparison of Reflective Photoplethysmography for Detection of Heart Rate, Blood Oxygen Saturation, and Respiration Rate at Various Anatomical Locations. Sensors 2019, 19, 1874. [Google Scholar] [CrossRef] [Green Version]
- Buekers, J.; Theunis, J.; De Boever, P.; Vaes, A.W.; Koopman, M.; Janssen, E.V.; Wouters, E.F.; Spruit, M.A.; Aerts, J.M. Wearable Finger Pulse Oximetry for Continuous Oxygen Saturation Measurements During Daily Home Routines of Patients With Chronic Obstructive Pulmonary Disease (COPD) Over One Week: Observational Study. JMIR Mhealth Uhealth 2019, 7, e12866. [Google Scholar] [CrossRef]
- Casalino, G.; Castellano, G.; Zaza, G. A mHealth solution for contact-less self-monitoring of blood oxygen saturation. In Proceedings of the 2020 IEEE Symposium on Computers and Communications (ISCC), Rennes, France, 7–10 July 2020; pp. 1–7. [Google Scholar] [CrossRef]
- Fioravanti, A.; Fico, G.; Arredondo, M.; Salvi, D.; Villalar, J. Integration of heterogeneous biomedical sensors into an ISO/IEEE 11073 compliant application. In Proceedings of the 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology, Buenos Aires, Argentina, 31 August–4 September 2010; pp. 1049–1052. [Google Scholar]
- Sharma, P.; Hui, X.; Zhou, J.; Conroy, T.B.; Kan, E.C. Wearable radio-frequency sensing of respiratory rate, respiratory volume, and heart rate. NPJ Digit. Med. 2020, 3, 1–10. [Google Scholar] [CrossRef]
- Ragot, M.; Martin, N.; Em, S.; Pallamin, N.; Diverrez, J.M. Emotion recognition using physiological signals: Laboratory vs. wearable sensors. Adv. Intell. Syst. Comput. 2018, 608, 15–22. [Google Scholar] [CrossRef] [Green Version]
- Liu, F.; Esparza, M.; Pavlovskaia, M.; Kaufman, G.; Dabbish, L.; Monroy-Hernández, A. Animo: Sharing Biosignals on a Smartwatch for Lightweight Social Connection. In Proceedings of the ACM on Interactive, Mobile, Wearable and Ubiquitous Technologies; Association for Computing Machinery: New York, NY, USA, 2019; Volume 3, pp. 1–19. [Google Scholar] [CrossRef] [Green Version]
- Pham, N.; Dinh, T.; Raghebi, Z.; Kim, T.; Bui, N.; Nguyen, P.; Truong, H.; Banaei-Kashani, F.; Halbower, A.; Dinh, T.; et al. WAKE: A Behind-the-ear Wearable System for Microsleep Detection. In Proceedings of the 18th International Conference on Mobile Systems, Applications, and Services, Toronto, ON, Canada, 15–19 June 2020; Volume 20, p. 15. [Google Scholar] [CrossRef]
- Wan, C.; Chen, D.; Huang1and, Z.; Luo, X. A Wearable Head Mounted Display Bio-Signals Pad System for Emotion Recognition. Sensors 2021, 22, 142. [Google Scholar] [CrossRef]
- Menghini, L.; Gianfranchi, E.; Cellini, N.; Patron, E.; Tagliabue, M.; Sarlo, M. Stressing the accuracy: Wrist-worn wearable sensor validation over different conditions. Psychophysiology 2019, 56, e13441. [Google Scholar] [CrossRef]
- Sano, A.; Picard, R.W. Stress recognition using wearable sensors and mobile phones. In Proceedings of the 2013 Humaine Association Conference on Affective Computing and Intelligent Interaction, ACII 2013, Geneva, Switzerland, 2–5 September 2013; pp. 671–676. [Google Scholar] [CrossRef]
- Ollander, S.; Godin, C.; Campagne, A.; Charbonnier, S. A comparison of wearable and stationary sensors for stress detection. In Proceedings of the 2016 IEEE International Conference on Systems, Man, and Cybernetics, SMC 2016, Budapest, Hungary, 9–12 October 2017; pp. 4362–4366. [Google Scholar] [CrossRef]
- McDuff, D.; Rowan, K.; Choudhury, P.; Wolk, J.; Pham, T.; Czerwinski, M. A Multimodal Emotion Sensing Platform for Building Emotion-Aware Applications. arXiv 2019, arXiv:1903.12133. [Google Scholar]
- Castro-Garcia, J.A.; Molina-Cantero, A.J.; Merino-Monge, M.; Gomez-Gonzalez, I.M. An Open-Source Hardware Acquisition Platform for Physiological Measurements. IEEE Sens. J. 2019, 19, 11526–11534. [Google Scholar] [CrossRef]
- Al-khafajiy, M.; Baker, T.; Chalmers, C.; Asim, M.; Kolivand, H.; Fahim, M.; Waraich, A. Remote health monitoring of elderly through wearable sensors. Multimed. Tools Appl. 2019, 78, 24681–24706. [Google Scholar] [CrossRef] [Green Version]
- Akhbarifar, S.; Javadi, H.H.S.; Rahmani, A.M.; Hosseinzadeh, M. A secure remote health monitoring model for early disease diagnosis in cloud-based IoT environment. In Personal and Ubiquitous Computing; Springer: Berlin/Heidelberg, Germany, 2020; pp. 1–17. [Google Scholar]
- Lam, L.K.; Szypula, A.J. Wearable emotion sensor on flexible substrate for mobile health applications. In Proceedings of the 2018 IEEE Sensors Applications Symposium (SAS), Seoul, Korea, 12–14 March 2018; pp. 1–5. [Google Scholar] [CrossRef]
- Sundaravadivel, P.; Kougianos, E.; Mohanty, S.P.; Ganapathiraju, M.K. Everything You Wanted to Know about Smart Health Care: Evaluating the Different Technologies and Components of the Internet of Things for Better Health. IEEE Consum. Electron. Mag. 2018, 7, 18–28. [Google Scholar] [CrossRef]
- GlobalData. Wearable Technology in Healthcare; GlobalData: London, UK, 2019; pp. 1–83. Available online: https://www.tauli.cat/institut/wp-content/uploads/2020/06/wearables-GlobalData_WearableTechnologyinHealthcare_220819.pdf (accessed on 14 June 2022).
- Ahad, A.; Tahir, M.; Sheikh, M.A.; Ahmed, K.I.; Mughees, A.; Numani, A. Technologies trend towards 5g network for smart health-care using iot: A review. Sensors 2020, 20, 4047. [Google Scholar] [CrossRef] [PubMed]
- Wang, Y.; Yang, B.; Hua, Z.; Zhang, J.; Guo, P.; Hao, D.; Gao, Y.; Huang, J. Recent advancements in flexible and wearable sensors for biomedical and healthcare applications. J. Phys. D Appl. Phys. 2021, 55, 134001. [Google Scholar] [CrossRef]
- Olmedo Aguirre, J.O.; Reyes Campos, J.; Alor Hernandez, G.; Machorro Cano, I.; Rodriguez Mazahua, L.; Sanchez Cervantes, J.L. Remote Healthcare for Elderly People Using Wearables: A Review. Biosensors 2022, 12, 73. [Google Scholar] [CrossRef]
- Kalovrektis, K.; Xenakis, A.A.; Stamoulis, G.; Xenakis, A.A.; Gotsinas, A.; Korinthios, I. A Customizable Wireless Sensor System for Estimating and Transmitting Human Stress Levels Biosignals with IoT Platform. J. Wirel. Netw. Commun. 2020, 2020, 9–16. [Google Scholar] [CrossRef]
- Maaoui, C.; Pruski, A. Emotion Recognition through Physiological Signals for Human-Machine Communication; InTech: London, UK, 2010; pp. 317–332. [Google Scholar]
- Cowley, B.; Filetti, M.; Lukander, K.; Torniainen, J.; Henelius, A.; Ahonen, L.; Barral, O.; Kosunen, I.; Valtonen, T.; Huotilainen, M.; et al. The psychophysiology primer: A guide to methods and a broad review with a focus on human-computer interaction. Found. Trends Hum.-Comput. Interact. 2016, 9, 151–308. [Google Scholar] [CrossRef] [Green Version]
- Patel, S.; Park, H.; Bonato, P.; Chan, L.; Rodgers, M. A review of wearable sensors and systems with application in rehabilitation. J. NeuroEng. Rehabil. 2012, 9, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Giannakakis, G.; Grigoriadis, D.; Giannakaki, K.; Simantiraki, O.; Roniotis, A.; Tsiknakis, M. Review on psychological stress detection using biosignals. IEEE Trans. Affect. Comput. 2019, 13, 440–460. [Google Scholar] [CrossRef]
- Hosseini, S.A.; Khalilzadeh, M.A. Emotional Stress Recognition System Using EEG and Psychophysiological Signals: Using New Labelling Process of EEG Signals in Emotional Stress State. In Proceedings of the 2010 International Conference on Biomedical Engineering and Computer Science, Wuhan, China, 23–25 April 2010; pp. 1–6. [Google Scholar] [CrossRef]
- Fletcher, R.R.; Amemori, K.I.; Goodwin, M.; Graybiel, A.M. Wearable wireless sensor platform for studying autonomic activity and social behavior in non-human primates. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS, San Diego, CA, USA, 28 August–1 September 2012; pp. 4046–4049. [Google Scholar] [CrossRef]
- Trung, T.Q.; Lee, N.E. Flexible and Stretchable Physical Sensor Integrated Platforms for Wearable Human-Activity Monitoring and Personal Healthcare. Adv. Mater. 2016, 28, 4338–4372. [Google Scholar] [CrossRef] [PubMed]
- Girardi, D.; Lanubile, F.; Novielli, N. Emotion detection using noninvasive low cost sensors. In Proceedings of the 2017 7th International Conference on Affective Computing and Intelligent Interaction, ACII 2017, San Antonio, TX, USA, 23–26 October 2018; pp. 125–130. [Google Scholar] [CrossRef] [Green Version]
- Nie, J.; Hu, Y.; Wang, Y.; Xia, S.; Jiang, X. SPIDERS: Low-cost wireless glasses for continuous in-situ bio-signal acquisition and emotion recognition. In Proceedings of the 5th ACM/IEEE Conference on Internet of Things Design and Implementation, IoTDI 2020, Sydney, Australia, 21–24 April 2020; pp. 27–39. [Google Scholar] [CrossRef]
- Domínguez-Jiménez, J.A.; Campo-Landines, K.C.; Martínez-Santos, J.C.; Delahoz, E.J.; Contreras-Ortiz, S.H. A machine learning model for emotion recognition from physiological signals. Biomed. Signal Process. Control 2020, 55, 101646. [Google Scholar] [CrossRef]
- Cegar, D.D.; Barreda-Ángeles, M.; Kukolj, D.; Callet, P.L. Modelling effects of S3D visual discomfort in human emotional state using data mining techniques. Multimed. Tools Appl. 2020, 79, 19803–19829. [Google Scholar] [CrossRef]
- Tomé, A.M.; Hidalgo-Muñoz, A.R.; López, M.M.; Teixeira, A.R.; Santos, I.M.; Pereira, A.T.; Vázquez-Marrufo, M.; Lang, E.W. Feature extraction and classification of biosignals: Emotion valence detection from EEG signals. In Proceedings of the International Conference on Bio-Inspired Systems and Signal Processing, BIOSIGNALS 2013, Barcelona, Spain, 11–14 February 2013; pp. 54–60. [Google Scholar] [CrossRef]
- Basu, S.; Jana, N.; Bag, A.; Mahadevappa, M.; Mukherjee, J.; Kumar, S.; Guha, R. Emotion Recognition Based on Physiological Signals using Valence-Arousal Model. In Proceedings of the 2015 Third International Conference on Image Infonnation Processing Emotion, Waknaghat, India, 21–24 December 2015; pp. 50–55. [Google Scholar]
- Kołodziej, M.; Majkowski, A.; Rak, R.J.; Tarnowski, P.; Pielaszkiewicz, T. Analysis of facial features for the use of emotion recognition. In Proceedings of the 19th International Conference Computational Problems of Electrical Engineering, Banska Stiavnica, Slovakia, 9–12 September 2018; pp. 1–4. [Google Scholar] [CrossRef]
- Masai, K.; Sugiura, Y.; Ogata, M.; Suzuki, K.; Nakamura, F.; Shimamura, S.; Kunze, K.; Inami, M.; Sugimoto, M. AffectiveWear: Toward recognizing facial expression. In Proceedings of the ACM SIGGRAPH 2015 Emerging Technologies, SIGGRAPH 2015, Los Angeles, CA, USA, 9–13 August 2015. [Google Scholar] [CrossRef]
- Masai, K.; Sugiura, Y.; Ogata, M.; Kunze, K.; Inami, M.; Sugimoto, M. Facial expression recognition in daily life by embedded photo reflective sensors on smart eyewear. In Proceedings of the International Conference on Intelligent User Interfaces, Proceedings IUI, Sonoma, CA, USA, 7–10 March 2016; pp. 317–326. [Google Scholar] [CrossRef]
- Sevil, M.; Rashid, M.; Hajizadeh, I.; Park, M.; Quinn, L.; Cinar, A. Physical activity and psychological stress detection and assessment of their effects on glucose concentration predictions in diabetes management. IEEE Trans. Biomed. Eng. 2021, 68, 2251–2260. [Google Scholar] [CrossRef]
- Schmidt, P.; Reiss, A.; Dürichen, R.; Laerhoven, K.V. Wearable-Based Affect Recognition-A Review. Sensors 2018, 19, 4079. [Google Scholar] [CrossRef] [Green Version]
- Dzedzickis, A.; Kaklauskas, A.; Bucinskas, V. Human Emotion Recognition: Review of Sensors and Methods. Sensors 2020, 20, 592. [Google Scholar] [CrossRef] [Green Version]
- Imotions. Galvanic Skin Response. The Complete Pocket Guide. Imotions–Biometric Research, Simplified. 2017. Available online: https://imotions.com/blog/galvanic-skin-response/ (accessed on 14 June 2022).
- Schuurmans, A.A.; de Looff, P.; Nijhof, K.S.; Rosada, C.; Scholte, R.H.; Popma, A.; Otten, R. Validity of the Empatica E4 wristband to measure heart rate variability (HRV) parameters: A comparison to electrocardiography (ECG). J. Med. Syst. 2020, 44, 1–11. [Google Scholar] [CrossRef]
Figure 1.
Final design of the platform and fastening system. (a) Side view of the platform. (b) Electrodes of the EDA/bioimpedance subsystem. (c) View with the electrodes placed in their final position and fixed by means of the EVA foam. (d) PCB with SpO2 and HR measure systems and fastening system in pink. (1) Microcontroller, (2) camera, (3) LIPO battery, (4) EDA electrode, (5) pulse oximeter, (6) temperature sensor.
Figure 1.
Final design of the platform and fastening system. (a) Side view of the platform. (b) Electrodes of the EDA/bioimpedance subsystem. (c) View with the electrodes placed in their final position and fixed by means of the EVA foam. (d) PCB with SpO2 and HR measure systems and fastening system in pink. (1) Microcontroller, (2) camera, (3) LIPO battery, (4) EDA electrode, (5) pulse oximeter, (6) temperature sensor.
Figure 2.
PCBs of the platform. (a) Top view of the main PCB, attached to the front of the device by means of a lightweight plastic bracket. It contains the power and communications control system (RS232 and Bluetooth, plus a USB port), the main microcontroller, the chip for EDA, the gyroscope and accelerometers (IMU) and the LIPO battery (including its charger) that powers all the electronics. (b) Bottom view of the same PCB, where the camera connection port and its power sources are located (left and center, respectively). The size of this PCB is 38 mm × 124 mm. (c) Top view of the pulse oximeter PCB, where the microcontroller and most of the components associated with the device can be seen. (d) Bottom view of the pulse oximeter PCB, where the corresponding microprocessor stands out. The size of this PCB is 18 mm × 20 mm.
Figure 2.
PCBs of the platform. (a) Top view of the main PCB, attached to the front of the device by means of a lightweight plastic bracket. It contains the power and communications control system (RS232 and Bluetooth, plus a USB port), the main microcontroller, the chip for EDA, the gyroscope and accelerometers (IMU) and the LIPO battery (including its charger) that powers all the electronics. (b) Bottom view of the same PCB, where the camera connection port and its power sources are located (left and center, respectively). The size of this PCB is 38 mm × 124 mm. (c) Top view of the pulse oximeter PCB, where the microcontroller and most of the components associated with the device can be seen. (d) Bottom view of the pulse oximeter PCB, where the corresponding microprocessor stands out. The size of this PCB is 18 mm × 20 mm.
Figure 3.
Block diagram of the complete sensoring and processing system.
Figure 3.
Block diagram of the complete sensoring and processing system.
Figure 4.
Circuit diagram for the ST sensor. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). The names in green are the names of the Integrated Circuits. In red are the pin numbers of the ICs.
Figure 4.
Circuit diagram for the ST sensor. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). The names in green are the names of the Integrated Circuits. In red are the pin numbers of the ICs.
Figure 5.
Block diagram of the blood oxygen saturation and heart rate monitoring subsystem.
Figure 5.
Block diagram of the blood oxygen saturation and heart rate monitoring subsystem.
Figure 6.
Circuit diagram for the SpO2 and HR monitoring subsystem. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). The names in green are the names of the Integrated Circuits. In red are the pin numbers of the ICs.
Figure 6.
Circuit diagram for the SpO2 and HR monitoring subsystem. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). The names in green are the names of the Integrated Circuits. In red are the pin numbers of the ICs.
Figure 7.
Circuit diagram for the IMU sensor. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). The names in green are the names of the Integrated Circuits. In red are the pin numbers of the ICs.
Figure 7.
Circuit diagram for the IMU sensor. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). The names in green are the names of the Integrated Circuits. In red are the pin numbers of the ICs.
Figure 8.
Connection of the electrodes for the three measurement modes, 2, 3 and 4 electrodes. Limiting/isolation resistors are included. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). See the text for an explanation of the two, three and four wire connection.
Figure 8.
Connection of the electrodes for the three measurement modes, 2, 3 and 4 electrodes. Limiting/isolation resistors are included. This figure has been generated with KiCAD (
https://www.kicad.org, accessed on 14 June 2022). See the text for an explanation of the two, three and four wire connection.
Figure 9.
Block diagram for the measurement system. Observe the 3-electrode multiplexed configuration (see text).
Figure 9.
Block diagram for the measurement system. Observe the 3-electrode multiplexed configuration (see text).
Figure 10.
EDA measurement process. A complete cycle, including multiplexing and reading of the left (blue) and right (orange) channels, is completed in 250 ms–150 ms of which, approximately, the system remains hibernated. See text for further details. The interrupt signals correspond to the reaching of the FIFO memory threshold and the time at which the MCU receives the stored data.
Figure 10.
EDA measurement process. A complete cycle, including multiplexing and reading of the left (blue) and right (orange) channels, is completed in 250 ms–150 ms of which, approximately, the system remains hibernated. See text for further details. The interrupt signals correspond to the reaching of the FIFO memory threshold and the time at which the MCU receives the stored data.
Figure 11.
Circuit diagram for the EDA measurement unit, AD5941.
Figure 11.
Circuit diagram for the EDA measurement unit, AD5941.
Figure 12.
General description of the software and firmware of the developed platform.
Figure 12.
General description of the software and firmware of the developed platform.
Figure 13.
Step response of the body temperature measurement systems. (a) Empatica E4 wristband. (b) MsWH. s. s. Both subfigures have been obtained using the Matlab System Identification Toolbox.
Figure 13.
Step response of the body temperature measurement systems. (a) Empatica E4 wristband. (b) MsWH. s. s. Both subfigures have been obtained using the Matlab System Identification Toolbox.
Figure 14.
Time response comparison of the body temperature measurement systems. It can be clearly seen that the time response of the Empatica E4 wristband is much faster than in the case of the MsHW system. Horizontal axis: time in seconds. Vertical axis: Amplitude. Blue line: MsWH. Red line: Empatica E4 wristband.
Figure 14.
Time response comparison of the body temperature measurement systems. It can be clearly seen that the time response of the Empatica E4 wristband is much faster than in the case of the MsHW system. Horizontal axis: time in seconds. Vertical axis: Amplitude. Blue line: MsWH. Red line: Empatica E4 wristband.
Figure 15.
Representation of the raw data of the RH measurement. Horizontal axis: frame number. Vertical axis: Amplitude. Blue line: MsWH. Red line: Empatica E4 wristband.
Figure 15.
Representation of the raw data of the RH measurement. Horizontal axis: frame number. Vertical axis: Amplitude. Blue line: MsWH. Red line: Empatica E4 wristband.
Figure 16.
HR synchronized signals. Horizontal axis: frame number. Vertical axis: Amplitude. Blue line: MsWH. Red line: Empatica E4 wristband.
Figure 16.
HR synchronized signals. Horizontal axis: frame number. Vertical axis: Amplitude. Blue line: MsWH. Red line: Empatica E4 wristband.
Figure 17.
Difference between MsWH and Empatica E4 HR signals. Histogram. Horizontal axis: error (%). Vertical axis: Number of occurrences. Red line: Gaussian fit. The maximum is located at −0.14%.
Figure 17.
Difference between MsWH and Empatica E4 HR signals. Histogram. Horizontal axis: error (%). Vertical axis: Number of occurrences. Red line: Gaussian fit. The maximum is located at −0.14%.
Figure 18.
Conductance calibration. Experimental results. MsWH and E4.
Figure 18.
Conductance calibration. Experimental results. MsWH and E4.
Figure 19.
MsWH conductance measurement. Raw data. Left axis: module (S). Right axis: phase (degrees).
Figure 19.
MsWH conductance measurement. Raw data. Left axis: module (S). Right axis: phase (degrees).
Figure 20.
Comparison of measurement of skin response conductance using MsWH and E4 over the same subject. Blue line (left axis): E4 measurement (S). Red line (right axis): MsWH measurement (S).
Figure 20.
Comparison of measurement of skin response conductance using MsWH and E4 over the same subject. Blue line (left axis): E4 measurement (S). Red line (right axis): MsWH measurement (S).
Table 1.
Sensors and main characteristics. N.L.→Nonlinearity. S.R.→Sample rate. Acc.→Accuracy. The sample rate of the AD5941 is 800 ksps. Data are filtered and downsampled in order to obtain a final data sequence at only 4 Hz, which is enough to register EDA activity.
Table 1.
Sensors and main characteristics. N.L.→Nonlinearity. S.R.→Sample rate. Acc.→Accuracy. The sample rate of the AD5941 is 800 ksps. Data are filtered and downsampled in order to obtain a final data sequence at only 4 Hz, which is enough to register EDA activity.
Sensor | Device/Measure | Range | Resolution (#bits) | Bandwidth (Hz) | Other |
---|
| Accelerometer | ±2 g4 g ±8 g16 g | 14 | 8, 16, 31, 63 125, 250, 500, 1000 | N.L.: 0.5%FS |
BNO055 | Gyroscope | ±125, ±250, ±500 ±1000, ±2000 (°/s) | 16 | 12, 23, 32, 64 47, 116, 230, 523 | N.L.: ±0.05%FS |
| Magnetometer | ±1300 T X, Y axis ±2500 T Z axis | 13/13/15 | 2, 6, 8, 10 15, 20, 25, 30 | N.L.: 1% |
MAX30205 | ST | 0 °C to +50 °C | 16 | 20 (max.) | Acc.: 0.1 °C (37–39 °C) |
MAX30101 | SpO2 + HR | 640–980 nm | 18 | 50–3200 | - |
AD5941 | EDA | 100 k–10 M 10–0.1 S | 16 | 4 | S.R.: 800 ksps, Acc.: 1% |
Table 2.
Sensors and main power consumption characteristics.
Table 2.
Sensors and main power consumption characteristics.
Device/Sensor | Theoretical | Measured | Sample Rate |
---|
| | (mA) | (fps) |
---|
MAX30205 (ST) | 600 A | 0.604 | 4 |
MAX30101 + MAX32664 (SpO2 + HR) | - | 9.8 to 19.5 | 100 |
AD5941 (EDA) | - | 1.4 | 4 |
ESP32-WROVER-IE | 500 mA | 100 | - |
OV2640 | 42.7mA | 53.2 | 12.5 |
BNO055 | 12.3 mA | 12.3 | 20 |
Table 3.
EDA measurement. Conductance calibration. Experimental results. First column: Theoretical value of the resistance set. Second: related conductance values. Third: high-precision measurement of resistances. Fourth: related high-precision conductance values. Fifth and sixth: experimental conductance measurements through the E4 wristband and MsWH, respectively. Seventh and nineth: absolute and relative errors for the E4 wristband. Eighth and tenth: absolute and relative errors for the MsWH.
Table 3.
EDA measurement. Conductance calibration. Experimental results. First column: Theoretical value of the resistance set. Second: related conductance values. Third: high-precision measurement of resistances. Fourth: related high-precision conductance values. Fifth and sixth: experimental conductance measurements through the E4 wristband and MsWH, respectively. Seventh and nineth: absolute and relative errors for the E4 wristband. Eighth and tenth: absolute and relative errors for the MsWH.
() | () | () | () | () | () | () | () | (%) | (%) |
---|
250,000 | 4.0 | 251,000 | 3.984064 | 3.88 | 3.987 | 0.1042 | −0.0031 | 2.615 | −0.079 |
500,000 | 2.0 | 500,001 | 1.999996 | 1.85 | 1.999 | 0.1484 | 0.00072 | 7.4188 | 0.036 |
1,000,000 | 1.0 | 1,006,000 | 0.994036 | 0.889 | 0.9968 | 0.1047 | −0.0028 | 10.531 | −0.276 |
2,000,000 | 0.5 | 2,000,500 | 0.499875 | 0.5188 | 0.499 | −0.01897 | 0.00089 | −3.795 | 0.178 |
3,000,000 | 0.33 | 3,001,100 | 0.333211 | 0.342 | 0.332 | −0.0088 | 0.0012 | −2.654 | 0.353 |
4,000,000 | 0.25 | 4,001,000 | 0.249938 | 0.2523 | 0.249 | −0.0024 | 0.00096 | −0.956 | 0.386 |
5,000,000 | 0.2 | 5,001,200 | 0.199952 | 0.1985 | 0.1988 | 0.0014 | 0.0011 | 0.689 | 0.563 |
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).