plot_sgy.m

来自「对sgy格式的地震数据」· M 代码 · 共 28 行

function plot_sgy(filename)
%PLOT_SGY Plot sgy  data without reel
%
%filename:实际上是包括路径和文件名的字符串
%example:plot_sgy('d:\data\wyt_pc.sgy');

%From CUP
%Author: geoer,2007-7-11.21-07
   FID = fopen(filename,'r','n');
status = fseek(FID,114,'bof');
    ns = fread(FID,1,'short');%得到每道采样点数
  total=60+ns;%每道所占四字节个数
  frewind(FID);
  fseek(FID, 0, 'eof');
position=ftell(FID);
trace_num=position/(total*4);
 frewind(FID);
for k=1:trace_num
position=(total*(k-1)+60)*4;
status=fseek(FID,position,'bof');
if status==0
D(:,k)=fread(FID,ns,'float');
else
return;
end
end
fclose(FID);
wigb(D);