l_gc.m

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

function [curve,info]=l_gc(wlog,mnem,depths)
% Function extracts curve with mnemonic "mnem" from log data set "wlog"
% If global variable S4M.case_sensitive is set to false, the case of the curve
% mnemonic is disregarded
% Written by: E. R., December 28, 2000
% Last updated: July 27, 2003; second output argument
%
%           [curve,info]=l_gc(wlog,mnem,depths)
% INPUT
% wlog      log structure
% mnem      curve mnemonic
% depths    vector of depth values for which the curve values are required;
%           these values need not coincide with depth values in "wlog".
%           if not given, all samples of curve are output
% OUTPUT
% curve     column vector with curve
% info      relevant row of "wlog.curve_info" with mnemonic name, units of 
%           measurement, and description

global S4M

mnems=wlog.curve_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
   if nargin == 2
      curve=wlog.curves(:,idx);
   else
      if wlog.step > 0
         curve=interp1(wlog.curves(:,1),wlog.curves(:,idx),depths,'*linear');
      else
         curve=interp1q(wlog.curves(:,1),wlog.curves(:,idx),depths);
      end
   end
   if nargout == 2
      info=wlog.curve_info(idx,:);
   end
   return
end

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

error(' Abnormal termination')