bshp - BSU program which shapes one seismic data set to another data set (Levinson)
bshp [ -h | infile1 | infile2 | iswopt | iswch | tmin | tmax | npf | sf | infil3 ]
Basic Seismic Utilities (BSU) to shape signals. Determines least squares shaping filter using the Levinson algoritm to solve the normal equations. Subroutines which do the shaping were written by E. A. Robinson (1967). The output is the application of the filters to a selected data set. In one option, one designs filters to apply to file1 (infile1*filters=infile2). An alternative option is to apply the filters to a third data set (design of filters is same, but application is to a file3. If infile3 contains impulses, one can capture the filters in a bsegy data set. The sample interval, number of samples, and number of traces must be identical for files 1, 2, and 3. A file of impulse functions can be generated with program bdum.
Options
|
-h |
Online help giving details on command line arguments |
infile1
Input file name (file to shape to match infile2)
infile2
Input file name (file which is
target for infile1 to match)
(design: infile1*filter=infile2)
|
iswopt |
Switch to select type of computation |
0=profile mode (one average filter for entire data set)
1=trace mode (each trace has its own filter)
|
iswch |
Switch to direct filter application |
0=apply filters to infile3 (design is infile1-->infile2)
1=apply filters to infile1 (design: infile1*filter=infile2)
|
tmin |
Start of design gate on infile1 (seconds) | ||
|
tmax |
End of design gate on infile1 (seconds) | ||
|
npf |
Number of points in filter (samples) | ||
|
sf |
Stability factor 0<sf<.01 typical (fraction of zero lag autocorrelation) |
infile3
Input file name (only required if iswch=0)
NOTE:
If invoked with no options, will prompt user for input
parameters.
EXAMPLE:
bshp f001.seg f002.seg 0 1 0. 0.1 51 .01
Shape f001.seg to match waveforms in f002.seg using profile mode (one filter for all the traces in f001.seg), apply the filter to f001.seg. The autocorrelation will be computed on data from 0 to 0.1 seconds in f001.seg, and only 51 samples of the autocorrelation will be computed (also the number of samples in the filter). A stability factor of 0.01xC(0) will be added to C(0), the zero lag autocorrelation (bias the spectrum).
bshpxxxx.seg
Named according to convention (first 4char bshp, the next 4char are the first 4char of the input file name, suffix .seg)
standard output
produces a progress bar
bshpxxxx.lst
Echo check of input parameters. The values "ase" listed are the average least squares error of the resulting match between the actual and desired output. The smaller the "ase" values, the better the result.
bhelp(1), bdum(1), shape(3), fold(3), eureka(3), dot(3), cross(3), tlist(3)
No known bugs
Copyright © 2017 by Paul Michaels
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
P. Michaels, PE. <pm@cgiss.boisestate.edu>