header_index1.m

三谱混合相位子波估计！！ 这是我的一篇文章所涉及到的matlab 源代码

function [index,ier]=header_index1(seismic,mnem)
% Function outputs index of header in seismic structure "seismic".
% If no header mnemonic is found matching "mnem" then "index" is []
%
%         [index,ier]=header_index1(seismic,mnem)
% INPUT
% seismic  seismic structure whose headers are requested
% mnem   header mnemonics
% OUTPUT
% index  index of header with mnemonic "mnem"
% ier    error indicator; If only one output argument is given, the function will abort
%        with an error message if no header or more than one header is found
%        Otherwise: ier=0 if no error was detected and ier=1 if an error was found     

global S4M

if S4M.case_sensitive
   index=find(ismember(seismic.header_info(:,1),mnem));
else
   index=find(ismember(lower(seismic.header_info(:,1)),lower(mnem)));
end

%       Check for errors
if nargout < 2 
  if isempty(index)     % Print error message
    disp([' Header with mnemonic  "',mnem,'" not found'])
    disp(' The following header mnemonics exist: ')
    disp(seismic.header_info(:,1)');
    error(' Abnormal termination')

  elseif length(index) > 1     % Print error message
    disp([' More than one header with mnemonic "',mnem,'" found.'])
    disp(' The following header mnemonics exist: ')
    disp(seismic.header_info(:,1)');
    error(' Abnormal termination')

  else
    return
  end

else
   if isempty(index) | length(index) > 1
      ier=1;
   else
      ier=0;
   end
end