-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathRobustACC.m
106 lines (92 loc) · 2.86 KB
/
RobustACC.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
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
function Robustacc = RobustACC(InputMartix)
%RobustACC - Description
%
% Syntax: Robustacc = RobustACC(InputMartix)
%
% This function is used to calculte ROBUST by each point
%
% ROBUST = mean(ACC(i-1,j),ACC(i+1,j),ACC(i,j-1),ACC(i,j+1));
SIZE = size(InputMartix);
Robustacc = zeros(SIZE(1),SIZE(2));
for i = 2,SIZE(1) - 1
for j = 2,SIZE(2) -1
r = corrcoef(InputMartix(i,j,:),InputMartix(i-1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i+1,j,:));
ACC2 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC3 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j+1,:));
ACC4 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2 ACC3 ACC4]);
end
end
j = 1;
for i = 2,SIZE(1) - 1
r = corrcoef(InputMartix(i,j,:),InputMartix(i-1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i+1,j,:));
ACC2 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j+1,:));
ACC3 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2 ACC3]);
end
j = SIZE(2);
for i = 2,SIZE(1) - 1
r = corrcoef(InputMartix(i,j,:),InputMartix(i-1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i+1,j,:));
ACC2 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC3 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2 ACC3]);
end
i = 1;
for j = 2,SIZE(2) - 1
r = corrcoef(InputMartix(i,j,:),InputMartix(i+1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC2 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC3 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2 ACC3]);
end
i = SIZE(1);
for j = 2,SIZE(2) - 1
r = corrcoef(InputMartix(i,j,:),InputMartix(i-1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC2 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC3 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2 ACC3]);
end
i = 1;
j = 1;
r = corrcoef(InputMartix(i,j,:),InputMartix(i+1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j+1,:));
ACC2 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2]);
i = 1;
j = SIZE(2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i+1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC2 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2]);
i = SIZE(1);
j = 1;
r = corrcoef(InputMartix(i,j,:),InputMartix(i-1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j+1,:));
ACC2 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2]);
i = SIZE(1);
j = SIZE(2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i-1,j,:));
ACC1 = r(1,2);
r = corrcoef(InputMartix(i,j,:),InputMartix(i,j-1,:));
ACC2 = r(1,2);
Robustacc(i,j) = mean([ACC1 ACC2]);
end