Identification of blood pressure reflecting personalized traits using bilateral photoplethysmography

1.Introduction

Blood pressure (BP) is the most common method for diagnosing hypertension. BP refers to pressure exerted on the walls of blood vessels as blood flows along the vessels, where the pressure during heart contraction is systolic BP (SBP) and the pressure during heart expansion is called diastolic BP (DBP) [1]. BP may adapt according to the condition or surrounding situations, increasing due to increased vascular resistance as the blood vessel walls become stiffer with aging [2]. Even when measured under the same conditions, BP may vary between morning and evening. Moreover, due to the so-called white coat effect, people may become more tense or nervous in front of others wearing a white coat, such as a doctor or nurse, and their BP appears higher than usual. Furthermore, BP may appear even higher in hypertensive patients [3]. Therefore, BP is sensitive to individual traits.

Cuff-based BP measurement methods have been widely used to date. However, due to limitations in continuous measurements and poor portability, various cuffless BP measurement methods have recently been developed [4, 5, 6, 7, 8, 9, 10, 11]. The trend among these studies show the measurement of pulse transit time [4, 5, 6, 7] and pulse arrival time [8] of the pulse wave from the heart to the tips of fingers, with radiofrequency radar [4], phonocardiography [5], and photoplethysmography (PPG) [5, 6, 9, 10, 11] as the primary measuring instruments. Technique using PPG alone [10] and finger-pressing method using PPG and force sensor [11] have recently been developed. However, such PPG-based BP measurement methods are also suspected to have difficulty accurately reflecting individualized traits due to their sensitivity.

Methods that use PPG to measure bio-signals are non-invasive and capable of measuring various signals such as respiratory rate, pulse rate, pulse wave, and oxygen saturation (SpO2) using a single probe [12, 13, 14, 15]. PPG, which uses light to measure the volume inside the body, is suspected to be a very effective method for measuring even very small variance in blood flow in the minor artery and capillaries distributed throughout tissues based on the contraction-relaxation period of the heart through variance in transmitted light, even in case the absorptivity of the physical substance is not familiar.

The present study analyzed waveform patterns and signal cycle characteristics of PPG and presents a simple technique for discriminating BP status that reflects individualized traits. With a short measurement time of 60–120 s, convenience of being a non-invasive measurement method on the left or right fingers, and use as a measurement device attachable to portable smartphones, PPG is expected to become a technique that can conveniently measure and manage BP anywhere at any time. In particular, data measured from the index finger of both hands could be used to extract BP values that reflect capillary elasticity and is expected to overcome the limitations of sensitivity of existing SBP/DBP-based measurement methods through two-dimensional spatial domain recognition based on BP pattern analyses.

2.Methods and materials

PPG waveform amplitude has been reported to accurately express vascular distensibility, even under situations with large fluctuations in cardiac output, such as under anesthesia [16]. As a method for discriminating BP status that reflects individualized traits in PPG signals, the present study extracted the traits in two stages using PPG waveform amplitude. The first stage comprised data collection and preprocessing. After preprocessing, amplitudes were generated as normalized pulse wave signals that allow comparison of amplitude size between subjects. In the second stage, traits were extracted as amplitudes from the normalized pulse wave signals. The extracted traits were categorized into three BP groups (normotensive, prehypertensive, and hypertensive) using linear discriminant analysis (LDA). Figure 1 shows the flowchart for the technique proposed in the present study.

Figure 1.

Flowchart of the proposed technique.

2.1Data collection

The study population consisted of members of the general population of various ages that have agreed to participate in the study. The present study recruited candidates after receiving approval from the Institutional Review Board of Chosun University. Sitting BP was measured from members of the general population of various ages using oscillometric and PPG methods. Oscillometric BP measurements were used as reference data for the PPG BP measurements. People who had difficulty with the sitting BP measurement were excluded from the study population. A total of 105 people aged 20–79 years participated in the study. BP was categorized into normotensive, prehypertensive, and hypertensive. The hypertension stages 1 and 2 group (14.3%) consisted of those with an SBP 140–179 mmHg/DBP 90–110 mmHg; prehypertension group (29.5%) included those with an SBP 121–139 mmHg/DBP 81–89 mmHg; and the normal group comprised individuals with an SBP ⩽ 120 mmHg/DBP ⩽ 80 mmHg. PPG measurements were performed after the description of the PPG measurement method is supplied and each participant provided written informed consent.

The measuring instrument was INNO-MEDU 100 (INNOTEMS Co., Ltd., Korea), a development kit for a bio-signal measurement system using medical sensors. The measurement method involved connecting the SpO2 measurement sensor and placing the PPG probes on both index fingers with the participant in a sitting pose. Subsequently, PPG signals were measured for 120 s each and the measured signals were sent wirelessly to a PC for later analysis. Figure 2a shows the measurement equipment used in the present study and Fig.  2b shows the screenshot of the PPG measurement software.

Figure 2.

a. PPG signal measuring equipment; b. Screenshot showing PPG measurement results.

2.2Preprocessing

The comparison of waveform amplitude size between the subjects from PPG pulse wave signals measured from each subject required resampling of a certain number of waveform data points within a set period of time. In the present study, two resampling stages were performed. In the first stage, resampling was performed at a sampling rate of 1 kHz to extract 66,000 data points over 60 s from the initially extracted waveform data. The resampled data were designed to make it easier to extract heart rate within the waveforms. In the second stage, resampling was performed to eliminate overlapping waveforms present within a single heart rate among the resampled data. The sampling interval was adjusted so the pulse wave would contain the average heart rate per minute in the overlapping waveforms. The sampling was applied 23 intervals. Figure 3 shows the pulse wave signals after the two-stage resampling process.

Figure 3.

PPG signals after the two-stage resampling process.

2.4Identification

The present study applied a LDA algorithm [17] to optimally classify the BP distribution into three categories (normotensive, prehypertensive, and hypertensive) based on the MLA features. In short, LDA reduces the dimension of the feature vector for the data in such a way as to maximize the ratio of between-class scatter and within-class scatter.

If the number of subjects who participated in the experiment is M=∑i=1cCi, then M represents the PPG samples measured from all subjects and C represents the class to be classified. In the present experiment, The PPG samples of 105 subjects were classified into three groups. If the training set is X={Xi∗}i=1C, each class in C classes is formed by Xi∗={Xi⁢j∗}j=1Ci. Here, Xi⁢j∗ represents the jth feature value of class i, while SW and SB denote the within-class scatter and between-class scatter in the training set, respectively, as follows:

(2)

SW=∑i=1c∑j=1ci(Xi⁢j∗-μi)⁢(Xi⁢j∗-μi)T,

(3)

SB=∑i=1c(μi-μ)⁢(μi-μ)T

(4)

μi=1Mi⁢∑x∗∈CiX∗,μ=1M⁢∑∀X∗CX∗=1C⁢∑i=1cμi

Here, μi represents the average of each class (ci) to be classified and μ represents the average of all classes (c). LDA-based feature representation is obtained by a linear projection, W=ψT⁢X∗. At this time, ψ means the axes of the transformed space acquired by maximizing the ratio of between-class scatter (SB) to the within-class scatter (SW) in the following Eq. (5). A set of N feature basis vectors {ψn}n=1N in Eq. (5) is computed as the first N(⩽M) eigenvectors with the largest eigenvalues of SW-1⁢SB.

(5)

ψ=arg⁡max⁡|ψT⁢SB⁢ψ||ψT⁢SW⁢ψ|,ψ={ψ1,…,ψN}

Figure 4.

Result of classifying BP into three categories (normotensive, prehypertensive, and hypertensive) using PPG.

Figure 5.

a. Analysis of first-eigenvector and second-eigenvector of hypertension stages 1 and 2; b. MLA waveform distribution of hypertension stage 2; c. MLA waveform distribution of hypertension stage 1.

Figure 6.

Distribution of BP in subjects with a DBP ⩾ 100 mmHg.

Figure 7.

a. Analysis of first-eigenvector and second-eigenvector of prehypertension; b. MLA waveform distribution of prehypertension 2 (high-level); c. MLA waveform distribution of prehypertension 1 (low-level).

Figure 8.

a. Analysis of first-eigenvector and second-eigenvector of normal BP; b. MLA waveform distribution of normal 3 (high-level); c. MLA waveform distribution of normal 2 (middle-level); d. MLA waveform distribution of normal 1 (low-level).

Figure 9.

BP region discrimination chart using PPG.

3.Results

The classification of BP into three categories (normotensive, prehypertensive, and hypertensive) was determined using the first-eigenvector (1-EV) and second-eigenvector (2-EV) of the LDA algorithm using MLA feature values. The experimental data consisted of MLA data obtained from the second finger of each right and left hand of a total of 105 participants aged 20–79 years. MLA feature values consisted of a total of 120 data points, 60 from each hand. Figure 4 shows the results of classifying BP into three categories using 1-EV and 2-EV in the LDA algorithm by applying the proposed feature extraction method.

Hypertension stages 1 and 2 showed a 1-EV value within the range of -0.02 to +0.02 and a 2-EV value of ⩾ 0. In particular, the 2-EV value of +0.02 to 0.059 represented hypertension stage 1 and ⩾+0.06 represented hypertension stage 2; a higher positive 2-EV value indicated a higher BP (Fig. 5a). Figure 5b shows the entire waveforms from MLA performed on left and right hand of subjects who belong to the region with 2-EV value ⩾ 0.06 (hypertension stage 2), while Fig. 5c shows that of those who belong to the category of 2-EV value 0.02–0.059 (hypertension stage 1). In addition, subjects with a DBP ⩾ 100 showed a value of ⩽ 0.0010 as the difference between their individual mean value and the overall mean of the MLA results from all hypertension stages 1 and 2 subjects (Fig. 6). It is suspected that this phenomenon was due to the reduced flexibility in vascular elasticity being reflected in the MLA feature values using PPG.

Prehypertensive BP was distributed mostly in the region in which 1-EV was negative and 2-EV was ⩽+0.02 (Fig. 7a). In particular, middle-aged subjects (mean age, 60.5 ± 12 years) showed BP patterns in prehypertension 2, a high-level region with 1-EV -0.02–0 and 2-EV ⩽+0.02 (Fig. 7b). The younger age group (mean age, 39.2 ±14.5 years) showed BP patterns in prehypertension 1, a low-level region with 1-EV <-0.02 and 2-EV ⩽+0.02 (Fig. 7c). On the other hand, subjects belonging to the prehypertensive group with an SBP/DBP of 122–135 mmHg/74–85 mmHg rarely appeared in normal 2, a middle-level region with 1-EV 0 to +0.02 and 2-EV 0 to +0.02; however, BP patterns of those aged 20–39 years appeared in this region. Subjects S1, S2, and S3 in Fig. 7a represent such cases.

A normal BP was distributed mostly in the region with 1-EV -0.02 to a positive number and 2-EV ⩽+0.02. In particular, BP patterns for normal 3, a high-level region consisting of age group (mean age, 46.3 ± 17.4 years) with 1-EV -0.02 to 0 and 2-EV ⩽+0.02 showed the distribution of a pattern that overlapped with that of prehypertension 2, a high-level region. Subjects S1–S13 in Fig. 8a belonged to this category (Fig. 8b). The age group (mean age, 50.2 ± 14.4 years) with 1-EV 0–0.02 and 2-EV ⩽+0.02 appeared in normal 2, a middle-level region (Fig. 8c), whereas the age group (mean age, 47.4 ± 18.5 years) with 1-EV > 0.02 and 2-EV ⩽+0.02 appeared in normal 1, a low-level region (Fig. 8d). Based on combined analysis of the BP classification into three categories using PPG, “BP region discrimination chart using PPG” was constructed (Fig. 9).

4.Discussion

PPG-based BP measurement offered a method for overcoming the limitations of existing oscillometric methods, including wearing a cuff and being a non-continuous measurement method. However, despite these advantages, existing SBP/DBP-based assessments have limited ability to reflect individualized traits due to their sensitivity. The findings in the present study showed that the method of BP range recognition using waveform pattern analysis of PPG signals can overcome the sensitivity issue of SBP/DBP-based methods and could reflect individualized traits. The findings in the present study could be summarized as follows.

First, an important feature in the discrimination of BP into three categories (normotensive, prehypertensive, and hypertensive) reflecting individualized traits appeared as MLA features in the PPG signal. The BP discrimination results showed that hypertension stages 1 and 2 showed distribution of 2-EV values ⩾ 0. In particular, when the distribution of 1-EV values was within the range of -0.02 to +0.02, BP appeared higher when the positive value of 1-EV was larger and the positive value of 2-EV was larger. The distribution of prehypertension appeared mostly in the region with the 1-EV value being negative and 2-EV being ⩽+0.02. In particular, the BP pattern appearing in the region with a 1-EV 0 to -0.02 and a 2-EV 0 to +0.02 was the BP pattern of middle-aged subjects aged 40 years or older. On the other hand, the BP pattern that rarely appeared in the region with a 1-EV 0 to +0.02 and a 2-EV 0 to +0.02 was the BP pattern of subjects aged 20–39 years. The distribution of a normal BP appeared mostly in the region with a 1-EV -0.02 to a positive number and a 2-EV ⩽+0.02. In particular, the BP pattern that appeared in the 1-EV 0 to -0.02 region was the BP pattern of middle-aged subjects aged 40 years or older, which included regions overlapping with distribution of the prehypertension pattern.

Subjects aged 20–39 years who were diagnosed as prehypertensive by oscillometric BP measurement but appeared in the region overlapping with normal 2, the middle-level region, on the “BP region discrimination chart using PPG” and middle-aged subjects aged 40 years or older diagnosed with normal BP by oscillometric BP measurement but appeared in the region overlapping with prehypertension 2, the high-level region, indicated a tendency of BP pattern served as an example of how well the individualized traits reflected the MLA features. Thus, we believe that MLA features are key parameters that can provide a good explanation for individualized traits in the discrimination of three BP categories.

Second, the “BP region discrimination chart using PPG” was assessed as a meaningful BP assessment method for discriminating among the three BP categories using the MLA features of PPG signals. Middle-aged subjects aged 40 years or older who showed a normal BP on the oscillometric BP measurement showed a tendency of a BP pattern that overlapped partially with the prehypertension region in the BP region discrimination method. We suspect that the BP pattern appeared in the prehypertension region due to reduced elasticity of the capillaries caused by vascular aging and vascular calcification despite their having a normal BP level from taking antihypertensive medication. In other words, this proves that BP level may temporarily appear normal due to medication, but the fundamental treatment of vascular elasticity is difficult. Moreover, subjects aged 20–39 years who belonged to the prehypertension group based on SBP/DBP (122–135 mmHg/74–85 mmHg) in oscillometric BP measurement rarely showed a normal 2, middle-level region in the BP region discrimination method are considered cases of sensitivity of the BP to the external environment, such as the white-coat effect, that could be improved.

Third, in the discrimination of BP into three categories using the LDA classification algorithm by applying PPG signals, successful recognition was possible using the 1-EV and 2-EV components. Hypertension stages 1 and 2 showed distribution when the 2-EV value was ⩾ 0 and they were distributed within the range of 1-EV -0.02 to +0.02. On the other hand, prehypertension was distributed mostly in the region with 1-EV being a negative number and 2-EV ⩽+0.02. Normal BP was distributed mostly in the region with 1-EV -0.02 to a positive number and 2-EV ⩽+0.02. Therefore, 1-EV and 2-EV components were determined to be significant parameters for discriminating among the three BP categories. In particular, the overlapping regions were interpreted as regions that well reflected the individualized traits.

5.Conclusions

The present study described a simple technique for discriminating among three BP categories by reflecting individualized traits using waveform pattern analysis on non-invasive and simple-to-measure PPG signals. The findings of this study showed that MLA of PPG signals measured from both index fingers is an important parameter for discriminating among three BP categories (normotensive, prehypertensive, and hypertensive). Moreover, the two-dimensional region assessment method by 1-EV and 2-EV components for discriminating among three BP categories well reflected individualized traits through the overlapping regions; thus, it is expected to be a new BP assessment technique that can accurately assess BP by improving upon the sensitivity of existing SBP/DBP-based BP assessments. Based on the findings of the present study, it is believed that cuffless PPG could become a technique that can conveniently measure and manage BP anywhere and at any time and be attached to a portable smartphone used as the main device.