def main():
opts, ifiles = getOptions()
ccpara = CCConfig()
gsac = loadData(ifiles, opts, ccpara)
opts.ccpara = ccpara
ccpara.twcorr = opts.twcorr
ccpara.cchdrs = [opts.ipick, opts.wpick]
# check data coverage, initialize quality factors
checkCoverage(gsac, opts)
initQual(gsac.saclist, opts.ccpara.hdrsel, opts.ccpara.qheaders)
if opts.auto_on:
autoiccs(gsac, opts)
elif opts.auto_on_all:
print ('Selecting all seismograms..')
hdrsel = opts.ccpara.hdrsel
for sacdh in gsac.saclist:
sacdh.selected = True
sacdh.sethdr(hdrsel, 'True ')
autoiccs(gsac, opts)
else:
stkdh, stkdata, quas = ccWeightStack(gsac.saclist, opts)
gsac.stkdh = stkdh
saveData(gsac, opts)
if opts.plotiter:
plotiter(stkdata)
def main():
opts, ifiles = getOptions()
ccpara = CCConfig()
gsac = loadData(ifiles, opts, ccpara)
opts.ccpara = ccpara
ccpara.twcorr = opts.twcorr
ccpara.cchdrs = [opts.ipick, opts.wpick]
# check data coverage, initialize quality factors
checkCoverage(gsac, opts)
initQual(gsac.saclist, opts.ccpara.hdrsel, opts.ccpara.qheaders)
if opts.auto_on:
autoiccs(gsac, opts)
elif opts.auto_on_all:
print('Selecting all seismograms..')
hdrsel = opts.ccpara.hdrsel
for sacdh in gsac.saclist:
sacdh.selected = True
sacdh.sethdr(hdrsel, 'True ')
autoiccs(gsac, opts)
else:
stkdh, stkdata, quas = ccWeightStack(gsac.saclist, opts)
gsac.stkdh = stkdh
saveData(gsac, opts)
if opts.plotiter:
plotiter(stkdata)
def getDataOpts():
'Get SAC Data and Options'
opts, ifiles = getOptions()
pppara = PPConfig()
gsac = loadData(ifiles, opts, pppara)
opts.pppara = pppara
opts.qheaders = pppara.qheaders
opts.qfactors = pppara.qfactors
opts.qweights = pppara.qweights
opts.hdrsel = pppara.hdrsel
opts.pick_on = True
initQual(gsac.saclist, opts.hdrsel, opts.qheaders)
sortSeis(gsac, opts)
return gsac, opts
def getDataOpts(): 'Get SAC Data and Options' opts, ifiles = getOptions() pppara = PPConfig() gsac = loadData(ifiles, opts, pppara) opts.pppara = pppara opts.qheaders = pppara.qheaders opts.qfactors = pppara.qfactors opts.qweights = pppara.qweights opts.hdrsel = pppara.hdrsel opts.pick_on = True initQual(gsac.saclist, opts.hdrsel, opts.qheaders) sortSeis(gsac, opts) return gsac, opts
def getDataOpts():
'Get SAC Data and Options'
opts, ifiles = getOptions()
pppara = PPConfig()
qcpara = QCConfig()
ccpara = CCConfig()
mcpara = MCConfig()
gsac = loadData(ifiles, opts, pppara)
mcpara.delta = opts.delta
opts.qheaders = qcpara.qheaders
opts.qweights = qcpara.qweights
opts.qfactors = qcpara.qfactors
opts.hdrsel = qcpara.hdrsel
opts.fstack = ccpara.fstack
ccpara.qqhdrs = qcpara.qheaders
ccpara.twcorr = opts.twcorr
# find phase:
if opts.phase is None:
phase = findPhase(ifiles[0])
print ('Found phase to be: ' + phase + '\n')
else:
phase = opts.phase
mcpara.phase = phase
opts.qcpara = qcpara
opts.ccpara = ccpara
opts.mcpara = mcpara
opts.pppara = pppara
# more options:
opts.upylim_on = False
opts.twin_on = True
opts.sort_on = True
opts.pick_on = True
opts.zero_on = False
opts.nlab_on = True
opts.ynormtwin_on = True
#checkCoverage(gsac, opts)
initQual(gsac.saclist, opts.hdrsel, opts.qheaders)
return gsac, opts
def main():
opts, ifiles = getOptions()
mcpara = MCConfig()
gsac = loadData(ifiles, opts, mcpara)
if opts.phase is None:
phase = findPhase(ifiles[0])
print('Found phase to be: ' + phase + '\n')
mcpara.phase = phase
opts.fstack = mcpara.fstack
if opts.filemode == 'sac' and os.path.isfile(opts.fstack):
print('Read array stack file: ' + opts.fstack)
gsac.stkdh = SacDataHdrs(opts.fstack, opts.delta)
getParams(gsac, mcpara, opts)
if opts.allseis_on:
solution = mccc(gsac, mcpara)
gsac.selist, gsac.delist = gsac.saclist, []
else:
initQual(gsac.saclist, mcpara.hdrsel, [])
gsac.selist, gsac.delist = seleSeis(gsac.saclist)
solution = mccc(gsac, mcpara)
saveData(gsac, opts)
def main():
opts, ifiles = getOptions()
mcpara = MCConfig()
gsac = loadData(ifiles, opts, mcpara)
if opts.phase is None:
phase = findPhase(ifiles[0])
print ('Found phase to be: ' + phase + '\n')
mcpara.phase = phase
opts.fstack = mcpara.fstack
if opts.filemode == 'sac' and os.path.isfile(opts.fstack):
print ('Read array stack file: ' + opts.fstack )
gsac.stkdh = SacDataHdrs(opts.fstack, opts.delta)
getParams(gsac, mcpara, opts)
if opts.allseis_on:
solution = mccc(gsac, mcpara)
gsac.selist, gsac.delist = gsac.saclist, []
else:
initQual(gsac.saclist, mcpara.hdrsel, [])
gsac.selist, gsac.delist = seleSeis(gsac.saclist)
solution = mccc(gsac, mcpara)
saveData(gsac, opts)
def getDataOpts():
'Get SAC Data and Options'
opts, ifiles = getOptions()
pppara = PPConfig()
qcpara = QCConfig()
ccpara = CCConfig()
mcpara = MCConfig()
gsac = loadData(ifiles, opts, pppara)
# set defaults
filterParameters = {}
filterParameters['apply'] = False
filterParameters['advance'] = False
filterParameters['band'] = 'bandpass'
filterParameters['lowFreq'] = 0.05
filterParameters['highFreq'] = 0.25
filterParameters['order'] = 2
filterParameters['reversepass'] = False
opts.filterParameters = filterParameters
# override defaults if already set in SAC files
firstSacdh = gsac.saclist[0]
if hasattr(firstSacdh, 'user0'):
filterParameters['lowFreq'] = firstSacdh.user0
if hasattr(firstSacdh, 'user1'):
filterParameters['highFreq'] = firstSacdh.user1
if hasattr(firstSacdh, 'kuser0'):
filterParameters['band'] = firstSacdh.kuser0
if hasattr(firstSacdh, 'kuser1'):
filterParameters['order'] = int(firstSacdh.kuser1)
mcpara.delta = opts.delta
opts.qheaders = qcpara.qheaders
opts.qweights = qcpara.qweights
opts.qfactors = qcpara.qfactors
opts.hdrsel = qcpara.hdrsel
opts.fstack = ccpara.fstack
ccpara.qqhdrs = qcpara.qheaders
ccpara.twcorr = opts.twcorr
# find phase:
if opts.phase is None:
phase = findPhase(ifiles[0])
print('Found phase to be: ' + phase + '\n')
else:
phase = opts.phase
mcpara.phase = phase
opts.qcpara = qcpara
opts.ccpara = ccpara
opts.mcpara = mcpara
opts.pppara = pppara
# more options:
opts.upylim_on = False
opts.twin_on = True
opts.sort_on = True
opts.pick_on = True
opts.zero_on = False
opts.nlab_on = True
opts.ynormtwin_on = True
#checkCoverage(gsac, opts)
initQual(gsac.saclist, opts.hdrsel, opts.qheaders)
return gsac, opts
def getDataOpts():
'Get SAC Data and Options'
opts, ifiles = getOptions()
pppara = PPConfig()
qcpara = QCConfig()
ccpara = CCConfig()
mcpara = MCConfig()
gsac = loadData(ifiles, opts, pppara)
# set defaults
filterParameters = {}
filterParameters['apply'] = False
filterParameters['advance'] = False
filterParameters['band'] = 'bandpass'
filterParameters['lowFreq'] = 0.05
filterParameters['highFreq'] = 0.25
filterParameters['order'] = 2
filterParameters['reversepass'] = False
opts.filterParameters = filterParameters
# override defaults if already set in SAC files
firstSacdh = gsac.saclist[0]
if hasattr(firstSacdh, 'user0'):
filterParameters['lowFreq'] = firstSacdh.user0
if hasattr(firstSacdh, 'user1'):
filterParameters['highFreq'] = firstSacdh.user1
if hasattr(firstSacdh, 'kuser0'):
filterParameters['band'] = firstSacdh.kuser0
if hasattr(firstSacdh, 'kuser1'):
filterParameters['order'] = int(firstSacdh.kuser1)
mcpara.delta = opts.delta
opts.qheaders = qcpara.qheaders
opts.qweights = qcpara.qweights
opts.qfactors = qcpara.qfactors
opts.hdrsel = qcpara.hdrsel
opts.fstack = ccpara.fstack
ccpara.qqhdrs = qcpara.qheaders
ccpara.twcorr = opts.twcorr
# find phase:
if opts.phase is None:
phase = findPhase(ifiles[0])
print ('Found phase to be: ' + phase + '\n')
else:
phase = opts.phase
mcpara.phase = phase
opts.qcpara = qcpara
opts.ccpara = ccpara
opts.mcpara = mcpara
opts.pppara = pppara
# more options:
opts.upylim_on = False
opts.twin_on = True
opts.sort_on = True
opts.pick_on = True
opts.zero_on = False
opts.nlab_on = True
opts.ynormtwin_on = True
#checkCoverage(gsac, opts)
initQual(gsac.saclist, opts.hdrsel, opts.qheaders)
return gsac, opts