# 大脑提取每一个体素26领域的matlab代码

%-------------- 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


posted @ 2014-06-19 09:24  二郎那个三郎  阅读(770)  评论(0编辑  收藏  举报