MNL-MotionCorrection/metric_gui.m at master · davidbichara/MNL-MotionCorrection · GitHub

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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
%Uses sumstatisics.m to create a metric and tag which gives an estimate of
%movement as well as where/what direction it occured

move_score = 0;
iso_score = 0;
consis_score = 0;

lower_bound = 0;
upper_bound = 0;

%integral compared to first rep threshold
if sum_stat(1) > 50000
    move_score = move_score + 1;
end

%integral compared to prev rep threshold
if sum_stat(2,1) > 10000
    move_score = move_score + 1;
end

%compares integral of sections of reps to full integral
index_count = [8 9 10 11];
index_string = [1 6 12 18 24];
average_area = sum_stat(2,1)/4;
for ii = 1:4
    if sum_stat(index_count(ii),1) > average_area
        lower_bound = index_string(ii);
        upper_bound = index_string(ii+1);
        if lower_bound ~= 0
            figure;
            p1 = plot(find(final_movement1), final_movement1, 'Color', 'red');
            hold on
            p2 = plot(find(final_movement2), final_movement2, 'Color', 'Blue');
            hold off
            xlabel('Repetition', 'FontSize', 22);
            ylabel('Amount of Voxels Changed', 'FontSize', 22);
            axis([lower_bound upper_bound 0 9000]);
            p1.LineWidth = 7;
            p2.LineWidth = 7;
            set(gca,'YTick',0:1000:5000)
            set(gca,'fontsize', 18)
            ay = gca;
            ay.YRuler.Exponent = 0;
            lower_bound = 0;
            upper_bound = 0;
            x = input('Y/N: Does this set look irregular:', 's');
            if x == 'Y'
                move_score = move_score + 5;
            end
            if x == 'N'
                return
            end
        end
    end
end


%ratio of max-mean compared to first rep threshold
if sum_stat(12,1) > 2
    move_score = move_score + 1;
end


%adjusted/non-adjusted range compared to prev rep threshold
if sum_stat(15,1) < 70
    move_score = move_score + 1;
end

if sum_stat(15,1) < 50
    iso_score = iso_score +1;
end


%adjusted/non-adjusted mean compared to prev rep threshold
if sum_stat(17,1) > 70 && sum_stat(17,1) < 87
    consis_score = consis_score +1;
    move_score = move_score + 1;
end

if sum_stat(17,1) <= 70
    move_score = move_score + 2;
    iso_score = iso_score +1;
end


%adjusted/non-adjusted mean compared to first rep threshold
if sum_stat(19,1) < 90
    disp('AP movement may have occured')
    move_score = move_score + .5;
end

%establishing feeback thresholds for types of movement
if consis_score ~= 0
    disp("Consistent movement may have occured")
    iso_score = 0;
end


if iso_score ~= 0
    disp("Isolated movement may have occured")
end

%Activation functions for movement score

%move_metric = atan(move_score);

%e_move_metric = exp(move_score);