segy_struct.m

对sgy格式的地震数据

function [segy] = segy_struct;
%segy格式道头每个字的数据类型,利用了结构体数组来构建.数据类型是根据fread,fwrite的数据类型来确定的.
      segy.tracl= 'int';     %一条测线中的道顺序号
      segy.tracr= 'int';     %在本卷磁带中的道顺序号
       segy.fldr= 'int';     %原始的野外记录号
      segy.tracf= 'int';     %原始野外记录中的道号
         segy.ep= 'int';     %震源点号
        segy.cdp= 'int';     %cdp号
       segy.cdpt= 'int';     %cdp道集中的道号
       segy.trid= 'short';   %道识别号
        segy.nva= 'short';   %产生这一道的垂直叠加道数
        segy.nhs= 'short';   %产生这一道的水平叠加道数
       segy.duse= 'short';   %数据类型
     segy.offset= 'int';     %炮检距
      segy.gelev= 'int';     %接受点高程
      segy.selev= 'int';     %炮点高程
     segy.sdepth= 'int';     %炮点低于地面的深度
       segy.gdel= 'int';     %接受点的基准面高程
       segy.sdel= 'int';     %炮点的基准面高程
      segy.swdep= 'int';     %炮点水深
      segy.gwdep= 'int';     %接受点的水深
     segy.scalel= 'short';   %41-68字节的值根据这里给出的因子得到真值
     segy.scalco= 'short';   %73-88字节的值根据这里给出的因子得到真值
         segy.sx= 'int';     %炮点x坐标
         segy.sy= 'int';     %炮点y坐标
         segy.gx= 'int';     %检波点x坐标
         segy.gy= 'int';     %检波点y坐标
     segy.counit= 'short';   %坐标单位
      segy.wevel= 'short';   %风化层速度
     segy.swevel= 'short';   %风化层下的速度
        segy.sut= 'short';   %震源点处的井口时间
        segy.gut= 'short';   %接收点处的井口时间
      segy.sstat= 'short';   %炮点静校正
      segy.gstat= 'short';   %接收点的静校正
      segy.tstat= 'short';   %应用的总静校正量
       segy.laga= 'short';   %延迟时间A
       segy.lagb= 'short';   %延迟时间B
      segy.delrt= 'short';   %时间延迟时间
       segy.muts= 'short';   %起始切除时间
       segy.mute= 'short';   %结束切除时间
         segy.ns= 'unsigned short';   %本道采样点数-------------------------
         segy.dt= 'unsigned short';   %本道的采样间隔-----------------------
       segy.gain= 'short';   %野外仪器的增益类型
        segy.igc= 'short';   %仪器增益常数
        segy.igi= 'short';   %仪器起始增益
       segy.corr= 'short';   %相关码
        segy.sfs= 'short';   %起始扫描频率
        segy.sfe= 'short';   %结束扫描频率
       segy.slen= 'short';   %扫描长度
       segy.styp= 'short';   %扫描类型
       segy.stas= 'short';   %扫描道起始斜坡长度
       segy.stae= 'short';   %扫描道终止斜坡长度
      segy.tatyp= 'short';   %斜坡类型
      segy.afilf= 'short';   %滤假频的频率
      segy.afils= 'short';   %滤假频的陡度
     segy.nofilf= 'short';   %陷波频率
     segy.nofils= 'short';   %陷波陡度
        segy.lcf= 'short';   %低截频率
        segy.hcf= 'short';   %高截频率
        segy.lcs= 'short';   %低截频率陡度
        segy.hcs= 'short';   %高截频率陡度
       segy.year= 'short';   %数据记录的年
        segy.day= 'short';   %数据记录的日
       segy.hour= 'short';   %数据记录的小时
     segy.minute= 'short';   %数据记录的分
        segy.sec= 'short';   %数据记录的秒
     segy.timbas= 'short';   %时间代码
       segy.trwf= 'short';   %道加权因子
     segy.grnors= 'short';   %覆盖开关位置1的检波器道号
     segy.grnofr= 'short';   %在野外原始记录中道号1的检波器号
     segy.grnlof= 'short';   %在野外原始记录中最后一道的检波器号
       segy.gaps= 'short';   %缺口大小
      segy.otrav= 'short';   %在测线的开始或结束处的斜坡位置
         segy.d1= 'float';   %181---240为自定义,以下不知道含义
         segy.f1= 'float';   %
         segy.d2= 'float';   %
         segy.f2= 'float';   %
     segy.ungpow= 'float';   %
    segy.unscale= 'float';   %
        segy.ntr= 'int';     %
       segy.mark= 'short';   %
      segy.unass= 'short';   %
      segy.trace= 'float';   %读取整道时，这个数据类型对么？