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