1

.sac ファイルの情報を使用して、グリッドに対してプロットする必要があります。さまざまな ObsPy 関数を使用して地震記録をプロットできることは知っていますst.plot()が、グリッドに対して取得できないようです。また、ここに示されている例「Python でプロットにグリッドを描画するにはどうすればよいですか?」 に従ってみましたが、UTCDatetime を使用するように x 軸を構成しようとすると問題が発生します。私はこの種のPythonとプログラミングが初めてなので、アドバイス/ヘルプをいただければ幸いです。

使用されるさまざまなリソース: " http://docs.obspy.org/tutorial/code_snippets/reading_seismograms.html "

" http://docs.obspy.org/packages/autogen/obspy.core.stream.Stream.plot.html#obspy.core.stream.Stream.plot "

4

1 に答える 1

1

のメソッドは、実際にはグリッドを自動的に生成します。たとえば、デフォルトの例を次のようにプロットするとStream、plot()

from obspy.core import read
st = read()  # without filename an example file is loaded
tr = st[0]   # we will use only the first channel
tr.plot()

ここに画像の説明を入力

http://docs.obspy.org/packages/autogen/obspy.core.stream.Stream.plot.htmlで指摘されているようにnumber_of_ticks、tick_formatおよびパラメータで遊んでみるとよいでしょう。tick_rotation

ただし、より詳細な制御が必要な場合は、matplotlib図を入力パラメーターとしてplot()メソッドに渡すことができます。

from obspy.core import read
import matplotlib.pyplot as plt

fig = plt.figure()

st = read('/path/to/file.sac')
st.plot(fig=fig)

# at this point do whatever you want with your figure, e.g.
fig.gca().set_axis_off()

# finally display your figure
fig.show()

それが役に立てば幸い。

于 2013-09-02T00:07:54.390 に答える