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

Popular posts from this blog

Why does Ruby on Rails generate add a blank line to the end of a file? -

keyboard - Smiles and long press feature in Android -

node.js - Bad Request - node js ajax post -