如何进行数据体处理

llllxxxxtttt

用Matlab进行三维数据的处理，将三维数据看成一个数据体M*N*K，每次取出N*K个数据，共取M次(相当于一个面一个面的取)，将整个数据体取完；然后将每个N*K点的图像，取其第 i （0<i<=M）列，循环完，将会得到一个第i列的N*K图像，然后处理这个图像，一直将M个N*K图像都取完，处理完；最后在把处理完的数据按照原来的顺序，还原成M*N*K数据体。下附我的程序，我的程序出现的问题是处理后排序和原来不一样，往大师指点，看看错误处在哪里。
fp=fopen('af_line.dat','r');
fp1=fopen('middle_af_crossline.dat','w');
img_rec=zeros(pt,line);
for j_trc=1:trc
    clc;disp('第二遍横测线（联络测线）处理中...');
    disp(j_trc);
    for i_line=1:line
        img=fread(fp,[pt trc],'float32');
        img_rec(:,i_line)=img(:,j_trc);
    end
    frewind(fp);
    img_rec=whole_domain_SVD(img,2,Energy_rate2);%img;
    fwrite(fp1,img_rec,'float32');

end
fclose(fp);
fclose(fp1);
%%%将数据体处理后的数据进行整理，按原始的数据体排列顺序展示%%%
fp=fopen('middle_af_crossline.dat','r');
fp1=fopen('af_crossline.dat','w');
% disp('数据整理中...');
for j_trc=1:trc
    clc;disp('数据整理中（please wait）...');
%     disp(j_trc);
    for i_line=1:line
        img=fread(fp,[pt trc],'float32');
        img_rec(:,i_line)=img(:,j_trc);
    end
    frewind(fp);
    fwrite(fp1,img_rec,'float32');
end
fclose(fp);
fclose(fp1);
disp('end2');

llllxxxxtttt

问题已解决！