- %这个程序用来绘制一个放在n行m列的数组中的数据,***原始数据顺序***是1列,2列,......m列,文件名为filename
- clear;%清除所有工程中变量
- filename='seismicZ.H@';%waveform_record_cal_1 waveform_record_obs
- n=401;%n行
- m=46;%m列
- noshot=1;
- fid=fopen(filename,'r');%打开data文件,把地址付给变量fid,原始数据顺序是1行,2行,......n行
- for i=1:noshot-1;
- nmdata=fread(fid,[n,m],'float32');%把数据(二进制格式)付给m行n列数组
- end;
- nmdata=fread(fid,[n,m],'float32');
- imagesc(nmdata);