Similar documents
Introduction to Wavelet Transform. Chapter 7 Instructor: Hossein Pourghassem

Fourier and Wavelets

FEASIBILITY STUDY OF PHOTOPLETHYSMOGRAPHIC SIGNALS FOR BIOMETRIC IDENTIFICATION. Petros Spachos, Jiexin Gao and Dimitrios Hatzinakos

EE216B: VLSI Signal Processing. Wavelets. Prof. Dejan Marković Shortcomings of the Fourier Transform (FT)

Sensor, Signal and Information Processing (SenSIP) Center and NSF Industry Consortium (I/UCRC)

Unlock with Your Heart: Heartbeat-based Authentication on Commercial Mobile Phones

Exam in 1TT850, 1E275. Modulation, Demodulation and Coding course

Introduction to Wavelets Michael Phipps Vallary Bhopatkar

RC Circuit Activity. Retrieve a power cord and a voltage sensor from the wire rack hanging on the wall in the lab room.

ECE 476/ECE 501C/CS Wireless Communication Systems Winter Lecture 6: Fading

Next Generation Biometric Sensing in Wearable Devices

Objectives. Presentation Outline. Digital Modulation Revision

Virtual Spirometry and Activity Monitoring using Multi-channel Electrical Impedance Plethysmographs in Ambulatory Settings

Validation of the Happify Breather Biofeedback Exercise to Track Heart Rate Variability Using an Optical Sensor

Fingerprint Segmentation using the Phase of Multiscale Gabor Wavelets

Harmonic Analysis of Power System Waveforms Based on Chaari Complex Mother Wavelet

Objectives. Presentation Outline. Digital Modulation Lecture 03

ECE 476/ECE 501C/CS Wireless Communication Systems Winter Lecture 6: Fading

ECE 476/ECE 501C/CS Wireless Communication Systems Winter Lecture 6: Fading

Feature Extraction Techniques for Dorsal Hand Vein Pattern

Amplitude Frequency Phase

ECG Biometrics using Intuitive Bases and Support Vector Machines. Peter Sam Raj

Practical Applications of the Wavelet Analysis

Automobile Independent Fault Detection based on Acoustic Emission Using FFT

Detection, localization, and classification of power quality disturbances using discrete wavelet transform technique

CO-CHANNEL SPEECH DETECTION APPROACHES USING CYCLOSTATIONARITY OR WAVELET TRANSFORM

PLL EXERCISE. R3 16k C3. 2π π 0 π 2π

Wavelet-based image compression

EENG473 Mobile Communications Module 3 : Week # (12) Mobile Radio Propagation: Small-Scale Path Loss

WAVELET TRANSFORM BASED METHOD FOR EDDY CURRENT TESTING OF CLADDING TUBES

Voice Activity Detection

Transient-based Identification of Wireless Sensor Nodes

Design Considerations for Wrist- Wearable Heart Rate Monitors

Study on the Algorithm of Vibration Source Identification Based on the Optical Fiber Vibration Pre-Warning System

TRANSFORMS / WAVELETS

Chapter 9. Carrier Acquisition and Tracking. March 5, 2008

II. Random Processes Review

Narrow- and wideband channels

Design of IIR Filter Using Model Order Reduction. Techniques

Effects of multipath propagation on design and operation of line-of-sight digital radio-relay systems

Wheel Health Monitoring Using Onboard Sensors

Elham Torabi Supervisor: Dr. Robert Schober

1 ONE- and TWO-DIMENSIONAL HARMONIC OSCIL- LATIONS

Biomedical Signals. Signals and Images in Medicine Dr Nabeel Anwar

Analysis of Visible Light Communication Using Wireless Technology

Throughput-optimal number of relays in delaybounded multi-hop ALOHA networks

DSP First. Laboratory Exercise #7. Everyday Sinusoidal Signals

Use of Matched Filter to reduce the noise in Radar Pulse Signal

FEASIBILITY OF SINGLE-ARM SINGLE-LEAD ECG BIOMETRICS. Peter Sam Raj, Dimitrios Hatzinakos

3 Department of Computer science and Application, Kurukshetra University, Kurukshetra, India

Satellite Navigation Principle and performance of GPS receivers

Gait Recognition Using WiFi Signals

SUPPLEMENT TO THE PAPER TESTING EQUALITY OF SPECTRAL DENSITIES USING RANDOMIZATION TECHNIQUES

Digital Signal Processing PW1 Signals, Correlation functions and Spectra

Chapter 5. Array of Star Spirals

Evoked Potentials (EPs)

Mechanics of Materials and Structures

Signal Processing Techniques for Damage Detection with Piezoelectric Wafer Active Sensors and Embedded Ultrasonic Structural Radar ABSTRACT Keywords

Time-Frequency Analysis of Millimeter-Wave Radar Micro-Doppler Data from Small UAVs

Digital modulations (part 1)

Chapter 4. Part 2(a) Digital Modulation Techniques

SAR AUTOFOCUS AND PHASE CORRECTION TECHNIQUES

ME 365 FINAL EXAM. Monday, April 29, :30 pm-5:30 pm LILY Problem Score

COMPRESSIVE SENSING BASED ECG MONITORING WITH EFFECTIVE AF DETECTION. Hung Chi Kuo, Yu Min Lin and An Yeu (Andy) Wu

WIND SPEED ESTIMATION AND WIND-INDUCED NOISE REDUCTION USING A 2-CHANNEL SMALL MICROPHONE ARRAY

Distinguishing Identical Twins by Face Recognition

Automatic Peak Picking Using Wavelet De-noised Spectra in Automated Struture Determination.

International Journal of Scientific & Engineering Research, Volume 7, Issue 12, December ISSN IJSER

Evaluation of Biometric Systems. Christophe Rosenberger

Class XII Chapter 7 Alternating Current Physics

Practical Application of Wavelet to Power Quality Analysis. Norman Tse

Chapter 5. Signal Analysis. 5.1 Denoising fiber optic sensor signal

A Novel Image Deblurring Method to Improve Iris Recognition Accuracy

In this lecture we consider four important properties of time series analysis. 1. Determination of the oscillation phase.

Real Time Detection and Classification of Single and Multiple Power Quality Disturbance Based on Embedded S- Transform Algorithm in Labview

PERFORMANCE TESTING EVALUATION REPORT OF RESULTS

Segmentation of Fingerprint Images

Time-Slotted Round-Trip Carrier Synchronization

Mobile Radio Propagation Channel Models

\quad. ond e r fulpr ert e i s h f

Heart Rate Monitoring using Adaptive Noise Cancellation

Biometric: EEG brainwaves

, answer the next six questions.

Radar signal detection using wavelet thresholding

BEAT DETECTION BY DYNAMIC PROGRAMMING. Racquel Ivy Awuor

Lab 1: Basic RL and RC DC Circuits

ECE 201: Introduction to Signal Analysis

FAULT DETECTION OF FLIGHT CRITICAL SYSTEMS

DSP First Lab 08: Frequency Response: Bandpass and Nulling Filters

Muhammad Ali Jinnah University, Islamabad Campus, Pakistan. Fading Channel. Base Station

WAVELETS: BEYOND COMPARISON - D. L. FUGAL

Written Exam Channel Modeling for Wireless Communications - ETIN10

Ranging detection algorithm for indoor UWB channels and research activities relating to a UWB-RFID localization system

Bicorrelation and random noise attenuation

A simple output voltage control scheme for single phase wavelet modulated inverters

Chapter 5. Frequency Domain Analysis

Empirical Evidence for Correct Iris Match Score Degradation with Increased Time-Lapse between Gallery and Probe Matches

WAVELET ANALYSIS TO DETECT THE KNOCK ON INTERNAL COMBUSTION ENGINES

ANALYSIS OF PARTIAL DISCHARGE SIGNALS USING DIGITAL SIGNAL PROCESSING TECHNIQUES

IDENTICAL AND FRATERNAL TWIN RECOGNITION USING PHOTOPLETHYSMOGRAM SIGNALS

Transcription:

. 36.9

d = 4 p = 5

.8 1.5

(1 : N) 1 : N 1 : (N + 1) (N + 1) th G I (nt P ) (nt N) (nf P ) (nf N) F AR = nf P I F RR = nf N G T AR = nt P G T RR = nt N I G + I = N T AR = 1 F RR T RR = 1 F AR

nt P + nt N Accuracy = G + I N G << I G = 1 I = 1 99 + nt P = G << I G I.1 1.1

1 2

Amplitude 4 Vasoconstriction PPG signal with its characteristic points Systolic Peak 3 2 1 Diastolic Peak -1 Inter Pulse Interval (IPI) Diacrotic Notch Diastolic Point 8 8.5 9 9.5 1 1.5 11 Time(s)

23 94 17

94 82.3 5.29 9 9.44 k = 1 94.44 k = 1 87.22

k = 1 4.2 3.7

5.29 + + + + π 96 9.53

+ +

52 19 35 1

Figure 3.1: Plux pulse sensor and black Velcro strip used in this study Figure 3.2: OpenSignal Desktop Application Figure 3.3: Bitalino MCU and Placement of pulse sensor on fingertip with black Velcro strip on the go anywhere via Bluetooth. Hence, these devices indeed emulate practical recording devices such as the smart-watch. Followed protocol before each recording: An open call for participation was circulated to invite as many people as possible. Common protocol followed during each recording is as follows. Prior to recording sessions, all participants were asked to sit down on a chair in a suitable position and relax for a while to achieve a normal heart rate. Participants were briefed about objective of this study and its importance. Any questions from them were answered and an informed consent form was signed for data collection. Then, Plux sensor was placed on either their left or right index finger-tip, irrespective of which fingertip was used in earlier sessions in multi-session configurations. Proper sensor contact was ensured by adjusting Velcro strip. It was found that a tight grip with Velcro strip resulted in poor signals. Therefore in all cases, it was made sure that Velcro strip provided moderate grip. Unconstrained subject: To minimize motion artifacts, subjects were requested to avoid any hand 21

Amplitude Amplitude Amplitude Amplitude Filtered Relax signal for subject ID 22 Filtered Relax signal for subject ID 1 2 1 4 2.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered Relax signal for subject ID 59.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered Relax signal for subject ID 37 1 6 5 4 2-5.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) -2.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) 84 4 4

Amplitude Amplitude Amplitude Amplitude Filtered Relax signal for subject ID 39 Filtered Relax signal for subject ID 7 6 4 4 2-2 -4.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered after exercise signal for subject ID 39 2-2 3.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered after exercise signal for subject ID 7 5 2 1-5.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) -1.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) 3

Amplitude Amplitude Amplitude Amplitude Amplitude Amplitude Amplitude Amplitude Amplitude Amplitude 15 1 5 Filtered signal from session 1 for subject ID 37 8 6 4 2-2 Filtered signal from session 1 for subject ID 1.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered signal from session 2 for subject ID 37.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered signal from session 2 for subject ID 1 2 15 1 5 8 6 4 2-2.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s).5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) 6 2 4-2 2 1-1 6 4 2 Filtered signal from sesson 1 for subject ID 5.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered signal from session 2 for subject ID 5.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered signal from session 3 for subject ID 5.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) 6 4 2-2 4 2 5-5 Filtered signal from session 1 for subject ID 3.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered signal from session 2 for subject ID 3.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s) Filtered signal from session 3 for subject ID 3.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 Time (s)

Days Days 1 9 8 7 Time difference between First and Second Session 97 89 77 69 67 63 63 14 6 5 4 29 29 32 3 29 29 29 29 3 31 3 3 27 28 2 1 29 29 23 22 3 29 21 22 19 18 17 17 15 14 14 14 1 2 3 4 5 6 7 8 9 1 11 12 13 14 15 16 17 18 19 2 21 22 23 24 25 26 27 28 29 3 31 32 33 34 35 36 37 Subject ID Time difference between first, second and third sessions 238 2 197 197 23 22 21 21 25 Time difference between first and second session Time difference between second and third session Time difference between first and third session 168 168 173 175 17 173 171 172 175 15 1 18 97 77 1 67 19 112 63 63 111 11 18 16 13 9 91 9 89 89 5 29 29 3 27 31 28 23 21 19 18 17 14 1 2 3 4 5 6 7 8 9 1 11 12 13 14 15 16 Subject ID

37 14 36 37 16 128 6

± ± ± ±

4 th.5 4 x norm = x min(x) max(x) min(x)

Data Acquisition Filtering Peak Detection Segmentation Pre-Processing Outlier Segments Removal PPG Segments PPG Segments Feature Extraction Feature Reduction PPG Feature Vector Generation PPG Feature Vectors PPG Feature Vectors System Database Train a Classifier Model Trained Classifier Model PPG Feature Vectors Predict with Trained Classifier Model Yes No Score > Threshold? Verification Phase Accept Reject Enrollment Phase.2.7.4.5.3

P1(f) P1(f) 1 2 X: Y: 518.4 Frquency plot of Raw PPG signal for subjectid 67 X: 1.417 Y: 12.47 1 1 1 Frquency plot of Filtered PPG signal for subjectid 67 X: 1.417 Y: 12.56 1 1-1 1-2 1-2 1-3 X: Y:.5761 1-4 2 4 6 8 112141618222242628332343638442444648 f (Hz) 1-4 2 4 6 8 112141618222242628332343638442444648 f (Hz)

Amplitude Amplitude Amplitude Amplitude x w f s l = 1 + k w s i i th [ ] i 1 l s i s i s 2 i s i.2.7.4 f s <.4 f s [ ].3 f s 1.3 f s 1 length(sysp eaks) (j) j th x(z minscan) x(z) 1 Peak detection on filtered signal from session 2 1 Peak detection on filtered signal from session 29.8.8.6.6.4.4.2.2 1 2 3 4 5 6 Time (s) 1 2 3 4 5 6 Time (s) 1 Peak detection on filtered signal from session 5 1 Peak detection on filtered signal from session 17.8.8.6.6.4.4.2.2 1 2 3 4 5 6 Time (s) 1 2 3 4 5 6 Time (s)

Amplitude Amplitude.5 PPG Segments Waterfall diagram of PPG segments 1.5 1.8 1.6.4.5.2 1.5.1.2.3.4.5.6.7.8.9 1 Time (s) 8 6.8 Segments 4.6 2.4.2 Time (s) 1 1.5.5.2

e jwt F (ω) = x(t).e jwt F (ω) ω ω ST F T (τ, ω) = x(t).w(t τ).e jwt x(t) ψ(t) a > W T x (a, b) W T x (a, b) = 1 a ( ) t b x(t). ψ a a a = 2, 2 2, 2 4... τ

Amplitude Amplitude Daubechies 5 Mother Wavelet 1 PPG Segment -1-1 1-1 -1 1-1 1 2 3 4 5 6 7 8 9 1 3 Samples.9.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 Time (s) 32 1 = 32

CWT Coefficient CWT of a PPG Segment 1.5 1.5 -.5-1 3-1.5 1 2 3 1 4 5 6 7 8 Samples 9 1 2 Scales

d = 4 p = 5 p BP (x(t)) = sign(x(t + i p d + 1) x(t)).2 i + sign(x(t + i + d) x(t)).2 i+p i= 1, x + ϵ, sign(x) =, x(t) d p x(t) p d x(t) k t + i + d > k t + i < p + d ϵ

n x(t) w n n n n x(t) 1s 1 k = 1 p = (4, 4) d = (1, 1) w = (5, 4) = (1, 2) p = 4 4 d = 1 x(t) 2 2p = 256 255 256 w = 5 = 1 l = 1 + k w = 6 1536 w = 4 = 2 124 256

ˆR xx [m] = N m 1 i= x[i] x[i + m] ˆR xx [] x[i] x[i + m] m =, 1, 2...M 1 N M << N Z {Z j } K j=1 Z k {z ki } N k i=1 k W

J(W ) W = arg W J(W ) = arg W ( W T ) S b W W T S w W S b S w S b = S w = K N k (µ k µ)(µ k µ) T k=1 K N k (z ki µ k )(z ki µ k ) T k=1 i=1 µ k 1 N k Nk i=1 z ki Z k µ J(W ) W W max(w T S b W ) s.t. W T S w W = 1 L(W, λ) = W T S b W λ(w T S w W 1) L(W, λ) W = = S bw = λs w W S w W m = K 1 (S w ) 1 (S b ) S b K 1 S w (S w ) 1 (S b ) K 1 K 1 W z = W T z z W S w 32 S w 32 32 = 1 S w

S w S b S w a S w a = S b a a T S b a/a T S w a a S b a T S b a = S b S w S b W T S b W = Λ W T S w W = I S b S b S w N N K S b = N k (µ k µ)(µ k µ) T = Φ b Φ T b (N N) k=1 Φ b = [ N 1 (µ 1 µ), N 2 (µ 2 µ), N 3 (µ 3 µ)...](n K) K N N Φ b Φ T b Φ T b Φ bv i = λ i v i Φ b Φ b Φ T b Φ bv i = λ i Φ b v i S b = Φ b Φ T b u i = Φ b v i S b u i = λ i u i S b Φ T b Φ b v i S b K K Φ T b Φ b

S w S w S w W = arg W J(W ) = arg W ( W T S b W ) η(w T S b W ) + (W T S w W ) η W η η > S w k(, ) = exp η σ ) ( 2 σ 2

y i { 1, 1} i R d i = 1...N 1,b,ξ i 2 2 + C N i=1, y i ( T i b) 1 ξ i i = 1, 2,..., N. ξ i ξ i ξ i 1 ξ i > 1 i ξ i C >

a i b N y( ) = a i y i i. + b y( ) = i=1 N a i y i k( i, ) + b k( i, ) i=1 a i a i > ( 2 ) k(, ) = exp σ σ 2 97.56

k k ν ν ν C N ξ i = C + ξ i + C i i I + i I ξ i C + C I + I C + C

C + n + = C n C + C = n n + n + n d d 51 GDI = 1 i p 2 i i p i i

S w 2 6 η σ k( i, j ) = exp ) ( i j 2 σ 2 81 81 65 σ σ σ k( i, j ) i j σ k( i, j ) 1 σ σ σ σ σ η η η σ 3 65 1.9 i 32 η σ σ [1 3, 1 17 ] η =.1 σ

[1 5, 1 7 ] σ = 1 5, 1 6, 1 7 η = 1 3, 1 4, 1 5 σ 84 C σ σ σ σ C σ σ C σ σ C + C + σ a i

51

1 : 4 1 : 8 1 : 16 1 : 32 1 : 64 1 : 128 1 : 2 : 16 1 : 2 : 32 1 : 2 : 64 1 : 2 : 128, 1 : 4 : 16 1 : 4 : 32 1 : 4 : 64 1 : 4 : 128 1 : 4 : 136

σ η =.1

EER (%) 5 4 3 3.91 29.15 36.81 DLDA & R-KDA Performance comparison with CWT 45.4 4.69 44.33 33.15 31.6 39.3 38.53 36.96 35.94 3.16 34.69 CWT with DLDA CWT with R-KDA CWT without reduction 37.5 37.8 31.25 31.25 2 1 7.82 3.94 3.57 Single Session Short Time Lapse Across Exercise Two Sessions Long Time Lapse Three Sessions, Train on 1st, Test on 2nd BioSec.Lab PPG Database Configurations Three Sessions, Train on 1st, Test on 3rd Three Sessions, Train on 2nd, Test on 3rd σ = 1 5, 1 6, 1 7 η = 1 3, 1 4, 1 5 η = 1 4 3.57 29.19

4.69 31.6 ϵ

1 st 2 nd 1 st 3 rd 2 nd 3 rd

η =.1 σ

EER (%) 45 4 35 3 25 LBP 4 points with DLDA LBP 4 points with R-KDA LBP 4 points without reduction LBP 6 points with DLDA LBP 6 points with R-KDA LBP 6 points without reduction 31. 29.9 31.68 27.27 3.62 29.9 DLDA & R-KDA Performance Comparison with 1DMRLBP 47.28 41.76 4.59 41. 4.37 39.82 39.17 38.86 38.5138.89 37.76 36.92 37.29 36.8936.7337.9 36.9 34.8 42.41 41.36 39.91 39.93 34.18 32.5 33.11 31.82 29.728.82 26.56 25.63 2 15 1 5 6.5 6.57 6.38 5.2 5.76 5.64 Single Session Short Time Lapse Across Exercise Two Sessions Long Time Lapse Three sessions, Train on 1st, Test on 2nd Three Sessions, Train on 1st, Test on 2nd Three Sessions, Train on 2nd, Test on 3rd BioSec.Lab PPG Database Configurations p = (4, 4), d = (1, 2) p = 4, d = 1 p = 4, d = 2 n n 1 ϵ =.1 ϵ =.1 p = (4, 4, 4, 4) p = (4, 4) p = (4, 4, 4, 4) p = (4, 4, 4, 4, 2, 2) η = 1 3, 1 4, 1 5 σ = 1 5, 1 6, 1 7 η = 1 4 η =.1 σ 4 6

p = (4, 4, 4, 4) p = (4, 4, 4, 4, 2, 2) 5.2 27.27, 36.9, 34.8 7 8

σ η = 1e 5

EER (%) 5 AC with DLDA AC with R-KDA 4 AC without reduction 34.86 DLDA & R-KDA Performance Comparison with Auto-Correlation 48.84 47.42 47.32 47.24 45.95 45.63 45.42 44.41 43.64 43.64 43.61 43.41 42. 4.13 4.54 39.58 38.24 41.81 3 23.43 23.81 2 1 Single Session Short Time Lapse Across Exercise Two Sessions Long Time Lapse Three Sessions Train on 1st, Test on 2nd BioSec.Lab PPG Database Configurations Three Sessions Train on 1st, Test on 3rd Three Sessions Train on 2nd, Test on 3rd 23.42 η = 1 5 η 43.64 43.16 4.13

EER (%) 5 Performance Comparison of Best Results from all Features 5. 49.17 4 CWT 1DMRLBP AC DWT 45.45 43.64 43.61 44.23 4.69 36.9 44.89 4.13 34.8 44.41 38.51 35.94 4.54 38.24 36.67 3 23.43 29.15 27.27 31.6 32.5 3.16 31.25 25.63 2 14.41 1 5.2 3.57 Single Session Short Time Lapse Across Exercise Two Sessions Long Time Lapse Three Sessions Long Time Lapse, Train on 1st, Test on 2nd BioSec.Lab Database Configurations Three Sessions Long Time Lapse, Train on 1st, Test on 3rd Three Sessions Long Time Lapse, Train on 2nd, Test on 3rd 1

FRR FRR FRR FRR FRR FRR FRR 1.9.8.7 ROC Curves Comparison for Single Session CWT LBP AC DWT 1.9.8.7 ROC Curves Comparison for Two Sessions Long Time Lapse CWT LBP AC DWT.6.6.5.5.4.4.3.3.2.2.1.1.1.2.3.4.5.6.7.8.9 1 FAR.1.2.3.4.5.6.7.8.9 1 FAR 1.9.8 ROC Curves Comparison for Across Exercise CWT LBP AC DWT 1.9.8 ROC Curves Comparison for Short Time Lapse CWT LBP AC DWT.7.7.6.6.5.5.4.4.3.3.2.2.1.1 1.9.8.1.2.3.4.5.6.7.8.9 1 FAR ROC Curves Comparison for Three Sessions Long Time Lapse, Train on 1st, Test on 2nd CWT LBP AC DWT 1.9.8 ROC Curves Comparison for Three Sessions Long Time Lapse, Train on 1st, Test on 3rd.1.2.3.4.5.6.7.8.9 1 FAR CWT LBP AC DWT 1.9.8 ROC Curves Comparison for Three Sessions Long Time Lapse, Train on 2nd, Test on 3rd CWT LBP AC DWT.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR

3.57 27.27 84 27.27 35.94 27.27 31.6 3 164 36.9

1 5 7 1 5 2 5 4 1 78 1 1 1 15 p = (4, 4, 4, 4) p = (4, 4, 4, 4, 2, 2)

EER (%) 4 Capnobase Database Results with different Feature 35 3 DLDA + Random Forest DLDA + SVM No Reduction + Random Forest No Reduction + SVM 34. 35.71 25 23.81 2 15 1 9.4 5...28. 2.3 2.38 2.15 1.17.99 1.11 1.34.35 CWT LBP 4 LBP 6 AC Features..35. 4.3 1

FRR FRR ROC Curves Comparison for Capnobase Database.5 CWT LBP4 LBP6.4 AC DWT.3 1.9.8.7.6 ROC Curves Comparison for DEAP Database CWT LBP4 LBP6 AC DWT.2.1.5.4.3.2.1.1.2.3.4.5 FAR.1.2.3.4.5.6.7.8.9 1 FAR 128 8 2 31 4 = 124 4 12.79

EER (%) 35 3 DLDA + Random Forest DLDA + SVM No Reduction + Random Forest No Reduction + SVM DEAP Database Results with different Features 3.45 25 22.73 2 2.86 18.56 18.31 19.94 19.64 15 13.5 16.4 12.79 16.31 15.21 14.83 14.86 14.61 13.48 1 5 CWT LBP 4 LBP 6 AC Features 9.53 1 14 1 35 1

44

3.57. 27.27 35.94 164 36.9 12.79 1 5 7 14 164

FRR FRR FRR FRR FRR FRR FRR ROC Curve Comparison for Single Session with CWT.3 DLDA R-KDA.25 Without Reduction.2.15.1.5 1.9.8.7.6.5.4.3.2.1 ROC Curve Comparison for Two Sessions Long Time Lapse with CWT DLDA R-KDA Without Reduction.5.1.15.2.25.3 FAR ROC Curve Comparison for Across Exercise with CWT 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Short Time Lapse with CWT 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Three Sessions Long Time Lapse Train on 1st, Test on 2nd with CWT 1 DLDA.9 R-KDA Without Reduction.8 ROC Curve Comparison for Three Sessions Long Time Lapse Train on 1st, Test on 3rd with CWT 1 DLDA.9 R-KDA Without Reduction.8.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Three Sessions Long Time Lapse Train on 2nd, Test on 3rd with CWT 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR

FRR FRR FRR FRR FRR FRR FRR ROC Curve Comparison for Single Session with 1DMRLBP.3 LBP4 DLDA LBP4 R-KDA.25 LBP4 without Reduction LBP6 DLDA LBP6 R-KDA LBP6 without reduction.2.15.1.5 1.9.8.7.6.5.4.3.2.1 ROC Curve Comparison for Two Sessions Long Time Lapse with 1DMRLBP LBP4 DLDA LBP4 R-KDA LBP4 without Reduction LBP6 DLDA LBP6 R-KDA LBP6 without reduction.5.1.15.2.25.3 FAR ROC Curve Comparison for Across Exercise with 1DMRLBP 1 LBP4 DLDA.9 LBP4 R-KDA LBP4 without Reduction.8 LBP6 DLDA LBP6 R-KDA.7 LBP6 without reduction.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for shorttimelapse with 1DMRLBP 1.9.8.7.6.5.4.3.2.1 LBP4 DLDA LBP4 R-KDA LBP4 without Reduction LBP6 DLDA LBP6 R-KDA LBP6 without reduction.6.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Three Sessions Long Time Lapse, Train on 1st, Test on 2nd with 1DMRLBP 1 LBP4 DLDA.9 LBP4 R-KDA LBP4 without Reduction.8 LBP6 DLDA LBP6 R-KDA.7 LBP6 without reduction ROC Curve Comparison for Three Sessions Long Time Lapse, Train on 1st, Test on 3rd with 1DMRLBP 1 LBP4 DLDA.9 LBP4 R-KDA LBP4 without Reduction.8 LBP6 DLDA LBP6 R-KDA.7 LBP6 without reduction.6.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Three Sessions Long Time Lapse, Train on 2nd, Test on 3rd with 1DMRLBP 1 LBP4 DLDA.9 LBP4 R-KDA LBP4 without Reduction.8 LBP6 DLDA LBP6 R-KDA.7 LBP6 without reduction.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR

FRR FRR FRR FRR FRR FRR FRR ROC Curve Comparison for Single Session with AC 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1 1.9.8.7.6.5.4.3.2.1 ROC Curve Comparison for Two Sessions Long Time Lapse with AC DLDA R-KDA Without Reduction.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Across Exercise with AC 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Short Time Lapse with AC 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Three Sessions Long Time Lapse, Train on 1st, Test on 2nd with AC 1 DLDA.9 R-KDA Without Reduction.8 ROC Curve Comparison for Three Sessions Long Time Lapse, Train on 1st, Test on 3rd with AC 1 DLDA.9 R-KDA Without Reduction.8.1.2.3.4.5.6.7.8.9 1 FAR ROC Curve Comparison for Three Sessions Long Time Lapse, Train on 2nd, Test on 3rd with AC 1 DLDA.9 R-KDA Without Reduction.8.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR.7.6.5.4.3.2.1.1.2.3.4.5.6.7.8.9 1 FAR