In this example we will create a dataset for the sulphur_peak sequence, which occurred near the idyllic town of Soda Springs Idaho.
In order to proceed you need to first download the cookiecutter and twine python packages. (we assume you have already installed ObsPlus):
pip install cookiecutter
pip install twineRun cookiecutter to fill in the template.
cookiecutter gh:niosh-mining/opsdataFollow the prompts when asked for name, email, etc.
The values entered for the example are:
- author_name: Derrick Chambers
- email: [email protected]
- dataset_name: sulphur_peak,
- project_url: github.com/d-chambers/opsdata_sulphur_peak,
- dataset_description: A 2017 earthquake sequence occurring in southern Idaho.
- version: 0.1.0
The created directory looks as follows:
- opsdata_sulphur_peak
- opsdata_sulphur_peak
- core.py
- __init.py__
- sulphur_peak
- example_data_file.txt
- version.py
- LICENSE
- README.md
- test_data
- requirements.txt
- setup.py
- tests
- conftest.py
- test_sulphur_peak.py
The next thing to do is fill out the README.md file. It contains
information about the dataset, including a summary, authors, references, etc.
Feel free to add additional headings and customize this in any way you wish.
The Sulphur Peak's page is shown below, but you may wish to add a bit more detail:
# Obsplus Dataset: sulphur_peak
A 2017 southern Idaho (USA) earthquake sequence. This dataset was created as an example
for plubishing ObsPlus datasets.
# More info
This dataset was collected by the University of Utah Seismograph Stations (UUSS) as well a
temporary network installed by UUSS and the United States Geological Survey (USGS). More info
can be found at the
[sequence's webpage](https://quake.utah.edu/monitoring-research/sulphur-peak-earthquake).
This dataset only includes events available from USGS with minimum magnitudes >= 4.0 (29 events)
and stations within 0.5 degrees of the approximate cluster.
# References
Koper, K. D., Pankow, K. L., Pechmann, J. C., Hale, J. M., Burlacu, R., Yeck, W. L., ... &
Shearer, P. M. (2018). Afterslip enhanced aftershock activity during the 2017 earthquake
sequence near Sulphur Peak, Idaho. Geophysical Research Letters, 45(11), 5352-5361.
Shokrgozar, A., & Mashal, M. (2019, April). The Mw 5.3 Sulphur Peak Earthquake in Soda
Springs, Idaho: Perspectives from Earthquake Engineering. In Structures Congress (pp. 306-320).
You may also want to ensure the default licence (BSD) is adequate for you data.
You can include any data files you wish in the distribution, but these
should be small in size. Included files should be added to the source directory,
on the same level as core.py which shares a name with the dataset. In this case
it is opsdata_sulphur_peak/opsdata_sulphur_peak/sulphur_peak.
As a rule of thumb, you should only include files you cannot get from a datacenter or other online resource. For example, if you had reprocessed some events it would be appropriate to include the quakeml (if it is small) with the new data.
Next we need to create the code which will download the data. This dataset logic
lives in the core.py file. In the example this is in relative path
opsdata_sulphur_peak/opsdata_sulphur_peak/core.py. This class is well commented
and shows which methods need to be defined. The first few lines look like this:
Paths to included files can then be accessed with source_path / 'file_name.txt'.
"""
ObsPlus instructions for downloading dataset.
"""
from pathlib import Path
from obsplus import DataSet
from opsdata_sulphur_peak import __version__, source_path
class Sulphur_peak(DataSet):
"""
A 2017 earthquake sequence occurring in southern Idaho.
"""
name = "sulphur_peak"
base_path = Path(__file__).parent
version = __version__
# --- functions used to specify how data are downloaded
def download_events(self):
""" download event data and store them in self.event_path """
Path(self.event_path).mkdir(exist_ok=True, parents=True)
def download_waveforms(self):
""" download waveform data and store them in self.waveform_path """
Path(self.waveform_path).mkdir(exist_ok=True, parents=True)
def download_stations(self):
""" download station data and store them in self.station_path """
Path(self.station_path).mkdir(exist_ok=True, parents=True)We need to define how event, station, and waveform data are downloaded separately. If your dataset doesn't have one of these types of data simply don't modify the method.
Adding the download event logic, the code now looks like this:
from pathlib import Path
import obsplus
import obspy
from obspy.clients.fdsn import Client
from opsdata_sulphur_peak import __version__, source_path
class Sulphur_peak(obsplus.DataSet):
"""
A 2017 earthquake sequence occurring in southern Idaho.
"""
name = "sulphur_peak"
base_path = Path(__file__).parent
version = __version__
longitudes = (-111.6, -111.25)
latitudes = (42.5, 42.7)
starttime = obspy.UTCDateTime('2017-01-01')
endtime = obspy.UTCDateTime('2019-01-01')
min_magnitude = 4.0
station_search_degrees = 0.5
# --- functions used to specify how data are downloaded
def download_events(self):
""" download event data and store them in self.event_path """
client = Client("IRIS")
# Define inputs to event query
event_params = dict(
minlongitude=self.longitudes[0],
maxlongitude=self.longitudes[1],
minlatitude=self.latitudes[0],
maxlatitude=self.latitudes[1],
starttime=self.starttime,
endtime=self.endtime,
minmagnitude=self.min_magnitude, # minimum magnitude
)
# Query iris for data and plot
cat = client.get_events(**event_params)
# save event to self.event_path
obsplus.EventBank(self.event_path).put_events(cat)The download_stations method for the example is:
def download_stations(self):
""" download station data and store them in self.station_path """
station_path = Path(self.station_path)
station_path.mkdir(exist_ok=True, parents=True)
client = Client("IRIS")
station_params = dict(
starttime=self.starttime,
endtime=self.endtime,
latitude=sum(self.latitudes) / 2.,
longitude=sum(self.longitudes) / 2.,
maxradius=self.station_search_degrees, # radius in degrees
channel='H*', # only include channels that start with H
level='response', # indicates we want instrument responses
)
inv = client.get_stations(**station_params)
inv.write(str(station_path / 'inventory.xml'), 'stationxml')The download_stations method for the example is:
def download_waveforms(self):
""" download waveform data and store them in self.waveform_path """
client = Client("IRIS")
NSLC = ['network', 'station', 'location', 'channel']
# Times before/after origin to include in waveforms
time_before = 5 # time before origin to download
time_after = 95 # time after origin to download
# Since we need to use the inventory and catalog to determine which
# channels/times to get data for, make sure those are downloaded first.
if self.events_need_downloading:
self.download_events()
if self.stations_need_downloading:
self.download_stations()
# Load inventory and catalog into dataframes.
inv_df = obsplus.stations_to_df(self.station_path)
cat_df = obsplus.events_to_df(self.event_path)
# Create bulk requests for each event and put into a wave bank.
wbank = obsplus.WaveBank(self.waveform_path)
for _, event_series in cat_df.iterrows():
event_df = inv_df[NSLC]
time = obsplus.utils.to_utc(event_series['time'])
event_df['starttime'] = time - time_before
event_df['endtime'] = time + time_after
bulk = obsplus.utils.pd.get_waveforms_bulk_args(event_df)
st = client.get_waveforms_bulk(bulk)
wbank.put_waveforms(st)In order to test the nascent dataset we need to first install it in editable (development) mode.
Now we are ready to run some simple tests. This is done using pytest like so:
pytest . -x --pdbThe "x" flag tells pytest to stop if any tests fail and the 'pdb' flag will automatically start a debugging session on the failed test (or unhandled exception).
The first time the tests run the data will be downloaded and a json file with the file checksums will be created in the source directory.
By default the tests will only ensure the dataset was downloaded, meaning the code added in step 5 ran without errors. You are certainly welcome (and encouraged) to add data-specific tests when it makes sense.
If the tests pass you are ready to publish the package!
We recommend you publish your dataset to both github (or some other comparable service) and pypi. Github allows users to file issues and ask questions, while pypi makes it simple for datasets to be downloaded and discovered.
The Sulphur Peak dataset was published to https://github.com/d-chambers/opsdata_sulphur_peak. After creating a new repo I did the following:
git init
git add remote origin https://github.com/d-chambers/opsdata_sulphur_peak
git push origin masterIn order to upload to pypi you need to first create a pypi account. Then, from the dataset directory, you can simply run the following commands:
python setup.py sdist bdist_wheel
twine upload dist/*Enter your pypi username and password as prompted.
Congragulations, you created and published an ObsPlus dataset!