for てごし

seis.py

import obspy
from obspy import read
from matplotlib import pylab as plt
import numpy as np
import function

samplerate = 100 #サンプリング周波数

st = read('1912011200_N.MMOH_U.s', debug_headers=True)#このファイルはメールで送るます
data = st[0].data

# tは時間方向での離散値．stの長さをサンプルレートで切っているので，サンプル数分の配列になってます．
t = np.arange(0, len(data))/samplerate

plt.plot(t,data)
plt.show()