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