getcshot.m

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

function [shot]=getcshot(shotnum)
if(shotnum>0&shotnum<86)
  s=num2str(shotnum);
  s1=s;
  if(shotnum<10)
    s1='0';
    s1(2:2)=s;
  end
  filen='cshot00.dat';
  filen(6:7) = s1;
  fid1=fopen(filen,'r','b');
  shot1=fread(fid1,512*48,'float');
  shot=reshape(shot1,512,48);
  clear shot1;
  status=fclose(fid1);
else
  'Wrong shot number'
   break
end;