makefile2

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

	< bei23.H Transp plane=23 > $@

elfnmod.H: elf.H ${DATADIR}/elfvel.H ${BINDIR}/NMO.x 
	< elf.H ${BINDIR}/NMO.x vel=${DATADIR}/elfvel.H > $@

elfhole.H elfmask.H: elfnmod.H ${BINDIR}/Shotholes.x 
	< elfnmod.H ${BINDIR}/Shotholes.x perc=0.6 > elfhole.H mask=junke.H
	< junke.H ${TRANSP} > elfmask.H

elfshot.H: elfnmod.H ${BINDIR}/Cmp2shot.x
	< elfnmod.H ${BINDIR}/Cmp2shot.x > $@

elfshot3.H: elfshot.H
	< elfshot.H Window n3=3 f3=800 > $@

elfsimp.H: elfshot3.H
	< elfshot3.H Window n3=1 f3=0 > junk1.H
	< elfshot3.H Window n3=1 f3=1 > junk2.H
	< elfshot3.H Window n3=1 f3=2 > junk3.H
	Add junk1.H junk3.H | Scale dscale=0.5 > junk4.H
	Add scale=1,-1 junk4.H junk2.H > junk5.H
	Cat junk4.H junk5.H > $@

marmsimp.H: marmshot3.H
	< marmshot3.H Window n3=1 f3=0 > junk1.H
	< marmshot3.H Window n3=1 f3=1 > junk2.H
	< marmshot3.H Window n3=1 f3=2 > junk3.H
	Add junk1.H junk3.H | Scale dscale=0.5 > junk4.H
	Add scale=1,-1 junk4.H junk2.H > junk5.H
	Cat junk4.H junk5.H > $@

simp.H: shot.cup.H
	< shot.cup.H Window n3=1 f3=78 > junk1.H
	< shot.cup.H Window n3=1 f3=79 > junk2.H
	< shot.cup.H Window n3=1 f3=80 > junk3.H
	Add junk1.H junk3.H | Scale dscale=0.5 > junk4.H
	Add scale=1,-1 junk4.H junk2.H > junk5.H
	Cat junk4.H junk5.H > $@

synnmod.H: ${DATADIR}/synth.H ${DATADIR}/synvel.H \
	${BINDIR}/NMO.x ${BINDIR}/Mutter.x
	< ${DATADIR}/synth.H Transp plane=23 > syn23.H
	echo "d2=0.032" >> syn23.H
	< syn23.H ${BINDIR}/Mutter.x slope0=0.7 mute=no > synmut.H
	< synmut.H ${BINDIR}/NMO.x vel=${DATADIR}/synvel.H > syn23.H
	< syn23.H Transp plane=23 > $@

marmshots.H: marmnmod.H ${BINDIR}/Cmp2shot.x
	< marmnmod.H ${BINDIR}/Cmp2shot.x > junk.H
	< junk.H Window f3=103 j3=2 n3=240 f2=8 > $@

marmffts.H: marmshots.H ${BINDIR}/Stretch.x ${BINDIR}/FFT1.x
	< marmshots.H Window f1=10 | ${BINDIR}/Stretch.x nout=4096 > junkms.H
	< junkms.H ${BINDIR}/FFT1.x > ffts.H
	< ffts.H ${TRANSP} > junkms.H
	< junkms.H ${TRANSP} plane=23 > $@
	Rm junkms.H

marmfills.H: marmffts.H ${BINDIR}/Shotfill.x
	< marmffts.H ${BINDIR}/Shotfill.x > $@

marmfilds.H: marmfills.H ${BINDIR}/Stretch.x ${BINDIR}/FFT2.x marmshots.H
	< marmfills.H ${TRANSP} plane=23 > junkms.H
	< junkms.H ${TRANSP} > ffts.H
	< ffts.H ${BINDIR}/FFT2.x > junkms.H
	< junkms.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< marmshots.H Window n1=10 > junk2.H
	Cat junk2.H junk2.H > junk3.H
	Cat axis=1 junk3.H junk1.H > $@

marmshot3.H: marmshots.H ${BINDIR}/Mutter.x
	< marmshots.H Window n3=3 f3=63 | \
	${BINDIR}/Mutter.x slope0=0.0012 > $@

${RESDIR}/marmshot3.v ${RESDIR}/marmshot3.v3: marmshot3.H
	< marmshot3.H Window n3=1 f3=0 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 1" out=shot1.v
	< marmshot3.H Window n3=1 f3=1 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 2" out=shot2.v
	< marmshot3.H Window n3=1 f3=2 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 3" out=shot3.v
	vp_Movie shot1.v  shot2.v shot3.v > ${RESDIR}/marmshot3.v3
	vp_OverUnderAniso shot1.v shot3.v > ${RESDIR}/marmshot3.v

marmstrd.H: marmnmod.H ${BINDIR}/Stretch.x
	< marmnmod.H Window f1=10 | ${BINDIR}/Stretch.x nout=4096 > $@

beistrd.H: beinmod.H ${BINDIR}/Stretch.x
	< beinmod.H Window f1=15 | \
	${BINDIR}/Stretch.x nout=2048 > $@

synstrd.H: synnmod.H ${BINDIR}/Stretch.x
	< synnmod.H Window f1=50 | ${BINDIR}/Stretch.x nout=512 > $@

elfstrd.H: elfhole.H ${BINDIR}/Stretch.x
	< elfhole.H Window f1=20 | ${TRANSP} plane=23 | \
	${BINDIR}/Stretch.x nout=2048 > $@

marmstrd3.H: marmshot3.H ${BINDIR}/Stretch.x
	< marmshot3.H Window f1=10 | ${BINDIR}/Stretch.x nout=4096 > $@

synstrd3.H: synshot3.H ${BINDIR}/Stretch.x
	< synshot3.H Window f1=50 | ${BINDIR}/Stretch.x nout=512 > $@

elfstrd3.H: elfshot3.H ${BINDIR}/Stretch.x
	< elfshot3.H Window f1=20 | ${BINDIR}/Stretch.x nout=2048 > $@

marmfftd.H: ${BINDIR}/FFT1.x marmstrd.H
	< marmstrd.H ${BINDIR}/FFT1.x > junkm.H
	< junkm.H ${TRANSP} plane=31 > $@
	Rm junkm.H

beifftd.H: ${BINDIR}/FFT1.x beistrd.H
	< beistrd.H ${BINDIR}/FFT1.x > junkb.H
	< junkb.H ${TRANSP} > junkb1.H
	Rm junkb.H
	< junkb1.H ${TRANSP} plane=23 > $@
	Rm junkb1.H

elffftd.H: ${BINDIR}/FFT1.x elfstrd.H
	< elfstrd.H Window n2=500 | ${BINDIR}/FFT1.x > junke1.H
	< elfstrd.H Window f2=500 | ${BINDIR}/FFT1.x > junke2.H
	< junke1.H ${TRANSP} > junke3.H
	< junke3.H ${TRANSP} plane=23 > junke1.H
	< junke2.H ${TRANSP} > junke3.H
	< junke3.H ${TRANSP} plane=23 > junke2.H
	Rm junke3.H
	Cat axis=1 junke1.H junke2.H > elffftd.H

elffill.H: elffftd.H elfmask.H ${BINDIR}/Ofilp.x
	< elffftd.H ${BINDIR}/Ofilp.x known=elfmask.H niter=200 > $@

elfhill.H: elffftd.H elfmask.H ${BINDIR}/Hfilp.x
	< elffftd.H ${BINDIR}/Hfilp.x known=elfmask.H niter=200 > $@

elffinv.H: ${BINDIR}/FFT2.x elffill.H
	< elffill.H Window n1=500 > junke3.H
	< junke3.H ${TRANSP} plane=23 > junke1.H
	< junke1.H ${TRANSP} > junke3.H
	< junke3.H ${BINDIR}/FFT2.x > junke1.H
	< elffill.H Window f1=500 > junke3.H
	< junke3.H ${TRANSP} plane=23 > junke2.H
	< junke2.H ${TRANSP} > junke3.H
	< junke3.H ${BINDIR}/FFT2.x > junke2.H
	Rm junke3.H
	Cat axis=2 junke1.H junke2.H > $@

elfhinv.H: ${BINDIR}/FFT2.x elfhill.H
	< elfhill.H Window n1=500 > junkg3.H
	< junkg3.H ${TRANSP} plane=23 > junkg1.H
	< junkg1.H ${TRANSP} > junkg3.H
	< junkg3.H ${BINDIR}/FFT2.x > junkg1.H
	< elfhill.H Window f1=500 > junkg3.H
	< junkg3.H ${TRANSP} plane=23 > junkg2.H
	< junkg2.H ${TRANSP} > junkg3.H
	< junkg3.H ${BINDIR}/FFT2.x > junkg2.H
	Rm junkg3.H
	Cat axis=2 junkg1.H junkg2.H > $@

elffild.H: elffinv.H ${BINDIR}/Stretch.x elfhole.H
	< elffinv.H ${BINDIR}/Stretch.x inv=1 | \
	${TRANSP} plane=23 > junke2.H
	< elfhole.H Window n1=20 > junke1.H 
	Cat axis=1 junke1.H junke2.H > $@

elfhild.H: elfhinv.H ${BINDIR}/Stretch.x elfhole.H
	< elfhinv.H ${BINDIR}/Stretch.x inv=1 | \
	${TRANSP} plane=23 > junkg2.H
	< elfhole.H Window n1=20 > junkg1.H 
	Cat axis=1 junkg1.H junkg2.H > $@

beistrm.H: beistrd.H
	< beistrd.H Reverse which=4 | Pad n3=26 | Reverse which=4 > $@

beimask.H: beistrm.H ${BINDIR}/Holes.x
	<beistrm.H ${BINDIR}/Holes.x perc=0 >junkb.H mask=$@

beifftm.H: ${BINDIR}/FFT1.x beistrm.H
	< beistrm.H ${BINDIR}/FFT1.x > junkb.H
	< junkb.H ${TRANSP} > junkb1.H
	Rm junkb.H
	< junkb1.H ${TRANSP} plane=23 > $@
	Rm junkb1.H

synfftd.H: ${BINDIR}/FFT1.x ${BINDIR}/Shift.x synstrd.H
	< synstrd.H ${BINDIR}/FFT1.x > junks.H
	< junks.H ${TRANSP} > junks1.H
	Rm junks.H
	< junks1.H ${TRANSP} plane=23 > $@
	Rm junks1.H

beistck.H: beifftd.H ${BINDIR}/Finstack.x
	< beifftd.H ${BINDIR}/Finstack.x > $@

elfstck.H: elffftd.H ${BINDIR}/Finstack.x
	< elffftd.H ${BINDIR}/Finstack.x > $@

beistcm.H: beifftm.H ${BINDIR}/Ofilm.x beimask.H
	< beifftm.H ${BINDIR}/Ofilm.x \
	niter=100 known=beimask.H ns=24 zmin=3 > $@

synstck.H: synfftd.H ${BINDIR}/Finstack.x
	< synfftd.H ${BINDIR}/Finstack.x > $@

beistup.H: beinmod.H
	< beinmod.H ${TRANSP} plane=23 | Stack norm=no > $@

synstup.H: synnmod.H
	< synnmod.H ${TRANSP} plane=23 | Stack norm=no > $@

elfstup.H: elfnmod.H
	< elfnmod.H Stack norm=no > $@

marmfftd3.H: ${BINDIR}/FFT1.x ${BINDIR}/Shift.x marmstrd3.H
	< marmstrd3.H ${BINDIR}/FFT1.x > junkm0.H
	< junkm0.H ${BINDIR}/Shift.x > junkm1.H
	< junkm1.H ${TRANSP} | ${TRANSP} plane=23 > $@

synfftd3.H: ${BINDIR}/FFT1.x synstrd3.H
	< synstrd3.H ${BINDIR}/FFT1.x > junks1.H
	< junks1.H ${TRANSP} | ${TRANSP} plane=23 > $@

elffftd3.H: ${BINDIR}/FFT1.x elfstrd3.H
	< elfstrd3.H ${BINDIR}/FFT1.x > junke1.H
	< junke1.H ${TRANSP} | ${TRANSP} plane=23 > $@

marmshin%.H: marmfftd3.H ${BINDIR}/Infill.x
	< marmfftd3.H ${BINDIR}/Infill.x eps=$* > $@

synshin%.H: synfftd3.H ${BINDIR}/Infill.x
	< synfftd3.H ${BINDIR}/Infill.x eps=$* > $@

elfshin%.H: elffftd3.H ${BINDIR}/Infill.x
	< elffftd3.H ${BINDIR}/Infill.x eps=$* > $@

marmshnd%.H: marmshot3.H marmshin%.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< marmshin$*.H Transp | ${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< marmshot3.H Window n1=10 n3=1 f3=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

synshnd%.H: synshot3.H synshin%.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< synshin$*.H Transp | ${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< synshot3.H Window n1=50 n3=1 f3=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

elfshnd%.H: elfshot3.H elfshin%.H \
	${BINDIR}/Stretch.x ${BINDIR}/FFT2.x ${BINDIR}/Shift.x 
	< elfshin$*.H Transp | ${BINDIR}/Shift.x inv=y> transp.H
	< transp.H ${BINDIR}/FFT2.x > finv.H
	< finv.H ${BINDIR}/Stretch.x inv=1 > junk1.H
	< elfshot3.H Window n1=20 n3=1 f3=1 > junk2.H
	Cat axis=1 junk2.H junk1.H > $@

marmbench.H: marmshot3.H
	< marmshot3.H Window n3=1 f3=1 > $@

marmdiff%.H: marmshnd%.H marmbench.H
	Add scale=1,-1 marmshnd$*.H marmbench.H > junk.H
	Cat marmshnd$*.H junk.H > $@

synbench.H: synshot3.H
	< synshot3.H Window n3=1 f3=1 > $@

syndiff%.H: synshnd%.H synbench.H
	Add scale=1,-1 synshnd$*.H synbench.H > junk.H
	Cat synshnd$*.H junk.H > $@

elfbench.H: elfshot3.H
	< elfshot3.H Window n3=1 f3=1 > $@

elfdiff%.H: elfshnd%.H elfbench.H
	Add scale=1,-1 elfshnd$*.H elfbench.H > junk.H
	Cat elfshnd$*.H junk.H > $@

${RESDIR}/synshot3.v ${RESDIR}/synshot3.v3: synshot3.H
	< synshot3.H Window n3=1 f3=0 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 1" out=shot1.v
	< synshot3.H Window n3=1 f3=1 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 2" out=shot2.v
	< synshot3.H Window n3=1 f3=2 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 3" out=shot3.v
	vp_Movie shot1.v shot2.v shot3.v > ${RESDIR}/synshot3.v3
	vppen vpstyle=n big=n gridnum=2,1 \
	< shot1.v shot3.v > ${RESDIR}/synshot3.v

${RESDIR}/synshotin.v: syndiff0.H syndiff2.H
	< syndiff0.H Byte pclip=99.5 > diff1.A
	< syndiff2.H Byte pclip=99.5 > diff2.A
	< diff1.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="interpolated shot, no reg" out=shot1.v
	< diff1.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error, no reg" out=diff1.v
	< diff2.A Window n3=1 f3=0 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="interpolated shot, reg" out=shot2.v
	< diff2.A Window n3=1 f3=1 | Ta2vplot >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error, reg" out=diff2.v
	vppen vpstyle=n big=n gridnum=2,2 < \
	shot1.v diff1.v shot2.v diff2.v > $@

${RESDIR}/elfshot3.v ${RESDIR}/elfshot3.v3: elfshot3.H
	< elfshot3.H Window n3=1 f3=0 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 1" out=shot1.v
	< elfshot3.H Window n3=1 f3=1 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 2" out=shot2.v
	< elfshot3.H Window n3=1 f3=2 | \
	${OGREY} label2=offset crowd1=0.85 title="Shot 3" out=shot3.v
	vp_Movie shot1.v shot2.v shot3.v > ${RESDIR}/elfshot3.v3
	vppen vpstyle=n big=n gridnum=2,1 \
	< shot1.v shot3.v > ${RESDIR}/elfshot3.v

elfsimp.v: elfsimp.H
	< elfsimp.H Window f1=400 | Byte pclip=99.5 > 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,1 < diff1.v diff2.v > $@

${RESDIR}/elfshotin.v: elfdiff1.3.H elfsimp.H
	< elfdiff1.3.H Window f1=400 n3=1 f3=0 | Grey clip=8973 >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="interpolated shot" out=shot.v
	< elfdiff1.3.H Window f1=400 n3=1 f3=1 | Grey clip=8973 >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="error" out=diff.v
	< elfsimp.H Window f1=400 n3=1 f3=0 | Grey clip=8973 >/dev/null \
	label1=time label2=offset crowd1=0.85 \
	title="(shot1+shot3)/2" out=diff1.v
	< elfsimp.H Window f1=400 n3=1 f3=1 | Grey clip=8973 >/dev/null \