AgentSkillsCN

obspy

利用ObsPy进行地震学数据处理。帮助用户读取地震波形，对时间序列进行滤波与处理，从FDSN服务获取数据，并开展地震分析。当Claude需要执行以下操作时使用：(1) 读取地震数据格式（MiniSEED、SAC、GSE2、SEGY）；(2) 对波形进行滤波或处理；(3) 从IRIS/USGS/FDSN服务获取数据；(4) 按震级或位置搜索地震事件；(5) 绘制地震记录或频谱图；(6) 去除仪器响应；(7) 分析台站元数据。

SKILL.md
--- frontmatter
name: obspy
description: |
  Seismology data processing with ObsPy. Helps with reading seismic waveforms,
  filtering/processing time series, fetching data from FDSN services, and
  earthquake analysis. Use when Claude needs to: (1) Read seismic data formats
  (MiniSEED, SAC, GSE2, SEGY), (2) Filter or process waveforms, (3) Fetch data
  from IRIS/USGS/FDSN services, (4) Search for earthquakes by magnitude/location,
  (5) Plot seismograms or spectrograms, (6) Remove instrument response,
  (7) Analyze station metadata.

ObsPy - Seismology Data Processing

Quick Reference

python
from obspy import read, UTCDateTime
from obspy.clients.fdsn import Client

# Read local file (MiniSEED, SAC, etc.)
st = read("data.mseed")
tr = st[0]                          # First trace
print(tr.stats)                     # Metadata

# Fetch from FDSN
client = Client("IRIS")
t = UTCDateTime("2023-02-06T01:17:00")
st = client.get_waveforms("IU", "ANMO", "00", "LHZ", t, t + 3600)
st.plot()

Key Classes

ClassPurpose
StreamContainer for multiple Trace objects
TraceSingle waveform with data + metadata
UTCDateTimePrecise time handling
InventoryStation/channel metadata
CatalogEarthquake event information

Essential Operations

Read and Inspect

python
st = read("data.mseed")              # Auto-detect format
tr = st[0]
print(tr.stats.station, tr.stats.channel, tr.stats.sampling_rate)

Filter and Process

python
st.detrend("demean")                 # Remove mean
st.detrend("linear")                 # Remove trend
st.taper(max_percentage=0.05)        # Taper edges
st.filter("bandpass", freqmin=0.1, freqmax=10.0)

Fetch Waveforms from FDSN

python
client = Client("IRIS")
t1 = UTCDateTime("2023-02-06T01:17:00")
st = client.get_waveforms("IU", "ANMO", "00", "LHZ", t1, t1 + 3600)
st.write("output.mseed", format="MSEED")

Search Earthquakes

python
client = Client("USGS")
cat = client.get_events(
    starttime=UTCDateTime("2023-01-01"),
    endtime=UTCDateTime("2023-12-31"),
    minmagnitude=7.0
)
event = cat[0]
print(f"M{event.magnitudes[0].mag} at {event.origins[0].latitude}")

Get Station Metadata

python
inv = client.get_stations(
    network="IU", station="ANMO",
    level="response"                 # Required for response removal
)

Remove Instrument Response

python
st = client.get_waveforms("IU", "ANMO", "00", "LHZ", t, t + 3600)
inv = client.get_stations(network="IU", station="ANMO", level="response")
st.remove_response(inventory=inv, output="VEL")  # VEL, DISP, or ACC

Trim and Select

python
st.trim(UTCDateTime("2023-01-01"), UTCDateTime("2023-01-01T01:00:00"))
st_z = st.select(channel="*Z")       # Vertical only
st_bh = st.select(channel="BH*")     # BH channels

Merge and Handle Gaps

python
st.print_gaps()                      # Check for gaps
st.merge(method=1, fill_value="interpolate")

Writing Data

python
st.write("output.mseed", format="MSEED")
st.write("output.sac", format="SAC")

Error Handling

python
from obspy.clients.fdsn.header import FDSNNoDataException

try:
    st = client.get_waveforms("IU", "ANMO", "00", "LHZ", t1, t2)
except FDSNNoDataException:
    print("No data available")

Common Tips

  1. •Always detrend and taper before filtering to avoid artifacts
  2. •Use level="response" when fetching stations for instrument correction
  3. •Check for gaps before processing with st.print_gaps()
  4. •Wildcards work in queries: station="A*", channel="BH?"
  5. •UTCDateTime accepts ISO strings, timestamps, datetime objects

References

Scripts