Evaluation of PPG Biometrics for Authentication in different states
Umang Yadav
Sherif N. Abbas
Dimitrios Hatzinakos
The Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto
10 Kings College Road, Toronto, ON, Canada, M5S 3G4.

Abstract

0.8

Amongst all medical biometric traits, Photoplethysmograph (PPG) is the easiest to acquire. PPG records the
blood volume change with just combination of Light Emitting Diode and Photodiode from any part of the body. With
IoT and smart homes’ penetration, PPG recording can easily be integrated with other vital wearable devices. PPG
represents peculiarity of hemodynamics and cardiovascular system for each individual. This paper presents nonfiducial method for PPG based biometric authentication.
Being a physiological signal, PPG signal alters with physical/mental stress and time. For robustness, these variations
cannot be ignored. While, most of the previous works focused only on single session, this paper demonstrates extensive performance evaluation of PPG biometrics against single session data, different emotions, physical exercise and
time-lapse using Continuous Wavelet Transform (CWT) and
Direct Linear Discriminant Analysis (DLDA). When evaluated on different states and datasets, equal error rate (EER)
of 0.5%-6% was achieved for 45-60s average training time.
Our CWT/DLDA based technique outperformed all other
dimensionality reduction techniques and previous work.

0.6

1. Introduction
Risk of impersonation and violation in security of any
system is ever lasting threat. This risk can cause financial ruins and, is also life threatening in today’s Internet of
Things (IoT) and Smart Home Era. Hence, need of low
cost and non-invasive authentication system has never been
greater. This is why in last two decades, biometrics research has emerged as very important area of research. Photoplethysmograph (PPG) is a non-invasive electro-optical
method which measures the volume of the blood flowing
through the body part under testing. It reflects the pulsative
actions of the arteries through the interaction of the oxygenized hemoglobin and photons. It is believed that every person has unique hemodynamics and cardiovascular system.
Since, PPG captures this unique characteristics, this paper

©

Systolic
Peak
Amplitude

arXiv:1712.08583v1 [cs.CR] 22 Dec 2017

{umang,sseha,dimitris}@ece.utoronto.ca

Diastolic Peak

0.4

0.2

Figure 1: PPG recording
with LED and PhotoDiode (PD)

Diastolic
Point

Inter
Pulse
Interval

Dicrotic
Notch

0

Time

Figure 2: Sample PPG Signal

presents method to leverage it for biometric authentication.
PPG being a biological signal, it is harder to steal or replicate. It has advantages of inherent anti-spoofing and liveness detection over traditional biometrics modalities. Also,
as shown in Fig.1, PPG can be recorded with just combination of LED and Photo-Diode (PD) from any part of the
body, which provides greater flexibility for systems design.
Since, PPG recording only requires LED and PD, it is very
cost effective compared to other biometric traits. In context
of medical biometrics, PPG recording doesn’t require any
kind of gel (EEG), external stimulus (TEOAE) or multiple
electrodes (ECG) and can be conveniently recorded from
virtually any part of body. All such factors give PPG biometrics extra edge on other medical biometrics traits.
Apart from authentication, PPG can be used for many
clinical applications such as measuring oxygen saturation,
blood pressure, detecting peripheral vascular diseases etc.
Wearable device companies such as FitbitTM are already using PPG signals for Heart Rate(HR) measurements. PPG
based biometric system can be integrated into such devices
and smart homes’ IoT. Thus, PPG has inherent advantage of
being portable and can give user unobtrusive experience.
In practice, PPG signal is generally impaired by many
common noise sources during acquisition such as motion
artifacts, sensor movements, respiration, premature ventricular contraction (PVC) and ambient light. In addition, since
PPG embodies proper functioning of heart, under physi-

2018 IEEE. Accepted at 11th IAPR/IEEE International conference on Biometrics, 2018.
This work has been supported by NSERC.

cal and emotional stress, it changes along with Heart Rate
(HR). Empirical comparison of PPG signal under rest condition and exercise shows considerable change in shape of
PPG signal and its spectrum. Although indiscernible, it is
also important to assess system performance under different
emotions as they influence functioning of autonomic nervous system and heart. Furthermore, like any other biometric traits, it is crucial to take into account long-term behavior
with time-lapse. In this paper, system performance is evaluated against all such variations using 3 different datasets.
For PPG based authentication, many fiducial and nonfiducial based approaches have been proposed in past. Fiducial based methods rely on detecting fiducial points of PPG
signals shown in Fig.2 such as systolic peak, diastolic peak,
dicrotic notch, inter-pulse interval, amplitudes of peaks, etc.
Given variability in PPG shape in different states, fiducial
detection on raw PPG signal might be unsuccessful or incorrect. Therefore, many researchers extended the idea to first
derivative (FD) and second derivatives (SD) of raw PPG signals and used similar points on FD and SD as features for
authentication. For example, Kavsaoglu et.al. extracted 40
features from raw PPG, its FD and SD with kNN showed
94.44% accuracy [6]. On non-fiducial side, Spachos et.al.
showed feasibility of PPG as biometrics identifier reporting
0.5% and 25% EER on two different datasets using temporal features and LDA [10]. In most recent study, Karimian
et.al. showed non-fiducial approach with DWT and kNN,
reported 1.31% mean EER for 42 subjects [4]. Sarkar et.al.
evaluated variations with emotions using DEAP dataset by
approximating PPG beat with sum of gaussians and showed
accuracy of 90.53% for cross-emotions evaluation [9].
Many of the previous methods have been focused on
fiducial approach. But fiducial points detection in every
condition is error-prone. Therefore, in this paper nonfiducial approach is adopted. Also, most of the prior work
has been focused on single session evaluation. For robustness, it is important to check variability in different states.
We present Continous Wavelet Transform (CWT) and Direct Linear Discriminant Analysis (DLDA) based method
for authentication. CWT reveals unique time-frequency behaviour of each individual while DLDA selects CWT coefficients that give maximum class separabillity. Such feature
selection is necessary for better results against variations.
Additionally, taking more pragmatic approach and training
of model was limited to 45 − 60 seconds of signal and system performance was evaluated for different length of authentication time.
Section 2 and 3 present methodology and detailed description of system comprising pre-processing, template
creation and template matching. Section 4 presents experimental results in different emotional states, exercise, rest
condition and under time lapse. Section 5 concludes the
paper with observations and future work.

Data
Acquisition

Filtering

Peak
Detection

False Peak
Removal

Segmentation

Pre-processing

CWT
Coefficients

Direct-LDA based
Dimensionality
reduction

Template Creation

Pearson s
Distance
Decision
System
Database
Template
Matching

Figure 3: System Flow Diagram

2. PPG Authentication System and
Methodology
System flow is similar to many biometrics authentication system. As shown in Fig.3 system consists of the three
main blocks; pre-processing, template creation and template matching.

2.1. Pre-processing
• Filtering: Butterworth IIR Band pass filter of cutoff frequencies 0.5-5 Hz and order 38 has been applied on whole
signal for the removal of the power line interference, motion artifacts, baseline wanders. After filtering amplitude
was normalized to have dynamic range of one.
• Peak Detection: Successful application of pattern recognition techniques depends on reliable peak detection. Unlike ECG, there is no standard proven algorithm for systolic peak detection in PPG. In practice, most of the local maxima finding algorithms work reasonably well with
PPG. To find peak locations, signal is first squared to emphasize amplitude difference between local maxima and
other points. Thereafter, peaks were detected based on
amplitude prominence.
• False Peak Removal: Many times motion artifacts, respiration, ambient light falls into passband of filter. In such
cases, filtering doesn’t help. Also, filtering degrades the
amplitude of PPG peaks and causes false peak detections.
To tackle this issue, simple false peak removal technique
based on distance was employed. For any healthy subject
based on his/her age and physical activity state if heart
rate (HR) found using distance between two peaks was
outside of normal range of HR, then those peaks were removed.
• Segmentation: For the features extraction, signal is segmented around detected peaks. Length of each segments
is set using median HR found from signal. For given peak
location i, segment can be represented as follows:
S(j, 1.....2r + 1) = SIG(i − 2r....i + 2r)

(1)

where r is median HR calculated from signal, j is peak
number, SIG is signal and S is segment. Since, each pulse

segment covers, 3.5 to 4 pulses, it is larger in length and
gives better frequency resolution with CWT. Hence, this
choice of length was found to be producing good results.
Furthermore, each two consecutive segments in S were
averaged. This was necessary to improve SNR, damping
any unwanted random noise as segment size was larger.

2.2. Template Creation
• Features Extraction: CWT has been selected for the features extraction. Unlike discrete wavelet transform, CWT
gives finer resolution without skipping samples. CWT reveals unique details of time-frequency variations in a single person. In time domain, Wavelet transform of signal
is given by
1
W Tx (a, b) = √
a

Z∞
x(t). ψ

∗



t−b
a


dt

(2)

−∞

In our case, x(t) is a PPG segment. ψ(t) is mother
wavelet which is basically bandpass filter that changes
with scale factor a and translation b. By such translation and scaling, mother wavelet separates mixed frequency components of given PPG segment. Since analytic wavelets are more suitable for oscillatory signal
and morse wavelet can approximate many other analytic
wavelets by selecting appropriate parameters, in experiments ‘Analytic Morse’ wavelet was used as mother
wavelet.
• Dimensionality Reduction: Linear discriminant analysis (LDA) based dimensionality reduction was applied on
features vectors generated in previous stage. CWT along
with LDA emphasizes time-frequency behavior in such a
way that it increases the inter subject variability and reduces the intra subject variability. Such features selection
is also important because PPG being physiological signal
changes with time. Prior to dimensionality reduction, features were normalized on time axis using zero padding.
K
Given a training set Z = {Zj }j=1 containing K classes
k
each having Zk ={zki }N
i=1 where Nk is number of feature
vectors for each class k, LDA finds projection weight W
that maximizes the fisher criterion function J(W ),

J(W ) = arg max
W

|W T Sb W |
|W T Sw W |

(3)

here Sb and Sw are the between class and within class
scatter matrix respectively defined as,
Sb =

K
X

Nk (µk − µ)(µk − µ)T

(4)

k=1

Sw =

Nk
K X
X
k=1 i=1

(zki − µk )(zki − µk )T

(5)

PNk
where µk = N1k i=1
zki is mean of class Zk and µ is
overall mean vector. LDA finds W as m=K-1 most significant eignevectors of (Sw )−1 (Sb ) that corresponds to
first m largest eigenvalues. After obtaining weight W, all
the input training PPG segments are subjected to linear
projection of z̃ = W T z. In practice, number of features
might be larger than number of training samples. In that
case, LDA suffers from small sample size problem. We
used Direct-LDA (DLDA) to address this issue[11]. Projected training vectors z̃ together with weight W is saved
in gallery for template matching.

2.3. Template Matching
Prior to template matching of test data, it is passed
through all the blocks i.e. pre-processing and template creation. Weights of the LDA is pre-saved from the training
session and test feature vectors are projected in similar way.
Template matching is carried out using Pearson’s distance
between training templates and test vectors. Pearson’s distance Sp between two vectors a and b is defined as:
Sp (a, b) = 1 − p

cov(a, b)
cov(a, a).cov(b, b)

(6)

where cov(a, b) is the covariance matrix of a and b. Test
vector is accepted if Sp ≤ threshold.

3. Other Implemented Methods
Apart from DLDA, few other techniques were also implemented to reproduce and compare the results.
• DWT/kNN: It is based on method presented in [4] by
Karimian et.al. It is a non-fiducial method which uses
DWT coefficients extracted using coiflet wavelet from
PPG segment as features. It implements two steps process based on kolmogrov smirnov based correlation filter
(ksCBF) and Kernel PCA (KPCA) to remove correlated
features and to reduce dimensionality. Classification is
carried out using local density factor based kNN.
• AC/LDA [3]: It is widely used in ECG based recognition
system. It doesn’t require peak detection in PPG. Instead,
signal is blindly segmented into overlapping windows
of predefined length. Then, normalized auto-correlation
(AC) of each window is calculated as following:
Rˆxx [m] =

PN −|m|−1
i=0

x[i] x[i + m]
ˆ
Rxx [0]

(7)

where, x[i] represents windowed PPG singal, x[i+m] delayed PPG singal with time lag of m = 0, 1, 2.....M − 1
and M << N . These AC windows are projected on
lower dimensional space using LDA. Finally template
matching is done with euclidean distance.

Table 1: Capnobase dataset results with DLDA when evaluated with different number of randomly selected consecutive test
segments (nTest) and trained using only first 45 seconds.
nTest
mean EER
std. EER

2
1.12%
0.44%

5
1.09%
0.45%

10
0.96%
0.38%

20
0.85%
0.29%

30
0.80%
0.27%

40
0.80%
0.25%

50
0.74%
0.21%

100
0.58%
0.14%

All
0.46%
0.00%

Table 2: Capnobase dataset results with different methods when evaluated with All test segments. Here std.EER is zero.
Method
EER

CWT/DLDA CWT/KDDA Openset CWT/KPCA CWT/LDA CWT/PCA Karimian et.al. AC/LDA
0.46%
2.32%
2.50%
2.38%
2.32%
4.01%
1.51±0.48%
4.82%

• Openset Validation: It follows same system as decribed in Fig.3 except doing dimensionality reduction. In
openset validation, CWT PPG feature segments are stored
in gallery without subjecting them to dimensionality reduction. During template matching, same process is followed and CWT features of claimed person is compared
against CWT feature segments stored in gallery without
projecting them to any lower dimensional space.
• Other Subspace Learning Techniques: Apart from
DLDA other subspace learning techniques such as PCA,
LDA, KPCA and Kernel Direct Discriminant Analysis
(KDDA) [8] were also implemented. For all of these techniques same system as Fig.3 was followed with DLDA dimensionality block replaced with different technique. For
both KPCA and KDDA, Gaussian kernel was used. Reduced dimensionality of all these techniques was set to
N − 1, where N is number of subjects in dataset.

4. Experimental Results and Discussion
4.1. Performance Metric
The performance of all techniques was evaluated in verification mode (1 to 1 matching). In this configuration,
subject can claim an enrolled ID and based on matching
score system rejects or accepts the claim. In this setting,
error can be characterized by False Rejection Rate (FRR)
and False Acceptance Rate (FAR). ROC (Receiver Operating Characteristic) curve plots FRR vs FAR for different
operating points. Operating points are matching score or
threshold values. Trade off appears between FAR and FRR.
EER (Equal Error Rate) is an operating point on ROC where
FAR=FRR. For experiments, EER is chosen as a metric to
compare different techniques. In practice, choice of operating point is application dependent.

4.2. Results
In experiments, evaluation strategy was designed such
that it addressed the issues of robustness against different
stress conditions, against time and speed. To assess the
system performance against all such conditions, three different datasets were used. All of three datasets contained

data recorded in single session, in different emotional stress,
physical stress and data recorded after time-lapse. In experiments, we ran multiple simulations varying different parameters such as CWT scale, segment length, distance metric etc. It was also found that using larger segment length
gave better results in general, since CWT gave better frequency resolution. Segment length was set based on median HR such that it covered 3.5-4 cycles of PPG. In addition, two consecutive time segments were averaged. This
reduced variance due to noise in larger segment and across
whole period for one subject. Experimental results also supported this hypothesis. One more engineering problem was
to select number of scales from CWT to maximize discrimination while keeping feature vector small. It was found
that choosing coefficients from only one scale, which filters
signal between 1-2 Hz gave better results. 1-2 Hz is also
most prominent frequency range of PPG signal. However
due to diverse nature and different sampling frequency of
each of three dataset, best CWT scale for classification was
different for each of the three dataset.
Another challenge is of speed. PPG is slower signal
compared to traditional biometric traits and requires larger
training and testing time. In experiments taking more pragmatic approach, average training time was limited between
45 − 60s. On the other side to assess effect of limited authentication time, different number of test segments were
selected to vary available test time. Experimental results
for each of three dataset under various configurations is presented in subsequent discussions.
Capnobase dataset was used for single session evaluation
[5]. It consists of 8 minute single session data recorded in
relax condition from 42 subjects. Training time was fixed
to first 45s of signal, while different number of consecutive test segments (nTest) were selected randomly over 50
iterations. EER was calculated for each iteration. To select
best scale, numerous experiments were carried out varying
training time, nTest and CWT Scale. Based on these experiments 3rd highest scale was chosen for further experiments. Similar experiments were also carried out to select best scale for other two datasets. Results in Tab.1 and
2 present mean EER and standard deviation in EER (std.

ROC curves comparisons for across emotion evaluation, DEAP dataset

60

50

False Rejection Rate (%)

Table 3: Across session evaluation of DEAP Dataset using
DLDA, nTest is number of selected test segments

OpenSet
Direct LDA
Kernel PCA
LDA
PCA
DWT/kNN
AC/LDA
Kernel Direct DA

nTest
Mean EER
std. EER

40

2
2.61%
1.14%

5
2.30%
0.97%

10
2.13%
0.87%

20
2.00%
0.82%

All
2.11%
0.87%

30

20

10

0
0

10

20

30

40

50

60

False Acceptance Rate(%)

Figure 4: DEAP dataset ROC Curves
EER) over 50 iterations for training time of 45 seconds only.
CWT/DLDA performed better than every other technique.
During the experiments, it was found the increasing training time decreased EER. However, drop wasn’t significant.
Also, increasing testing time (nTest) also decreased EER.
Nevertheless, EER of 1.12 ± 0.44% was achieved using 45s
of training and just 2 segments which required 6-7s in time.
Next, system was tested for emotional robustness using DEAP Dataset [7]. DEAP dataset consists of data
recorded in 40 different kinds of emotions spread over valence arousal plane for 32 subjects. Data was recorded on
same day with short breaks and baseline emotional stimuli in between. For each subject, model was trained using data recorded in one type of emotion, that is total one
minute of training of data and tested it against rest of the
dataset. For each subject, dataset consisted of 39 genuine
trials and 1240 (40 data samples from rest of the 31 subjects) imposter trials. EER was calculated over 40 iteration, each time model trained with different emotion. Here,
nT est were not selected randomly, as test data were from
different sessions. Tab.3 presents mean EER and std. EER
for DLDA. Remarkably, EER of 2.61±1.14 % was achieved
using 60s of training time and 2 test segments or 6-7s of test
time. Fig.4 presents comparison of different techniques using ROC curves. It can be noticed that all other techniques
performed poorly compared to DLDA. Also CWT/DLDA
produced better result than [9] on same dataset.
Apart from emotions, PPG signal is susceptible to variations due to exercise because of heart rate change. For robust system, these variations can not be overlooked. In addition, as shown in Fig.5, PPG being physiological signal,
changes with time. Therefore, consequences of ignoring it
can be disastrous. Like all other biometrics traits, perma-

nence test on PPG is important. Since, no large dataset is
available in public to estimate the effect of physical stress
and time-lapse, BioSec.Lab PPG dataset was developed at
University of Toronto [1]. PPG signals were recorded from
fingertip using Plux Sensor over two sessions [2]. First session was conducted in two parts. First, PPG was recorded
in relax condition for 3 minutes. Then subjects were asked
to perform some form of intense exercise such as climbing up the stairs very fast to increase HR. PPG signal was
again recorded for 3 minutes just after exercise. In total 41
subjects participated in first session. In second session, signals were only recorded in relax condition for 3 minutes.
Two sessions were separated atleast 2 weeks apart in time.
However only 34 subjects out of 41 participated in second
session. For evaluation, Model was trained using only initial 45 seconds of relax data from first session. To assess,
effect of physical stress, model was then tested on exercise
data with different nTest selected from the initial signal nonrandomly. To estimate time-robustness, again model was
trained using 45 seconds of relax data from first session and
then tested on session 2 data non-randomly. Tab.4 presents
results for each case where single session performance on
45 seconds of training data is included for comparison.
From Tab.4, it can be observed that EER increased by
5% because of time-lapse and exercise when tested with
whole signal. However, from Fig.6 it can be seen that,
CWT/DLDA performed far better than any other technique. Surprisingly, DWT/kNN method, which performed
very badly for DEAP dataset, had decreasing EER for
BioSec.Lab dataset, across different sessions. This can be
due to that fact that correlation based filter empolyed in this
method does not take into account robustness of features
across sessions while removing them. Hence it is possible that because of reduced number of features, DWT/kNN
gave decreasing error across 2 sessions for BioSec.Lab
dataset, but on average when tested across 40 sessions in
DEAP dataset, it performed worst. Also, ideally EER
should decrease with increase in nTest. This trend is visible in Tab.1. But in Tab.3 and Tab.4, it is not consistent.
This is due to the quality of samples. For example in Tab.4
across exercise, nTest= 20 has lower EER compared to All
test segments. Because, it is possible that whole test signal (all segments) would have more deviated samples under
physical or mental stress compared to only subset (e.g. 20
segments) of test signal.

Data from subject 25 in session1

1

Comparison of performance for different methods for BioSec.Lab PPG dataset

Data from subject 5 in session1

1

50
45

0.5

0.5

40
0

50

100

150

200

0

250

subject 25 after doing exercise

1

100

150

200

250

subject 5 after doing exercise

1

0.5

50

0.5

0

0
0

50

100

150

200

250

0

50

150

200

250

35
30
25
20
15

subject 5 in session2

subject 25 in session2

1

100

mean EER (%)

0

0

1

0.5

10
5

0.5

0
0

DLDA

0
0

50

100

150

200

250

0

50

100

150

200

KPCA

PCA

Single Session

Figure 5: Variations in PPG signal with exercise and timelapse for BioSec.Lab PPG dataset
Table 4:
BioSec.Lab PPG Dataset results using
CWT/DLDA using only 45s of training data, here
nTest is number of selected test segments.
nTest
Mean EER
std. EER
EER
EER

2
5
10
20
Single Session Evaluation
1.05% 1.00% 0.97% 0.86%
0.22% 0.18% 0.17% 0.02%
Across Exercise Evaluation
2.50% 2.50% 2.95% 2.50%
Across Session Evaluation
6.86% 8.65% 8.82% 8.82%

LDA

KDDA

AC/LDA DWT/kNN Openset

250

All
0.86%
0.02%
5.00%
5.88%

5. Conclusion
In this paper, CWT/DLDA based method was presented
for PPG based authentication. By nature, CWT reveals idiosyncratic time-frequency behavior of PPG for each individual. To get holistic understanding, method was evaluated for different datasets under different conditions. It was
found that, system performed better with larger training and
testing time. However by considering more practical scenario training time was fixed to only 45-60s. Better results
of our method compared to previous works is attributed to
the fact that, analytic wavelet are more suitable for oscillatory signals compared to non-analytic and DLDA boosts
system performance by only choosing features that maximizes discriminality in class specific way. It was found that,
physical stress and time-lapse has adverse effect on system.
In future, we plan to address these issues with more robust feature selection, continuous authentication and training model under many different conditions.

References
[1] Biosec.lab ppg dataset.
http://www.comm.
utoronto.ca/˜biometrics/PPG_Dataset/.

After Exercise

After Time-Lapse

Figure 6: Comparison of EER in different cases

[2] Pulse sensor.
https://store.plux.info/
bitalino-sensors/42-pulsesensor.html.
Last accessed on 21 Dec 2017.
[3] F. Agrafioti and D. Hatzinakos. Ecg based recognition using
second order statistics. In 6th Annual Communication Networks and Services Research Conference (cnsr 2008), pages
82–87, May 2008.
[4] N. Karimian, M. Tehranipoor, and D. Forte. Non-fiducial
ppg-based authentication for healthcare application. In 2017
IEEE EMBS International Conference on Biomedical Health
Informatics (BHI), pages 429–432, Feb 2017.
[5] W. Karlen, S. Raman, J. M. Ansermino, and G. A. Dumont.
Multiparameter respiratory rate estimation from the photoplethysmogram. IEEE Transactions on Biomedical Engineering, 60(7):1946–1953, July 2013.
[6] A. R. Kavsaolu, K. Polat, and M. R. Bozkurt. A novel feature ranking algorithm for biometric recognition with ppg
signals. Computers in Biology and Medicine, 49:1 – 14,
2014.
[7] S. Koelstra, C. Muhl, M. Soleymani, J. S. Lee, A. Yazdani,
T. Ebrahimi, T. Pun, A. Nijholt, and I. Patras. Deap: A
database for emotion analysis using physiological signals.
IEEE Transactions on Affective Computing, 3(1):18–31, Jan
2012.
[8] J. Lu, K. N. Plataniotis, and A. N. Venetsanopoulos. Face
recognition using kernel direct discriminant analysis algorithms. IEEE Transactions on Neural Networks, 14(1):117–
126, Jan 2003.
[9] A. Sarkar, A. L. Abbott, and Z. Doerzaph. Biometric authentication using photoplethysmography signals. In 2016 IEEE
8th International Conference on Biometrics Theory, Applications and Systems (BTAS), pages 1–7, Sept 2016.
[10] P. Spachos, J. Gao, and D. Hatzinakos. Feasibility study
of photoplethysmographic signals for biometric identification. In 2011 17th International Conference on Digital Signal Processing (DSP), pages 1–5, July 2011.
[11] H. Yu and J. Yang. A direct lda algorithm for highdimensional data with application to face recognition. Pattern Recognition, 34(10):2067 – 2070, 2001.

