如何进行数据体处理
本帖最后由 llllxxxxtttt 于 2010-11-30 16:05 编辑用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,,'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,,'float32');
img_rec(:,i_line)=img(:,j_trc);
end
frewind(fp);
fwrite(fp1,img_rec,'float32');
end
fclose(fp);
fclose(fp1);
disp('end2'); 问题已解决!
页:
[1]