wavematlab.txt

Matlab地震波显示程序 Matlab的矩阵数据类型与地震波数据类似

代码都很简单，比如第一个就是matlab中低级i/o操作，如下示例，读取无卷头的sgy文件。有卷头的，加上3600即可。欢迎讨论弹性波有限元模拟等内容。
%主程序
global file;
clf reset;
set(gcf,'MenuBar','none','NumberTitle','off','Name','地震数据显示工具');
uimenu(gcf,'Label','选择地震数据','CallBack',['[FileName,PathName]=uigetfile(''*.sgy'',''读入地震数据'');', ...
'file=strcat(PathName,FileName);']);
uimenu(gcf,'Label','绘制剖面','CallBack',@plot_sgy);
uimenu(gcf,'Label','关闭','CallBack','close');
--------------------------------------------------------------------------------------------------------------------------------------------
function plot_sgy(src,event)
%PLOT_SGY Plot sgy  data without reel
%
%file:实际上是包括路径和文件名的字符串
%

%From CUP
%Author: geoer,2007-7-11
global file;
   FID = fopen(file,'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);