Code covered by the BSD License  

Highlights from
SeisLab 3.01

image thumbnail

SeisLab 3.01


Eike Rietsch (view profile)


24 Jul 2007 (Updated )

Release 10.0301 of a collection of functions for analysis/display of seismic and well-log data.

% Examples4SeismicWigglePlots
% The script illustrates the use of some of the parameters available to control
% a wiggle-trace plot

keep WF
global S4M    %#ok

%%     Get data for the examples

%% 1  All kewords are set to defaults
mytitle('Example 1: Default settings')

%% 2  Create a high-quality plot with gray background 
%    (setting the quality parameter to 'high' prevents vertical zero-deflection
%    lines which show up under some circumstances --- in particular when 
%    saving or printing the plot; it is much slower than the draft-quality plot)
s_wplot(seismic,{'quality','high'},{'background','gray'}, ...
                {'title','Example 2: high-quality and gray background'})

%% 3  Use the "quality" keyword to create a spike display 
%    (intended to display reflection coefficients) and background color defind
%    via an RGB triplet; the thickness of the spikes can be controlled via 
%    keyword "wiggle_width".
s_wplot(seismic,{'quality','spikes'},{'background',[0.25,0.25,0.25]}, ...
                {'title','Example 3: Spike display with dark-gray background'})

%% 4  Trace annotation by CDP number, but only annotate traces 1, 3, 7,
%     8, and 12; this assumes that the header "CDP" exists
s_wplot(seismic,{'annotation','CDP'},{'aindex',[1 3 7 8 12]}, ...
                {'title','Example 4: Trace annotation by CDP number'})
%% 5   No wiggle but peak filled read and trough filled gray
s_wplot(seismic,{'peak_fill','red'},{'trough_fill','gray'},{'wiggle_color',[]}, ...
                {'title','Example 5: Peak and trough fill but no wiggle_color'})

%% 6  Top plot: traces plotted from right to left (default is left to right);
%    it uses a trace deflection of 0.9 traces spacings; The font
%    size of the axis annotation etc. is 13 points.
%    Bottom plot: traces plotted from left to right (default); 
%    the plot uses the output of the first plot function to display the 
%    the traces with the same trace deflection; it is wiggle-only and 
%    uses the default 11-point font.
pfigure        % Plot figure in portrait format
   aux=s_wplot(seismic,{'direction','r2l'},{'deflection',0.9},{'figure','old'}, ...
   s_wplot(seismic,{'direction','l2r'},{'scale',aux.scale},{'figure','old'}, ...
mysuptitle('Example 6: trace deflection is 0.9 trace spacings',{'factor',0.93})

%% 7a Plot the 5-th and the 8-th trace in red (brute force approach)
bool=false(1,12);                % The seismic dataset has 12 traces
temp.traces(:,bool)=NaN;         % Set traces 5 and 7 to NaN
temp.null=NaN;                   % Dataset has null values (this is only done to avoid warnings)
s_wplot(temp,{'deflection',1})   % Plot the datset
hold on
temp.traces(:,~bool)=NaN;        % Set all traces except 5 and 7 to NaN
temp.null=NaN;                   % Dataset has null values (this is only done to avoid warnings)
s_wplot(temp,{'deflection',1},{'wiggle_color','red'},{'peak_fill','red'}, ...
             {'figure','old'})   % Plot this dataset on top of the one 
                                                      % plotted before
mytitle('Example 7a: Traces plotted in black and red')

%% 7b Plot the 5-th and the 8-th trace in different colors (handle manipulation)
aux=s_wplot(seismic,{'trough_fill','white'},{'deflection',1});   % Plot the datset.
                                  % Notice that the trough-fill color is specified
set(aux.handles([5,8],1),'FaceColor','green')          % Peak fill (green)
set(aux.handles([5,8],2),'FaceColor',[0.6,0.6,0.6])    % Trough fill (gray)
set(aux.handles([5,8],3),'Color','red','LineWidth',1)  % Wiggle width and color (red)
mytitle('Example 7b: Traces plotted in black and other colors')

%% 9 Wiggle trace over color
hold on
s_wplot(seismic,{'figure','old'},{'peak_fill',''},{'wiggle_width',1.5}, ...
mytitle('Example 9: Wiggle trace over color')

%% 10 Custom wiggle plot with color-coded troughs
s_xplot(seismic,{'title','Example 10: Custom wiggle plot with color-coded troughs'})

%% 11 Use plot scale factor of one dataset for another dataset
seismic2=0.5*seismic1; % Create a dataset with twice the amplitudes
   aux=s_wplot(seismic1,{'figure','old'},{'title', 'Dataset'});
   s_wplot(seismic2,{'figure','old'},{'scale',aux.scale},{'title', '0.5*Dataset (same scale)'})
%   mytitle(' ')
   mysuptitle('Example 11: Plot of two equally scaled datasets',{'color','blue'})

Contact us