toolbox.agc(shot, None, **params)
params['primary'] = 'fldr'
params['secondary'] = 'tracf'
params['wiggle'] = True
toolbox.display(shot, **params)
#fk plot
params['dx'] = 33.5 #m
#~ toolbox.fk_view(shot, **params)
#~ #fk filter design
params['fkVelocity'] = 2000
params['fkSmooth'] = 20
params['fkFilter'] = toolbox.fk_design(shot, **params)
shot = toolbox.fk_filter(shot, None, **params)
toolbox.display(shot, **params)
##############end of testing
#~ data, nparams = toolbox.initialise("prepro.su")
#~ toolbox.agc(data, None, **params)
#~ data = toolbox.fk_filter(data, None, **params)
#~ #nmo
#~ params['vels'] = np.fromfile('vels_full.bin').reshape(-1, params['ns'])
#~ params['smute'] = 150
#~ toolbox.nmo(data, None, **params)
#~ data.tofile("fk_nmo_gathers.su")
#~ toolbox.agc(data, None, **params)
#~ #stack
#~ stack = toolbox.stack(data, None, **params)
import toolbox
import numpy as np
import pylab
stack, params = toolbox.initialise("fk_stack.su")
stack['fldr'] = 1
params['dx'] = 33.5/2.0 #m
params['fkVelocity'] = 6000
params['fkSmooth'] = 20
params['fkFilter'] = toolbox.fk_design(stack, **params)
stack = toolbox.fk_filter(stack, None, **params)
#bandpass
params['lowcut'] = 10.0
params['highcut'] = 100.0
toolbox.bandpass(stack, None, **params)
stack.tofile("model_filtered.su")
#display
#~ params['primary'] = None
#~ params['secondary'] = 'cdp'
#~ toolbox.display(stack, **params)
#~ pylab.show()