读取sg2地震数据程序

clc clear all fid01 = fopen('20001.sg2', 'r');% Open the seg2 file for reading. if fid01 == -1 error('Unable to open file.') end id= fread(fid01,1,'int16');%文件描述块ID char_id = dec2hex (id) %3A55 bb=fread(fid01,1,'int16');%所使用得SEG标准版本号 M = fread(fid01,1,'uint16');%从文件描述块中读取道指针子块大小 N = fread(fid01,1,'uint16');%从文件描述块中读取采样道数N fread(fid01,1,'uint8'); %为01h或02h，这主要取决于第9，10字节所存放的字符串结束符个数 fread(fid01,1,'uint16');%字符串结束符，其将文件描述块和道描述快中得多个字符串分开 fread(fid01,1,'uint8');%行结束符得大小（01h或02h） fread(fid01,1,'uint16');%行结束符 fread(fid01,9,'uint16');%SEG02 保留部分 sym Pointers for i=1:N Pointers(i)=fread(fid01,1,'int32'); end for i=1:N fread(fid01);%读数指针拨至文件开始位置 fseek(fid01,Pointers(i),'bof');%将读书指针拨至第i道得道描述块起始位置 temp= fread(fid01,1,'int16'); %读取道描述块得0，1字节 char_temp= dec2hex(temp); %4422 B_zise= fread(fid01, 1 ,'uint15'); %读取道描述块得大小 D_zise= fread(fid01, 1 ,'ulong'); %读取数据块字节的大小 S_number= fread(fid01, 1 ,'ulong'); %读取数据块中采样点数（S number） D_format= fread(fid01, 1 ,'uchar'); %读取数据存储格式代码 chx = dec2hex(D_format); fread(fid01,19,'uchar'); fseek(fid01, B_zise-32,'cof');%将读数指针拨至数据块的起始位置 for j= 1: S_number %for j= 1: 1 Data(i,j)= fread (fid01,1,'float');%该数组为数据存储数组 end end fclose (fid01); %文件存储部分，几种方式，根据需要选择 %波形变面积显示地震记录 dx=10; %道间距 dt= 0.0025;%采样间隔 zy=0.5; %增益 x= 0: dx: (N-1)*dt; z= 0: dt: (S_number-1)*dt; %wigb(Data,zy,x,z); %子函数 %------------------------------------- %