s_gh.m

基于Matlab的地震数据处理显示和测井数据显示于处理的小程序

function [header,info]=s_gh(seismic,mnem)
% Function gets header with mnmeonic "mnem" from seismic data set "seismic"
% If global variable "S4M.case_sensitive" is set to false, the case of the 
% header mnemonic is disregarded
%
% Written by: E. R., December 27, 2000
% Last updated: July 27, 2003; second output argument
%
%            [header,info]=s_gh(seismic,mnem)
% INPUT
% seismic    seismic data set
% mnem       header mnemonic
% OUTPUT
% header     row vector with header
% info       relevant row of "seismic.header_info" with mnemonic name, units of 
%            measurement, and description

global S4M

if iscell(mnem)
   mnem=mnem{1};
end

if strcmpi(mnem,'trace_no')       % Implied header "trace_no"
   header=1:size(seismic.traces,2);
   if nargout > 1
      info={'trace_no','','Trace number'};
   end
   return
end

mnems=seismic.header_info(:,1);

if S4M.case_sensitive
   idx=find(ismember(mnems,mnem));
else 
   idx=find(ismember(lower(mnems),lower(mnem)));
end

if ~isempty(idx) & length(idx) == 1
   header=seismic.headers(idx,:);
   if nargout == 2
      info=seismic.header_info(idx,:);
   end
   return
end

% Handle error condition
if isempty(idx)
   disp([' Requested header "',mnem,'" has not found. Available headers are:'])
   disp(mnems')
else
   disp([' More than one header found: ',cell2str(mnems(idx),', ')])
end

error(' Abnormal termination')