def main2(): data = IPOC(xcorr_append='/1bit', use_local_LVC=True) t1 = UTC('2010-01-01') stream0_1 = data.getRawStream(t1, 'PB01', component='Z') stream0_2 = data.getRawStream(t1, 'PB02', component='Z') stream2_1 = data.getStream(t1, 'PB01', component='Z') stream2_2 = data.getStream(t1, 'PB02', component='Z') plot_streams([stream0_1, stream0_2, stream2_1, stream2_2], [None, None, 0.1, 0.1]) plotPSD([stream0_1, stream0_2, stream2_1, stream2_2], 4096) ipshell()
from sito import read from sito.imaging import plotTrace stream = read(args.file_station) plotTrace(stream, **kwargs) else: from sito.imaging import plotTrace2 station = args.file_station if station.startswith('PB') or station == 'LVC': from sito.data import IPOC data = IPOC(xcorr_append=args.xcorr_append) elif station == 'PKD': from sito.data import Parkfield data = Parkfield(xcorr_append=args.xcorr_append) else: raise argparse.ArgumentError('Not a valid station name') day = UTCDateTime(args.date) if args.xcorr_append is None: stream = data.getRawStream(day, station, component=args.component) else: stream = data.getStream(day, station, component=args.component) if stream[0].stats.is_fft: stream.ifft() plotTrace(stream, component=args.component, **kwargs) if args.save is None: from matplotlib.pyplot import show show()
if args.date is None: from sito import read from sito.imaging import plotTrace stream = read(args.file_station) plotTrace(stream, **kwargs) else: from sito.imaging import plotTrace2 station = args.file_station if station.startswith('PB') or station == 'LVC': from sito.data import IPOC data = IPOC(xcorr_append=args.xcorr_append) elif station == 'PKD': from sito.data import Parkfield data = Parkfield(xcorr_append=args.xcorr_append) else: raise argparse.ArgumentError('Not a valid station name') day = UTCDateTime(args.date) if args.xcorr_append is None: stream = data.getRawStream(day, station, component=args.component) else: stream = data.getStream(day, station, component=args.component) if stream[0].stats.is_fft: stream.ifft() plotTrace(stream, component=args.component, **kwargs) if args.save is None: from matplotlib.pyplot import show show()