BSU comes with some example bash scripts. The scripts xplot and psplot are examples of how we can combine trace equalization with plotting, and default a number of bplt parameters. As can be seen from the listing below, the script prompts the user for some values, then scales the data with either bscl or bequ before plotting.
#!/bin/sh # $Id: psplot,v 1.2 2009-05-05 19:43:43 pm Exp $ # # Script to plot seismic data in Postscript format using BSU program bplt $ #Author: P. Michaels <pm@cgiss.boisestate.edu> #set -x XDIM=5 YDIM=2 if test "$1" = "-h" then echo "USAGE: xplot filename tmax scaling" echo 'Scaling Choices:' echo '1= Peak Absolute Value of profile' echo '2= L2 Norm of profile' echo '3= Trace by trace L2 Norm' else if test "$1" = '' then echo 'Enter input file name' read FILEN else FILEN=$1 fi if test "$2" = '' then echo 'Enter tmax' read TMAX else TMAX=$2 fi if test "$3" = '' then echo 'Enter Scaling Choice' echo '1= Peak Absolute Value of profile' echo '2= L2 Norm of profile' echo '3= Trace by trace L2 Norm' read SCL else SCL=$3 fi NAME=`basename $FILEN .seg` NAME4=`echo $NAME | gawk -F "" '{print $1$2$3$4}' ` case $SCL in 1) bscl $FILEN 1 500 3 AMP=`gawk '/Peak Absolute Value/ {print $4}' bscl$NAME4.lst` rm -f bscl$NAME4.* PFILEN=$FILEN bplt $PFILEN 1 1 0 1 500 0 $TMAX 1 $AMP 200 $XDIM $YDIM ;; 2) bscl $FILEN 1 500 1 AMP=`gawk '/L2 Norm of Data Set=/ {print $6}' bscl$NAME4.lst` rm -f bscl$NAME4.* PFILEN=$FILEN bplt $PFILEN 1 1 0 1 500 0 $TMAX 1 $AMP 200 $XDIM $YDIM ;; 3) bequ $FILEN 0 $TMAX PFILEN=bequ$NAME4.seg AMP=4 bplt $PFILEN 1 1 0 1 500 0 $TMAX 1 $AMP 200 $XDIM $YDIM rm -f bequ$NAME4.* ;; esac rm -f bplt*.lst fi
There are two ways to run this script. One can simply type psplot and then be prompted for the data set name, the maximum time to plot, and the type of scaling to use. Alternatively, one can type all the arguments on the command line:
psplot twav.seg .2 3
which produces a plot as shown in Figure 3.
|
Another useful script is xplot. It plots the data using the X11 x-window device, and is very similar to the psplot script.
pm 2018-04-08