大脑提取每一个体素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