A collection of scripts to steer HypoDD by Felix Waldhauser.
Clone or download
krischer Merge pull request #8 from krischer/revert-6-bump_version_for_obspy_1.1
Revert "bump version number for obspy 1.1.0 compatibility"
Latest commit ed733e9 Oct 3, 2017

README.md

hypoDDpy

This a collection of tools to run HypoDD by Felix Waldhauser.

It takes event files in the QuakeML format, station data in the SEED format and waveform data in any format ObsPy can read and does all the rest.

The output is one QuakeML file with the relocated events having one additional Origin node. The events that could not be relocated will not be changed.

Flowchart 1

Did It Help You in Your Research?

If yes, consider citing this software! Just click on this button here to get all the required information:

DOI

Requirements

  • Python 2.6 or 2.7
  • NumPy
  • matplotlib
  • progressbar
  • ObsPy (Tested on 1.0.2)

Installation

hypoDDpy currently works with HypoDD 2.1b which you will have to acquire separately from Felix Waldhauser.

Put the archive here:

hypoddpy/src/HYPODD_2.1b.tar.gz

The src directory will likely not exist.

Then run either of the following two commands, depending on which Python module installer you prefer:

pip install -v -e .
python setup.py develop

The in-place install is a good idea because there is a chance that you will have to adjust the source code.

Running it

It is steered via a Python script that you will have to create. It should be rather self-explanatory.

After you created it, simply run it to perform the relocation.

import glob
from hypoddpy import HypoDDRelocator


# Init the relocator with the working directory and some necessary
# configuration values.
#
# The working dir is where all the working files and some output files will be
# stored.
# All the other attributes are related to the cross correlation and should be
# self-explanatory.
relocator = HypoDDRelocator(working_dir="relocator_working_dir",
    cc_time_before=0.05,
    cc_time_after=0.2,
    cc_maxlag=0.1,
    cc_filter_min_freq=1.0,
    cc_filter_max_freq=20.0,
    cc_p_phase_weighting={"Z": 1.0},
    cc_s_phase_weighting={"Z": 1.0, "E": 1.0, "N": 1.0},
    cc_min_allowed_cross_corr_coeff=0.4)

# Add the necessary files. Call a function multiple times if necessary.
relocator.add_event_files(glob.glob("events/*.xml"))
relocator.add_waveform_files(glob.glob("waveform/*.mseed"))
relocator.add_station_files(glob.glob("station/*.xml"))

# Setup the velocity model. This is just a constant velocity model.
relocator.setup_velocity_model(
    model_type="layered_p_velocity_with_constant_vp_vs_ratio",
    layer_tops=[(-10000, 5.8)],
    vp_vs_ratio=1.73)

# Start the relocation with the desired output file.
relocator.start_relocation(output_event_file="relocated_events.xml")