#~ data, params = toolbox.initialise("geometries.su")
#trim data
#~ params['ns'] = 1500
#~ data = toolbox.slice(data, None, **params)
#~ data.tofile("geom_short.su")
#initialise dataset
data, params = toolbox.initialise("geom_short.su")
#agc
#~ toolbox.agc(data, None, **params)
params['gamma'] = 1.5
toolbox.tar(data, None, **params)
kills = [270, 300, 374, 614] #fldr
mask = toolbox.build_mask(data['fldr'], kills)
data = data[mask]
data.tofile("prepro.su")
#display
#~ params['primary'] = 'fldr'
#~ params['secondary'] = 'tracf'
#~ params['wiggle'] = True
#~ toolbox.display(data, **params)
#~ pylab.show()
#AGC
toolbox.agc(data, None, **params)
#select the velocities we want to check
params['velocities'] = np.arange(1000, 6000, 100)
params['smoother'] = 25
params['smute'] = 150
#calculate our analysis locations
cdps = np.unique(data['cdp'])
sg_centers = np.arange(800, 3000, 200)
#loop over the supergathers
for sg in sg_centers:
gathers = range(sg - 5, sg + 6)
mask = toolbox.build_mask(data['cdp'], gathers)
supergather = data[mask]
supergather.sort(order=['offset'])
#display it
toolbox.semb(supergather, **params)
#this is out holding dictionary - we'll put our velocity solutions in it.
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[800] = (2456.0, 0.153), (2772.1, 0.413), (3003.2, 0.612), (
#trim data
#~ params['ns'] = 1500
#~ data = toolbox.slice(data, None, **params)
#~ data.tofile("geom_short.su")
#initialise dataset
data, params = toolbox.initialise("geom_short.su")
#agc
#~ toolbox.agc(data, None, **params)
params['gamma'] = 1.5
toolbox.tar(data, None, **params)
kills = [270, 300, 374, 614] #fldr
mask = toolbox.build_mask(data['fldr'], kills)
data = data[mask]
data.tofile("prepro.su")
#display
#~ params['primary'] = 'fldr'
#~ params['secondary'] = 'tracf'
#~ params['wiggle'] = True
#~ toolbox.display(data, **params)
#~ pylab.show()
toolbox.agc(data, None, **params)
#select the velocities we want to check
params['velocities'] = np.arange(1000, 6000, 100)
params['smoother'] = 25
params['smute'] = 150
#calculate our analysis locations
cdps = np.unique(data['cdp'])
sg_centers = np.arange(800, 3000, 200)
#loop over the supergathers
for sg in sg_centers:
gathers = range(sg-5, sg+6)
mask = toolbox.build_mask(data['cdp'], gathers)
supergather = data[mask]
supergather.sort(order=['offset'])
#display it
toolbox.semb(supergather, **params)
#this is out holding dictionary - we'll put our velocity solutions in it.
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[800]= (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[1000]= (2443.9, 0.122), (2662.7, 0.413), (2942.4, 0.980), (3027.5, 1.668), (3890.8, 5.893),
vels[1200]= (2370.9, 0.138), (3051.8, 0.765), (3367.9, 1.209), (3720.5, 1.730), (4012.3, 2.602), (4292.0, 5.862), (4158.3, 3.689), (2711.4, 0.459),
vels[1400]= (2419.5, 0.092), (2820.8, 0.520), (3112.6, 1.026), (3392.2, 1.668), (3562.5, 2.128), (3720.5, 2.755), (3830.0, 3.367), (4085.3, 5.939),
vels[1600]= (2285.8, 0.107), (2662.7, 0.520), (3173.4, 1.087), (3501.7, 1.730), (3647.6, 2.526), (3781.3, 3.184), (4097.5, 5.862),
vels[1800]= (2322.3, 0.153), (2638.4, 0.566), (3015.3, 1.133), (3355.8, 1.806), (3538.2, 2.556), (3708.4, 3.260), (4000.2, 5.893),
