Graph connection Laplacian allows for enhanced outcomes of consumer camera based photoplethysmography imaging ============================================================================================================= * Stefan Borik * Hau-Tieng Wu * Kirk H. Shelley * Aymen A. Alian ## Abstract **Object** This work aims to introduce a novel method to mitigate the global phase deviation inherent in photoplethysmography imaging (PPGI) due to hemodynamics. **Method** We model the facial vascular network captured by a consumer camera as a two-dimensional manifold, where the complex dynamics of the vascular tree leads to intricate phase variations across skin sites. Utilizing PPGI, we sample the vector field on the facial manifold encoding these intricate phase variations over different skin sites resulting from blood volume modulations. We propose leveraging the graph connection Laplacian (GCL) technique to quantify the global phase deviation, with the hypothesis that correcting this deviation can improve the quality of the PPGI signal and that the phase deviation encodes valuable anatomical and physiological information. **Result** The proposed algorithm yields a higher-quality global PPGI signal by correcting the global phase deviation estimated by GCL, emphasizing waveform features such as the dicrotic notch. The perfusion map, with the global phase deviation estimated by GCL as intensity, reflects skin perfusion dynamics influenced by varying travel distances and anatomical structures. **Conclusion** This algorithm enhances the quality of the global PPGI signal, facilitating the analysis of morphological parameters and showing promise for advancing PPGI applications in scientific research and clinical practice. ## 1. Introduction The landscape of cardiovascular health monitoring has undergone a notable transformation with the introduction of smartphone cameras. This technology allows for extraction and analysis of photoplethysmography (PPG) signals. Camera-based PPG imaging (PPGI) has evolved from the original contact-based methods which traces its origin back to the early 20th century [1]. Initially developed to measure blood volume changes in the microvascular bed, PPG has evolved over time to incorporate imaging modalities, providing a non-invasive means of assessing physiological parameters [2]. The advent of consumer electronics and advanced imaging technologies (like CCD earlier [3] and now CMOS sensors [4]) has propelled PPGI into the spotlight, particularly in the last two decades, starting with pioneering works of Blazek et al. group [3], [5]. The PPGI typically involves capturing variations in light absorption or reflection at the skin surface, primarily driven by pulsatile blood flow changes. This technique has found diverse applications, ranging from pulse oximetry [6], [7], [8], [9] and cardiovascular monitoring to more recent innovations in facial PPG for health assessments [10], [11]. The ability of PPGI to provide contactless measurements has garnered interest also in another fields such as such as vital sign monitoring [12], [13], [14], [15], [16], wound assessment [17], [18], driver state estimation [19], [20], [21] or pain evaluation [22]. The continuous refinement of PPGI methodologies and the exploration of novel applications underscore its significance in modern healthcare and wellness monitoring. Researchers and practitioners are increasingly turning to additional techniques to broaden the comparative analysis and push the boundaries of PPGI capabilities. Starting from traditional spatial averaging [12], which primarily aims to increase the signal-to-noise ratio, different methods are used, differing from each other in their complexity or approach to signal extraction. At this point, we can mention the R-G algorithm [23] combining the red and green channels in the case of PPGI via an RGB camera, or the POS (Plane-Orthogonal-to-Skin) algorithm [24], which works with all three layers, or others such as, also in this case, the historically sorted G [12], which works only with the green layer or greyscale data, the aforementioned G-R [23], PCA [25], ICA [26], CHROM [27], PBV [28], 2SR [29], nonlinear-type time-frequency analysis like synchrosqueezing transform [30], the aforementioned POS [24] or Face2PPG [31]. Another approach is deep or machine learning methods [32], [33], which bring additional potential for feature extraction, classification and understanding of complex links within PPG signals. Another approach is using a lock in amplification method with the reference PPG signal either from external device or by averaging of skin pixels from whole recorded area [34], [35]. There is unavoidable global phase (or temporal) deviation across a body of tissue caused by the properties of the vascular tree, which is a system of interconnecting and branching elastic vessels where the pulse waves (i.e. direct and backward) propagate in different directions and velocities from the heart towards the periphery and back. This is further complicated by the interaction of this vascular tree with the autonomic system causing both localized and global changes in vascular tone. The main novelty of this paper is the use of synchronization based on the graph connection Laplacian (GCL) [36], [37] to quantify and modify the global phase deviation. Already as mentioned above in a physiological context, phase deviations in PPGI signals at different skin sites are intricately linked to the dynamic blood flow pattern in the vascular network. Our hypothesis is that this hemodynamic-induced phase variation embodies valuable anatomical information and potentially other relevant details of interest. We model the vascular network underlying the facial pixels, which in this case are recorded using a consumer camera, as a two-dimensional manifold that encodes the phase of localized blood flow. The goal is to use PPGI over each pixel to reconstruct this phase information and thus obtain an integral and improved global PPGI curve that could be used to analyze its morphological parameters and at the same time improve the spatial resolution of PPGI in mapping skin perfusion dynamics. ## 2. Materials and Methods In this section, we provide a detailed description of our algorithm (the flow chart is shown in Fig. 1), which aims to extract tissue perfusion information from video recording acquired by consumer electronics (a smartphone camera) and improve the properties of the detected PPG curve in terms of its morphology and the detection of spatial dynamics across the recorded region. ![Fig. 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/09/2024.03.07.24303826/F1.medium.gif) [Fig. 1:](http://medrxiv.org/content/early/2024/03/09/2024.03.07.24303826/F1) Fig. 1: Algorithm flowchart. ### 2.1. Experimental setup The experimental setup is based on a study conducted on volunteers undergoing a lower body negative pressure (LBNP) procedure [38], with their faces being recorded using a consumer-quality camera, namely the iPhone 11 (Apple, Cupertino, CA, USA). During recording, automatic settings such as autofocus were turned off while exposure was locked to prevent changes in the image that could reduce the quality of the perfusion information. Video was recorded in H264 compression format with a frame rate of 30 Hz. The subject’s face in supine position was illuminated with white diffuse LED light. As part of the comprehensive LBNP protocol, we selected only a 30 s segment from the baseline portion of the recording. The study protocol was approved by the Institutional Review Board of Yale University (No. 2000031899). Informed consent was obtained from all subjects prior to the study. All methods were performed in accordance with the study protocol and with the Declaration of Helsinki. To verify the functionality and assess the performance of the algorithm, we selected 5 subjects (2 females and 3 males) aged 36.2±11.39. ### 2.2. Photoplethysmography signal extraction from video The first step was the extraction of PPG signals from the video recording. In this respect, the camera sensor can be thought of as a PPG sensor matrix, where in theory each pixel of the image sequence can carry spatial information about the blood supply of a given location. The SNR of consumer cameras is usually not sufficient to extract the PPG signal from only one pixel, and different techniques have been used to increase the SNR. In most cases, the first step is spatial averaging of pixels either in the form of spatial down sampling [9] or through a moving kernel with the desired overlap [39]. In our case, we average the pixel values in a box of size 30×30 px2 with 10 px overlap which represents cumulating pixel information from the areas of approx. 1 cm2 at given camera distance and its field of view. We call an overlapping region of size 30×30 pixels a channel and suppose we obtain *L* channels. This step will not only increase the SNR of the extracted PPGI signals, but also reduce the number of signal vectors to work with while still maintain a reasonable spatial resolution. The result of this process (see Fig. 2) is a raw PPG signal that we decompose into individual color layers, including red (R), green (G), and blue (B), and store on disk for further analysis. Next, we work only with the green (G) channel, which contains the most information about perfusion. This fact is related to the technical realization of the color camera (Bayer mask), the wavelength of the green light and its penetration into the tissue, as well as the absorption and scattering properties of the skin. The main luminophore in this case is blood hemoglobin [39], [40], [41]. Denote ![Graphic][1], where *N* = 900 is the number of sampling points over the 30 s segment with the frame rate of 30 Hz, to be the raw PPGI signal collected from the *i*th pixel. ![Fig. 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/09/2024.03.07.24303826/F2.medium.gif) [Fig. 2:](http://medrxiv.org/content/early/2024/03/09/2024.03.07.24303826/F2) Fig. 2: Spatial averaging and obtaining PPGI signals ### 2.3. Signal pre-processing and filtering Further signal processing proceeds by segmenting the skin using a Gaussian Mixture Model (GMM) to ensure that we continue to work only with living skin pixels where we can expect PPG signals to be present (see the dark blue region in Figure 1(a)). The segmentation is semi-automatic and allows the user to control the detected clusters and add them as needed. The next step involves the application of a discrete wavelet transform (DWT) for each PPGI channel to select a band with the range 0.4 – 4 Hz that we can assume cardiac activity. For their extraction, we use the Daubechies 4 wavelet (db4), with the number of DWT levels/octaves set to 5. Denote the resulting signals ![Graphic][2]. After the DWT, we select a segment with minimal motion artifact in the following way. First, find all channels associated with the subject’s forehead region (see the red box/ROI in Figure 1(a)) and obtain average of the associated signals, which is denoted as ![Graphic][3]. From ![Graphic][4], we manually extract a high quality ten-second interval that is less impacted by motion artifacts. Then, restrict ![Graphic][5] on the selected ten-second interval, and denote the resulting signals ![Graphic][6]. The process continues by extracting the first and second harmonics of each PPGI signal ![Graphic][7] based on the continuous wavelet transform (CWT) using the time-frequency ridge method, followed by the inverse CWT and interpolating from the original 30 Hz to 900 Hz in order to retrieve subtle phase shifts across the signals, which is denoted as *x**i* ∈ ℝ9,000. Since the PPGI signals are generally noisy and we use relatively small averaging regions to keep the spatial information, the second harmonic component in the CWT spectrum is usually overlaid with noise. To extract it, we use the frequency ridge by simply considering doubling the extracted ridge of the first harmonic component, and then extract the second harmonic component. Denote the first and second harmonics of the *i*th channel as *h*1,*i* ∈ ℝ9,000 and *h*2,*i* ∈ ℝ9,000 respectively. Next, the signal is trimmed to the length of only 3 consecutive heartbeats, taking into account both the phase alignment dependence of the data length (influenced by heart rate variability) and the exact selection starting from the peak or valley of the PPG curve. A Hilbert transform is then performed. ### 2.4. Synchronization of the extracted PPG signals Due to the inevitable global phase deviation caused by hemodynamics, the main novelty of this paper is applying the graph connection Laplacian [36], [37] to adjust the global phase deviation. In a physiological context, the phase variations in PPGI across different pixels are intricately linked to the dynamic progression of blood flow within the vascular network. Our hypothesis posits that this phase deviation induced by hemodynamics encapsulates valuable anatomical and physiological information and potentially other pertinent details of interest. We model the vascular network underlying the facial pixels as a 2-dim manifold and the geometric structure encodes the phase of blood flow changes. The goal is to use the PPGI sensors over each pixel to recreate this phase information and hence the vascular network structure. We first estimate the pairwise phase deviation between any two PPGI channels *x**i* and *x**j* using their 1st harmonics *h*1,*i* and *h*1,*j*. Compute: ![Formula][8] where *H* converts the real-form signal into its analytic companion, defined as ![Graphic][9] *Hilbert*(*h*1,*i*), with *Hilbert* the Hilbert transform, and ⟨, ⟩is the inner product of two complex signals. The entry *c**i,j* estimates the global phase shift of *x**i* and *x**j*. See [38] for more discussion about the nonlinear relationship between harmonic phase and PPG morphology and Discussion for more technical details. Since we have *L*PPGI channels, we obtain a Hermitian function *C* of size *L*×*L*, called the graph connection Laplacian (GCL), so that the (*i, j*)th entry of *C*is *c**i,j* if the ith channel and jth channel are within 5 cm distance, and 0 otherwise. We find the eigenvector of *C*, denoted as *w*1, *w*2 …, with the associated eigenvalues ranked from large to small. Then, the global phase deviation of the *i*th PPGI channel is then corrected by: ![Formula][10] where *Re* means taking the real part and the superscript * means taking the complex conjugation. The philosophy underlying this correction is that *w*1(*i*) recovers the blood flow phase underlying the *i*th channel. See discussion section for more technical details of this GCL algorithm. ### 2.5. Construction of the final non-contact PPG signal We construct several non-contact PPG signals. The first one is by a simple averaging of all unaligned PPGI channels; that is, ![Graphic][11]. The second one is averaging all aligned PPGI channels; that is ![Graphic][12]. ### 2.6. Construction of perfusion maps We propose the following approach to visualize perfusion patterns in the face. We generate images encoding the energy of the first harmonic, denoted as *I**E*,1, which is defined as the energy of the 1st harmonic of each channel. We call these maps *energy maps*. The next set of mappings is generated from the top eigenvector of GCL, *w*1. The phase encoded in ![Graphic][13], reflects the phase shift of the *l*th channel over the face. Generate a colored image on the face region with the color in each pixel reflecting the entry of *ϑ* 1. Denote the resulting images as *I**phase*. We call these maps *phase maps*. For each image, to avoid the impact of oversaturation, we replace all values larger than 95% percentile by the 95% percentile before plotting. For phase maps, we also trim colormap based on the range of median phase ±1rad of entire skin area in order to highlight subtle difference in phase shifts across the PPGI channels. ## 3. Statistical analysis To quantify the performance of phase synchronization for PPGI signal ![Graphic][14], we consider the following metric. For the raw unaligned signals, we apply the Hilbert transform and obtain the complex form ![Graphic][15]. We then convert the amplitude and phase into the phasor form at the center point: ![Formula][16] where *c* denotes the center point that we extract the phasor and ![Graphic][17]. For the aligned signals, which are the output of (2), we apply the Hilbert transform to obtain the complex form ![Graphic][18]. We then convert the amplitude and phase into the phasor form at the center point: ![Formula][19] where *c* denotes the center point that we extract the phasor and ![Graphic][20]. Denote the median of ![Graphic][21] as ![Graphic][22] and denote the median of ![Graphic][23] as ![Graphic][24]. To check if ![Graphic][25] and ![Graphic][26] have the same distribution, we apply the Kolmogorov-Smirnov test with p<0.05 viewed as statistically significant. Also, apply F-test to check if ![Graphic][27] and ![Graphic][28] have the same variance with p<0.05 viewed as statistically significant. We also performed Wilcoxon signed-rank test to investigate if the phase shift on the cheeks and forehead, where the phase deviation of the *i*th channel is quantified by the angle of ![Graphic][29], are different with p<0.05 viewed as statistically significant. Bonferroni correction is applied to correct the multiple tests. ## 4. Results The first step to verify the quality of the phase alignment of the PPGI curves was to plot the phasors separately for 1st harmonic component. The distribution of phasors of all PPGIs channels evaluated by (6) of all subjects is shown in Fig. 3, where we plot ![Graphic][30] as a vector, with the magnitude ![Graphic][31] and angle *φ**i*(*c*). It can also be noticed in Fig. 3 that the unaligned phasors (blue arrows) tend to be directed and arranged in two or more clusters for some subjects, while aligned phasors are in general well clustered in one direction. The Kolmogorov-Smirnov test showed that the two distributions are significantly different with *p* < 10−4, and the F-test showed that the two distributions have significantly different variance with *p* < 10−4. On the other hand, we shall emphasize the observation that dominant direction ![Graphic][32] of aligned PPGI signals differs from subject to subject. This is the consequence of the nonuniqueness of GCL approach inherited from the freedom of eigendecomposition. ![Fig. 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/09/2024.03.07.24303826/F3.medium.gif) [Fig. 3:](http://medrxiv.org/content/early/2024/03/09/2024.03.07.24303826/F3) Fig. 3: Visualizaton of phase alignment. Five compass plots in the top row depict arrows representing aligned (red) and unaligned (blue) states across subjects, where columns represent different subjects. In the bottom row, the alignment is carried out using the second harmonic. After aligning the PPGI signals and verifying the quality of their alignment, we obtain various final non-contact PPG signals, including *x*(*ua*−*mean*), *x*(*a*−*mean*)shown in Fig. 4. One of the interesting effects of phase alignment is the appearance of the dicrotic notch, which is most pronounced for Subject 2 in Fig. 4. Another observation is that the average of phase aligned PPGIs contributes to maximizing the amplitude of the PPGI signal. ![Fig. 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/09/2024.03.07.24303826/F4.medium.gif) [Fig. 4:](http://medrxiv.org/content/early/2024/03/09/2024.03.07.24303826/F4) Fig. 4: Comparison of averaging phase unaligned and phase aligned PPGI waveforms from all channels to generate the final non-contact PPG signal. Solid blue curves are from the unaligned PPGI signals, while dashed red curves are from the aligned PPGI signals. From left to right: results from subjects 1 to 5 respectively. Next, we demonstrate the perfusion maps derived from PPGI. See Fig. 5 and Fig. 6 for our 5 different subjects, who were selected based on specific characteristics such as skin tone, facial hair, facial mask or extensive makeup. We will pay specific attention to the difference between cheeks and forehead, the two dominant areas on the face that are supplied by different arteries. Fig. 5: Results of the two perfusion maps for the first three subjects with diverse characteristics, including skin tone, facial hair, mask. For each subject, the left subfigure is the energy map, and the right subfigure is the phase map. Figures not shown due to the Medrxiv’s policy. Fig. 6: Results of the two perfusion maps for the last two subjects with diverse characteristics, including skin tone and makeup. For each subject, the left subfigure is the energy map, and the right subfigure is the phase map. Figures not shown due to the Medrxiv’s policy. For subject 1 we can observe different captured amplitudes in the cheeks and glabella region in the energy maps, which might be interpreted as ‘‘cheeks are more perfused’’. Another interesting observation for subject 1 is its phase map associated with *w*1, where we can see clustering at specific locations in the image. There is a slightly different phase shift in the right side of the face. The forehead area also exhibits different phase shift compared to other parts of the face which could be related to different arterial supply. Subject 2 is specific to his facial hair (chin and beard), which may be a source of motion artifacts that may deceptively influence the detected PPGI signal amplitude. Again, the cheeks have more dominant energy compared to forehead also in this case. In the phase map, while it is less obvious, the phase in the forehead is different from that in the cheek. Subject 3 is specific in that he is wearing a face shroud and also has goggles on his eyes. For this case, it is not possible to obtain information over cheek, so a comparison between cheeks and forehead is not possible. For subject 4 (Fig. 6), we notice a distinct vertical line across the middle of the forehead, which may be a wrinkle or an ongoing supratrochlear vein. This area is even more pronounced in the phase map. Again, we can see a phase difference between the forehead and cheeks in the phase map. We can also see an interesting result for subject 5 (Fig. 6), where the glowing eye surroundings can be tracked in the energy maps. In this case, this is the result of different makeup on the face and around the eyes. In this case, this is a consequence of the different makeup on the face and around the eyes. In this case, the phase difference between the forehead and cheeks is less clear in the phase map. Wilcoxon signed-rank test showed that the phase deviation is significantly different in forehead and cheeks for subject 1, 2 and 4 (with *p* < 10−4). We did not test subject 2 because we cannot access the cheeks due to the facial mask. We also omitted subject 5 due the extensive makeup. ## 5. Discussion and Conclusion Our aim is to advance the signal processing techniques for PPGI further in two directions, both of which can be applied in different areas. The first attempt is applying the GCL-based synchronization algorithm to construct a better-quality global PPG, which contains details that are more challenging to obtain using traditional approaches, like dicrotic notch and its precisely localized position. The novelty is viewing GCL as a coupler for all PPGI signals so that we can determine the phase discrepancy arises from the pulse wave travel through the vascular tree, and hence synchronize all PPGI signals. The results show that this unique method proposed by us can enhance the features of the PPG waveform (see Fig. 4) and boost its amplitude. The second attempt is introducing novel images called perfusion images, which encode hemodynamics extracted from GCL. The novelty is using the phase information hidden in the eigenvectors of GCL, which directly encodes the global phase deviation caused by the blood flow in the complicated vascular network at different locations of the imaged area. Physiologically, the supraorbital and supratrochlear arteries dominate the blood supply to the forehead, while the facial artery predominantly supplies blood to the cheeks. The supraorbital and supratrochlear arteries are supplied by the internal carotid while the facial artery originates from the external carotid artery. Due to the difference in the travel distance and anatomical structure, it is reasonable to hypothesize that the PPG phases are different in these two regions. Our findings support this hypothesis. We shall elaborate some technical details. In general, PPG signal is not sinusoidal and can be represented by multiple harmonics [42]; that is, we can model a clean PPG signal as ![Graphic][33], where *A* is a positive smooth function describing the amplitude modulation, *s* is a 1-periodic function describing the morphology of the nonsinusoidal oscillation, *φ*(*t*) is a strictly monotonically increasing function describing the phase of the PPG so that its derivative describes the instantaneous heart rate, *L*> 1 is the number of harmonics, and α*l* > 0 and *β**l* ∈ [0,2*π*) come from the Fourier series of *s*. In this case, recovering the phase *φ*(*t*) from the Hilbert transform of *x*(*t*) is inaccurate due to the existence of harmonics. This situation is worsened by the existence of noise. Instead, we could better recover the phase from the fundamental component (the first harmonic) *h*1(*t*) ≔ *A*(*t*)α1cos (2*πφ*(*t*) + *β*1), thanks to the Nuttall theorem [43] and the assumption that *A*(*t*) in general oscillates slowly compared with the cos (2*πφ*(*t*) + *β*1). This is the main reason we design the algorithm with the first harmonic in (1). On the other hand, note that since the PPGI signals *x*1 and *x**j* are usually noisy, *c**i,j* is a noisy estimate of the global phase shift of *x**i* and *x**j*. The capability of recovering phase from the usually noisy PPGI signal with the GCL-based synchronization algorithm is supported by the robustness of GCL to noise [37], [44]. In addition to the robustness of GCL, we shall also discuss the geometric meaning of GCL. Consider a simplest and ideal situation that PPGI is clean and faithfully representative of the blood flow changes, and there is no geometric constraint among pixels. In this case, *C* can be constructed without setting 0; that is, the (*i, j*)-th entry of *C*is simply *c**i,j* for any *i, j*. In this case, *C* is a rank one matrix as *C*= *θθ*∗, where *θ* ∈ ℂ*L* so that *θ*(*i*) encodes the phase of the blood flow underlying the *i*th channel. In this simplest case, the top eigenvector of *C*gives the full information of phase in each channel. However, in practice the PPGI is noisy, and there are geometric constraints imposed by the vascular network. We thus only trust the phase information when two channels are close by, which leads to the construction of *C* in our algorithm. Geometrically, the top eigenvector of *C* gives us an estimate of the most synchronized phase information over the vascular network on the face. In this paper, we focus only on the top eigenvector. The higher eigenvectors of *C* encode more geometric information associated with the underlying geometric structure associated with the blood vessel distribution. We will explore this topic in our future work. It is noteworthy to mention that the phase of PPG signals (in our case, it is captured by the first eigenvector of GCL) has been studied in different contents. In [12], the authors applied Discrete Fourier Transform to possibly reveal carotid artery or differentiating normal skin with comparison with port wine stains. Another approach can be found in [45], where the lock-in amplification is used for phase estimation. It is worth noting that the phase shift can be calculated across different frequency bands across PPGI signals (in our case, we focus on the frequency band associated with the first harmonic), e.g. low-frequency oscillation phase distribution which shows as promising tool to study autonomic nerve system responses to external stimuli [10]. Studying and mapping the phase shift across imaged area may reveal the locations that the propagating pulse wave reached at the same time, or the positions of the subcutaneous tissue perforators. It also encodes the locations where motion or ballistocardiographic (BCG) artifacts may be present, or we may have encountered a counter-phase of the original PPGI signal [46], [47], [48]. To the best of our knowledge, perfusion maps generated using the phases of the first eigenvector of GCL have not yet been published and may provide important robust phase estimation tool which is immune to inevitable noise with theoretical support and can be used in different PPGI scenarios. Examples of applications of our proposed signal processing technique could be used in different fields, e.g. the body’s response to external stimuli, such as temperature (cold [10] or heat [49]) or mental stress [50]. While contact PPG has been widely applied in practice, there are situations that PPGI is needed. For example, it is not possible to place a contact sensor on the patient’s skin when the patient’s skin is damaged (injury, burns), or when monitoring the vital signs of premature newborns. Other applications are the extraction of physiological parameters such as HR or blood pressure [47], [48]. Given that the dicrotic notch and its determined position play an important role in studying the condition of the arterial tree, particularly with regard to its stiffness (as indicated by augmentation or stiffness index, etc.), our proposed algorithm shows promise in effectively capturing these features. As such, it holds potential for application in this direction. In addition, as some studies have already pointed, the PPGI has still untapped potential. Namely the challenge is to focus on truly spatial mapping of physiological parameters, especially in the context of ANS regulation or allergic reactions [10], [11], [51], [52]. There are also approaches in which the respiratory waveform is extracted from the PPGI recording or respiratory activity is detected using a method such as optical flow [11], [53] where the proposed algorithm could help too. Another option is to deploy this method in extracting information about SpO2, and thus in evaluation of the oxygenation level of arterial or venous or mixed blood [8], [9]. Finally, we want to summarize what our paper offers and how to use the tools. To obtain a higher quality global PPGI signal when integrating all PPGI signals from the sensed area, we use GCL for their phase alignment. This results in higher amplitude and highlighted details such as dicrotic notch. If we want to work with spatial information and analyze the distribution of perfusion parameters across the imaged area, a good step here is to analyze higher eigenvectors of *C*. If we are interested in the phase relations between the PPGI waveforms, we can again use the GCL, which conveys information about the global phase of a given PPGI waveform at a given location in the image. Here we would like to point out the limitations of this work. Although this is a small sample size, we would like to stress that this is mainly a methodological work using examples of subjects with a wide variety of characteristics. Even if the GCL is immune to noise, the quality of the PPGI signals may not always be sufficient. Instead of simple spatial averaging, other raw signal processing methods presented in the Materials and Methods should also be considered, which may offer higher SNR or robustness to motion artifacts in particular. Moreover, capturing the subtle differences between the phase shifts given by the pulse wave propagation velocity and the geometrical dimensions of the analyzed region is challenging, especially in terms of the demands on the sampling rate (interpolation may no longer reveal and reconstruct all details) and the size of the processed data. We may need a ground truth of the perfusion dynamics on the face in some form to further validate the proposed algorithm. In conclusion, despite the shortcomings, we believe that this novel approach will open the door for new applications of PPGI and its use not only in science but also in clinical practice. ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Conflicts of Interests/Financial Disclosures None ## Author contributions SB: idea, literature review, data analysis and write-up. HTW: idea, literature review, data analysis and write-up. KHS: idea, literature review, data collection, and write-up. AAA: idea, literature review, data collection, and write-up. ## Acknowledgments The work of Stefan Borik was supported by the National Scholarship Programme of the Slovak Republic (Application No. 47518). Heartfelt thanks also go to Prof. Vladimir Blazek from RWTH Aachen, a native of the former Czechoslovakia, who is considered as the pioneer of photoplethysmography imaging and who contributed greatly to the field and the scientific and personal growth of Stefan Borik. * Received March 7, 2024. * Revision received March 7, 2024. * Accepted March 9, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## Reference 1. [1]. A. B. Hertzman, “Photoelectric plethysmography of the fingers and toes in man,” Proceedings of the Society for Experimental Biology and Medicine, vol. 37, no. 3, pp. 529–534, 1937. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3181/00379727-37-9630&link_type=DOI) 2. [2]. A. A. Alian and K. H. Shelley, “Photoplethysmography,” Best Practice & Research Clinical Anaesthesiology, vol. 28, no. 4, pp. 395–406, Dec. 2014, doi: 10.1016/j.bpa.2014.08.006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bpa.2014.08.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25480769&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F09%2F2024.03.07.24303826.atom) 3. [3]. T. Wu, V. Blazek, and H. J. Schmitt, “Photoplethysmography imaging: a new noninvasive and noncontact method for mapping of the dermal perfusion changes,” presented at the Optical Techniques and Instrumentation for the Measurement of Blood Composition, Structure, and Dynamics, International Society for Optics and Photonics, 2000, pp. 62–70. 4. [4]. Y. Sun and N. Thakor, “Photoplethysmography Revisited: From Contact to Noncontact, From Point to Imaging,” IEEE Trans. Biomed. Eng., vol. 63, no. 3, pp. 463–477, Mar. 2016, doi: 10.1109/TBME.2015.2476337. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TBME.2015.2476337&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F09%2F2024.03.07.24303826.atom) 5. [5]. O. Such, S. Acker, and V. Blazek, “Mapped hemodynamic data acquisition by near infrared CCD imaging,” in Proceedings of the 19th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. “Magnificent Milestones and Emerging Opportunities in Medical Engineering” (Cat. No.97CH36136), Chicago, IL, USA: IEEE, 1997, pp. 637–639. doi: 10.1109/IEMBS.1997.757694. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/IEMBS.1997.757694&link_type=DOI) 6. [6]. W. Verkruysse, M. Bartula, E. Bresch, M. Rocque, M. Meftah, and I. Kirenko, “Calibration of Contactless Pulse Oximetry,” Anesthesia & Analgesia, vol. 124, no. 1, pp. 136–145, Jan. 2017, doi: 10.1213/ANE.0000000000001381. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1213/ANE.0000000000001381&link_type=DOI) 7. [7]. A. Moço and W. Verkruysse, “Pulse oximetry based on photoplethysmography imaging with red and green light: Calibratability and challenges,” J Clin Monit Comput, vol. 35, no. 1, pp. 123–133, Feb. 2021, doi: 10.1007/s10877-019-00449-y. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10877-019-00449-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F09%2F2024.03.07.24303826.atom) 8. [8]. M. Van Gastel, W. Wang, and W. Verkruysse, “Reducing the effects of parallax in camera-based pulse-oximetry,” Biomed. Opt. Express, vol. 12, no. 5, p. 2813, May 2021, doi: 10.1364/BOE.419199. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1364/BOE.419199&link_type=DOI) 9. [9]. M. Van Gastel and W. Verkruysse, “Contactless SpO 2 with an RGB camera: experimental proof of calibrated SpO 2,” Biomed. Opt. Express, vol. 13, no. 12, p. 6791, Dec. 2022, doi: 10.1364/BOE.471332. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1364/BOE.471332&link_type=DOI) 10. [10]. S. Borik, S. Lyra, V. Perlitz, M. Keller, S. Leonhardt, and V. Blazek, “On the spatial phase distribution of cutaneous low-frequency perfusion oscillations,” Sci Rep, vol. 12, no. 1, p. 5997, Dec. 2022, doi: 10.1038/s41598-022-09762-0. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-022-09762-0&link_type=DOI) 11. [11]. S. Borik et al., “On the cardiorespiratory coordination assessed by the photoplethysmography imaging technique,” Sci Rep, vol. 13, no. 1, p. 14645, Sep. 2023, doi: 10.1038/s41598-023-41828-5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-023-41828-5&link_type=DOI) 12. [12]. W. Verkruysse, L. O. Svaasand, and J. S. Nelson, “Remote plethysmographic imaging using ambient light.,” Optics express, vol. 16, no. 26, pp. 21434–21445, 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1364/OE.16.021434&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19104573&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F09%2F2024.03.07.24303826.atom) 13. [13]. L. Tarassenko, M. Villarroel, A. Guazzi, J. Jorge, D. A. Clifton, and C. Pugh, “Non-contact video-based vital sign monitoring using ambient light and auto-regressive models,” Physiological measurement, vol. 35, no. 5, p. 807, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/0967-3334/35/5/807&link_type=DOI) 14. [14]. M. Kumar, A. Veeraraghavan, and A. Sabharwal, “DistancePPG: Robust non-contact vital signs monitoring using a camera,” Biomedical Optics Express, vol. 6, no. 5, p. 1565, May 2015, doi: 10.1364/boe.6.001565. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1364/boe.6.001565&link_type=DOI) 15. [15]. N. Blanik et al., “Remote vital parameter monitoring in neonatology–robust, unobtrusive heart rate detection in a realistic clinical scenario,” Biomedical Engineering/Biomedizinische Technik, vol. 61, no. 6, pp. 631–643, 2016. 16. [16]. W. Wang, A. C. den Brinker, S. Stuijk, and G. de Haan, “Robust heart rate from fitness videos,” Physiological measurement, vol. 38, no. 6, p. 1023, 2017. 17. [17].1. A. V. Clough and 2. C.-T. Chen M. Huelsbusch and V. Blazek, “Contactless mapping of rhythmical phenomena in tissue perfusion using PPGI,” A. V. Clough and C.-T. Chen, Eds., Apr. 2002, p. 110. doi: 10.1117/12.463573. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1117/12.463573&link_type=DOI) 18. [18]. J. E. Thatcher et al., “Multispectral and photoplethysmography optical imaging techniques identify important tissue characteristics in an animal model of tangential burn excision,” Journal of Burn Care & Research, vol. 37, no. 1, pp. 38–52, 2016. 19. [19]. J. Kuo, S. Koppel, J. L. Charlton, and C. M. Rudin-Brown, “Evaluation of a video-based measure of driver heart rate,” Journal of Safety Research, vol. 54, p. 55.e29-59, Sep. 2015, doi: 10.1016/j.jsr.2015.06.009. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jsr.2015.06.009&link_type=DOI) 20. [20]. C. S. Pilz, S. Zaunseder, U. Canzler, and J. Krajewski, “Heart rate from face videos under realistic conditions for advanced driver monitoring,” Current Directions in Biomedical Engineering, vol. 3, no. 2, pp. 483–487, Sep. 2017, doi: 10.1515/cdbme-2017-0102. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1515/cdbme-2017-0102&link_type=DOI) 21. [21]. S. Leonhardt, L. Leicht, and D. Teichmann, “Unobtrusive Vital Sign Monitoring in Automotive Environments—A Review,” Sensors, vol. 18, no. 9, p. 3080, 2018. 22. [22]. M. Koeny et al., “Using photoplethysmography imaging for objective contactless pain assessment,” Acta Polytechnica, vol. 54, no. 4, pp. 275–280, 2014. 23. [23]. M. Hülsbusch, “An image-based functional method for opto-electronic detection of skinperfusion,” RWTH Aachen (in German), 2008. 24. [24]. W. Wang, A. C. den Brinker, S. Stuijk, and G. de Haan, “Algorithmic Principles of Remote PPG,” IEEE Trans. Biomed. Eng., vol. 64, no. 7, pp. 1479–1491, Jul. 2017, doi: 10.1109/TBME.2016.2609282. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TBME.2016.2609282&link_type=DOI) 25. [25]. M. Lewandowska and J. Nowak, “Measuring Pulse Rate with a Webcam,” j med imaging hlth inform, vol. 2, no. 1, pp. 87–92, Mar. 2012, doi: 10.1166/jmihi.2012.1064. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1166/jmihi.2012.1064&link_type=DOI) 26. [26]. M.-Z. Poh, D. J. McDuff, and R. W. Picard, “Advancements in noncontact, multiparameter physiological measurements using a webcam,” IEEE transactions on biomedical engineering, vol. 58, no. 1, pp. 7–11, 2010. 27. [27]. G. de Haan and V. Jeanne, “Robust Pulse Rate From Chrominance-Based rPPG,” IEEE Trans. Biomed. Eng., vol. 60, no. 10, pp. 2878–2886, Oct. 2013, doi: 10.1109/TBME.2013.2266196. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TBME.2013.2266196&link_type=DOI) 28. [28]. G. De Haan and A. Van Leest, “Improved motion robustness of remote-PPG by using the blood volume pulse signature,” Physiological measurement, vol. 35, no. 9, p. 1913, 2014. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F09%2F2024.03.07.24303826.atom) 29. [29]. W. Wang, S. Stuijk, and G. de Haan, “A Novel Algorithm for Remote Photoplethysmography: Spatial Subspace Rotation,” IEEE Trans. Biomed. Eng., vol. 63, no. 9, pp. 1974–1984, Sep. 2016, doi: 10.1109/TBME.2015.2508602. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TBME.2015.2508602&link_type=DOI) 30. [30]. H.-T. Wu, G. Lewis, M. Davila, I. Daubechies, and S. Porges, “Optimizing Estimates of Instantaneous Heart Rate from Pulse Wave Signals with the Synchrosqueezing Transform,” Methods Inf Med, vol. 55, no. 05, pp. 463–472, 2016, doi: 10.3414/ME16-01-0026. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3414/ME16-01-0026&link_type=DOI) 31. [31]. C. Á. Casado and M. B. López, “Face2PPG: An Unsupervised Pipeline for Blood Volume Pulse Extraction From Faces,” IEEE J. Biomed. Health Inform., vol. 27, no. 11, pp. 5530–5541, Nov. 2023, doi: 10.1109/JBHI.2023.3307942. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/JBHI.2023.3307942&link_type=DOI) 32. [32]. L. C. Lampier, C. T. Valadão, L. A. Silva, D. Delisle-Rodríguez, E. M. D. O. Caldeira, and T. F. Bastos-Filho, “A deep learning approach to estimate pulse rate by remote photoplethysmography,” Physiol. Meas., vol. 43, no. 7, p. 075012, Jul. 2022, doi: 10.1088/1361-6579/ac7b0b. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/1361-6579/ac7b0b&link_type=DOI) 33. [33]. M. Rong and K. Li, “A Blood Pressure Prediction Method Based on Imaging Photoplethysmography in combination with Machine Learning,” Biomedical Signal Processing and Control, vol. 64, p. 102328, Feb. 2021, doi: 10.1016/j.bspc.2020.102328. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bspc.2020.102328&link_type=DOI) 34. [34]. S. Borik, P. Procka, J. Kubicek, and C. H. Antink, “Skin tissue perfusion mapping triggered by an audio-(de) modulated reference signal,” Biomedical Optics Express, vol. 13, no. 7, pp. 4058–4070, 2022. 35. [35]. A. A. Kamshilin, S. Miridonov, V. Teplov, R. Saarenheimo, and E. Nippolainen, “Photoplethysmographic imaging of high spatial resolution,” Biomedical optics express, vol. 2, no. 4, pp. 996–1006, 2011. 36. [36]. A. Singer and H. -T. Wu, “Vector diffusion maps and the connection Laplacian,” Comm Pure Appl Math, vol. 65, no. 8, pp. 1067–1144, Aug. 2012, doi: 10.1002/cpa.21395. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cpa.21395&link_type=DOI) 37. [37]. N. El Karoui and H.-T. Wu, “Graph connection Laplacian methods can be made robust to noise,” Ann. Statist., vol. 44, no. 1, Feb. 2016, doi: 10.1214/14-AOS1275. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1214/14-AOS1275&link_type=DOI) 38. [38]. A. Alian, K. Shelley, and H.-T. Wu, “Amplitude and phase measurements from harmonic analysis may lead to new physiologic insights: lower body negative pressure photoplethysmographic waveforms as an example,” J Clin Monit Comput, vol. 37, no. 1, pp. 127–137, Feb. 2023, doi: 10.1007/s10877-022-00866-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10877-022-00866-6&link_type=DOI) 39. [39]. S. Borik, S. Lyra, M. Paul, C. H. Antink, S. Leonhardt, and V. Blazek, “Photoplethysmography imaging: camera performance evaluation by means of an optoelectronic skin perfusion phantom,” Physiological measurement, vol. 41, no. 5, p. 054001, 2020. 40. [40]. Jihyoung Lee, K. Matsumura, K. Yamakoshi, P. Rolfe, S. Tanaka, and T. Yamakoshi, “Comparison between red, green and blue light reflection photoplethysmography for heart rate monitoring during motion,” in 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka: IEEE, Jul. 2013, pp. 1724–1727. doi: 10.1109/EMBC.2013.6609852. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/EMBC.2013.6609852&link_type=DOI) 41. [41]. Y. Maeda, M. Sekine, T. Tamura, A. Moriya, T. Suzuki, and K. Kameyama, “Comparison of reflected green light and infrared photoplethysmography,” in 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Vancouver, BC: IEEE, Aug. 2008, pp. 2270–2272. doi: 10.1109/IEMBS.2008.4649649. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/IEMBS.2008.4649649&link_type=DOI) 42. [42]. A. Alian, Y.-L. Lo, K. Shelley, and H.-T. Wu, “Reconsider phase reconstruction in signals with dynamic periodicity from the modern signal processing perspective,” FoDS, vol. 4, no. 3, p. 355, 2022, doi: 10.3934/fods.2022010. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3934/fods.2022010&link_type=DOI) 43. [43]. A. H. Nuttall and E. Bedrosian, “On the quadrature approximation to the Hilbert transform of modulated signals,” Proc. IEEE, vol. 54, no. 10, pp. 1458–1459, 1966, doi: 10.1109/PROC.1966.5138. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/PROC.1966.5138&link_type=DOI) 44. [44]. V. Huroyan, G. Lerman, and H.-T. Wu, “Solving Jigsaw Puzzles by the Graph Connection Laplacian,” SIAM J. Imaging Sci., vol. 13, no. 4, pp. 1717–1753, Jan. 2020, doi: 10.1137/19M1290760. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/19M1290760&link_type=DOI) 45. [45]. X. Xue, Y. Dong, and Y. Pan, “A Synchronous Detection Algorithm for Quasi-Periodic Signal Components and Its Application in Photoplethysmographic Imaging,” IEEE Trans. Instrum. Meas., vol. 71, pp. 1–13, 2022, doi: 10.1109/TIM.2022.3214492. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TIM.2022.3214492&link_type=DOI) 46. [46]. A. V. Moço, S. Stuijk, and G. de Haan, “New insights into the origin of remote PPG signals in visible light and infrared,” Sci Rep, vol. 8, no. 1, p. 8501, Dec. 2018, doi: 10.1038/s41598-018-26068-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-26068-2&link_type=DOI) 47. [47]. A. V. Moco, S. Stuijk, and G. De Haan, “Ballistocardiographic Artifacts in PPG Imaging,” IEEE Trans. Biomed. Eng., vol. 63, no. 9, pp. 1804–1811, Sep. 2016, doi: 10.1109/TBME.2015.2502398. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TBME.2015.2502398&link_type=DOI) 48. [48]. A. V. Moço, S. Stuijk, and G. De Haan, “Motion robust PPG-imaging through color channel mapping,” Biomed. Opt. Express, vol. 7, no. 5, p. 1737, May 2016, doi: 10.1364/BOE.7.001737. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1364/BOE.7.001737&link_type=DOI) 49. [49]. N. P. Podolyan et al., “Imaging photoplethysmography quantifies endothelial dysfunction in patients with risk factors for cardiovascular complications,” Biomedical Signal Processing and Control, vol. 86, p. 105168, Sep. 2023, doi: 10.1016/j.bspc.2023.105168. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bspc.2023.105168&link_type=DOI) 50. [50]. D. McDuff et al., “Non-contact imaging of peripheral hemodynamics during cognitive and psychological stressors,” Scientific Reports, vol. 10, no. 1, pp. 1–13, 2020. 51. [51].1. V. Blazek, 2. J. Kumar 3. V. S. Leonhardt, and 4. M. Mukunda Rao C. Blazek and V. Blazek, “Selected Clinical Applications of Functional PPGI Perfusion Mapping in Dermatology,” in Studies in Skin Perfusion Dynamics, V. Blazek, J. Kumar V. S. Leonhardt, and M. Mukunda Rao, Eds., in Biological and Medical Physics, Biomedical Engineering., Singapore: Springer Singapore, 2021, pp. 207–221. doi: 10.1007/978-981-15-5449-0_12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-981-15-5449-0_12&link_type=DOI) 52. [52]. S. Leonhardt, “Concluding remarks and new horizons in skin perfusion studies,” in Studies in Skin Perfusion Dynamics, Springer, 2021, pp. 223–232. 53. [53]. W. Wang and A. C. den Brinker, “Algorithmic insights of camera-based respiratory motion extraction,” Physiol. Meas., vol. 43, no. 7, p. 075004, Jul. 2022, doi: 10.1088/1361-6579/ac5b49. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/1361-6579/ac5b49&link_type=DOI) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: /embed/inline-graphic-6.gif [7]: /embed/inline-graphic-7.gif [8]: /embed/graphic-3.gif [9]: /embed/inline-graphic-8.gif [10]: /embed/graphic-4.gif [11]: /embed/inline-graphic-9.gif [12]: /embed/inline-graphic-10.gif [13]: /embed/inline-graphic-11.gif [14]: /embed/inline-graphic-12.gif [15]: /embed/inline-graphic-13.gif [16]: /embed/graphic-5.gif [17]: /embed/inline-graphic-14.gif [18]: /embed/inline-graphic-15.gif [19]: /embed/graphic-6.gif [20]: /embed/inline-graphic-16.gif [21]: /embed/inline-graphic-17.gif [22]: /embed/inline-graphic-18.gif [23]: /embed/inline-graphic-19.gif [24]: /embed/inline-graphic-20.gif [25]: /embed/inline-graphic-21.gif [26]: /embed/inline-graphic-22.gif [27]: /embed/inline-graphic-23.gif [28]: /embed/inline-graphic-24.gif [29]: /embed/inline-graphic-25.gif [30]: /embed/inline-graphic-26.gif [31]: /embed/inline-graphic-27.gif [32]: /embed/inline-graphic-28.gif [33]: /embed/inline-graphic-29.gif