#apply tar
params['gamma'] = 5
toolbox.tar(workspace, None, **params)
#copy vels from previous exercise
vels = {}
vels[225] = (0.06, 1537.38) , (0.28, 2876.21) , (0.87, 4608.10)
vels[270] = (0.05, 1525.09) , (0.18, 2483.16) , (0.36, 3171.00) , (0.66, 4079.93) , (0.98, 4816.90)
vels[315] = (0.04, 1365.42) , (0.14, 2728.82) , (0.22, 3134.15) , (0.57, 4116.78) , (0.74, 4571.25) , (0.97, 5013.43)
vels[360] = (0.04, 1697.05) , (0.10, 2520.01) , (0.21, 2937.62) , (0.43, 3244.70) , (0.64, 3981.67) , (0.98, 4239.61)
vels[405] = (0.06, 1439.11) , (0.27, 2753.38) , (0.49, 3957.10) , (0.97, 5381.92)
vels[450] = (0.06, 1340.85) , (0.41, 2741.10) , (0.52, 3625.47) , (0.02, 1144.32) , (0.29, 3060.45) , (0.54, 3711.45) , (0.97, 4313.31)
vels[495] = (0.04, 1611.07) , (0.11, 3072.74) , (0.23, 3318.39) , (0.35, 3772.86) , (0.48, 3981.67) , (0.94, 5099.41)
vels[539] = (0.04, 2028.69) , (0.11, 3072.74) , (0.32, 3883.41) , (0.51, 4485.27) , (0.96, 5222.24)
vels[584] = (0.06, 1623.36) , (0.20, 2495.44) , (0.32, 3121.87) , (0.95, 4411.57)
#build vels
vels = toolbox.build_vels(vels, **params)
params['primary'] = None
params['smute'] = 30
params['vels'] = vels
v100 = toolbox.co_nmo(workspace, None, **params)
toolbox.agc(v100, None, **params)
section100 = toolbox.stack(v100, None, **params)
toolbox.cp(section100, 'stack100.su', None)
toolbox.display(section100, None, **params)
pylab.show()
import toolbox
import pylab
import numpy as np
model, mparams = toolbox.initialise("model_filtered.su")
gathers, params = toolbox.initialise("fk_nmo_gathers.su")
params['model'] = model
params['gate'] = (.4, 1.0) #seconds
params['maxshift'] = 4 #samples
toolbox.trim(gathers, None, **params)
toolbox.apply_statics(gathers, None, **params)
stack = toolbox.stack(gathers, None, **params)
params['gamma'] = -1
toolbox.tar(stack, None, **params)
stack.tofile("trim_stack2.su")
cdps = np.unique(data['cdp'])
#recreate original velocity field
#~ vels = {}
#~ vels[753]= (2456.0, 0.153), (2772.1, 0.413), (3003.2, 0.612), (3076.1, 0.704), (3270.7, 1.056), (3367.9, 1.668), (3538.2, 2.204), (3671.9, 3.566), (3915.1, 5.908),
#~ vels[3056]=(2456.0, 0.153), (2772.1, 0.413), (3003.2, 0.612), (3076.1, 0.704), (3270.7, 1.056), (3367.9, 1.668), (3538.2, 2.204), (3671.9, 3.566), (3915.1, 5.908),
#~ params['cdp'] = cdps
#~ params['vels'] = toolbox.build_vels(vels, **params)
#~ np.array(params['vels']).tofile('vels_initial.bin')
params['vels'] = np.fromfile('vels_initial.bin').reshape(-1, params['ns'])
#~ pylab.imshow(params['vels'].T, aspect='auto')
#~ pylab.colorbar()
#~ pylab.show()
#agc and stack
toolbox.agc(data, None, **params)
params['smute'] = 200
toolbox.nmo(data, None, **params)
stack = toolbox.stack(data, None, **params)
params['gamma'] = -1
toolbox.tar(stack, None, **params)
stack.tofile('field_stack.su')
#display
params['primary'] = None
params['secondary'] = 'cdp'
toolbox.display(stack, **params)
pylab.show()
import toolbox
import numpy as np
import pylab
#~ for mute in range(110,200, 10):
#initialise
data, params = toolbox.initialise('prepro.su')
#load vels
params['vels'] = np.fromfile('vels_full.bin').reshape(-1, params['ns'])
params['smute'] = 150 #mute
#normal moveout correction
toolbox.nmo(data, None, **params)
#AGC
toolbox.agc(data, None, **params)
#stack
stack = toolbox.stack(data, None, **params)
params['gamma'] = -1
toolbox.tar(stack, None, **params)
data['tstat'] /= 2
toolbox.apply_statics(stack, None, **params)
#~ stack.tofile("smute_%d.su" %mute)
stack.tofile("1st_vels_stack_elev.su")
# i think I like 150
vels[584] = (0.06, 1623.36) , (0.20, 2495.44) , (0.32, 3121.87) , (0.95, 4411.57)
#build vels
vels = toolbox.build_vels(vels, **params)
params['primary'] = None
params['highcut'] = 100
params['lowcut'] = 30
params['smute'] = 30
params['vels'] = vels
v100 = toolbox.co_nmo(workspace, None, **params)
toolbox.agc(v100, None, **params)
section100 = toolbox.stack(v100, None, **params)
toolbox.bandpass(section100, None, **params)
toolbox.display(section100, None, **params)
params['vels'] = vels * .9
v90 = toolbox.co_nmo(workspace, None, **params)
toolbox.agc(v90, None, **params)
section90 = toolbox.stack(v90, None, **params)
toolbox.bandpass(section90, None, **params)
toolbox.display(section90, None, **params)
params['vels'] = vels *1.1
v110 = toolbox.co_nmo(workspace, None, **params)
toolbox.agc(v110, None, **params)
section110 = toolbox.stack(v110, None, **params)
#build vels
vels = toolbox.build_vels(vels, **params)
#
params['vels'] = vels
params['smute'] = 30
#normal moveout
v100 = toolbox.co_nmo(workspace, None, **params)
#agc
toolbox.agc(v100, None, **params)
#trace mix
params['mix'] = 5
toolbox.trace_mix(workspace, None, **params)
#stack
section100 = toolbox.stack(v100, None, **params)
#bandpass filter
params['lowcut'] = 30
params['highcut'] = 100
toolbox.bandpass(section100, None, **params)
toolbox.cp(section100, 'final_stack.su', **params)
#display
toolbox.display(section100, None, **params)
pylab.show()
import toolbox
import numpy as np
import pylab
#initialise
data, params = toolbox.initialise('prepro.su')
#load vels
params['vels'] = np.fromfile('vels_full.bin').reshape(-1, params['ns'])
params['smute'] = 200
#normal moveout correction
nmo = toolbox.nmo(data, None, **params)
#AGC
toolbox.agc(nmo, None, **params)
#stack
stack = toolbox.stack(nmo, None, **params)
params['gamma'] = -1
toolbox.tar(stack, None, **params)
stack.tofile("field_stack2.su")
params['primary'] = None
params['secondary'] = 'cdp'
toolbox.display(stack, **params)
pylab.show()
