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.

% Seismic_examples2
%       Examples of the usage of seismic functions 

keep WF
global S4M   %#ok

%     Create a minimum-phase wavelet

%     Create constant-Q filters for Q = 150,100,80,60 and 1 sec travel time

%     Convolve the constant-Q filters with the wavelet

%     Prepend the original wavelet (since the original wavelet is shorter 
%     than the filtered ones the missing samples are, by default, set to 
%     NaN's; this will give a warning in R 14, if the data set is plotted}

%     Give the data set "wavelets" a descriptive name'Original wavelet and its Q-filtered versions';

%     Plot the result with the following conditions:
%         plot only the first 200 ms;
%         annotate traces 2 to 5 with the Q-value of the filter (the Q-value
%             of the first trace is NaN as it is the original wavelet); 
%         fill the wiggle troughs with light gray
s_wplot(wavelets,{'times',0,200},{'aindex',2:5},{'annotation','Q'}, ...

%     Display real and imaginary parts of the Fourier transform of the wavelets
%     for frequencies from 0 to 80 Hz
mytitle('Real (black) and imaginary (red) parts of the Fourier transform of the wavelets')

Contact us