地震数据segy读取的MATLAB函数

由Thomas Mejer Hansen编写，读取Segy数据，提取道头和卷头信息，以及GUI包

% SegyTraceHeader=CheckSegyTraceHeader（SegyTraceHeader）;
%
% Checks that all fields of the SegyTraceHeader is set. 
% If not they are initialized.
%

%
% （C） 2004 Thomas Mejer Hansen tmh@gfy.ku.dk/thomas@cultpenguin.com
%
%
%
%    This program is free software; you can redistribute it and/or modify
%    it under the terms of the GNU General Public License as published by
%    the Free Software Foundation; either version 2 of the License or
%    （at your option） any later version.
%
%    This program is distributed in the hope that it will be useful
%    but WITHOUT ANY WARRANTY; without even the implied warranty of
%    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%    GNU General Public License for more details.
%
%    You should have received a copy of the GNU General Public License
%    along with this program; if not write to the Free Software
%    Foundation Inc. 59 Temple Place Suite 330 Boston MA  02111-1307  USA
%
function SegyTraceHeader=CheckSegyTraceHeader（SegyTraceHeader）;
%% SET UP SegyTraceHeader
  NA=0;
  
  if ~isfield（SegyTraceHeader‘TraceSequenceLine‘） SegyTraceHeader.TraceSequenceLine=NA;end
  if ~isfield（SegyTraceHeader‘TraceSequenceFile‘） SegyTraceHeader.TraceSequenceFile=NA;end
  if ~isfield（SegyTraceHeader‘FieldRecord‘） SegyTraceHeader.FieldRecord=NA;end
  if ~isfield（SegyTraceHeader‘TraceNumber‘） SegyTraceHeader.TraceNumber=NA;end
  if ~isfield（SegyTraceHeader‘EnergySourcePoint‘） SegyTraceHeader.EnergySourcePoint=NA;end
  if ~isfield（SegyTraceHeader‘cdp‘） SegyTraceHeader.cdp=NA;end
  if ~isfield（SegyTraceHeader‘cdpTrace‘） SegyTraceHeader.cdpTrace=NA;end
  if ~isfield（SegyTraceHeader‘TraceIdenitifactionCode‘） SegyTraceHeader.TraceIdenitifactionCode=NA;end
  if ~isfield（SegyTraceHeader‘NSummedTraces‘） SegyTraceHeader.NSummedTraces=NA;end
  if ~isfield（SegyTraceHeader‘NStackedTraces‘） SegyTraceHeader.NStackedTraces=NA;end
  if ~isfield（SegyTraceHeader‘DataUse‘） SegyTraceHeader.DataUse=NA;end
  if ~isfield（SegyTraceHeader‘offset‘） SegyTraceHeader.offset=NA;end
  if ~isfield（SegyTraceHeader‘ReceiverGroupElevation‘） SegyTraceHeader.ReceiverGroupElevation=NA;end
  if ~isfield（SegyTraceHeader‘SourceSurfaceElevation‘） SegyTraceHeader.SourceSurfaceElevation=NA;end
  if ~isfield（SegyTraceHeader‘SourceDepth‘） SegyTraceHeader.SourceDepth=NA;end
  if ~isfield（SegyTraceHeader‘ReceiverDatumElevation‘） SegyTraceHeader.ReceiverDatumElevation=NA;end
  if ~isfield（SegyTraceHeader‘SourceDatumElevation‘） SegyTraceHeader.SourceDatumElevation=NA;end
  if ~isfield（SegyTraceHeader‘SourceWaterDepth‘） SegyTraceHeader.SourceWaterDepth=NA;end
  if ~isfield（SegyTraceHeader‘GroupWaterDepth‘） SegyTraceHeader.GroupWaterDepth=NA;end
  if ~isfield（SegyTraceHeader‘ElevationScalar‘） SegyTraceHeader.ElevationScalar=NA;end
  % Multiply/divide next number for following 4 values
  if ~isfield（SegyTraceHeader‘SourceGroupScalar‘） SegyTraceHeader.SourceGroupScalar=NA;end
  if ~isfield（SegyTraceHead

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-09-17 16:48  SegyMAT\
     文件        9418  2004-07-26 15:15  SegyMAT\CheckSegyTraceHeader.m
     文件        2900  2004-07-26 15:15  SegyMAT\Contents.m
     文件       17656  2002-05-22 11:40  SegyMAT\DIAAboutSegymat.fig
     文件        3129  2002-05-06 10:40  SegyMAT\DIAAboutSegymat.m
     文件       36968  2002-06-16 12:04  SegyMAT\DIAOpenSEGY.fig
     文件       10779  2002-09-17 16:27  SegyMAT\DIAOpenSEGY.m
     目录           0  2016-09-17 16:48  SegyMAT\GUI\
     文件          19  2004-01-11 16:44  SegyMAT\GUI\compile_segymat.m
     文件       44896  2003-10-24 20:40  SegyMAT\GUI\segymat.fig
     文件       24149  2004-02-19 16:15  SegyMAT\GUI\segymat.m
     文件       28264  2003-10-24 16:52  SegyMAT\GUI\segymat_about.fig
     文件        8270  2003-10-24 16:52  SegyMAT\GUI\segymat_about.m
     文件        9360  2005-08-27 20:54  SegyMAT\GUIEditSegyHeader.fig
     文件       13424  2003-03-03 16:09  SegyMAT\GUIEditSegyHeader.m
     文件       31344  2002-06-05 11:43  SegyMAT\GUIEditSegyTraceHeader.fig
     文件       12093  2004-12-03 00:40  SegyMAT\GUIEditSegyTraceHeader.m
     文件        9856  2002-09-17 16:28  SegyMAT\GUIEditTextualFileHeader.fig
     文件        7200  2003-03-03 16:09  SegyMAT\GUIEditTextualFileHeader.m
     文件       18648  2002-06-16 11:44  SegyMAT\GUIPlotXY.fig
     文件        4418  2004-07-26 15:15  SegyMAT\GUIPlotXY.m
     文件       59528  2002-06-12 15:22  SegyMAT\GUISegyMAT.fig
     文件       16830  2004-07-26 15:15  SegyMAT\GUISegyMAT.m
     文件        7824  2006-09-18 14:54  SegyMAT\GetSegyHeader.m
     文件       10069  2006-09-18 14:54  SegyMAT\GetSegyHeaderBasics.m
     文件        1451  2004-07-26 15:15  SegyMAT\GetSegyTrace.m
     文件        2761  2004-07-26 15:15  SegyMAT\GetSegyTraceData.m
     文件        9427  2004-07-26 15:15  SegyMAT\GetSegyTraceHeader.m
     文件        1279  2004-07-26 15:15  SegyMAT\GetSegyTraceHeaderInfo.m
     文件        4840  2004-07-26 15:15  SegyMAT\InitSegyTraceHeader.m
     文件       24389  2005-08-27 20:54  SegyMAT\LICENSE
............此处省略40个文件信息