s_cplot.m

实现地震勘探中

if iscell(param.limits)
   param.limits=cat(2,param.limits{:});
end
  
ma=max(yi(~isnan(yi)));
if isempty(ma)
   iname=inputname(1);
   if strcmpi(iname,'')
      iname='data set';
   end
   alert([' All elements of ',iname,' are null values']);
   return
end
mi=min(yi(~isnan(yi)));
if ma*mi < 0
   ma=(ma-mi)*0.5;
%   ma=max(ma,-mi);
   mi=-ma;
elseif ma == mi
   ma=ma+10000*max(ma*eps,eps);
   mi=mi-10000*max(mi*eps,eps);
end
dmami=ma-mi;

if isempty(param.limits)
   if iscell(param.clim)
      param.clim=cat(2,param.clim{:});
   end
   param.limits(2)=ma-dmami*param.clim(2);
   param.limits(1)=mi+dmami*param.clim(1);
   if param.limits(1) >= param.limits(2)
      alert(' clim(1) + clim(2) must be less than 1; present values ignored')
      param.limits(1)=mi;
      param.limits(2)=ma;
   end
else
   if param.limits(1) >= param.limits(2)
      alert(' limits(1) must be less than limits(2); present values ignored')
      param.limits(1)=mi;
      param.limits(2)=ma;
   end
end

%       Define colormap to use
try
   cm=mycolormap(param.colormap);
catch %#ok
   disp(' Reqested colormap not found; seismic colormap used instead.')
   cm=mycolormap('seismic');
end

%       Flip color matrix if requested
if isyes(param.flip_colormap)
   cm=flipud(cm);
end 

%    Plot data
if isempty(param.shading)
   cplot_no1(yi,seismic.first,dti,xi,param.direction,cm, ...
         param.limits,ha)
else
   cplot_no2(yi,seismic.first,dti,xi,param.direction,cm, ...
         param.limits,ha,param.shading)
end


%  Title
if ~isempty(param.title)
   if iscell(param.title)		% Handle multi-line titles
      mytitle(param.title{1})
   else
      mytitle(param.title)
   end
end

%	Add annotation of horizontal axis
xtext=s_gd(seismic,param.annotation);
hunits=s_gu(seismic,param.annotation);

if ~isempty(xtext) && ~strcmpi(xtext,'not available')
   if ~isempty(hunits) && ~strcmpi(hunits,'n/a')
      xtext=[xtext,' (',hunits,')'];
   end
   xlabel([xtext,text2append])
end


if ~isempty(param.time_lines)
   if ~iscell(param.time_lines)
      param.time_lines={param.time_lines};
   end
   v=axis;
   width=1;
   for ii=1:length(param.time_lines)
      t1=ceil(v(3)/param.time_lines{ii})*param.time_lines{ii};
      temp=t1:param.time_lines{ii}:v(4);
      mygrid(temp,'h','k',width);
      width=width*1.5;
   end
end

box on
grid on 
zoom on
set(ha,'gridlinestyle','-','xgrid','off','Layer','top')

%	Button for figure-export menu
% figure_export_menu(figure_handle);

%	Button for scroll bars
cseismic_scrollbar_menu(figure_handle,seismic,param.direction)

%	Button for cursor tracking
if (isempty(param.tracking) && strcmpi(param.figure,'new'))  ||  ...
    strcmpi(param.tracking,'yes')       % Add cursor tracking
    
   [dummy,xinfo]=s_gh(seismic,param.annotation);   %#ok First output argument is not required
   y=linspace(seismic.first,seismic.last,nsamp);
   yinfo=info4time(seismic);
   initiate_3d_tracking(seismic.traces,xi,y,xinfo,yinfo,{'amplitude','','Amplitude'})
else
   yinfo=info4time(seismic);
   ylabel([yinfo{3},' (',yinfo{2},')'])
end

%	Button for menu to change image (color, etc. ...)
if strcmpi(param.imagemenu,'yes')
   myimagemenu  %	Create menu button to interactively change colors, etc.
end

%       Update colorbar
if isyes(param.colorbar)
   colorbar
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function cplot_no1(seis,stm,prm,xannot,plot_dir,cm,param_scale,ha)
% Function plot seismic data as image plot
%
%          cplot(seis,stm,prm,xannot,plot_dir,cm,param_scale,ha,colbar)
% INPUT
% seis   seismic traces
% stm    start time
% prm    sample interval
% xannot x-axis annotation
% plot_dir  plot direction ('l2r' or 'r2l'); default: 'l2r'
% cm     color matrix
% param_scale  a 2-element vector containing clow and chigh (see imagesc and colormap).
%
%        or by global variable param.clim which computes 
%        ma=max(seis(:)); 
%        if ma > 0
%           param_scale(2)=ma*param.clim(2);
%        else
%           param_scale(2)=ma/param.clim(2);
%        end
%        mi=min(seis(:)); 
%        if mi > 0
%           param_scale(1)=mi/param.clim(1);
%        else
%           param_scale(1)=mi*param.clim(1);
%        end
%        if param_scale(1) >= param_scale(2)
%           cc=mean(param_scale);
%           param_scale(1)=cc*(1-eps);
%           param_scale(2)=cc*(1+eps);
%        end
% ha     handle of axis

[nbin,ntr]=size(seis);

% 	Handle plot and CDP direction
if strcmp(plot_dir,'r2l') == 1 && xannot(1) < xannot(ntr), 
   flip_axis=1;
elseif strcmp(plot_dir,'r2l') ~= 1 && xannot(1) > xannot(ntr), 
   xannot=flipud(xannot(:));
   seis=fliplr(seis);
   flip_axis=1;
elseif strcmp(plot_dir,'r2l') == 1 && xannot(1) > xannot(ntr),
   xannot=flipud(xannot(:));
   flip_axis=0;
   seis=fliplr(seis);
else
   flip_axis=0;
end

%axis([xannot(1)-0.5,xannot(ntr)+0.5,stm,stm+(nbin-1)*prm])
if ~isempty(param_scale) && ~all(isconstant(seis))
   handle=imagesc(xannot,stm:prm:(nbin-1)*prm+stm,seis,param_scale);
else
   handle=imagesc(xannot,stm:prm:(nbin-1)*prm+stm,seis);
end 
set(handle,'Tag','image_displayed') % Create tag for image object which is  
                                    % used in "myimagemenu"
try
   colormap(cm),
catch  %#ok
   % Do nothing
end

set(ha,'XAxisLocation','top','YDir','Reverse');
if flip_axis == 1,
   set(ha,'XDir','Reverse');
end

drawnow

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function cplot_no2(seis,stm,prm,xannot,plot_dir,cm,param_scale,ha,param_shading)
% Function plot seismic data using "pcolor"
%
%          cplot(seis,stm,prm,xannot,plot_dir,cm,param_scale,ha,colbar)
% INPUT
% seis   seismic traces
% stm    start time
% prm    sample interval
% xannot x-axis annotation
% plot_dir  plot direction ('l2r' or 'r2l'); default: 'l2r'
% cm     color matrix
% param_scale  a 2-element vector containing clow and chigh (see imagesc and colormap).
%
%        or by global variable param.clim which computes 
%        ma=max(seis(:)); 
%        if ma > 0
%           param_scale(2)=ma*param.clim(2);
%        else
%           param_scale(2)=ma/param.clim(2);
%        end
%        mi=min(seis(:)); 
%        if mi > 0
%           param_scale(1)=mi/param.clim(1);
%        else
%           param_scale(1)=mi*param.clim(1);
%        end
%        if param_scale(1) >= param_scale(2)
%           cc=mean(param_scale);
%           param_scale(1)=cc*(1-eps);
%           param_scale(2)=cc*(1+eps);
%        end
% ha     handle of axis

[nbin,ntr]=size(seis);

% 	Handle plot and CDP direction
if strcmp(plot_dir,'r2l') == 1 && xannot(1) < xannot(ntr), 
   flip_axis=1;
   xannot=xannot(:);
elseif strcmp(plot_dir,'r2l') ~= 1 && xannot(1) > xannot(ntr), 
   xannot=flipud(xannot(:));
   seis=fliplr(seis);
   flip_axis=1;
elseif strcmp(plot_dir,'r2l') == 1 && xannot(1) > xannot(ntr),
   xannot=flipud(xannot(:));
   flip_axis=0;
   seis=fliplr(seis);
else
   flip_axis=0;
   xannot=xannot(:);
end

try
   colormap(cm),
catch  %#ok
   % Do nothing
end

%axis([xannot(1)-0.5,xannot(ntr)+0.5,stm,stm+(nbin-1)*prm])
seis=double(seis);
seis(end+1,end+1)=0;
seis(:,end)=seis(:,end-1);
% shading flat

handle=pcolor([0;xannot]+0.5,stm-prm/2:prm:(nbin-1)*prm+stm+prm/2,seis);
if ~isempty(param_scale)
   set(gca, 'CLim', param_scale);
end 
set(handle,'Tag','image_displayed') % Create tag for pplot object which is  
                                    % used in "myimagemenu"
try
   shading(param_shading)
catch  %#ok
   disp([' Unknown parameter for keyword "shading": ',param_shading])
end

set(ha,'XAxisLocation','top','YDir','Reverse');
if flip_axis == 1,
   set(ha,'XDir','Reverse');
end

drawnow