-
Notifications
You must be signed in to change notification settings - Fork 17
/
Inference_HRevaluation.m
66 lines (41 loc) · 1.68 KB
/
Inference_HRevaluation.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
clear
clc
load Inference_Physformer_TDC07_sharp2_hid96_head4_layer12_VIPL.mat
GT_list = importdata('VIPL_fold1_test1.txt');
total_samples = length(outputs_rPPG_concat);
HR_peaks = [];
HR_PSD = [];
HR_GT = [];
signal = double(outputs_rPPG_concat);
framerate = GT_list.data(1,2);
GT_HR = GT_list.data(1,3);
% signal_filtered = signal;
signal_filtered = bpfilter64(signal, framerate);
signal_filtered = (signal_filtered-mean(signal_filtered))/std(signal_filtered);
%% PSD for HR
% % Single long clip
% [Pg,f] = pwelch(signal_filtered,[],[],2^13,framerate);
% Frange = find(f>0.7&f<3); % consider the frequency within [0.7Hz, 4Hz].
% idxG = Pg == max(Pg(Frange));
% HR2 = f(idxG)*60;
% % Separate into three clips
signal_length = length(signal_filtered);
[Pg,f] = pwelch(signal_filtered(1:floor(signal_length/3)),[],[],2^13,framerate);
Frange = find(f>0.7&f<4); % consider the frequency within [0.7Hz, 4Hz].
idxG = Pg == max(Pg(Frange));
HR2_1 = f(idxG)*60;
[Pg,f] = pwelch(signal_filtered(floor(signal_length/3):2*floor(signal_length/3)),[],[],2^13,framerate);
Frange = find(f>0.7&f<4); % consider the frequency within [0.7Hz, 4Hz].
idxG = Pg == max(Pg(Frange));
HR2_2 = f(idxG)*60;
[Pg,f] = pwelch(signal_filtered(2*floor(signal_length/3):signal_length),[],[],2^13,framerate);
Frange = find(f>0.7&f<4); % consider the frequency within [0.7Hz, 4Hz].
idxG = Pg == max(Pg(Frange));
HR2_3 = f(idxG)*60;
HR2 = (HR2_1+HR2_2+HR2_3)/3;
HR_peaks = [HR_peaks; HR1];
HR_PSD = [HR_PSD; HR2];
HR_GT = [HR_GT; GT_HR];
%% calculate ErrorMean, ErrorSD, RMSE, R
Error_PSD = HR_PSD - HR_GT;
MAE = abs(Error_PSD)