l_append.m

地震、测井方面matlab代码,解释的比较详细

function wlog=l_append(wlog1,wlog2,varargin)
% Function appends one log structure to another. 
% The following cases are handled (or planned).
%    1. Both logs have the same depth values: 
%       curves on the second log which are identical to curves on the first log are 
%       not copied.
%       Only parameters of the second log with names different from those in the 
%       first log are copied to the output log.  
%    2. If this is not the case both logs are uniformly resampled. Unless given explicitely,
%       the sample interval is taken from the first log. If the first log is not sampled uniformly,
%       it is taken from the second log. If the second log is not sampled uniformly either, an error 
%       message is issued.
%       The string 'no_1' is appended to mnemonics of curves of the second log structure
%       which equal those of the first log structure. The same is true for parameters.
%       All other fields (such as wall name, country) are taken from the first log.
%       However, fields of the second log that do not exist in the first log are also
%       copied to the output log.
% Written by: E. R., December 26, 2000
% Last updated: March17, 2005: combine fields 'curve_types' of teh two well logs
%
%         wlog=l_append(wlog1,wlog2,varargin)
% INPUT
% wlog1    first input log 
% wlog2    log to be appended
% varargin    one or more cell arrays; the first element of each cell array is a keyword,
%             the other elements are parameters. Presently, keywords are:
%     'step'  stepsize of output log.
%             default: {'step',wlog1.step}  i.e. depth step of first log (must not be 0)
%     'units' depth units. Default: {'units',wlog1.curve_info{1,2}} i.e. units of first input log
%     
% OUTPUT 
% wlog  combination of the two input logs

global S4M

[nsamp1,ncurves1]=size(wlog1.curves);
[nsamp2,ncurves2]=size(wlog2.curves);

wlog.type='well_log';
wlog.tag='unspecified';
wlog.name=wlog1.name;

% 	First handle the case where the depth values und units of measurement of the two logs are identical

if nsamp1 == nsamp2 & ...
      strcmp(wlog1.curve_info{1,2},wlog1.curve_info{1,2}) &  ...
      max(abs(wlog1.curves(:,1)-wlog2.curves(:,2))) < 1.0e-5*(wlog1.last-wlog1.first)/nsamp1;        

  index=curve_indices(wlog2,wlog1.curve_info(2:end,1));
  if isempty(index)
     wlog.curves=[wlog1.curves,wlog2.curves];
     wlog.curve_info=[wlog1.curve_info;wlog2.curve_info];

  else
     nindex=find(~ismember(1:ncurves2,index));
     nindex=nindex(2:end);     
     wlog.curves=[wlog1.curves,wlog2.curves(:,nindex)];
     wlog.curve_info=[wlog1.curve_info;wlog2.curve_info(nindex,:)];
     for ii=index;
       mnem=wlog2.curve_info{ii,1};
       curve1=l_gc(wlog1,mnem);
       curve2=l_gc(wlog2,mnem);
       if max(abs(curve1-curve2)) <= 1.0e-6*max(curve1)
          alert([' Alert from "l_append": curve "',mnem,'" of second data set disregarded'])
          alert(['                        (same as curve "',mnem,'" of first data set'])
       else
          s_header(wlog,'add',[mnem,'_no_1'],curve2,wlog2.curve_info{ii,2},wlog2.curve_info{ii,3})
       end
     end
  end

%       Copy rest of fields
  wlog=copy_fields(wlog1,wlog);
  wlog=copy_fields(wlog2,wlog);
  return
end

% 	If depth values are different resample at least second data set

%     Set default values
append.step=wlog1.step;
append.units=wlog1.curve_info{1,2};

%       Decode and assign input arguments
append=assign_input(append,varargin);
if append.step == 0
  if wlog2.step == 0
     error(' Sample interval "step" must not be zero')
  end
  append.step=wlog2.step;
end

%       Check depth units
if ~strcmp(append.units,'m') & ~strcmp(append.units,'ft') & ~strcmp(append.units,'ms')
  error([' Unknown depth units: ',append.units])
end

%       Convert depth if necessary
if ~strcmp(append.units,wlog1.curve_info{1,2})
   if strcmp(wlog1.curve_info{1,2},'m')
      wlog1.curves(:,1)=wlog1.curves(:,1)/0.3048;
      wlog1.first=wlog1.first/0.3048;
      wlog1.last=wlog1.last/0.3048;
      wlog1.step=wlog1.step/0.3048;
   elseif strcmp(wlog1.curve_info{1,2},'ft')
      wlog1.curves(:,1)=wlog1.curves(:,1)*0.3048;
      wlog1.first=wlog1.first*0.3048;
      wlog1.last=wlog1.last*0.3048;
      wlog1.step=wlog1.step*0.3048;
   else
      error([' First input log has unknown depth units: ',wlog1.curve_info{1,2}])
   end
end

if ~strcmp(append.units,wlog2.curve_info{1,2})
   if strcmp(wlog2.curve_info{1,2},'m')
      wlog2.curves(:,1)=wlog2.curves(:,1)/0.3048;
      wlog2.first=wlog2.first/0.3048;
      wlog2.last=wlog2.last/0.3048;
      wlog2.step=wlog2.step/0.3048;
   elseif strcmp(wlog2.curve_info{1,2},'ft')
      wlog2.curves(:,1)=wlog2.curves(:,1)*0.3048;
      wlog2.first=wlog2.first*0.3048;
      wlog2.last=wlog2.last*0.3048;
      wlog2.step=wlog2.step*0.3048;
   else
      error([' Second input log has unknown depth units: ',wlog2.curve_info{1,2}])
   end
end

% 	Define new start and end depth
first=min([wlog1.first,wlog2.first]);
last=max([wlog1.last,wlog2.last]);
wlog.first=ceil(first/append.step)*append.step;
wlog.last=floor(last/append.step)*append.step;
wlog.step=append.step;
nsamp=round((wlog.last-wlog.first)/wlog.step)+1;
wlog.curves=zeros(nsamp,ncurves1+ncurves2-1);
depth=(0:nsamp-1)'*(wlog.last-wlog.first)/(nsamp-1)+wlog.first;

if wlog1.step > 0
   wlog.curves(:,2:ncurves1)=interp1(wlog1.curves(:,1),wlog1.curves(:,2:ncurves1),depth,'*linear');
else
   wlog.curves(:,2:ncurves1)=interp1q(wlog1.curves(:,1),wlog1.curves(:,2:ncurves1),depth);
end

if wlog2.step > 0
   wlog.curves(:,ncurves1+1:ncurves1+ncurves2-1) = ...
                       interp1(wlog2.curves(:,1),wlog2.curves(:,2:ncurves2),depth,'*linear');
else
   wlog.curves(:,ncurves1+1:ncurves1+ncurves2-1) = ...
                       interp1q(wlog2.curves(:,1),wlog2.curves(:,2:ncurves2),depth);
end
wlog.curves(:,1)=depth;

%       Handle curve information field
curve_info=wlog2.curve_info(2:end,:);

if S4M.case_sensitive
   idx=find(ismember(curve_info(:,1),wlog1.curve_info(:,1)));
else
   idx=find(ismember(lower(curve_info(:,1)),lower(wlog1.curve_info(:,1))));
end

if ~isempty(idx)
   for ii=idx(:)'
     curve_info{ii,1}=[curve_info{ii,1},'_no_1'];
   end
end

%       Combine 'curve_info'
wlog.curve_info=[wlog1.curve_info;curve_info];
wlog.curve_info{1,2}=append.units;


%       Handle 'curve_types'
if isfield(wlog1,'curve_types')
   if isfield(wlog2,'curve_types')
      curve_types=wlog2.curve_types;
      if S4M.case_sensitive
         idx=find(ismember(curve_types(:,1),wlog1.curve_types(:,1)));
      else
         idx=find(ismember(lower(curve_types(:,1)),lower(wlog1.curve_types(:,1))));
      end

      if ~isempty(idx)
         curve_types(idx,:)=[];
      end

%               Combine 'curve_types'
      wlog.curve_types=[wlog1.curve_types;curve_types];
      wlog.curve_info{1,2}=append.units;

   else
      wlog.curve_types=wlog1.curve_types;
   end

elseif  isfield(wlog2,'curve_types')
   wlog.curve_types=wlog2.curve_types;
end


%       Handle parameters
wlog=combine_parameters(wlog1,wlog2,wlog);

%       Copy rest of fields
wlog=copy_fields(wlog1,wlog);
wlog=copy_fields(wlog2,wlog);

		if 0

if isfield(wlog1,'parameter_info') & isfield(wlog2,'parameter_info')
  par1=wlog1.parameter_info(:,1);
  par2=wlog2.parameter_info(:,1);
  par2m=par2;
  par2_info=wlog2.parameter_info;
  idx=find(ismember(par2,par1));
  if ~isempty(idx)
     for ii=idx(:)'
        par2m{ii}=[par2m{ii},'_no_1'];
        par2_info{ii,1}=par2m{ii};
     end
  end
  wlog.parameter_info=[wlog.parameter_info;par2_info];

  for ii=1:length(par2)
 %   wlog.parameter_info=setfield(wlog.parameter_info,par2m{ii}, ...
 %                        getfield(wlog2.parameter_info,par2{ii}));
 %   wlog.parameter_info=[wlog.parameter_info;
    wlog=setfield(wlog,par2m{ii},getfield(wlog2,par2{ii}));
  end 
end
		end

%       Check for NaNs in the combined structure
if ~isfield(wlog,'null')
  index=find(isnan(wlog.curves(:,2:end)));
  if ~isempty(index)
     wlog.null=NaN;
  end
end  

wlog=l_rm_nulls(wlog,'all');