s_ispectrum.m

实现地震勘探中

function table=s_ispectrum(seismic,varargin)
% Function computes and displays the spectra of seismic data in interactively
% selected windows. Null values in the input data set are replaced by zeros
% prior to spectrum computation.
%
% Written by: E. R.: January 1, 2001
% Last updated: December 6, 2007: add time-range selection parameter; 
%                                 change parameter "taper" to "window"
% 
%           s_ispectrum{seismic,varargin}
% INPUT
% seismic   seismic data set
% varargin  one or more cell arrays; the first element of each cell array is a
%           keyword, the other elements are parameters. Presently, keywords are:
%     'annotation'  Header mnemonic used to annotate the horizontal axis
%                  Default: {'annotation','trace_no'}
%     'colors'     Colors to be used for consecutive curves.
%                  Possible values: any permissible colors and line styles
%                  Default: {'colors','r','g','b','m','k','c','y',...
%                            'r--','g--','b--','m--','k--','c--','y--' ...
%                            'r:','g:','b:','m:','k:','c:','y:'};
%     'frequencies'  Two positive numbers representing the range of frequencies to 
%                  display. The first number must be non-negative and smaller than 
%                  the second. If the second number is greater than the Nyquist 
%                  frequency of the data set with the smallest sample interval, it 
%                  is set to the Nyquist frequency.
%                  Default: {'frequencies',0,inf}
%     'legend'     Figure legend (curve annotation).
%                  Default: names of the seismic input data sets.
%     'lloc'       Location of figure legend. Possible values are 1,2,3,4,5;
%                  Default: {'lloc',5}
%     'linewidth'  Line width of curves. Default: {'linewidth',2}           
%                  Default: {'plot','amp'} ... plot amplitude spectrum only
%     'minamp'     Minimum amplitude in dB to display if keyword "scale" is
%                  set to 'log'. Must not be larger than -10. 
%                  Default: {'minamp',-50}
%     'normalize'  Establish if the amplitude spectra are to be normalized.
%                  Possible values: 'yes' and 'no'. Default: {'normalize','yes'}
%     'padding'    Traces with fewer than "padding" samples are padded with
%                  zeros. This parameter is ignored if the number of samples 
%                  per trace exceed "padding". Default:{'padding',256}
%     'plottype'   Type of seismic plot to use. Options are 'color' and 'wiggle'.
%                  Default: {'plottype','wiggle') for fewer than 
%                                                 "S4M.ntr_wiggle2color" traces
%                           {'plottype','color'}  for "S4M.ntr_wiggle2color"
%                                                 or more traces
%     'scale'      Set linear or logarithmic scale (dB) for amplitude spectrum.
%                  Possible values: 'linear', 'log'. Default: {'scale','linear'}
%     'times'      2-element or 3-element cell array; time range to plot
%                  {'times',vector of first and last time to plot} or ('times',first,last}. 
%                  Default: {'times',seismic.first,seismic.last}  which is
%                          equivalent to {'times',[seismic.first,seismic.last]}
%     'traces'     2-element or 3-element cell array. The second element can
%                  be an array of trace numbers or it can be a string.
%                  If it is a string it can be a header mnemonic or it 
%                  can contain a logical expression involving header values 
%                  to include. A "pseudo-header" 'trace_no' can also be used.
%                  If the second element is a string containing a header 
%                  mnemonic there must be a third element containing a 
%                  vector of values. (see "s_select")
%                  Default:  {'traces',[]}      which is equivalent to 
%                            {'traces',1:ntr}   where ntr denotes the 
%                                  number of traces in the input data set 
%                                  (ntr = size(seismic.traces,2))
%     'window'      Window applied to the data before computing the FFT
%           Possible names are :
%            	'Hamming', 'Hanning', 'Nuttall',  'Papoulis', 'Harris'
%            	'Rect',    'Triang',  'Bartlett', 'BartHann', 'Blackman'
%            	'Gauss',   'Parzen',  'Kaiser',   'Dolph',    'Hanna'
%            	'Nutbess', 'spline'
%                  Default: {'window','rect'}   (no window)
% OUTPUT
% table  table structure with corner points of windows used;
%        the table is set to [] if no windows are picked.
%        Also, in this case, the global variable ABORTED is set to "true";
%        it is set to false, otherwise.
%        If an output data set is specified the program will stop and only 
%        resume when picking is complete.
%
% EXAMPLES
%        seismic=s_data;
%        table=s_ispectrum(seismic,{'frequencies',0,80},{'padding',128},{'annotation','cdp'})

% UPDATE HISTORY
%        August 28, 2006: add trace selection parameter


global ABORTED S4M 

ABORTED=true;

%       Set defaults for input parameters
param.annotation='trace_no';
param.colors={'r','g','b','k','c','m','y', ...
              'r--','g--','b--','k--','c--','m--','y--', ...
              'r:','g:','b:','k:','c:','m:','y:'};
param.linewidth=3;
param.lloc=1;
param.minamp=-50;
param.normalize='yes';
param.frequencies={0,inf};
param.padding=128;
param.plot='amp';
param.plottype=[];
param.scale='linear';
param.taper=[];
param.times=[seismic.first,seismic.last];
param.traces=[];
param.window='rect';

%       Replace defaults by actual input arguments
param=assign_input(param,varargin);

%       Check deprecated parameter
if ~isempty(param.taper)
   alert('Parameter "taper" has been deprecated. Use "window" instead.')
   param.window=param.taper;
end

%       Select subset of the data (if requested)
if ~isempty(param.traces) || ~isempty(param.times)

   if iscell(param.traces)
      param.traces=cell2num(param.traces);
   end

   seismic=s_select(seismic,{'traces',param.traces},{'times',param.times(1),param.times(2)});

%   ntr=size(seismic.traces,2);
end

%       Determine if the seismic data wipp be represented by a wiggle plot or
%       by a color plot 
if isempty(param.plottype)
   if size(seismic.traces,2) <= S4M.ntr_wiggle2color
      param.plottype='wiggle';
   else
      param.plottype='color';
   end
else
   if ~ismember(param.plottype,{'wiggle','color'})
      error('Parameter "plottype" must be either "wiggle" or "color".')
   end
end

if param.minamp > -10
   alert('Parameter for keyword "minamp" must not be greater than -10.')
   param.minamp=-10;
end

if strcmpi(param.plottype,'color')
   aux=s_cplot(seismic,{'annotation',param.annotation},{'tracking','no'});
elseif strcmpi(param.plottype,'wiggle')
   aux=s_wplot(seismic,{'annotation',param.annotation},{'tracking','no'});
else
   error([' Unknown plot type: ',param.plottype])
end
SeisHandle=aux.figure_handle;
set(SeisHandle,'DoubleBuffer','on');

userdata=struct('SeisHandle',SeisHandle);
userdata.param=param;
userdata.seismic=seismic;
userdata.exit=false;
set(SeisHandle,'UserData',userdata)

%	Create a menu button that allows interactive window picking
menuHandle=menu2pick_frequency_windows(SeisHandle);

%	Create a help button
create_help_button('s_ispectrum')

if nargout > 0
   waitfor(menuHandle)
   if ishandle(SeisHandle)
      userdata=get(SeisHandle,'UserData');
   else
      userdata=[];
   end
end

if isfield(userdata,'exit')  &&  userdata.exit
   try2delete(SeisHandle)
end

if nargout > 0
   try 
      temp=userdata.SpecWindows;
   catch
      temp=[];
   end
   if ~isempty(temp)

      ABORTED=false;
      table.type='table';
      table.tag='unspecified';
      table.format='loose';
      table.name='Windows for spectrum estimation';
      table.column_info={'left',s_gu(seismic,param.annotation),param.annotation;
                      'right',s_gu(seismic,param.annotation),param.annotation;
                      'first',seismic.units,'Window start time';
                      'last',seismic.units,'Window end time'};
      if strcmpi(param.plottype,'wiggle')
         headers=s_gh(seismic,param.annotation);
	 table.left=headers(temp(:,1));
         table.right=headers(temp(:,2));
      else
         table.left=temp(:,1);
         table.right=temp(:,2);
      end
      table.first=temp(:,3);
      table.last=temp(:,4);
   else
      table=[];
   end
end