An automated P-phase Arrival Time Picker with SNR output

PPHASEPICKER is a powerful tool for automatically picking P-phase onsets.
1.5K 次下载
更新时间 2019/4/3

查看许可证

PPHASEPICKER is a powerful tool for automatically picking P-phase onsets with high precision without requiring detection interval or threshold settings. The algorithm detects P-phase onset in single-component acceleration or broadband velocity records using the histogram method. It also computes signal-to-noise ratio (SNR).

PPHASEPICKER has been integrated into the “Automated Processing and Review Interface for Strong Motion Data (PRISM)” software (https://earthquake.usgs.gov/research/software/prism/) of the U.S. Geological Survey in order to identify pre-event time-window for systematic and automated processing of large numbers of strong-motion data.

An example MatLAB code is provided in zip file to show how to run PPHASEPICKER using three sample waveforms, one from strong earthquake and others from micro seismic events.

For your intended application, please note that filter corner frequencies may need to be updated manually within the code based on the frequency content of the input record.

Reference:

Kalkan, E. (2016). “An Automatic P-phase Arrival Time Picker“, Bull. of Seis. Soc. of Am., 106(3): 971-986, doi: 10.1785/0120150111.

引用格式

Dr. Erol Kalkan, P.E. (2024). An automated P-phase Arrival Time Picker with SNR output (https://www.mathworks.com/matlabcentral/fileexchange/57729-an-automated-p-phase-arrival-time-picker-with-snr-output), MATLAB Central File Exchange. 检索时间: .

MATLAB 版本兼容性
创建方式 R2014a
兼容任何版本
平台兼容性
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
版本 已发布 发行说明
6.0.2

Description updated

6.0.1

Updated notes

6.0.0.0

Added signal-to-noise ratio computation
Updated demo example file
Updated the reference
Updated description

5.0.0.0

Code is modified to prevent numerical errors due to very low amplitude input signals

4.0.0.0

Added BSSA paper describing the picker
Fixed "nbins" for arbitrary sampling rates; tested for inhomogeneous (arbitrarily sampled) weak acceleration records from China.

3.0.0.0

Modified for low sampled input waveforms (lower than 100 samples-per-second, sps) and third example (waveform with 80 sps) is included in demo file.

2.0.0.0

Tested for MatLAB 2013a, 2014a and 2016a.

1.0.0.0

Text is revised
Typo has been fixed.
A new micro seismic earthquake example is included.
Histogram computation is added in order to eliminate need for digital signal processing toolbox.
"no bandpass filter" option is added.