-
Notifications
You must be signed in to change notification settings - Fork 0
/
mseedinfo.py
32 lines (24 loc) · 1.09 KB
/
mseedinfo.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
# Just display info on miniseed file
from obspy import read
from bokeh.plotting import figure, output_file, show
from bokeh.palettes import Set1_3
input_path = 'E:\project_amanah\S3\earthquake\dataset-avicenna\panyaungan2018\panyaungan2018-2.miniseed'
st = read(input_path)
print(st.__str__(extended=True))
print(st[0].stats)
print(st[3].stats)
print(st[0].data)
# output to static HTML file
output_file("waveform.html")
# create a new plot with a title and axis labels
p = figure(title="%s.%s %s Hz %s" % (st[0].stats.network, st[0].stats.station, st[0].stats.sampling_rate, st[0].stats.starttime),
x_axis_label='sample', y_axis_label='amplitude',
plot_height=120)
p.sizing_mode = 'scale_width'
# add a line renderer with legend and line thickness
#p.line(x, y, legend="Temp.", line_width=2)
p.line(range(len(st[0].data)), st[0].data, color=Set1_3[0], legend=st[0].stats.channel)
p.line(range(len(st[1].data)), st[1].data, color=Set1_3[1], legend=st[1].stats.channel)
p.line(range(len(st[2].data)), st[2].data, color=Set1_3[2], legend=st[2].stats.channel)
# show the results
show(p)