1 July 2011 Relations between ac-dc components and optical path length in photoplethysmography
Author Affiliations +
J. of Biomedical Optics, 16(7), 077012 (2011). doi:10.1117/1.3600769
Abstract
Photoplethysmography is used in various areas such as vital sign measurement, vascular characteristics analysis, and autonomic nervous system assessment. Photoplethysmographic signals are composed of ac and dc, but it is difficult to find research about the interaction of photoplethysmographic components. This study suggested a model equation combining two Lambert-Beer equations at the onset and peak points of photoplethysmography to evaluate ac characteristics, and verified the model equation through simulation and experiment. In the suggested equation, ac was dependent on dc and optical path length. In the simulation, dc was inversely proportionate to ac sensitivity (slope), and ac and optical path length were proportionate. When dc increased from 10% to 90%, stabilized ac decreased from 1 to 0.89 ± 0.21, and when optical path length increased from 10% to 90%, stabilized ac increased from 1 to 1.53 ± 0.40.
Lee, Sik Shin, and Lee: Relations between ac-dc components and optical path length in photoplethysmography

1.

Introduction

Photoplethysmography (PPG) is composed of a pulsatile component (ac) and nonpulsatile component (dc). ac is synchronized with the heart and related to arterial pulsation, while dc is related to light absorption in the tissue, vein, and diastolic arterial blood volume. dc has many physiological meanings and correlates with the ac waveform. 1, 2, 3, 4, 5, 6, 7, 8 According to the Lambert–Beer (or Beer–Lambert) law that explains mechanisms of PPG, PPG ac is generated by the optical path length (Δd) that changes with the increase of volume in blood vessels as they pulsate.9 Therefore, it can be said that PPG ac is in a consistent relationship generated by variations in dc and Δd.

In previous studies of ac, dc, and Δd, there were correlations between continuous ac and spontaneous variations of dc,10 but there has been no research on ac variation according to dc intensity. There have also been comparative studies of preop and postop blocking of sympathetic nerve systems or specific patients such as diabetics.3, 11, 12, 13 However, these studies only reported on experimental results and did not suggest a mathematical model that can explain the relationship between components, so they offer limited explanations of physiological results. Therefore, this study suggested an equation for the relationship between ac, dc, and components using the existing Lambert–Beer law to analyze the characteristics of PPG as a noninvasive vascular assessment tool. Results of the suggested equation were assessed through simple simulation and experiment.

2.

Material and Method

2.1.

Optical Model

Based on the Lambert–Beer law, the photodetector detects photons not absorbed into the blood. The waveform of general photoplethysmography used in many studies can be gained by inversing or subtracting the amount of light not absorbed at a constant (k).3 This is shown in Fig. 1, and the onset and peak of the photoplethysmography can be presented as Eqs. 1, 2.14

1

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} I_H = I_O {\rm e}^{ - b_{{\rm dc}} d_{{\rm dc}} } e^{ - b_{{\rm Hb} + {\rm HbO}_2 } d_{\min } } = \alpha \beta ^{d_{\min } } = k - I_{{\rm onset}}, \end{equation}\end{document} IH=IOebdcddcebHb+HbO2dmin=αβdmin=kIonset,

2

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} I_L = I_O {\rm e}^{ - b_{{\rm dc}} d_{{\rm dc}} } e^{ - b_{{\rm Hb} + {\rm HbO}_2 } d_{\max } } = \alpha \beta ^{d_{\max } } = k - I_{{\rm peak}}. \end{equation}\end{document} IL=IOebdcddcebHb+HbO2dmax=αβdmax=kIpeak.
Here, I L is the small value of reflected amount of light and the peak point of PPG. I H is the large value of the reflected amount of light and the onset point of photoplethysmography. b dc is a function of the absorption and the attenuation constants, and equation is b dc = εdc(λ)c dc. εdc(λ), c dc, d dc each represent extinction coefficient, concentration of absorbing substance, and optical path length, respectively, that form extra-vascular photoplethysmography dc. λ stands for wavelength, and Hb and HbO2 represent hemoglobin and oxyhemoglobin, respectively. d min is the basal diameter of arterial vessel before pulsation and d max is the maximum variation of arterial diameter during pulsation. The square of arterial diameter is proportionate to the amount of blood. If each ε and c are constant during a heartbeat, [TeX:] $\alpha = I_O {\rm e}^{ - b_{{\rm dc}} d_{{\rm dc}} }$ α=IOebdcddc and [TeX:] $\beta = e^{ - b_{{\rm Hb} + {\rm HbO}_2 } }$ β=ebHb+HbO2 can be substituted into Eqs. 1, 2.

Fig. 1

Typical photoplethysmographic signal. (a) A raw signal measured from a photodetector. (b) Final signal, constant k-reflected light intensity.

077012_1_1.jpg

PPG ac(I ac) is I LI H, the subtraction of onset from peak, and is shown in Eqs 3, 4, 5, and (6) using variations in optic path length Δd(= d maxd min).

3

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} I_{\rm ac} = I_L - I_H, \end{equation}\end{document} Iac=ILIH,

4

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} = \alpha \beta ^{d_{\min } } - \alpha \beta ^{d_{\max } }, \end{equation}\end{document} =αβdminαβdmax,

5

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} = \alpha \beta ^{d_{\min } } - \alpha \beta ^{\left({d_{\min } + \Delta d} \right)}, \end{equation}\end{document} =αβdminαβdmin+Δd,

6

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} I_{\rm ac} &=& \alpha ({\beta ^{d_{\min } } - \beta ^{d_{\min } + \Delta d} }) = \alpha \beta ^{d_{\min } } ({1 - \beta ^{\Delta d} })\nonumber\\ &=& (k - I_{{\rm onset}})\cdot({1 - \beta ^{\Delta d} }). \end{eqnarray}\end{document} Iac=α(βdminβdmin+Δd)=αβdmin(1βΔd)=(kIonset)·(1βΔd).
According to Eq. 6, I ac (a component of ac) is inversely proportionate to I onset, proportionate to arterial diameter variation Δd, and has nonlinear characteristics. I onset is a dc component of photoplethysmography and stands for nonpulsatile components and basal blood volume before pulsation. Δd is affected by arterial stiffness and generates I ac. When Eq. 6 is transformed like Eq. 7, Δd can be estimated.

7

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \Delta d = {\rm abs}\left[ {\ln \left({1 - \frac{{I_{{\rm ac}} }}{{k - I_{{\rm onset}} }}} \right)} \right]. \end{equation}\end{document} Δd=absln1IackIonset.
Numerical simulation was conducted to analyze the characteristics of ac components according to I onset and Δd. Constant k was set up as 10 because signals voltage range of acquisition in the most commercial product for measuring biomedical data is ±10 V.15, 16 I onset increased in units of 0.001 from 0 to 10. Δd is converted approximately 0.04 mm when photoplethysmography is generated,9 so it was increased from 0 to 1 mm in units of 1 μm. The size of PPG ac calculated from each combination of I onset and Δd was divided by the largest value, stabilized, and then analyzed.

2.2.

Subject

A total of 20 subjects (15 male, 5 female, age 23.1±3.4, BMI 22.4±2.0 kg/m2) who consented to the experimental objective and progress participated. Subjects did not have cardiovascular or respiratory disease and were prohibited from drinking alcohol, smoking cigarettes, and drinking coffee one day before the experiment.

2.3.

Experiment

In this experiment, respiration, ECG, and PPG were measured for 5 min in the supine position, as can be seen in Fig. 2. All signals were measured using ECG100C, PPG100C, and SKT100C (BIOPAC System, Inc., Goleta, California), with a 1 kHz sampling rate. PPG was measured on the left index finger, and using a TSD100 reflection-type sensor, it was set to 20× amplification in the 0.05 to 10 Hz bandwidth. ECG was used to identify the beat of PPG using a Lead II and amplified 2000× in the 0.5 to 35 Hz bandwidth.

Fig. 2

Experiment setting.

077012_1_2.jpg

Respiration was used to observe the rapid change in ac and dc due to deep inspiration, and it was measured via the Nasal method using a TSD202 sensor.17 The experiment setting was maintained at room temperature.

2.4.

Signal Processing and Analysis

The QRS of ECG was detected using the Pan-Tompkin algorithm.18 Constant 10 was added to PPG to make a relative baseline with no negative voltage, and the onset and peak point of PPG was detected using the adaptive mood detection algorithm.19 Noise was removed from respiration through an 11-point moving average filter. PPG beats with ECG beats, and no rapid changes in ac and dc due to deep inspiration were selected. Δd was derived from Eq. 5. dc and Δd were transformed from the 10th to the 90th percentile distribution for the characteristic analysis of ac, and the ac value was stabilized by dividing it with the mean ac value of dc 10th percentile. All signal processing was conducted using MATLAB (Mathworks, Natick, Massachusetts).

2.5.

Statistical Analysis

All data were expressed as mean ± standard deviation. SPSS (SPSS-IBM Inc, Chicago, Illinois) was used for statistical analysis. The ac characteristic curve was estimated using cubic regression (3rd polynomial regression). A p-value of <0.05 was considered statistically significant.

3.

Result

Figure 3 is a numerical simulation of stabilized ac characteristics with constant parameters (dc, Δd) in Eq. 6. Figure 3a shows characteristics of ac components according to Δd when dc is constant. As the dc value increased, the ac slope (sensitivity) toward increase in Δd decreased. Figure 3b shows the response of ac component characteristics according to dc when Δd is constant. When Δd increased, the slope of ac raised in dc increased.

Fig. 3

Computed ac characteristics of proposed model equation. (a) ac-Δd characteristics with fixed dc and (b) ac-dc characteristics with fixed Δd.

077012_1_3.jpg

Figure 4 shows the ac characteristics according to dc and Δd percentile distribution obtained from actual experiment. In Fig. 4a, when Δd was increased from the 10th percentile distribution to the 90th percentile distribution, stabilized ac increased from 1 ± 0 to 1.53 ± 0.40. The curve fitting equation was y = 0.000001368x 3 − 0.00022x 2 + 0.01621x + 0.8627(adj. R 2 = 0.9997, p < 0.0001). In Fig. 4b, when dc was increased from the 10th percentile distribution to the 90th percentile distribution, stabilized ac decreased from 1 ± 0 to 0.89 ± 0.21. The curve fitting equation was y = −0.3898x 3 + 0.3028x 2 − 0.07596x + 1.009(adj. R 2 = 0.9273, p < 0.0001).

Fig. 4

ac characteristics by increasing Δd and dc percentile distribution. (a) ac-Δd and (b) ac-dc.

077012_1_4.jpg

4.

Discussion

This study investigated ac characteristics according to the nonpulsatile component (dc) and optical path length (Δd) to analyze characteristics of PPG components.

In the equation suggested to analyze ac characteristics, ac was dependent on dc and Δd. dc was inversely proportionate to ac, and Δd was proportionate to ac. Nonpulsatile components are related to tissue, vein, and diastolic arterial blood volume. Light absorption of other tissue, such as skin pigment, dye, and bone, were disregarded because it does not change according to blood flow, which means weak changes in continuous PPG signals.

From nonpulsatile components, venous capacity during respiratory change is connected to respiratory-induced fluctuation in the dc component, which is highly correlated to central and periphery venous pressure.1, 6 Venous capacity corresponds with venous return that affects blood capacity of the heart in each contraction. If venous tone in consistent in the venous return curve, venous return decreases with the increase in central venous pressure and decreases stroke volume. At this point, the autonomous nervous system increases heartbeat and arterial pressure to maintain homeostasis (maintain cardiac output).2, 7 PPG ac is proportionate to stroke volume, so when PPG dc is increased due to increase in central venous pressure, this increases arterial pressure, leading to a decrease in PPG ac components.2, 7

The second nonpulsatile component in PPG is the arterial nonpulsatile component, which signifies the baseline blood capacity of the artery in contraction, and this baseline blood capacity is proportionate to diameter.20 Arterial nonpulsatile components increase as arterial capacity or diameter increases. According to the stress-strain curve of the artery, the diameter is nonlinearly proportionate to blood pressure, and high blood pressure has high stiffness.20 Therefore, an increase in arterial nonpulsatile components increases vascular stiffness, which then decreases ac. From the two cases above, ac decreases according to an increase in PPG dc shown in the proposed equation, and the experiment is in accord with existing physiological research.

ac characteristics according to dc and Δd generated transition characteristics starting from each median point. In experimental characteristics against dc, attenuation appeared differently at the 50th percentile, and as transition happened in Δd, ac increased. According to existing research,8, 21 on mechanical properties of normal subjects’ peripheral blood vessels, as blood pressure increases, compliance and distensibility of the artery exponentially decreased with an inflection point, and diameter exponentially increases. Therefore, it can be concluded that the transition characteristics of ac shown in attenuation and an increase of experimental results follow the mechanical properties of blood pressure-diameter.

In conclusion, an equation to explain ac characteristics was suggested, and its dependence on dc and Δd was verified through a simple simulation and experiment. Our results may be utilized as a basic analysis tool for noninvasive vascular assessment.

References

1. 

L. Nilsson, A. Johansson, and S. Kalman, “Respiratory variations in the reflection mode photoplethysmographic signal. Relationships to peripheral venous pressure,” Med. Biol. Eng. Comput. 41(3), 249–254 (2003).10.1007/BF02348428Google Scholar

2. 

L. Nilsson, A. Johansson, and S. Kalman, “Macrocirculation is not the sole determinant of respiratory induced variations in the reflection mode photoplethysmographic signal,” Physiol. Meas. 24(4) 925–937 (2003).10.1088/0967-3334/24/4/009Google Scholar

3. 

M. Nitzan, A. Babchenko, B. Khanokh, and D. Landau, “The variability of the photoplethysmographic signal—a potential method for the evaluation of the autonomic nervous system,” Physiol. Meas. 19(1), 93–102 (1998).10.1088/0967-3334/19/1/008Google Scholar

4. 

M. Nitzan, I. Faib, and H. Friedman, “Respiration-induced changes in tissue blood volume distal to occluded artery, measured by photoplethysmography,”J. Biomed. Opt. 11(4), 040506 (2006).10.1117/1.2236285Google Scholar

5. 

A. Reisner, P. Shaltis, D. McCombie, and H. Asada, “Utility of the photoplethysmogram in circulatory monitoring,” Anesthesiology 108(5), 950–958 (2008).10.1097/ALN.0b013e31816c89e1Google Scholar

6. 

K. H. Shelley, D. G. Silverman, A. J. Shelley, and R. G. Stout, “Method of assessing blood volume using photoelectric plethysmography,” U.S. Patent No. 2010/0016739 (2010).Google Scholar

7. 

K. H. Shelley, D. G. Silverman, and A. J. Shelley, “Volume status monitor: Peripheral venous pressure, hypervolemia and coherence analysis,” WIPO Patent No. 2010/045556 (2010).Google Scholar

8. 

Y. Tardy, J. J. Meister, F. Perret, H. R. Brunner, and M. Arditi, “Non-invasive estimate of the mechanical properties of peripheral arteries from ultrasonic and photoplethysmographic measurements,” Clin. Phys. Physiol. Meas. 12(1), 39–54 (1991).10.1088/0143-0815/12/1/003Google Scholar

9. 

T. K. Aldrich, M. Moosikasuwan, S. D. Shah, and K. S. Deshpande, “Length-normalized pulse photoplethysmography: a noninvasive method to measure blood hemoglobin,” Ann. Biomed. Eng. 30(10), 1291–1298 (2002).10.1114/1.1527046Google Scholar

10. 

M. Nitzan, A. Babchenko, and B. Khanokh, “Very low frequency variability in arterial blood pressure and blood volume pulse,” Med. Biol. Eng. Comput. 37(1), 54–58 (1999).10.1007/BF02513266Google Scholar

11. 

J. Zheng, S. Hu, S. Xin, and V. Crabtree, “Effect of postural changes on lower limb blood volume, detected with non-invasive photoplethysmography,” J. Med. Eng. Technol. 32(5), 358–364 (2008).10.1080/03091900701453990Google Scholar

12. 

X. F. Teng and Y. T. Zhang, “The effect of contacting force on photoplethysmographic signals,” Physiol. Meas. 25(5), 1323–1335 (2004).10.1088/0967-3334/25/5/020Google Scholar

13. 

G. Tanaka and Y. Sawada, “Examination of normalized pulse volume-blood volume relationship: toward a more valid estimation of the finger sympathetic tone,” Int. J. Psychophysiol. 48(3), 293–306 (2003).10.1016/S0167-8760(03)00056-4Google Scholar

14. 

J. Webster, Design of Pulse Oximeters, Taylor & Francis, London (1997).Google Scholar

15. 

National Instruments, http://www.ni.com.Google Scholar

16. 

17. 

J. Allen, J. R. Frame, and A. Murray, “Microvascular blood flow and skin temperature changes in the fingers following a deep nspiratory gasp,” Physiol. Meas. 23(2), 365–373 (2002).10.1088/0967-3334/23/2/312Google Scholar

18. 

J. Pan and W. J. Tompkins, “A real-time QRS detection algorithm,” IEEE Trans. Biomed. Eng. 32(3), 230–236 (1985).10.1109/TBME.1985.325532Google Scholar

19. 

H. Shin, C. Lee, and M. Lee, “Adaptive threshold method for the peak detection of photoplethysmographic waveform,” Comput. Biol. Med. 39(12), 1145–1152 (2009).10.1016/j.compbiomed.2009.10.006Google Scholar

20. 

W. W. Nichols, M. F. O’Rourke, and D. A. McDonald, McDonald's Blood Flow in Arteries: Theoretic, Experimental, and Cinical Principles, Oxford University Press, London (2005).Google Scholar

21. 

D. Hayoz, B. Rutschmann, F. Perret, M. Niederberger, Y. Tardy, V. Mooser, J. Nussberger, B. Waeber, and H. R. Brunner, “Conduit artery compliance and distensibility are not necessarily reduced in hypertension,” Hypertension 20(1), 1–6 (1992).Google Scholar

Chungkeun Lee, Hang Sik Shin, Myoungho Lee, "Relations between ac-dc components and optical path length in photoplethysmography," Journal of Biomedical Optics 16(7), 077012 (1 July 2011). http://dx.doi.org/10.1117/1.3600769
JOURNAL ARTICLE
5 PAGES


SHARE
KEYWORDS
Photoplethysmography

Blood

Biomedical optics

Electrocardiography

Tissue optics

Analytical research

Statistical analysis

Back to Top