def semb(workspace,**kwargs):
vels = kwargs['velocities']
nvels = vels.size
ns = kwargs['ns']
result = np.zeros((nvels,ns),'f')
for v in range(nvels):
panel = workspace.copy()
kwargs['vels'] = np.ones(kwargs['ns'], 'f') * vels[v]
nmo(panel, None, **kwargs)
result[v,:] += np.abs(_stack_gather(panel)['trace'])
pylab.imshow(result.T, aspect='auto', extent=(min(vels), max(vels),kwargs['ns']*kwargs['dt'],0.), cmap='gist_heat')
pylab.xlabel('velocity')
pylab.ylabel('time')
pylab.colorbar()
pylab.show()
def semb(workspace,**kwargs):
vels = kwargs['velocities']
nvels = vels.size
ns = kwargs['ns']
result = np.zeros((nvels,ns),'f')
for v in range(nvels):
panel = workspace.copy()
kwargs['vels'] = np.ones(kwargs['ns'], 'f') * vels[v]
nmo(panel, None, **kwargs)
result[v,:] += np.abs(_stack_gather(panel)['trace'])
pylab.imshow(result.T, aspect='auto', extent=(min(vels), max(vels),kwargs['ns']*kwargs['dt'],0.), cmap='gist_heat')
pylab.xlabel('velocity')
pylab.ylabel('time')
pylab.colorbar()
pylab.show()
#apply LMO
#~ print "applying lmo"
#~ params['lmo'] =1000.0
#~ lmo(workspace, None, **params)
#~ workspace['trace'][:,:30] *= 0
#~ workspace['trace'][:,1850:] *= 0
#~ params['lmo'] = -1000.0
#~ lmo(workspace, None, **params)
#apply our NMO
print "applying nmo"
params['smute'] = 30
v = [3000]
t = [0.5]
params['vels'] = toolbox.build_vels(t, v, ns=params['ns'])
nmo(workspace, None, **params)
#~ #apply AGC
toolbox.agc(workspace, None, **params)
#~ #stack
print "stacking"
stack(workspace, 'stack1.su', **params)
#view
params['primary'] = None
toolbox.display('stack1.su', None, **params)
pylab.show()
tar(workspace, None, **params)
#apply LMO
#~ print "applying lmo"
#~ params['lmo'] =1000.0
#~ lmo(workspace, None, **params)
#~ workspace['trace'][:,:30] *= 0
#~ workspace['trace'][:,1850:] *= 0
#~ params['lmo'] = -1000.0
#~ lmo(workspace, None, **params)
#apply our NMO
print "applying nmo"
params['smute'] = 30
v = [3000]
t = [0.5]
params['vels'] = toolbox.build_vels(t, v, ns=params['ns'])
nmo(workspace, None, **params)
#~ #apply AGC
toolbox.agc(workspace, None, **params)
#~ #stack
print "stacking"
stack(workspace, 'stack1.su', **params)
#view
params['primary'] = None
toolbox.display('stack1.su', None, **params)
pylab.show()