makefile2

国外免费地震资料处理软件包

	label1=time label2=offset crowd1=0.85 \
	title="error: (shot1+shot3)/2 - shot2" out=diff2.v
	vppen vpstyle=n big=n gridnum=2,2 < shot.v diff.v diff1.v diff2.v > $@

${RESDIR}/marmshotin.v: marmdiff2.H marmsimp.H 
	< marmdiff2.H Byte clip=2706 > diff.A
	< diff.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="interpolated shot" out=shot.v
	< diff.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error: interpolated shot - shot2" out=diff.v
	< marmsimp.H Byte clip=2706 > diff.A
	< diff.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="(shot1+shot3)/2" out=diff1.v
	< diff.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error: (shot1+shot3)/2 - shot2" out=diff2.v
	vppen vpstyle=n big=n gridnum=2,2 \
	 < shot.v diff.v diff1.v diff2.v > $@

fslice%.H: marmfftd.H ${BINDIR}/Average.x
	< marmfftd.H Window n3=1 f3=$* | ${BINDIR}/Average.x > $@

mfill%.H: marmfill.H
	< marmfill.H Window n3=1 f3=$* > $@

${RESDIR}/mslice.v ${RESDIR}/mslice.v3: \
	fslice2005.H fslice1962.H fslice1919.H \
	mfill2005.H mfill1962.H mfill1919.H
	< fslice2005.H Real | ${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=f1.v
	< fslice1962.H Real | ${OGREY} pclip=99 title="frequency=20" \
	transp=n label1=midpoint label2=offset out=f2.v
	< fslice1919.H Real | ${OGREY} pclip=99 title="frequency=30" \
	transp=n label1=midpoint label2=offset out=f3.v
	vp_OverUnderAniso f1.v f2.v f3.v > f.v
	< mfill2005.H Real | ${OGREY} pclip=99 title="frequency=10" \
	transp=n label1=midpoint label2=offset out=i1.v
	< mfill1962.H Real | ${OGREY} pclip=99 title="frequency=20" \
	transp=n label1=midpoint label2=offset out=i2.v
	< mfill1919.H Real | ${OGREY} pclip=99 title="frequency=30" \
	transp=n label1=midpoint label2=offset out=i3.v
	vp_OverUnderAniso i1.v i2.v i3.v > i.v
	vp_SideBySideAniso f.v i.v > junk.v
	vp_Movie f.v i.v > ${RESDIR}/mslice.v3
	< junk.v vppen vpstyle=n txscale=3 xsize=7.5 ysize=10 \
	> ${RESDIR}/mslice.v

marmfill.H: marmfftd.H marmmask.H ${BINDIR}/Ofilp.x 
	< marmfftd.H ${BINDIR}/Ofilp.x known=marmmask.H niter=200 >$@

marmfinv.H: marmfill.H ${BINDIR}/FFT2.x ${BINDIR}/Shift.x
	< marmfill.H ${TRANSP} plane=31 > junkm2.H
	< junkm2.H ${BINDIR}/Shift.x inv=y > junkm0.H
	Rm junkm2.H
	< junkm0.H ${BINDIR}/FFT2.x > $@
	Rm junkm0.H

beistkd.H: beistck.H beinmod.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< beistck.H ${TRANSP} > junkb0.H
	< junkb0.H ${BINDIR}/FFT2.x > junkb1.H
	< junkb1.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< beinmod.H Window n1=15 n3=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

elfstkd.H: elfstck.H elfnmod.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< elfstck.H ${TRANSP} > junke0.H
	< junke0.H ${BINDIR}/FFT2.x > junke1.H
	< junke1.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< elfnmod.H Window n1=20 n2=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

beistkm.H: beistcm.H beinmod.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< beistcm.H Window n2=1 f2=1 | ${TRANSP} > junkb0.H
	< junkb0.H ${BINDIR}/FFT2.x > junkb1.H
	< junkb1.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< beinmod.H Window n1=15 n3=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

synstkd.H: synstck.H synnmod.H ${BINDIR}/FFT2.x ${BINDIR}/Stretch.x
	< synstck.H ${TRANSP} > junks0.H
	< junks0.H ${BINDIR}/FFT2.x > junks1.H
	< junks1.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< synnmod.H Window n1=50 n3=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

beiagcd.H: beistkd.H
	< beistkd.H Smooth rect1=250 agc=1 > $@

beisgcd.H: beistup.H
	< beistup.H Smooth rect1=250 agc=1 > $@

${RESDIR}/beistack.v ${RESDIR}/beistack.v3: beistkd.H beistup.H
	< beistkd.H Tpow tpow=-0.5 | ${OGREY} title="OC Stack" out=ostack.v
	< beistup.H Tpow tpow=-0.5 | ${OGREY} title="NMO Stack" out=sstack.v
	vp_Movie sstack.v ostack.v > ${RESDIR}/beistack.v3
	vp_OverUnderAniso sstack.v ostack.v > ${RESDIR}/beistack.v

marmfild.H: marmfinv.H marmnmod.H \
	${BINDIR}/Stretch.x ${BINDIR}/NMO.x marmvpk.H 
	< marmfinv.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< marmnmod.H Window n1=10 > junk2.H
	Cat axis=1 junk2.H junk1.H > junk3.H
	< junk3.H ${BINDIR}/NMO.x vel=marmvpk.H inv=1 > $@

marmavg.H: marmnmod.H ${BINDIR}/Average2.x
	< marmnmod.H ${BINDIR}/Average2.x > $@

marmbak.H: marmavg.H ${BINDIR}/NMO.x marmvpk.H \
	${BINDIR}/Average2.x ${BINDIR}/Mutter.x
	< marmvpk.H ${BINDIR}/Average2.x > marmvbk.H
	< marmavg.H ${BINDIR}/NMO.x vel=marmvbk.H inv=y | \
	${BINDIR}/Mutter.x slope0=0.0014 > $@


#################### B-spline ##############################################

splan2_%.H: ${BINDIR}/Splineplane.x
	${BINDIR}/Splineplane.x p=$* eps=0.0001 nw=2 > $@

splan3_%.H: ${BINDIR}/Splineplane.x
	${BINDIR}/Splineplane.x p=$* eps=0.0001 nw=3 > $@

splan4_%.H: ${BINDIR}/Splineplane.x
	${BINDIR}/Splineplane.x p=$* eps=0.0001 nw=4 > $@

sthree%.H: splan%_0.7.H splan%_-1.H splan%_-0.1.H ${BINDIR}/Splinefilter.x \
	spike1.H spike2.H spike3.H
	< spike1.H Helicon div=1 filt=splan$*_0.7.H  > pl1.H
	< spike2.H Helicon div=1 filt=splan$*_-1.H   > pl2.H
	< spike3.H Helicon div=1 filt=splan$*_-0.1.H > pl3.H
	Add pl1.H pl2.H pl3.H | \
	${BINDIR}/Splinefilter.x post=y nw=-$* both=y | Window n2=75 > $@

${RESDIR}/sthree.v ${RESDIR}/sthree.v3: sthree2.H sthree4.H
	< sthree2.H ${SGREY} out=three2.v \
		title="Spline waves B-1 (linear)" titlesz=9
	< sthree4.H ${SGREY} out=three4.v \
		title="Spline waves B-3 (cubic)" titlesz=9
	vp_Movie three2.v three4.v > ${RESDIR}/sthree.v3
	vp_SideBySideAniso three2.v three4.v > ${RESDIR}/sthree.v

##################### 3-D real #############################################
vnmo.H: ${DATADIR}/vnmo.H
	< ${DATADIR}/vnmo.H ${TRANSP} | ${TRANSP} plane=23 > $@

head-win.H: ${DATADIR}/jupiter_win.H@@
	< ${DATADIR}/jupiter_win.H@@ Window n1=4 f1=11 > $@

offset-win.H: ${DATADIR}/jupiter_win.H@@
	< ${DATADIR}/jupiter_win.H@@ Window n1=1 f1=15 > $@

vel-win.H: vnmo.H ${BINDIR}/Velint.x head-win-transp.H
	< vnmo.H ${BINDIR}/Velint.x cmp=head-win-transp.H nw=4 > $@

nmo-win.H: ${DATADIR}/jupiter_win.H vel-win-transp.H \
	${BINDIR}/NMO3.x offset-win.H nmo.p
	< ${DATADIR}/jupiter_win.H ${BINDIR}/NMO3.x \
	offset=offset-win.H vel=vel-win-transp.H par=nmo.p > nmo-win.H
	echo hff=-1 >> $@
	echo gff=-1 >> $@

dat-win.H: nmo-win.H
	< nmo-win.H Window n1=250 f1=100 > $@

%-win-transp.H: %-win.H
	< $*-win.H ${TRANSP} > $@

smo2-win.H: dat-win-transp.H head-win-transp.H \
	grid.p ${BINDIR}/Invinter2.x gtens0.H
	< dat-win-transp.H ${BINDIR}/Invinter2.x coord=head-win-transp.H \
	filt=gtens0.H par=grid.p spline=n niter=50 > $@

%-dip.H: %.H ${BINDIR}/Smooth3.x
	< $*.H ${BINDIR}/Smooth3.x mask=$*.H nliter=10 > $@

%-pch.H: %-dip.H ${BINDIR}/Patch.x 
	< $*-dip.H ${BINDIR}/Patch.x pmax=4 np=400 > $@

bank%.H nh%.H lag%.H: ${BINDIR}/Splinebank.x
	${BINDIR}/Splinebank.x nw=$* pmax=4 np=400 nh=nh$*.H lag=lag$*.H \
	nt=100 niter=20 eps=0.0001 > bank$*.H

int2-win.H: dat-win-transp.H head-win-transp.H smo2-win-pch.H grid.p \
	${BINDIR}/Ninvinter2.x bank2.H nh2.H lag2.H 
	< dat-win-transp.H ${BINDIR}/Ninvinter2.x coord=head-win-transp.H \
	filt=bank2.H nh=nh2.H lag=lag2.H pch=smo2-win-pch.H \
	par=grid.p spline=n Xeps=1 niter=50 > $@


TITLE.bin = Binning
TITLE.int4 = B-spline+plane
TITLE.int2 = Linear+plane

${RESDIR}/int2-win.v: int2-win.H
	< int2-win.H ${TRANSP} plane=23 > junk.H
	< junk.H ${TRANSP} > transp.H
	< transp.H Byte gainpanel=all > int2.A
	< int2.A ${CUBE1} title=${TITLE.int2} out=$@

int4-win.H: dat-win-transp.H head-win-transp.H smo2-win-pch.H grid.p \
	${BINDIR}/Ninvinter2.x bank4.H nh4.H lag4.H ${BINDIR}/Splinefilter.x
	< dat-win-transp.H ${BINDIR}/Ninvinter2.x coord=head-win-transp.H \
	filt=bank4.H nh=nh4.H lag=lag4.H pch=smo2-win-pch.H \
	par=grid.p spline=y nw=4 eps=1 niter=50 | \
	${BINDIR}/Splinefilter.x nw=4 post=y both=y > $@

bin-win.H fold-win.H: grid.p dat-win-transp.H head-win-transp.H \
	${BINDIR}/Bin2.x
	${BINDIR}/Bin2.x par=grid.p hff=head-win-transp.H fold=fold-win.H \
	interp=1 <dat-win-transp.H > bin-win.H

${RESDIR}/fold-win.v: fold-win.H
	<fold-win.H Grey allpos=y pclip=100 \
	transp=n yreverse=n wantscalebar=1 \
	label1="In-line midpoint" label2="Cross-line midpoint" \
	wheretitle=t wherexlabel=b \
	title="Fold Map" out=$@ > /dev/null

${RESDIR}/%-win.v: %-win.H
	< $*-win.H ${TRANSP} plane=23 > junk.H
	< junk.H ${TRANSP} > transp.H
	< transp.H Byte gainpanel=all > $*.A
	< $*.A ${CUBE1} title=${TITLE.$*} out=$@

${RESDIR}/winslice.v: int4-win.H bin-win.H
	< bin-win.H Window n3=1 f3=132 | \
	Grey crowd1=0.8 title='t=0.928' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice1.v
	< bin-win.H Window n3=1 f3=164 | \
	Grey crowd1=0.8 title='t=1.056' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice3.v
	< bin-win.H Window n3=1 f3=196 | \
	Grey crowd1=0.8 title='t=1.184' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice2.v
	< bin-win.H Window n3=1 f3=228 | \
	Grey crowd1=0.8 title='t=1.312' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice4.v
	vp_OverUnderAniso slice1.v slice2.v slice3.v slice4.v > bin.v
	< int4-win.H Window n3=1 f3=132 | \
	Grey crowd1=0.8 title='t=0.928' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice1.v
	< int4-win.H Window n3=1 f3=164 | \
	Grey crowd1=0.8 title='t=1.056' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice3.v
	< int4-win.H Window n3=1 f3=196 | \
	Grey crowd1=0.8 title='t=1.184' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice2.v
	< int4-win.H Window n3=1 f3=228 | \
	Grey crowd1=0.8 title='t=1.312' >/dev/null \
	labelsz=18 titlesz=24 transp=n out=slice4.v
	vp_OverUnderAniso slice1.v slice2.v slice3.v slice4.v > int4.v
	vp_SideBySideAniso bin.v int4.v > slice.v
	< slice.v vppen vpstyle=n xsize=12 ysize=8 > $@

##################### 3-D vs splitting ############################
lintens1.H: ${BINDIR}/Spline0.x
	${BINDIR}/Spline0.x eps=0.001 tension=1 > $@

tens1.v: lintens1.H 
	Spike n1=60 n2=60 nsp=1 k1=30 k2=30 > inp.H
	<inp.H Helicon       div=1 filt=lintens1.H > div.H
	<div.H Helicon adj=1 div=1 filt=lintens1.H > div2.H
	<div2.H ${GREY} out=$@ pclip=99 title=factorization allpos=y

bob.H: bob.HH
	< bob.HH Dd esize=4 | Window > $@

bob1.v: bob.H 
	Spike n1=60 n2=60 nsp=1 k1=30 k2=30 > inp.H
	<inp.H Helicon       div=1 filt=bob.H > div.H
	<div.H Helicon adj=1 div=1 filt=bob.H > div2.H
	<div2.H ${GREY} out=$@ pclip=99 title=splitting allpos=y

${RESDIR}/bob.v: tens1.v bob1.v
	vp_SideBySideAniso tens1.v bob1.v > ${RESDIR}/bob.v

############ Interpolation beyond aliasing #############################

aliasp.H: ${BINDIR}/Aliasp.x
	${BINDIR}/Aliasp.x cycles=25 > $@

%-pad2.H %-mask4.H: %.H ${BINDIR}/LPad.x
	<$*.H ${BINDIR}/LPad.x jump=4 mask=$*-mask4.H > $*-pad2.H

aliasp-pq2.H: ${BINDIR}/Smooth2.x aliasp.H
	< aliasp.H ${BINDIR}/Smooth2.x nw=3 nj=4 eps=0.04 > $@

aliasp-p2.H: aliasp.H aliasp-pq2.H ${BINDIR}/Smooth.x 
	< aliasp.H ${BINDIR}/Smooth.x nw=3 nj=4 pfile=aliasp-pq2.H | \
	Transp > $@

aliasp-q2.H: aliasp-p2.H ${BINDIR}/Double.x 
	< aliasp-p2.H ${BINDIR}/Double.x nw=4 | \
		      ${BINDIR}/Double.x nw=4 | Transp > $@

%-mis4.H: %-pad2.H %-mask4.H %-q2.H ${BINDIR}/Miss1.x 
	< $*-pad2.H ${BINDIR}/Miss1.x nw=3 \
	dip=$*-q2.H mask=$*-mask4.H niter=100 > $@

${RESDIR}/aliasp2.v: aliasp.H aliasp-mis4.H
	<aliasp.H Window n1=575 | Grey title=Input \
	crowd=0.9 >/dev/null wantaxis=n out=inp.v
	<aliasp-mis4.H Window n1=575 | Grey title=Output \
	crowd=0.9 >/dev/null wantaxis=n out=out.v
	vp_SideBySideAniso inp.v out.v > $@

################ offcon impulse response ###############################
taper.H: spike.p ${BINDIR}/Dipfilter.x taper.p
	Spike  par=spike.p | \
	${BINDIR}/Dipfilter.x par=taper.p | \
	Bandpass flo=10 fhi=50 > $@

${RESDIR}/off-imp.v ${RESDIR}/off-imp.v3: taper.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT.x ${BINDIR}/Fincon.x \
	${BINDIR}/Prod.x ${BINDIR}/Oper.x
	<taper.H ${BINDIR}/Stretch.x dens=2 | Rtoc | \
	${BINDIR}/FFT.x sign1=1 sign2=1 > fft.H
	<fft.H ${BINDIR}/Oper.x h=0.5 | ${BINDIR}/Prod.x prod=fft.H > out.H
	< out.H ${BINDIR}/FFT.x sign1=-1 sign2=-1 | Real | \
	${BINDIR}/Stretch.x inv=1 dens=2 >imp.H
	<imp.H ${OGREY} title="Omega-K IDMO impulse response" out=imp1.v
	<taper.H ${BINDIR}/Stretch.x dens=2 | Rtoc | \
	${BINDIR}/FFT.x sign1=1 > fft.H
	<fft.H ${TRANSP} | \
	${BINDIR}/Fincon.x h0=0 dh=0.004 nh=125 | \
	${TRANSP} >out.H
	<out.H ${BINDIR}/FFT.x sign1=-1 | Real | \
	${BINDIR}/Stretch.x inv=1 dens=2 >imp.H
	<imp.H ${OGREY} title="FinDif OC impulse response" out=imp2.v
	vp_SideBySideAniso imp1.v imp2.v > ${RESDIR}/off-imp.v
	vp_Movie imp1.v imp2.v > ${RESDIR}/off-imp.v3

clean : jclean
	${RM} ${BINDIR}/*.x
	${RM} ${OBJDIR}/*.o

include ${SEPINC}/SEP.bottom