firstbreak

su 的源代码库

#! /bin/sh
# First break picking
# Author: Jack
# NOTE: Comment lines preceeding user input start with  #!#
#set -x

#!# Set parameters
input=cdpby100
cdp=601		# pick a cdp -- also used for naming output files
dxcdp=13.36
timefrac=2	# use first  nt/timefrac samples
fold=30 panelsize=6

### Label output according to cdp number and get the cdp
pickfile=fbpicks.$cdp
nt=`sugethw ns <$input | sed 1q | sed 's/.*ns=//'`
itmax=`bc -l <<-END
	$nt / $timefrac
END`
suwind <$input key=cdp min=$cdp max=$cdp count=$fold |
suwind itmax=$itmax >cdp.$cdp

### Do the picking
echo "Pick first breaks by moving mouse and typing 's', type 'Q' when done"
trace=1
while [ $trace -le $fold ]
do
	sushw <cdp.$cdp key=tracr a=1 b=1 |
	suwind key=tracr min=$trace count=$panelsize |
	suxwigb perc=99 xcur=0.25 nbpi=96 f2=$trace d2=1 wbox=1100 \
		label1="Time (sec)" label2="Velocity (m/sec)" \
		mpicks=mpicks.$trace

	trace=`bc -l <<-END
		$trace + $panelsize
	END`
done

set +x

### Combine the individual picks
echo "Editing pick files ..."
trace=1
while [ $trace -le $fold ]
do
	cat mpicks.$trace
	trace=`bc -l <<-END
		$trace + $panelsize
	END`
done |
awk '
	BEGIN	{dx='$dxcdp'}
		{print $1 "	" (NR-1)*dx}
' >$pickfile

a2b <$pickfile >temp outpar=picks.p

xgraph <temp par=picks.p nplot=1 \
	label1="Time (sec)" label2="Range (m)" \
	title="First Breaks" \
	grid1=solid grid2=solid \
	linecolor=2 style=seismic &

### Clean up
trace=1
while [ $trace -le $fold ]
do
	rm mpicks.$trace
	trace=`bc -l <<-END
		$trace + $panelsize
	END`
done
rm temp

echo "The first breaks are recorded in $pickfile"