OBSPY 中文教程》
[Pages:96]OBSPY
V 1.0 2020 / 04 / 12
Obspy python , python Obspy / Obspy Python obsp y
obspy
penghuiw@
ObsPy
1. UTCDateTime......................................................................................6 1) :..............................................................................................................................6 2) ...........................................................................................................................6 3) .......................................................................................................................6
2. Reading Seismograms ..............................................................................7 1) .......................................................................................................................7 2) ...................................................................................................................8 3) ...........................................................................................................................8
3. Waveform Plotting Tutorial .............................................................................8 1) ...........................................................................................................................9 2) .......................................................................................................................9 3) ...........................................................................................................................9 4) ...............................................................................................................9 5) One-Day ...........................................................................................................9 6) .................................................................................................................10 7) .............................................................................................................10 8) Matplotlib ..........................................................................................10
4. Retrieving Data from Data Centers .........................................10 1) FDSN ...............................................................................................................11 2) ArcLink.............................................................................................................................11 3) IRIS ..................................................................................................................12 4) Earthworm Wave ...................................................................................................12 5) NERIES ............................................................................................................12 6) NEIC.................................................................................................................................12 7) Syngine ...................................................................................................................12
5. Filtering Seismograms ............................................................................................13 6. Downsampling Seismograms .............................................................................13 7. Merging Seismograms............................................................................................14 8. Beamforming - FK Analysis FK ............................................................................15 9. Seismogram Envelopes ..................................................................................22 10. Plotting Spectrograms................................................................................24 11. Trigger/Picker Tutorial/...................................................................25
1) .................................................................................................................25 2) .........................................................................................................................26 3) .................................................................................................................27
a) STA/LTA ...........................................................................................................27 b) Z ......................................................................................................................28 c) STA/LTA ...........................................................................................................28 d) -STA-.........................................................................................................29 e) STA/LTA ...........................................................................................................30
4) .........................................................................................................31
5) .................................................................33
6) .....................................................................................................................35
a) Baer Picker...............................................................................................................35
b) AR Picker..................................................................................................................35
7) .........................................................................................................................36
12. Poles and Zeros, Frequency Response ........................................37
13. Seismometer Correction/Simulation..............................................39
1) .....................................................................................................39
a) StationXML Inventory ...............................................39
b) RESP .......................................................................................................40
c) / SEED XMLSEED ................................................42
2) PAZ .................................................................................................................43
14. Clone an Existing Dataless SEED File
SEED ............45
15. Export Seismograms to MATLAB MATLAB ...............................................46
16. Export Seismograms to ASCII ASCII .................................................47
1) Built-in ...........................................................................................................47
2) .....................................................................................................................48
17. Anything to MiniSEED MiniSEED ...............................................48
18. Beachball Plot.........................................................................................49
19. Basemap Plots........................................................................................................................51
1) Basemap plot..................................................................................51
2) Beachball Basemap plot ................................................................52
3) beachball Basemap ...................................................................................55
20. Interfacing R from Python python R ...........................................................56
21. Coordinate Conversions .................................................................................57
22. Hierarchical Clustering ...................................................................................57
23. Visualizing Probabilistic Power Spectral Densities ..........59
24. Array Response Function ......................................................................64
25. Continuous Wavelet Transform.............................................................65
1) obspy ......................................................................................................................65
2) MLPY.......................................................................................................................66
26. Time Frequency Misfit ..................................................................................68
1) .............................................................................................................68
2) .................................................................................................................69
3) .................................................................................................................71
4) .....................................................................................................................73
5) .....................................................................................................................74
27. Visualize Data Availability of Local Waveform Archive
..................................................................................................................................... 76
28. Travel Time and Ray Path Plotting..........................................76
1) .........................................................................................................................76
2) .............................................................................................................77
3) .................................................................................................................78
4) .............................................................................................................79
29. Cross Correlation Pick Correction ................................................82
30. Handling custom defined tags in QuakeML and the ObsPy Catalog/Event framework
QuakeML ObsPy /........................................85
31. Handling custom defined tags in StationXML with the Obspy Inventory Obspy I
nventory StationXML .....................................................................88
32. Creating a StationXML file from Scratch
Scratch StationXML ....92
33. Connecting to a SeedLink Server SeedLink ........................................94
1) create_client ...........................................................................................................95
a) INFO .......................................................................................95
b) .............................................................................................95
2) .............................................................................................96
1. UTCDateTime
ObsPy UTCDateTime POSIX Python datetime
1) :
>>> from obspy.core import UTCDateTime >>> UTCDateTime("2012-09-07T12:15:00") UTCDateTime(2012, 9, 7, 12, 15) >>> UTCDateTime(2012, 9, 7, 12, 15, 0) UTCDateTime(2012, 9, 7, 12, 15) >>> UTCDateTime(1347020100.0) UTCDateTime(2012, 9, 7, 12, 15)
>>> UTCDateTime("2012-09-07T12:15:00+02:00") UTCDateTime(2012, 9, 7, 10, 15)
2)
>>> time = UTCDateTime("2012-09-07T12:15:00") >>> time.year 2012 >>> time.julday 251 >>> time.timestamp 1347020100.0 >>> time.weekday 4
3)
>>> time = UTCDateTime("2012-09-07T12:15:00") >>> print(time + 3600) 2012-09-07T13:15:00.000000Z >>> time2 = UTCDateTime(2012, 1, 1) >>> print(time - time2) 21644100.0
2. Reading Seismograms
obspy SACMiniSEEDGSE2SEISANQ read Stream Stream Trace /
Trace data NumPy ndarray stats Stats Stats starttime endtime UTCDateTime
GSE2 ObsPy Stream Trace
>>> from obspy import read >>> st = read('') >>> print(st) 1 Trace(s) in Stream: .RJOB..Z | 2005-10-06T07:21:59.849998Z - 2005-1006T07:24:59.844998Z | 200.0 Hz, 36000 samples >>> len(st) 1 >>> tr = st[0] # assign first and only trace to new variable >>> print(tr) .RJOB..Z | 2005-10-06T07:21:59.849998Z - 2005-1006T07:24:59.844998Z | 200.0 Hz, 36000 samples
1)
Trace stats
>>> print(tr.stats) network: station: RJOB location: channel: Z starttime: 2005-10-06T07:21:59.849998Z endtime: 2005-10-06T07:24:59.844998Z sampling_rate: 200.0 delta: 0.005 npts: 36000 calib: 0.0948999971151 _format: GSE2 gse2: AttribDict({'instype': ' ', 'datatype': 'CM6', 'hang': -
1.0, 'auxid' >>> tr.stats.station 'RJOB' >>> tr.stats.gse2.datatype 'CM6'
2)
Trace data
>>> tr.data array([-38, 12, -4, ..., -14, -3, -9]) >>> tr.data[0:3] array([-38, 12, -4]) >>> len(tr) 36000
3)
Stream plot obspy.imaging
>>> st.plot()
3. Waveform Plotting Tutorial
obspy Stream Trace singlechannel threechannel
>>> from obspy.core import read >>> singlechannel = read('') >>> print(singlechannel) 1 Trace(s) in Stream: DK.COP..BHZ | 2009-02-19T00:00:00.025100Z - 2009-0219T23:59:59.975100Z | 20.0 Hz, 1728000 samples >>> threechannels = read('') >>> threechannels += read('') >>> threechannels += read('') >>> print(threechannels) 3 Trace(s) in Stream:
................
................
In order to avoid copyright disputes, this page is only a partial summary.
To fulfill the demand for quickly locating and searching documents.
It is intelligent file search solution for home and business.