obspy

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.
version	1.0.0
author	Geoscience Skills
license	MIT
tags	["Seismology","Waveforms","FDSN","Earthquake","Time Series","ObsPy","MiniSEED","Signal Processing"]
dependencies	["obspy>=1.4.0"]
complements	["segyio","disba"]
workflow_role	processing

ObsPy - Seismology Data Processing

Quick Reference

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

Class	Purpose
Stream	Container for multiple Trace objects
Trace	Single waveform with data + metadata
UTCDateTime	Precise time handling
Inventory	Station/channel metadata
Catalog	Earthquake event information

Essential Operations

Read and Inspect

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

Filter and Process

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