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()