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

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，如需转载请自行联系原作者