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