med_spec
Calculates the median spectra of waveforms from a seismic deployment
Science Score: 77.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
✓CITATION.cff file
Found CITATION.cff file -
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
✓DOI references
Found 5 DOI reference(s) in README -
✓Academic publication links
Links to: zenodo.org -
✓Committers with academic emails
2 of 4 committers (50.0%) from academic institutions -
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (11.7%) to scientific vocabulary
Repository
Calculates the median spectra of waveforms from a seismic deployment
Basic Info
- Host: GitHub
- Owner: tbartholomaus
- Language: Python
- Default Branch: master
- Size: 25.6 MB
Statistics
- Stars: 1
- Watchers: 0
- Forks: 2
- Open Issues: 0
- Releases: 1
Metadata Files
README.md
med_spec
Calculates the median spectrograms of waveforms from a seismic deployment. Method follows that described in Bartholomaus et al., 2015, GRL.
Please cite this software as:
Bartholomaus, T., & Terleth, Y. (2023). med_spec: Calculation of median seismic spectrograms, for purposes of quantifying tremor (Version 1.0.0) [Computer software]. https://doi.org/10.5281/zenodo.8102681
This repository consists of several key functions and a couple of helper functions to calculate the median power at a range of frequencies for seismic timeseries at single stations. These median spectra are less sensitive to brief, high power transient events such as earthquakes and icequakes, and thus are appropriate for use in noisy seismic environments, like glaciers.
The core of these scripts begins with medspecloop.py. This file loops over a long duration seismic deployment and calls getmedspectra_v1.py every time there's a snippet of data to calculate the median spectra of. The output of getmedspectrav1.py, which contains the ingredients of a median spectrogram, includes time, frequency and power. As of March 1, 2018, getmedspectrav1.py satisfies Parseval's theorem. Previous versions of data run with this function reported power that was too low due to windowing that reduced amplitudes by an unaccounted-for scalar. The file medspecloop_v3.py corrects the data using the instrument response and requires use of an appropriate RESP or station xml file.
As of version 3 of medspecloop (15 November 2018), neither medspecloop.py nor getmedspectrav1.py should require modification to run in a variety of configurations and with a variety of datasets. Key parameters and paths are defined in a parameter file, which includes all of these. An example parameter file is included in this repository as **medspec.par**.
At the command prompt, the successful execution of the medspec algorithm looks like:
```
python ./medspecloop.py medspec.par
```
A more feature-rich execution of this code could be
nohup python -u ./med_spec_loop.py med_spec.par LM BBWL > BBWL.log &
which is expected to run on a server with nohup preventing killing the process if the server is disconnected from. The -u flag forces the text output of the python script to stdout, so that it can then be redirected to a log file (> BBWL.log). Here network and station codes from the par file are also overridden by LM and BBWL. If present, these network, station, and channel codes (in that order) take precedence over whatever is in the par file.
Ancillary scripts include the following: plotmedspec.py: Read in the saved output from medspecloop.py and create spectrograms of the processed data.
The next file is loadallGHT.py: This script reads in a median spectrogram and then plots a time series of GHT amplitude by integrating the PSD over a specified range of frequencies.
Owner
- Name: Timothy Bartholomaus
- Login: tbartholomaus
- Kind: user
- Location: Moscow, ID
- Company: University of Idaho
- Website: www.tbartholomaus.org
- Twitter: timbartholomaus
- Repositories: 4
- Profile: https://github.com/tbartholomaus
Glaciologist and Associate Professor studying glacier and ice sheet change.
Citation (CITATION.cff)
cff-version: 1.2.0 message: "If you use this software, please cite it as below." authors: - family-names: "Bartholomaus" given-names: "Timothy" orcid: "https://orcid.org/0000-0002-1470-6720" - family-names: "Terleth" given-names: "Yoram" orcid: "https://orcid.org/0000-0002-6583-5258" title: "med_spec: Calculation of median seismic spectrograms, for purposes of quantifying tremor" version: 1.0.0 doi: 10.5281/zenodo.8102681 date-released: 2023-06-30 url: "https://github.com/tbartholomaus/med_spec"
GitHub Events
Total
Last Year
Committers
Last synced: 8 months ago
Top Committers
| Name | Commits | |
|---|---|---|
| Tim Bartholomaus | 2****s | 75 |
| Yoram Terleth | 8****h | 7 |
| Timothy Bartholomaus | t****s@k****u | 4 |
| Timothy Bartholomaus | t****s@c****u | 3 |
Committer Domains (Top 20 + Academic)
Issues and Pull Requests
Last synced: 8 months ago
All Time
- Total issues: 0
- Total pull requests: 7
- Average time to close issues: N/A
- Average time to close pull requests: about 2 months
- Total issue authors: 0
- Total pull request authors: 2
- Average comments per issue: 0
- Average comments per pull request: 0.14
- Merged pull requests: 7
- Bot issues: 0
- Bot pull requests: 0
Past Year
- Issues: 0
- Pull requests: 0
- Average time to close issues: N/A
- Average time to close pull requests: N/A
- Issue authors: 0
- Pull request authors: 0
- Average comments per issue: 0
- Average comments per pull request: 0
- Merged pull requests: 0
- Bot issues: 0
- Bot pull requests: 0
Top Authors
Issue Authors
Pull Request Authors
- tbartholomaus (5)
- yoramterleth (4)