Revisiting SEP tour with Madagascar and SCons - Revision history

Sfomel: /* SConstruct file */

2014-04-08T21:57:08Z

SConstruct file

← Older revision		Revision as of 21:57, 8 April 2014
Line 278:		Line 278:
	Here is a complete listing of the <tt>SConstruct</tt> file used in this		Here is a complete listing of the <tt>SConstruct</tt> file used in this
	example.		example.
	<python>		<syntaxhighlight lang="python">
	#########################################################		#########################################################
	# Setting up		# Setting up
Line 354:		Line 354:

	End()		End()
	</~~python~~>		</syntaxhighlight>

Sfomel: /* Advanced plotting */

2014-04-08T21:56:31Z

Advanced plotting

← Older revision		Revision as of 21:56, 8 April 2014
Line 222:		Line 222:
	Start by creating common plotting plotting arguments and plotting the		Start by creating common plotting plotting arguments and plotting the
	data in greyscale.		data in greyscale.
	<python>		<syntaxhighlight lang="python">
	plotpar = plotpar+' min1=.4 max1=.8 max2=1. min2=.05 poly=n'		plotpar = plotpar+' min1=.4 max1=.8 max2=1. min2=.05 poly=n'

	Plot('grey','windowed',		Plot('grey','windowed',
	'grey wheretitle=t wherexlabel=b' + plotpar)		'grey wheretitle=t wherexlabel=b' + plotpar)
	</~~python~~>		</syntaxhighlight>

	Next, plot the wiggle traces twice: the fist time, using thick black		Next, plot the wiggle traces twice: the fist time, using thick black
	lines (<tt>plotcol=0 plotfat=10</tt>), and the second time, using		lines (<tt>plotcol=0 plotfat=10</tt>), and the second time, using
	thinner white lines (<tt>plotcol=7 plotfat=5</tt>).		thinner white lines (<tt>plotcol=7 plotfat=5</tt>).
	<python>		<syntaxhighlight lang="python">
	Plot('wiggle1','windowed',		Plot('wiggle1','windowed',
	'wiggle plotcol=0 plotfat=10' + plotpar)		'wiggle plotcol=0 plotfat=10' + plotpar)
	Plot('wiggle2','windowed',		Plot('wiggle2','windowed',
	'wiggle plotcol=7 plotfat=3' + plotpar)		'wiggle plotcol=7 plotfat=3' + plotpar)
	</~~python~~>		</syntaxhighlight>

	The plots are combined by overlaying or by putting them side by side.		The plots are combined by overlaying or by putting them side by side.
	<python>		<syntaxhighlight lang="python">
	Result('overplot','grey wiggle1 wiggle2','Overlay')		Result('overplot','grey wiggle1 wiggle2','Overlay')
	Result('sidebyside','grey wiggle2','SideBySideIso')		Result('sidebyside','grey wiggle2','SideBySideIso')
	</~~python~~>		</syntaxhighlight>

	The resultant plots are shown in the figures.		The resultant plots are shown in the figures.

Sfomel: /* Normal Moveout */

2014-04-08T21:54:30Z

Normal Moveout

← Older revision		Revision as of 21:54, 8 April 2014
Line 194:		Line 194:
	windowed data		windowed data
	and pipes the result to a wiggle plotting command:		and pipes the result to a wiggle plotting command:
	<python>		<syntaxhighlight lang="python">
	Result('nmo','windowed',		Result('nmo','windowed',
	'''		'''
Line 200:		Line 200:
	wiggle pclip=100 max1=0.6 poly=y		wiggle pclip=100 max1=0.6 poly=y
	''')		''')
	</~~python~~>		</syntaxhighlight>

	Running <tt>scons -Q nmo.view</tt> produces		Running <tt>scons -Q nmo.view</tt> produces

Sfomel: /* Resampling */

2014-04-08T21:54:03Z

Resampling

← Older revision		Revision as of 21:54, 8 April 2014
Line 147:		Line 147:
	using Fourier interpolation.		using Fourier interpolation.
	Subsampling is accomplished with <tt>sfwindow</tt>.		Subsampling is accomplished with <tt>sfwindow</tt>.
	<python>		<syntaxhighlight lang="python">
	# decimate time axis by two		# decimate time axis by two
	Flow('subsampled','windowed','window j1=2')		Flow('subsampled','windowed','window j1=2')
	</~~python~~>		</syntaxhighlight>

	Running <tt>scons -Q subsampled.rsf</tt> produces		Running <tt>scons -Q subsampled.rsf</tt> produces
Line 168:		Line 168:
	#Inverse Fourier transform from frequency to time.		#Inverse Fourier transform from frequency to time.
	All three steps are conveniently combined into one using pipes.		All three steps are conveniently combined into one using pipes.
	<python>		<syntaxhighlight lang="python">
	# sinc interpolation in the Fourier domain		# sinc interpolation in the Fourier domain
	Flow('resampled','subsampled',		Flow('resampled','subsampled',
	'fft1 \| pad n1=102 \| fft1 inv=y opt=n \| window max1=0.8')		'fft1 \| pad n1=102 \| fft1 inv=y opt=n \| window max1=0.8')
	</~~python~~>		</syntaxhighlight>

	Why do we pad the Fourier domain to 102? The time length of the original data		Why do we pad the Fourier domain to 102? The time length of the original data
Line 181:		Line 181:
	the figure. Comparing this result with		the figure. Comparing this result with
	the previous plots, we can verify a fairly accurate reconstruction.		the previous plots, we can verify a fairly accurate reconstruction.
	<python>		<syntaxhighlight lang="python">
	Result('resampled','wiggle title=Resampled' + plotpar)		Result('resampled','wiggle title=Resampled' + plotpar)
	</~~python~~>		</syntaxhighlight>

	[[Image:resampled.png\|frame\|center\|To see this figure on your screen, run <tt>scons resampled.view</tt>]]		[[Image:resampled.png\|frame\|center\|To see this figure on your screen, run <tt>scons resampled.view</tt>]]

Sfomel: /* Windowing and plotting */

2014-04-08T21:53:03Z

Windowing and plotting

← Older revision		Revision as of 21:53, 8 April 2014
Line 81:		Line 81:
	Our next task is to window and plot a significant portion of the data. Add the		Our next task is to window and plot a significant portion of the data. Add the
	following line to the <tt>SConstruct</tt> file:		following line to the <tt>SConstruct</tt> file:
	<python>		<syntaxhighlight lang="python">
	Flow('windowed','Txx.HH','window n2=10 min1=0.4 max1=0.8')		Flow('windowed','Txx.HH','window n2=10 min1=0.4 max1=0.8')
	</~~python~~>		</syntaxhighlight>

	The window command selects the first ten traces and the time window between		The window command selects the first ten traces and the time window between
	0.4 and 0.8 seconds.		0.4 and 0.8 seconds.
	We will plot the windowed data with three different plotting programs.		We will plot the windowed data with three different plotting programs.
	<python>		<syntaxhighlight lang="python">
	plotpar = '''		plotpar = '''
	transp=y poly=y yreverse=y pclip=100 nc=100 allpos=n		transp=y poly=y yreverse=y pclip=100 nc=100 allpos=n
Line 95:		Line 95:
	for plot in ('wiggle','contour','grey'):		for plot in ('wiggle','contour','grey'):
	Result(plot,'windowed',plot + plotpar)		Result(plot,'windowed',plot + plotpar)
	</~~python~~>		</syntaxhighlight>

	For convenience, plotting parameters are put in a string called		For convenience, plotting parameters are put in a string called

Sfomel: /* Displaying the data */

2014-04-08T21:52:22Z

Displaying the data

← Older revision		Revision as of 21:52, 8 April 2014
Line 51:		Line 51:

	Add the following line to the <tt>SConstruct</tt> file:		Add the following line to the <tt>SConstruct</tt> file:
	<python>		<syntaxhighlight lang="python">
	Result('wiggle0','Txx.HH','wiggle')		Result('wiggle0','Txx.HH','wiggle')
	</~~python~~>		</syntaxhighlight>

Sfomel: /* Obtaining the test data */

2014-04-08T21:51:52Z

Obtaining the test data

← Older revision		Revision as of 21:51, 8 April 2014
Line 37:		Line 37:
	===Obtaining the test data===		===Obtaining the test data===
	Add a Fetch command as follows:		Add a Fetch command as follows:
	<python>		<syntaxhighlight lang="python">
	Fetch('Txx.HH','septour')		Fetch('Txx.HH','septour')
	</~~python~~>		</syntaxhighlight>

	Now, by running		Now, by running
Line 47:		Line 47:
	you can instruct SCons to connect to an anonymous data server and extract		you can instruct SCons to connect to an anonymous data server and extract
	(fetch) the data file "Txx.HH" from the "septour" directory.		(fetch) the data file "Txx.HH" from the "septour" directory.

	===Displaying the data===		===Displaying the data===

Sfomel: /* Setting up */

2014-04-08T21:51:17Z

Setting up

← Older revision		Revision as of 21:51, 8 April 2014
Line 25:		Line 25:
	Open a file named "SConstruct" in your favorite editor, start it with a		Open a file named "SConstruct" in your favorite editor, start it with a
	line		line
	<python>		<syntaxhighlight lang="python">
	from rsf.proj import *		from rsf.proj import *
	</~~python~~>		</syntaxhighlight>
	and end it with a line		and end it with a line
	<python>		<syntaxhighlight lang="python">
	End()		End()
	</~~python~~>		</syntaxhighlight>

	The first line tells Python to load the RSF project module.		The first line tells Python to load the RSF project module.

Sfomel at 20:30, 10 August 2013

2013-08-10T20:30:13Z

← Older revision		Revision as of 20:30, 10 August 2013
Line 1:		Line 1:
	<center><font size="-1">''This page was created from the LaTeX source in [http://~~rsf.svn.~~sourceforge.net/~~viewvc~~/rsf/trunk/book/rsf/rsf/tour.tex~~?view=markup~~ book/rsf/rsf/tour.tex] using [[latex2wiki]]''</font></center>		<center><font size="-1">''This page was created from the LaTeX source in [http://sourceforge.net/p/rsf/code/HEAD/tree/trunk/book/rsf/rsf/tour.tex book/rsf/rsf/tour.tex] using [[latex2wiki]]''</font></center>

	Many appreciative users were introduced to SEPlib		Many appreciative users were introduced to SEPlib

Sfomel: /* Acknowledgments */

2011-07-09T15:03:12Z

Acknowledgments

← Older revision		Revision as of 15:03, 9 July 2011
Line 266:		Line 266:
	Thanks to Joe Dellinger and Sándor Tálas for creating "SEP		Thanks to Joe Dellinger and Sándor Tálas for creating "SEP
	tour" and to James Rickett for updating it. Several generations of		tour" and to James Rickett for updating it. Several generations of
	SEP students contributed to SEPlib. We ~~try~~ to preserve all their good		SEP students contributed to SEPlib. We tried to preserve all their good
	ideas when refactoring SEPlib into ~~RSF~~.		ideas when refactoring SEPlib into Madagascar.

	The test dataset used in this paper is courtesy of Beltram Nolte and L. Neil Frazer.		The test dataset used in this paper is courtesy of Beltram Nolte and L. Neil Frazer.