Research Article  Open Access
RongChao Peng, XiaoLin Zhou, WanHua Lin, YuanTing Zhang, "Extraction of Heart Rate Variability from Smartphone Photoplethysmograms", Computational and Mathematical Methods in Medicine, vol. 2015, Article ID 516826, 11 pages, 2015. https://doi.org/10.1155/2015/516826
Extraction of Heart Rate Variability from Smartphone Photoplethysmograms
Abstract
Heart rate variability (HRV) is a useful clinical tool for autonomic function assessment and cardiovascular diseases diagnosis. It is traditionally calculated from a dedicated medical electrocardiograph (ECG). In this paper, we demonstrate that HRV can also be extracted from photoplethysmograms (PPG) obtained by the camera of a smartphone. Sixteen HRV parameters, including timedomain, frequencydomain, and nonlinear parameters, were calculated from PPG captured by a smartphone for 30 healthy subjects and were compared with those derived from ECG. The statistical results showed that 14 parameters (AVNN, SDNN, CV, RMSSD, SDSD, TP, VLF, LF, HF, LF/HF, nLF, nHF, SD1, and SD2) from PPG were highly correlated (, ) with those from ECG, and 7 parameters (AVNN, TP, VLF, LF, HF, nLF, and nHF) from PPG were in good agreement with those from ECG within the acceptable limits. In addition, five different algorithms to detect the characteristic points of PPG wave were also investigated: peak point (PP), valley point (VP), maximum first derivative (M1D), maximum second derivative (M2D), and tangent intersection (TI). The results showed that M2D and TI algorithms had the best performance. These results suggest that the smartphone might be used for HRV measurement.
1. Introduction
The heart rate (HR) of human is not constant but varies from one heartbeat to the next. Heart rate variability (HRV) is the physiological phenomenon of tiny fluctuations in the time intervals between heartbeats. It reflects the tenseness and the balance of the sympathetic and the vagus nerve activities and their effects on cardiovascular motion [1]. It is a noninvasive method for assessing the autonomic functions [1]. Numerous publications have validated that abnormal changes of HRV are related to several cardiological and noncardiological diseases like myocardial infarction, diabetic neuropathy, myocardial dysfunction, and tetraplegia [2]. Therefore, HRV is of significant importance and is widely used in clinical application.
HRV is traditionally determined by digital processing of electrocardiograms (ECG). The Rwave peaks of QRS complexes in ECGs are detected by computer algorithms and RtoR intervals (RRI) are calculated. Then, HRV parameters are computed using timedomain, frequencydomain, and nonlinear methods [2]. However, the traditional measurement of ECG has several limitations: (1) ECG instruments generally require three electrodes attached to specific anatomical positions, which limit the subjects’ activities and make them uncomfortable; (2) the electrodes may cause skin irritations for some special subjects with allergies; (3) the ECG instruments are usually operated by specially trained nurses in the hospital and are not suitable for daily use at home. Therefore, new technologies have been developed to measure HRV without ECG, such as photoplethysmography (PPG) [3–5], Finapres (continuous blood pressure monitoring) [6, 7], impedance plethysmography [8, 9], ballistocardiography [10], optical vibrocardiography [11], a microwave sensor [12], or a webcam [13, 14]. In this paper, we demonstrate that HRV can also be extracted from a smartphone.
Mobile phones have already shown promising applications in healthcare service [15]. As the new generation smartphones are becoming more powerful and more popular and with more builtin sensors, the smartphonebased healthcare applications are being rapidly developed [16]. In recent years, a new method was proposed to acquire PPG signals from the builtin camera of a smartphone, in which one only needs to press a finger on the camera lens and capture a video record with the builtin LED flash on [17]. This method is based on the principle that the intensity changes of video frames are associated with the variations of light absorption of blood. When the heart systoles and the capillaries in the fingertip are full of blood, more light is absorbed and the frame becomes darker; likewise, when the heart diastoles and the capillaries in the fingertip are full of less blood, less light is absorbed and the frame becomes brighter. PPG signals can thus be obtained by calculating the intensity changes of the frames, and then several physiological parameters such as HR [17–21], respiratory rate [19], pulse volume [21], and oxygen saturation [19] can be extracted from the PPG signals. With the addition of a microphone to detect the phonocardiogram signal, blood pressure can also be estimated [22]. This method is simple, lowcost, and easytouse, with great potential to be used in the healthcare service in the future.
However, to our best knowledge, the extraction of HRV from the smartphone PPG signals has not been well investigated, especially compared with ECG—the golden standard. Therefore, we comprehensively studied the extraction of HRV from smartphone PPG signals and compared the results with an ECG in order to assess the accuracy. Specifically, we used five algorithms to detect the characteristic points of the smartphone PPG signals: peak point (PP), valley point (VP), maximum first derivative (M1D), maximum second derivative (M2D), and tangent intersection (TI). The performances of these algorithms were also compared.
2. Methods
2.1. Data Acquisition
The experiment was approved by the Institutional Review Board of Shenzhen Institutes of Advanced Technology (registration number: SIATIRB140215H0040). Thirty subjects participated in the experiment (20 males and 10 females, age: 20–32 years, height: 150–183 cm, and weight: 40–90 kg). All the subjects were healthy and provided their informed consent. They were asked to refrain from caffeine, alcohol, cigarette, or strenuous exercise for 2 hours prior to the study.
In the experiment, all the subjects were instructed to lie in the supine position on a mattress and place their right index finger on the camera lens of an HTC S510e smart phone with the builtin LED flash turned on. A camera application (APP) in the smart phone was used to record the video of the fingertip with a resolution of 320 × 240 pixels at an unfixed sampling rate of 20–30 frames per second (fps). The sampling rate is unfixed due to the CPU processing load. ECG electrodes in the standard configuration were attached to the subjects to measure the ECG signals with a Finometer MIDI (Model II, Finapres Medical Systems B.V., The Netherlands). The ECG signals were digitalized at 200 Hz and automatically stored in the computer by BeatScope Easy software (Finapres Medical Systems B.V., The Netherlands). The experiment lasted at least 5 minutes for each subject and the subject was asked to keep still during this period.
2.2. Smartphone PPG Processing
All the data were processed offline. The 3GP format videos recorded by the HTC S510e smart phone were converted into AVI format using Pazera Free 3GP to AVI Converter 1.3 (). All further analysis was performed on the AVI videos in MatLab 7.0 (The Mathworks Inc., USA).
First, an 80 × 80 pixel region in the center of the video image was selected as the region of interesting (ROI). Then, the average intensity of the red channel in the ROI for each individual frame was calculated to generate a timeseries waveform (the raw PPG signal). The red channel was chosen because the intensity values of the green and blue channels were often tending to zero and contained no valuable information in most situations. As the smartphone PPG worked in the reflection mode, the generated waveform should be inverted to “normal” mode for further processing [23].
The raw PPG signals were often corrupted by random noise, baseline drifting, and baseline abrupt changes (increase/decrease). Baseline abrupt changes were possibly caused by sudden moves of the finger or sudden changes of the light illumination, or by other unknown reasons. They could not be completely removed by general digital filters as they contained wideband frequency components. We used a statistics method to solve this problem. First, we calculated the difference of the raw signal and then removed the outliers out of the range mean ± 5 × standard deviation (SD) and interpolated new values using cubic spline interpolation. At last, we reconstructed the new PPG signal by summation, the inverse of the difference. The range mentioned above was determined empirically, which meant that the probability of the outliers was 5.7330 × 10^{−7} if the difference of the PPG was normally distributed. It was the best range according to our data and could be adjusted if required.
The random noise and baseline drifting were reduced by a zerophase Butterworth lowpass filter with cutoff frequency of 10 Hz and a zerophase Butterworth highpass filter with cutoff frequency of 0.5 Hz, respectively. Zero phase filters were implemented by filtering the signal both forward and backward to eliminate phase distortion.
The PPG signals were then resampled to 800 Hz with cubic spline interpolation to increase the temporal resolution. For each cardiac circle, five algorithms were used to obtain the pulsetopulse interval (PPI) by detection of five different characteristic points, as illustrated in Figure 2.
(i) Maximum First Derivative (M1D). It is the location of the maximum value of the first derivative of the PPG signal [24]. The derivative is calculated by a fivepoint central difference equation (1).
(ii) Peak Point (PP). It is the location of the maximum amplitude in the PPG signal following the M1D point [24].
(iii) Valley Point (VP). It is the location of the minimum amplitude in the PPG signal preceding the M1D point [24].
(iv) Maximum Second Derivative (M2D). It is the location of the maximum value of the second derivative of the PPG signal [25]. The second derivative is calculated using a fivepoint central difference equation (1) and a subsequent sevenpoint central difference equation (2). Considerwhere is the input signal and is the output signal.
(v) Tangent Intersection (TI). It is the location of the intersection of the tangent line at the M1D point and the horizontal line passing the valley point [25]. The tangent line is the fitted line of five points centered at the M1D point.
2.3. ECG Processing
The ECG signals were first passed through a finite impulse response (FIR) lowpass filter with cutoff frequency of 11 Hz and then a FIR highpass filter with cutoff frequency of 5 Hz to reduce most of the noise and interference [26]. Thereafter, they were resampled to 800 Hz with cubic spline interpolation to increase the temporal resolution. Rwave peak detection was performed using Pan and Tompkins’ algorithm [26] and RRIs were obtained as the difference of successive Rwave peak locations. For both PPG and ECG signals, missed beats and false beats were manually identified and adjusted. An example of the obtained RRI and PPIs is shown in Figure 3.
2.4. HRV Parameters Calculation
Table 1 lists some commonly used HRV parameters. These parameters often refer to a professional term, NN intervals (“normaltonormal intervals”), which means that only regular heartbeats should be considered [2]. Therefore, ectopic beats of the RRI/PPI series were removed and replaced by cubic spline interpolation before HRV parameters calculation.

(i) TimeDomain Parameters. Seven parameters were calculated from RRI and PPI series in time domain: AVNN, CV, SDNN, RMSSD, SDSD, NN50, and pNN50, as described in Table 1.
(ii) FrequencyDomain Parameters. The RRI and PPI series were evenly resampled at 4 Hz using cubic spline interpolation, and the DC component was removed by subtracting the mean of the series. Then, a 16th order autoregressive (AR) model was employed to estimate the power spectral density. Seven parameters were calculated: TP, VLF, LF, HF, LF/HF, nLF, and nHF, as described in Table 1. The ranges of different frequency bands were in accordance with the standard definition [2].
(iii) Poincaré Parameters. The Poincaré plot is one of the most widely used methods for nonlinear HRV analysis. It is a plot of each RRI/PPI against its previous one. Two parameters were calculated from the Poincaré plot, SD1 and SD2, as described in Table 1.
2.5. Statistics Analysis
HRV parameters derived from smartphone PPG were compared with the corresponding parameters derived from ECG. The Pearson correlation coefficients were calculated and the linear regression equations were obtained. A value < 0.05 was considered statistically significant.
The agreement between the two devices (smartphone and ECG) was assessed using BlandAltman method [27]. The limit of agreement (LOA) was defined as bias ± 1.96 × SD (3)–(5) [27] and a BlandAltman ratio (BAR) was defined as the ratio of half the range of limits of agreement to the mean of the pairwise measurement means (6). Agreements were ranked as good (BAR < 10%), moderate (10% BAR < 20%), or insufficient (BAR 20%) [28]. Acceptable limit (AL) of agreement was defined as 20% of the mean of the pairwise measurement means (7), as there are wide interindividual variations for HRV measurement [2] and a limit greater than 20% is generally considered unacceptable. Consider where and are the HRV parameters derived from the ECG and the smartphone, respectively, and is the number of subjects.
3. Results
3.1. Agreement Analysis
Table 2 shows the Pearson correlation coefficients and linear regression equations between HRV parameters derived from the smartphone and the ECG. It was found that the correlation coefficients were >0.6 for all parameters except NN50 and pNN50. For timedomain parameters, the correlation coefficients of AVNN, CV, and SDNN were higher than those of RMSSD and SDSD. For frequencydomain parameters, all parameters showed strong correlations (, ). For nonlinear parameters, SD2 exhibited higher correlation (, ) than SD1 (, ). Nevertheless, a good correlation does not mean a good agreement as either adding a constant or multiplying a factor will still yield a good correlation. A good method for assessing agreement is BlandAltman analysis.
 
: HRV parameters derived from an electrocardiograph and : HRV parameters derived from a smartphone using five different algorithms. PP: peak point; VP: valley point; M1D: maximum first derivative; M2D: maximum second derivative, and TI: tangent intersection. A value <0.05 was considered statistically significant. HRV parameters are explained in Table 1. 
Table 3 shows the BlandAltman analysis of HRV parameters derived from the smartphone and the ECG. For the sake of simplicity, we speak of good/moderate agreement if three of the five algorithms were in good/moderate agreement. It was found that all the timedomain parameters showed insufficient agreements (BAR 20%), but the AVNN showed excellent agreement (BAR < 1%), indicating that smartphonederived HR can be a surrogate of ECGderived HR. This result was in line with Gregoski et al.’s [20] and Matsumura and Yamakoshi’s [21]. It was also found that all the frequencydomain parameters were in moderate agreement (BAR < 20%) except for TP and nLF which were in good agreement (BAR < 10%). TP, VLF, HF, and nHF were overestimated (bias > 0), while LF, LF/HF, and nLF were underestimated (bias < 0), implying that the smartphonederived HRV contains more noise, which can be observed in detail in Figure 1. For nonlinear parameters, SD2 showed good agreement (BAR < 10%) and SD1 showed insufficient agreement (BAR 20%).
 
Data are presented as bias ± 1.96 standard deviation (SD). ^{*}Bias ± 1.96 SD within the acceptable limits. BAR: BlandAltman ratio, PP: peak point, VP: valley point, M1D: maximum first derivative, M2D: maximum second derivative, and TI: tangent intersection. HRV parameters are explained in Table 1. 
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(f)
As shown in Table 3, a total number of 7 parameters (AVNN, TP, VLF, LF, HF, nLF, and nHF) were within the acceptable limits. Figure 4 shows the BlandAltman plots of different frequency components of HRV which are commonly used for assessing the autonomic functions. It was found that the limits of agreement for LF, HF, nLF, and nHF were all within their corresponding acceptable limits, meaning that the discrepancies between the smartphone PPG and the ECG for LF, HF, nLF, and nHF were not considerable. It was also found that the lower limit of agreement for LF/HF is out of the range of acceptable limits.
3.2. Algorithms Comparison
In terms of both bias and SD, we analyzed these data satisfying the condition BAR < 20% in Table 3 to evaluate the performance of the five algorithms mentioned above. For each HRV parameter, the best two algorithms with the least bias or SD would gain a star “^{*}” each. The overall performance was graded according to the number of total stars. As shown in Table 4, the M2D and the TI algorithms were better than the others, the PP and the VP were the worst, and the M1D was in the middle.
 
PP: peak point, VP: valley point, M1D: maximum first derivative, M2D: maximum second derivative, and TI: tangent intersection. SD: standard deviation. HRV parameters are explained in Table 1. 
4. Discussion
As previously mentioned, no researches have been reported to measure HRV with smartphone PPG. Nevertheless, many researches have been reported to measure HRV with traditional PPG (tPPG, i.e., a pulse oximeter), which can provide valuable information for our work. A good review of tPPGderived HRV can be found in [28], where the authors commented that it was controversial whether the tPPGderived HRV was a surrogate of the ECGderived HRV. Although a number of publications reported universally good agreement between tPPGderived HRV and ECGderived HRV for all parameters, many other studies found that shortterm parameters were more susceptible to disagreement between tPPG and ECG than longterm parameters. Our results are in accordance with the latter. As shown in Table 2, the coefficients of longterm parameters (AVNN, SDNN, and CV) are higher than the shortterm parameters (RMSSD, SDSD, NN50, and pNN50). As shown in Table 3, the lowfrequency parameters are in higher agreement than the highfrequency parameters, as the bias of LF is less than that of the HF and their BARs have no significant difference. The nonlinear parameter SD2 presenting the level of longterm HRV is in good agreement between ECG and smartphone PPG while SD1 presenting the level of shortterm HRV is in insufficient agreement.
The detection of the characteristic points is also an impact factor for the accuracy of smartphonederived HRV measurement [4]. As the morphology of pulse wave changes over time, an algorithm less sensitive to morphology variation will produce better accuracy. We graded these five algorithms: PP, VP, M1D, M2D, and TI. Our results indicate that M2D and TI algorithms are better than the others, which are in accordance with Chiu et al.’s [25], but they are slightly different from PosadaQuintero et al.’s [29]. PosadaQuintero et al. stated that the TI algorithm was better than the VP and M2D algorithms. This slight difference is possibly due to different algorithm details and different evaluation criteria. Another interesting research found that PP was more sensitive to waveform distortion than VP and M1D when the peripheral pulse was affected by cold temperature [24]. Overall, M2D and TI algorithms are better than PP, VP, and M1D. The fact that many researchers used the PP algorithm is worthy of attention [3, 4, 30, 31]. In addition, none of the algorithms is perfect or errorfree. A manual correction is usually required by visual inspection on the computer screen, which is timeconsuming. More efficiency algorithms are needed.
A possible consideration of the smartphonebased HRV analysis is the sampling rate. Our HTC smartphone has a sampling rate of approximately 20–30 Hz that may be considered not suitable for HRV analysis. In fact, the spectrum of the pulse signals has the vast majority power in the range of 0~10 Hz [32]. A sampling rate of 20 Hz is not less than the Nyquist rate and the temporal resolution can be improved by interpolation, which was confirmed by Sun et al.’s experiments in [14] where they compared different sampling rates of 200, 100, 50, and 20 fps and their results revealed no significant differences of these sampling rates. Moreover, as the new generation smartphones have more powerful CPUs and higher speed cameras, the sampling rate will be improved for better performance.
The color channel should be also considered. In the processing of the recorded video, most previous studies calculated the intensity of the green channel in RGB color model [17–19], because (oxy)haemoglobin absorbs more green light than red light and penetrates sufficiently deep into the skin as compared to blue light [33]. We chose the red channel over the green or blue channel, because we observed that the pixel values in green and blue channels were tending to zero and the changes of the red channel are more pronounced than the green and blue channels in most situations. Chandrasekaran et al. also observed that the prominent color was red and they used the red channel in their research [22]. Grimaldi et al. demonstrated that the distribution of the pixels in the green channel was not uniform for different models of the smartphones and the only channel that had similar characteristics was the red one [34]. This suggests that heterogeneous characteristics of different cameras in different smartphones should be taken into account and more robust algorithms are required. Some potential alternatives are principal component analysis (PCA), independent component analysis (ICA), or using other color models (e.g., CMY color model, HSI color model, and YUV color model).
Motion artifacts are another complicated problem and are tough to deal with. To our best knowledge, none of the reported studies have solved this problem very well. In our experiments, the subjects were instructed to lie on a mattress and keep their fingers as still as possible to minimize the motion artifacts. This is not practical in daily life as shorttime HRV testing usually takes 5 minutes that seems so long time for keeping still. Therefore, efficient motionresistant algorithms are required. Several motion artifacts detection algorithms in pulse oximeters could be applied in smartphonebased HRV analysis [35, 36].
5. Conclusion
Traditional ECG recordings require electrodes attached to body surface and are operated by specially trained nurses in the hospital. The new smartphonebased technology requires no more than placing a finger on the camera lens of a smartphone. It is lowcost and easytouse and can be used in daily life out of hospital. In the present study, we quantitatively investigated the measurement of HRV based on smartphone technology and compared the results with those derived from a standard ECG to assess the accuracy. The results suggest that the smartphone can be of potential use for HRV measurement at resting and would be applied in lowcost healthcare applications.
Conflict of Interests
The authors declare that they have no conflict of interests.
Acknowledgments
This work was supported in part by the National Basic Research Program 973 (no. 2010CB732606), the National Natural Science Foundation of China (no. 61401453), the Guangdong Innovation Research Team Fund for Lowcost Healthcare Technologies in China, the External Cooperation Program of the Chinese Academy of Sciences (no. GJHZ1212), the Key Lab for Health Informatics of Chinese Academy of Sciences, the Peacock Program to Attract Overseas HighCaliber Talents to Shenzhen, and the Shenzhen City Government (nos. KFJEWSTS097 and KFJEWSTS095).
References
 B. Pomeranz, R. J. Macaulay, M. A. Caudill et al., “Assessment of autonomic function in humans by heart rate spectral analysis,” The American Journal of Physiology, vol. 248, no. 1, pp. H151–H153, 1985. View at: Google Scholar
 A. J. Camm, M. Malik, J. T. Bigger et al., “Heart rate variability: standards of measurement, physiological interpretation, and clinical use,” Circulation, vol. 93, pp. 1043–1065, 1996. View at: Google Scholar
 K. Charlot, J. Cornolo, J. V. Brugniaux, J. P. Richalet, and A. Pichon, “Interchangeability between heart rate and photoplethysmography variabilities during sympathetic stimulations,” Physiological Measurement, vol. 30, no. 12, pp. 1357–1369, 2009. View at: Publisher Site  Google Scholar
 E. Gil, M. Orini, R. Bailón, J. M. Vergara, L. Mainardi, and P. Laguna, “Photoplethysmography pulse rate variability as a surrogate measurement of heart rate variability during nonstationary conditions,” Physiological Measurement, vol. 31, no. 9, pp. 1271–1290, 2010. View at: Publisher Site  Google Scholar
 S. Lu, H. Zhao, K. Ju et al., “Can photoplethysmography variability serve as an alternative approach to obtain heart rate variability information?” Journal of Clinical Monitoring and Computing, vol. 22, no. 1, pp. 23–29, 2008. View at: Publisher Site  Google Scholar
 P. S. McKinley, P. A. Shapiro, E. Bagiella et al., “Deriving heart period variability from blood pressure waveforms,” Journal of Applied Physiology, vol. 95, no. 4, pp. 1431–1438, 2003. View at: Google Scholar
 A. Suhrbier, R. Heringer, T. Walther, H. Malberg, and N. Wessel, “Comparison of three methods for beattobeatinterval extraction from continuous blood pressure and electrocardiogram with respect to heart rate variability analysis,” Biomedizinische Technik, vol. 51, no. 2, pp. 70–76, 2006. View at: Publisher Site  Google Scholar
 N. K. Kristiansen, J. Fleischer, M. S. Jensen, K. S. Andersen, and H. Nygaard, “Design and evaluation of a handheld impedance plethysmograph for measuring heart rate variability,” Medical and Biological Engineering and Computing, vol. 43, no. 4, pp. 516–521, 2005. View at: Publisher Site  Google Scholar
 E. F. Treo, M. C. Herrera, and M. E. Valentinuzzi, “Algorithm for identifying and separating beats from arterial pulse records,” BioMedical Engineering OnLine, vol. 4, article 48, 2005. View at: Publisher Site  Google Scholar
 J. H. Shin, S. H. Hwang, M. H. Chang, and K. S. Park, “Heart rate variability analysis using a ballistocardiogram during Valsalva manoeuvre and post exercise,” Physiological Measurement, vol. 32, no. 8, pp. 1239–1264, 2011. View at: Publisher Site  Google Scholar
 U. Morbiducci, L. Scalise, M. De Melis, and M. Grigioni, “Optical vibrocardiography: a novel tool for the optical monitoring of cardiac activity,” Annals of Biomedical Engineering, vol. 35, no. 1, pp. 45–58, 2007. View at: Publisher Site  Google Scholar
 G. Lu, F. Yang, Y. Tian, X. Jing, and J. Wang, “Contactfree measurement of heart rate variability via a microwave sensor,” Sensors, vol. 9, no. 12, pp. 9572–9581, 2009. View at: Publisher Site  Google Scholar
 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, 2011. View at: Publisher Site  Google Scholar
 Y. Sun, S. Hu, V. AzorinPeris, R. Kalawsky, and S. Greenwald, “Noncontact imaging photoplethysmography to effectively access pulse rate variability,” Journal of Biomedical Optics, vol. 18, no. 6, Article ID 61205, 2013. View at: Publisher Site  Google Scholar
 R. S. H. Istepanian, E. Jovanov, and Y. T. Zhang, “Introduction to the special section on mHealth: beyond seamless mobility and global wireless healthcare connectivity,” IEEE Transactions on Information Technology in Biomedicine, vol. 8, no. 4, pp. 405–414, 2004. View at: Publisher Site  Google Scholar
 M. N. K. Boulos, S. Wheeler, C. Tavares, and R. Jones, “How smartphones are changing the face of mobile and participatory healthcare: an overview, with example from eCAALYX,” BioMedical Engineering OnLine, vol. 10, article 24, 2011. View at: Publisher Site  Google Scholar
 E. Jonathan and M. Leahy, “Investigating a smartphone imaging unit for photoplethysmography,” Physiological Measurement, vol. 31, no. 11, pp. N79–N83, 2010. View at: Publisher Site  Google Scholar
 E. Jonathan and M. J. Leahy, “Cellular phonebased photoplethysmographic imaging,” Journal of Biophotonics, vol. 4, no. 5, pp. 293–296, 2011. View at: Publisher Site  Google Scholar
 C. G. Scully, J. Lee, J. Meyer et al., “Physiological parameter monitoring from optical recordings with a mobile phone,” IEEE Transactions on Biomedical Engineering, vol. 59, no. 2, pp. 303–306, 2012. View at: Publisher Site  Google Scholar
 M. J. Gregoski, M. Mueller, A. Vertegel et al., “Development and validation of a smartphone heart rate acquisition application for health promotion and wellness telehealth applications,” International Journal of Telemedicine and Applications, vol. 2012, Article ID 696324, 7 pages, 2012. View at: Publisher Site  Google Scholar
 K. Matsumura and T. Yamakoshi, “iPhysioMeter: a new approach for measuring heart rate and normalized pulse volume using only a smartphone,” Behavior Research Methods, vol. 45, no. 4, pp. 1272–1278, 2013. View at: Publisher Site  Google Scholar
 V. Chandrasekaran, R. Dantu, S. Jonnada, S. Thiyagaraja, and K. P. Subbu, “Cuffless differential blood pressure estimation using smart phones,” IEEE Transactions on Biomedical Engineering, vol. 60, no. 4, pp. 1080–1089, 2013. View at: Publisher Site  Google Scholar
 J. Allen, “Photoplethysmography and its application in clinical physiological measurement,” Physiological Measurement, vol. 28, no. 3, pp. R1–R39, 2007. View at: Publisher Site  Google Scholar
 X.Y. Zhang and Y.T. Zhang, “The effect of local mild cold exposure on pulse transit time,” Physiological Measurement, vol. 27, no. 7, pp. 649–660, 2006. View at: Publisher Site  Google Scholar
 Y. C. Chiu, P. W. Arand, S. G. Shroff, T. Feldman, and J. D. Carroll, “Determination of pulse wave velocities with computerized algorithms,” American Heart Journal, vol. 121, no. 5, pp. 1460–1470, 1991. View at: Publisher Site  Google Scholar
 J. Pan and W. J. Tompkins, “A realtime QRS detection algorithm,” IEEE Transactions on Biomedical Engineering, vol. 32, no. 3, pp. 230–236, 1985. View at: Google Scholar
 J. M. Bland and D. G. Altman, “Statistical methods for assessing agreement between two methods of clinical measurement,” The Lancet, vol. 1, no. 8476, pp. 307–310, 1986. View at: Google Scholar
 A. Schäfer and J. Vagedes, “How accurate is pulse rate variability as an estimate of heart rate variability? A review on studies comparing photoplethysmographic technology with an electrocardiogram,” International Journal of Cardiology, vol. 166, no. 1, pp. 15–29, 2013. View at: Publisher Site  Google Scholar
 H. F. PosadaQuintero, D. DelisleRodríguez, M. B. CuadraSanz, and R. R. F. de la VaraPrieto, “Evaluation of pulse rate variability obtained by the pulse onsets of the photoplethysmographic signal,” Physiological Measurement, vol. 34, no. 2, pp. 179–187, 2013. View at: Publisher Site  Google Scholar
 J.S. Wong, W.A. Lu, K.T. Wu, M. Liu, G.Y. Chen, and C.D. Kuo, “A comparative study of pulse rate variability and heart rate variability in healthy subjects,” Journal of Clinical Monitoring and Computing, vol. 26, no. 2, pp. 107–114, 2012. View at: Publisher Site  Google Scholar
 J. B. Bolkhovsky, C. G. Scully, and K. H. Chon, “Statistical analysis of heart rate and heart rate variability monitoring through the use of smart phone cameras,” in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, pp. 1610–1613, 2012. View at: Google Scholar
 C. T. Lee and L. Y. Wei, “Spectrum analysis of human pulse,” IEEE Transactions on Biomedical Engineering, vol. 30, no. 6, pp. 348–352, 1983. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 D. Grimaldi, Y. Kurylyak, F. Lamonaca, and A. Nastro, “Photoplethysmography detection by smartphone's videocamera,” in Proceedings of the 6th IEEE International Conference on Intelligent Data Acquisition and Advanced Computing Systems: Technology and Applications (IDAACS '11), pp. 488–491, Prague, Czech Republic, September 2011. View at: Publisher Site  Google Scholar
 M. J. Hayes and P. R. Smith, “Artifact reduction in photoplethysmography,” Applied Optics, vol. 37, no. 31, pp. 7437–7446, 1998. View at: Publisher Site  Google Scholar
 B. S. Kim and S. K. Yoo, “Motion artifact reduction in photoplethysmography using independent component analysis,” IEEE Transactions on Biomedical Engineering, vol. 53, no. 3, pp. 566–568, 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 RongChao Peng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.