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
|
%-------------- outer loop
for
x=
1
:
40
for
y =
1
:
48
for
z =
1
:
34
%----------inter loop
x=
20
;
y=
30
;
z=
15
;
k =
1
;
for
ztmp =-
1
:
1
z_ztmp = z+ztmp;
if
z_ztmp <
1
|| z_ztmp >
34
continue
;
end
for
ytmp =-
1
:
1
y_ytmp = y+ytmp;
if
y_ytmp <
1
|| y_ytmp >
48
continue
;
end
for
xtmp =-
1
:
1
if
xtmp ==
0
&& ytmp ==
0
&& ztmp ==
0
continue
;
end
x_xtmp = x+xtmp;
if
x_xtmp <
1
|| x_xtmp >
48
continue
;
end
index =
40
*
48
*(z+ztmp-
1
)+
40
*( y+ytmp-
1
)+(x+xtmp);
if
index <
1
|| index >
27524
continue
;
end
data(k,:) = full_nt_matrix(index , :);
k = k+
1
;
end
end
end
%% --------处理data
l = size(data,
1
);
if
l <=
1
continue
;
end
index =
40
*
48
*(z-
1
)+
40
*(y-
1
)+x;
time_series = full_nt_matrix(index,:);
j =
1
;
for
i=
1
:l
r = corr2(time_series,data(i,:));
if
r>
0.5
sim_data(j,:) = data(i,:);
j = j+
1
;
end
end
l = size(sim_data,
1
);
if
l>
0
result_time_series = sum( sim_data,
1
)/
9
;
end
tmp_V(x,y,z,:) = result_time_series;
%---------inter loop
end
end
end
|
本文转自二郎三郎博客园博客,原文链接:http://www.cnblogs.com/haore147/p/3795802.html,如需转载请自行联系原作者