plot_sgyr.m

对sgy格式的地震数据

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

%From CUP
%Author: geoer,2007-7-12.11-12
   FID = fopen(filename,'r','n');
status = fseek(FID,3714,'bof');
    ns = fread(FID,1,'short');%得到每道采样点数
  total=60+ns;%每道所占四字节个数
  frewind(FID);
  fseek(FID, 0, 'eof');
position=ftell(FID);
trace_num=(position-3600)/(total*4);
 frewind(FID);
for k=1:trace_num
position=3600+(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);