getoffset.m

来自「这是matlab在地球物理数据处理方面的源码」· M 代码 · 共 21 行

% *******************************************
% *** GetOffset                           ***
% *******************************************
%
      function [offset]=GetOffset(filename,trace1,trace2);
      nt=1250;
      nx=48; 
      fid=fopen(filename,'r','l');
      fseek(fid,3600+(trace1-1)*(nt+60)*4,'bof');
      m=0;
      for i=trace1:trace2
        m=m+1;
        fseek(fid,36,'cof');
        offset(m)=fread(fid,1,'int32');
        fseek(fid,200+nt*4 ,'cof');
      end
      fclose(fid)
      return
%
% *** end of program ***