You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
The input PSD file should have the same format as:
49
-
- output of the Noise Toolkit ntk_extractPsdHour.py script (see NTK PSD/PDF bundle, https://seiscode.iris.washington.edu/projects/noise-toolkit-pdf-psd-bundle)
49
+
- output of the Noise Toolkit ntk_extractPsdHour.py script (see NTK PSD/PDF bundle, https://github.com/iris-edu/noise-toolkit/wiki/Noise-Toolkit-PDF-PSD-bundle)
50
50
or
51
51
- PQLX's "exPSDhour script" (see http://pubs.usgs.gov/of/2010/1292/pdf/OF10-1292.pdf)
52
52
53
53
NOTES
54
54
55
55
- for more information:
56
56
. visit IRIS DMC Noise Toolkit Data Product web page at: http://ds.iris.edu/ds/products/noise-toolkit/
57
-
. visit product's Wiki page (https://seiscode.iris.washington.edu/projects/noise-toolkit-microseism-energy-me-bundle/wiki)
57
+
. visit product's Wiki page (https://github.com/iris-edu/noise-toolkit/wiki/Noise-Toolkit-Microseism-Energy-(ME)-bundle)
@@ -20,33 +20,33 @@ This repository is the new home for the Noise Toolkit on seiscode (https://seisc
20
20
**DESCRIPTION:**
21
21
22
22
The IRIS DMC Noise Toolkit data product is a collection of three open-source bundles that provide relatively simple
23
-
Python/ObsPy code bundles to:
23
+
Python 2.7 /ObsPy code bundles to:
24
24
25
25
- compute power spectral density (PSD) of station waveform data using customized parameters (Noise Toolkit PDF/PSD bundle)
26
26
- compute microseism energy from PSDs at different frequency bands (Noise Toolkit Microseism Energy bundle)
27
27
- perform frequency dependent polarization analysis of the waveform data (Noise Toolkit Polarization Analysis bundle)
28
28
29
29
**PDF/PSD bundle**
30
30
31
-
The PDF/PSD bundle provides three highly configurable Python scripts to calculate waveform spectra in Python/ObsPyP. Although this package takes advantage of FDSN Web service client for ObsPy to retrieve necessary waveform data, but it also allows users to process waveform data from their local files. This package provides PSD file collections similar to popular PQLX package (McNamara and Boaz, 2005) and therefore it is compatible with existing user programs (see the Wiki page). The scripts included in this bundle are:
31
+
The PDF/PSD bundle provides three highly configurable Python 2.7 scripts to calculate waveform spectra in Python/ObsPyP. Although this package takes advantage of FDSN Web service client for ObsPy to retrieve necessary waveform data, but it also allows users to process waveform data from their local files. This package provides PSD file collections similar to popular PQLX package (McNamara and Boaz, 2005) and therefore it is compatible with existing user programs (see the Wiki page). The scripts included in this bundle are:
32
32
33
-
- ntk_computePSD.py - a Python script to request waveforms and response data for given station(s) using the ObsPy FDSN client OR to read user's waveform data files (in SAC, MSEED, CSS, etc. format), compute PSDs and populate a file-based PSD database
34
-
- ntk_extractPsdHour.py - a Python script to extract PSDs for a given channel and bounding parameters from the PSD database. The output is similar to PQLX's exPSDhour script
35
-
- ntk_binPsdDay.py - a Python script to bin PSD's to daily files for a given channel and bounding parameters.
33
+
- ntk_computePSD.py - a Python 2.7 script to request waveforms and response data for given station(s) using the ObsPy FDSN client OR to read user's waveform data files (in SAC, MSEED, CSS, etc. format), compute PSDs and populate a file-based PSD database
34
+
- ntk_extractPsdHour.py - a Python 2.7 script to extract PSDs for a given channel and bounding parameters from the PSD database. The output is similar to PQLX's exPSDhour script
35
+
- ntk_binPsdDay.py - a Python 2.7 script to bin PSD's to daily files for a given channel and bounding parameters.
36
36
37
37
**Microseism Energy (ME) bundle**
38
38
39
-
Microseism Energy (ME) bundle is a collection of three Python scripts that are configurable and allow users to conveniently calculate and plot microseism energy temporal variations in the period band of interest using the available PSD values. By default, the bundle is configured to to calculate ME over 1-5 s band, targeting smaller local storms for coastal stations, 5-10 s for the secondary microseisms, 11-30 s for the primary microseisms and the 50-200 s band for the Earth hum using PSDs of three-component broadband seismic data (BH channels). The stored microseism energy values will be smoothed using a median sliding time window (e.g. 6 hours, 12 hours, 1 day, 4 days, and 16 days, etc.).
39
+
Microseism Energy (ME) bundle is a collection of three Python 2.7 scripts that are configurable and allow users to conveniently calculate and plot microseism energy temporal variations in the period band of interest using the available PSD values. By default, the bundle is configured to to calculate ME over 1-5 s band, targeting smaller local storms for coastal stations, 5-10 s for the secondary microseisms, 11-30 s for the primary microseisms and the 50-200 s band for the Earth hum using PSDs of three-component broadband seismic data (BH channels). The stored microseism energy values will be smoothed using a median sliding time window (e.g. 6 hours, 12 hours, 1 day, 4 days, and 16 days, etc.).
40
40
41
41
The scripts included in this bundle are:
42
42
43
-
- ntk_computePower.py – a Python script to calculate power of each PSD window (by default 1 hour) over selected bin period bands
44
-
- ntk_medianPower.py – a Python script to calculates median power from the PSD powers (by default hourly PSD power) using a sliding window of a given length (e.g., 12 hours)
45
-
- ntk_plotPower.py – a Python script to plot median powers computed from a series of PSD powers
43
+
- ntk_computePower.py – a Python 2.7 script to calculate power of each PSD window (by default 1 hour) over selected bin period bands
44
+
- ntk_medianPower.py – a Python 2.7 script to calculates median power from the PSD powers (by default hourly PSD power) using a sliding window of a given length (e.g., 12 hours)
45
+
- ntk_plotPower.py – a Python 2.7 script to plot median powers computed from a series of PSD powers
46
46
47
47
PSDs of seismic station waveform data needed for computation of the microseism energy are available from IRIS DMC via:
48
48
49
-
- Noise Toolkit PDF-PSD bundle — an open-source Python script bundle to compute PSDs
49
+
- Noise Toolkit PDF-PSD bundle — an open-source Python 2.7 script bundle to compute PSDs
50
50
- PSDs computed by IRIS’s MUSTANG noise-psd Web Service (http://service.iris.edu/mustang/noise-psd/1/)
51
51
52
52
**Polarization Attributes (POLAR) bundle**
@@ -60,11 +60,11 @@ The Polarization bundle of the Noise Toolkit is based on the eigen-decomposition
60
60
o phase difference between the vertical and principal horizontal components (Φ VH)
61
61
o phase difference between the horizontal components (Φ HH)
62
62
63
-
This bundle is composed of three Python scripts:
63
+
This bundle is composed of three Python 2.7 scripts:
64
64
65
-
- ntk_computePolarization.py – an ObsPy script to calculate polarization parameters for a given station and time window
66
-
- ntk_extractPolarHour.py – a Python script to extract polarization parameters for the given channels and bounding parameters
67
-
- ntk_binPolarDay.py – a Python script to bin polarization attributes to daily files for a given channel tag and bounding parameters
65
+
- ntk_computePolarization.py – an Python 2.7 script to calculate polarization parameters for a given station and time window
66
+
- ntk_extractPolarHour.py – a Python 2.7 script to extract polarization parameters for the given channels and bounding parameters
67
+
- ntk_binPolarDay.py – a Python 2.7 script to bin polarization attributes to daily files for a given channel tag and bounding parameters
0 commit comments