def process_energyratio (windowsize = 0.1):
writenow('PROCESS: Computing STA/LTA for traces.')
writenow(' \twindowsize = %f s'%windowsize)
params = dbase['params']
data = dbase['data_current'].reshape(params['flatdatadims'])
dt = params['dt']
windowsamps = int(windowsize/dt)
erdata = analysis.energyratio(data, windowsamps).reshape(params['datadims'])
if (keepdata):
dbase['data_energyratio'] = erdata
dbase['data_current'] = erdata
def process_energyratio(windowsize=0.1):
writenow('PROCESS: Computing STA/LTA for traces.')
writenow(' \twindowsize = %f s' % windowsize)
params = dbase['params']
data = dbase['data_current'].reshape(params['flatdatadims'])
dt = params['dt']
windowsamps = int(windowsize / dt)
erdata = analysis.energyratio(data,
windowsamps).reshape(params['datadims'])
if (keepdata):
dbase['data_energyratio'] = erdata
dbase['data_current'] = erdata
wlen = DEFAULT_WLEN
if len(sys.argv) >= 3:
wlen = int(sys.argv[2])
shift = -wlen
if len(sys.argv) >= 4:
shift = -int(sys.argv[3])
inifile = glob.glob('*.ini')[0]
projnm = inifile.split('.')[0]
ini = readini(inifile)
ns = ini['ns']
nr = ini['nr']
sf = SEGYFile(infile)
dt = sf.bhead['hdt'] * TUNIT
picks = (np.argmax(energyratio(sf[:], wlen), axis=1) + shift)
picks.shape = (ns, nr)
picks = picks * dt
with open('%s.picks'%(projnm,), 'w') as fp:
fp.writelines(['%d %d %f\n'%(isrc+1, irec+1, picks[isrc, irec]) for isrc in xrange(ns) for irec in xrange(nr)])