Python - How to transform counts in to m/s using the obspy module -
i have miniseed file singlechannel trace , assume data in counts (how can check units of trace?). need transform in m/s. checked obspy tutorial , main problem dont know how access poles , zeros , amplification factor miniseed file. also, need calibration file this?
here code:
from obspy.core import * st=read('/users/guilhermew/documents/projecto/dados sismicos 1 dia/2012_130_doc01.mseed') st.plot()
thanks in advance, guilherme
edit: understood how convert data. obspy has different ways achieve this, comes down removing instrument response waveform data. @robert barsch said, needed file instrument response metadata. came following code:
parser=parser("dir/parser/file") tr in stream_aux: stream_id=tr.stats.network+'.'+tr.stats.station+ '..' + tr.stats.channel paz=parser.getpaz(stream_id, tr.stats.starttime) df = tr.stats.sampling_rate tr.data = seissim(tr.data, df, paz_remove=paz)
im using seissim function convert data. problem output dosen't right (but cant seem post image)
this question should asked seismology community , not @ stackoverflow! how write obspy user mailinglist?
update: still feel answer he/she should ask directly @ obspy mailing list. however, in order give proper answer actual question: miniseed data format not contain meta information such poles , zeros or used unit. yes need file such resp, sac paz, dataless seed, full seed etc in order station specific meta data. apply seismometer correction read http://docs.obspy.org/tutorial/code_snippets/seismometer_correction_simulation.html
Comments
Post a Comment