Examples

The Pyrocko framework covers many aspects of seismological data handling and management. This section offers tutorials and real-world examples on how to utilize the framework.

Jupyter Notebooks

Example notebooks are available in the Pyrocko Notebooks Repository - Contributions to the collection are highly welcome!

Scripts collection

Annotated scripts and real-world use cases using Pyrocko.

• Seismic traces
  • Load, filter and save
  • Visual inspection of traces
  • Creating a trace object from scratch
  • Extracting part of a trace (trimming)
  • Time shifting a trace
  • Resampling a trace
  • Convert SAC to MiniSEED
  • Convert MiniSEED to ASCII
  • Finding the comparative misfits of mulitple traces
  • Restitute to displacement using poles and zeros
  • Restitute to displacement using SEED RESP response
• Squirrel-powered data access
  • Tutorial: writing a Squirrel based tool to calculate hourly RMS values
• Dataset management - The pile
  • Reorganizing a dataset into hour-files
  • Downsampling a whole dataset
  • Converting a dataset from Mini-SEED to SAC format
• ObsPy compatibility
  • Using the Snuffler on obspy.Stream
• Metadata read & write
  • StationXML import
  • Pyrocko stations to StationXML
  • Create a StationXML file with flat displacement responses
• GPS (GNSS) data handling
  • Loading a GNSS campaign from CSV files
  • Loading and mapping of GNSS campaign data from UNR
• Downloading seismic data (FDSN)
  • Seismic data from GEOFON
  • StationXML data manipulation
• Earthquake catalog
  • QuakeML import
  • Creating QuakeML from scratch
  • Searching the GlobalCMT catalog
  • Search for an event in GEOFON catalog
• Moment tensor conversions
  • Convert moment tensor components to strike, dip and rake
  • Strike, dip and rake to moment tensor
• Plotting functions
  • Generating topographic maps with automap
  • Plotting beachballs (focal mechanisms)
  • Hudson’s source type plot
  • Source radiation plot
  • Pseudo dynamic rupture - slip map, slip movie, source plots
• GmtPy - A Python interface to GMT 4/5
  • Introduction
  • Tutorial
  • Random remarks
• Snuffler markers
  • Read markers from file exported by Snuffler
• Traveltime calculation and raytracing
  • Calculate traveltimes in layered media
  • Calculate traveltimes in heterogeneous media
  • Traveltime table interpolation
• Forward modeling synthetic seismograms and displacements
  • Calculate synthetic seismograms from a local GF store
  • Calculate synthetic seismograms using the Pseudo Dynamic Rupture
  • Calculate Coulomb Failure Stress (CFS) changes using the Pseudo Dynamic Rupture
  • Calculate spatial surface displacement from a local GF store
  • Calculate spatial surface displacement using subfault dislocations
  • Calculate spatial surface displacement using the Pseudo Dynamic Rupture
  • Calculate spatial surface displacement and export Kite scenes
  • Calculate forward model of thrust faulting and display wrapped phase
  • Combining dislocation sources
  • Modelling viscoelastic static displacement
  • Creating a custom Source Time Function (STF)
• Kinematic/Dynamic source parameter modeling/inversion
  • Calculate subfault dislocations from tractions with Okada half-space equation
• Geodesic functions
  • Distance between points on earth
  • Azimuth and backazimuth
  • Conversion to carthesian coordinates
  • Relative carthesian coordinates to latitude and longitude
• Geographical datasets
  • Topography
  • GSHHG coastal database
  • Tectonic plates and boundaries (PB2003)
  • Global strain rate (GSRM1)
• Accessing crustal velocity databases
  • Crust 2.0 Database
  • Global Crustal Database
• Working with the guts package
  • Introduction to guts
  • guts and pyrocko
  • Implementing guts
  • Usage examples
  • Example guts YAML File: