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