feature/port-to-py3 #11
120
docs/gui.md
120
docs/gui.md
@ -86,26 +86,25 @@ Press right mouse button and click "View All" from the context menu to reset the
|
||||
|
||||
[//]: <> (Hack: We need these invisible spaces to add space to the first column, otherwise )
|
||||
|
||||
Icon | Description
|
||||
--- | ---
|
||||
<img src="../icons/newfile.png" alt="Create new project" width="64" height="64"> | Create a new project, for more information about projects see [Projects and Events](#projects-and-events).
|
||||
<img src="../icons/openproject.png" alt="Open project" width="64" height="64"> | Load a project file from disk.
|
||||
<img src="../icons/saveproject.png" alt="Save Project" width="64" height="64"> | Save all current events into an associated project file on disk. If there is no project file currently associated, you will be asked to create a new one.
|
||||
<img src="../icons/saveprojectas.png" alt="Save Project as" width="64" height="64"> | Save all current events into a new project file on disk. See [Saving projects](#saving-projects).
|
||||
<img src="../icons/add.png" alt="Add event data" width="64" height="64"> | Add event data by selecting directories containing waveforms. For more information see [Event folder structure](#event-folder-structure).
|
||||
<img src="../icons/openpick.png" alt="Load event information" width="64" height="64"> | Load picks/origins from disk into the currently displayed event. If a pick already exists for a station, the one from file will overwrite the existing one.
|
||||
<img src="../icons/openpicks.png" alt="Load information for all events" width="64" height="64"> | Load picks/origins for all events of the current project. PyLoT searches for files within the directory of the event and tries to load them for that event. For this function to work, the files containing picks/origins have to be named as described in [Event folder structure](#event-folder-structure). If a pick already exists for a station, the one from file will overwrite the existing one.
|
||||
<img src="../icons/savepicks.png" alt="Save picks" width="64" height="64"> | Save event information such as picks and origin to file. You will be asked to select a directory in which this information should be saved.
|
||||
<img src="../icons/openloc.png" alt="Load location information" width="64" height="64"> | Load location information from disk,
|
||||
<img src="../icons/Matlab_PILOT_icon.png" alt="Load legacy information" width="64" height="64"> | Load event information from a previous, MatLab based PILOT version.
|
||||
<img src="../icons/key_Z.png" alt="Display Z" width="64" height="64"> | Display Z component of streams in waveform plot.
|
||||
<img src="../icons/key_N.png" alt="Display N" width="64" height="64"> | Display N component of streams in waveform plot.
|
||||
<img src="../icons/key_E.png" alt="Display E" width="64" height="64"> | Display E component of streams in waveform plot.
|
||||
<img src="../icons/tune.png" alt="Tune Autopicker" width="64" height="64"> | Open the [Tune Autopicker window](#tuning).
|
||||
<img src="../icons/autopylot_button.png" alt="" width="64" height="64"> | Opens a window that allows starting the autopicker for all events ([Production run of the AutoPicker](#production-run-of-the-autopicker)).
|
||||
<img src="../icons/compare_button.png" alt="Comparison" width="64" height="64"> | Compare automatic and manual picks, only available if automatic and manual picks for an event exist. See [Comparison between automatic and manual picks](#comparison-between-automatic-and-manual-picks).
|
||||
<img src="../icons/locate_button.png" alt="Locate event" width="64" height="64"> | Run a location routine (NonLinLoc) as configured in the settings on the picks. See [Location determination](#location-determination).
|
||||
|
||||
| Icon | Description |
|
||||
|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| <img src="../icons/newfile.png" alt="Create new project" width="64" height="64"> | Create a new project, for more information about projects see [Projects and Events](#projects-and-events). |
|
||||
| <img src="../icons/openproject.png" alt="Open project" width="64" height="64"> | Load a project file from disk. |
|
||||
| <img src="../icons/saveproject.png" alt="Save Project" width="64" height="64"> | Save all current events into an associated project file on disk. If there is no project file currently associated, you will be asked to create a new one. |
|
||||
| <img src="../icons/saveprojectas.png" alt="Save Project as" width="64" height="64"> | Save all current events into a new project file on disk. See [Saving projects](#saving-projects). |
|
||||
| <img src="../icons/add.png" alt="Add event data" width="64" height="64"> | Add event data by selecting directories containing waveforms. For more information see [Event folder structure](#event-folder-structure). |
|
||||
| <img src="../icons/openpick.png" alt="Load event information" width="64" height="64"> | Load picks/origins from disk into the currently displayed event. If a pick already exists for a station, the one from file will overwrite the existing one. |
|
||||
| <img src="../icons/openpicks.png" alt="Load information for all events" width="64" height="64"> | Load picks/origins for all events of the current project. PyLoT searches for files within the directory of the event and tries to load them for that event. For this function to work, the files containing picks/origins have to be named as described in [Event folder structure](#event-folder-structure). If a pick already exists for a station, the one from file will overwrite the existing one. |
|
||||
| <img src="../icons/savepicks.png" alt="Save picks" width="64" height="64"> | Save event information such as picks and origin to file. You will be asked to select a directory in which this information should be saved. |
|
||||
| <img src="../icons/openloc.png" alt="Load location information" width="64" height="64"> | Load location information from disk, |
|
||||
| <img src="../icons/Matlab_PILOT_icon.png" alt="Load legacy information" width="64" height="64"> | Load event information from a previous, MatLab based PILOT version. |
|
||||
| <img src="../icons/key_Z.png" alt="Display Z" width="64" height="64"> | Display Z component of streams in waveform plot. |
|
||||
| <img src="../icons/key_N.png" alt="Display N" width="64" height="64"> | Display N component of streams in waveform plot. |
|
||||
| <img src="../icons/key_E.png" alt="Display E" width="64" height="64"> | Display E component of streams in waveform plot. |
|
||||
| <img src="../icons/tune.png" alt="Tune Autopicker" width="64" height="64"> | Open the [Tune Autopicker window](#tuning). |
|
||||
| <img src="../icons/autopylot_button.png" alt="" width="64" height="64"> | Opens a window that allows starting the autopicker for all events ([Production run of the AutoPicker](#production-run-of-the-autopicker)). |
|
||||
| <img src="../icons/compare_button.png" alt="Comparison" width="64" height="64"> | Compare automatic and manual picks, only available if automatic and manual picks for an event exist. See [Comparison between automatic and manual picks](#comparison-between-automatic-and-manual-picks). |
|
||||
| <img src="../icons/locate_button.png" alt="Locate event" width="64" height="64"> | Run a location routine (NonLinLoc) as configured in the settings on the picks. See [Location determination](#location-determination). |
|
||||
|
||||
### Array Map
|
||||
|
||||
@ -127,19 +126,19 @@ The eventlist displays event parameters. The displayed parameters are saved in t
|
||||
|
||||
<img src="images/gui/eventlist.png" alt="Eventlist">
|
||||
|
||||
Column | Description
|
||||
--- | ---
|
||||
Event | Full path to the events folder.
|
||||
Time | Time of event.
|
||||
Lat | Latitude in degrees of event location.
|
||||
Lon | Longitude in degrees of event location.
|
||||
Depth | Depth in km of event.
|
||||
Mag | Magnitude of event.
|
||||
[N] MP | Number of manual picks.
|
||||
[N] AP | Number of automatic picks.
|
||||
Tuning Set | Select whether this event is a Tuning event. See [Automatic Picking](#automatic-picking).
|
||||
Test Set | Select whether this event is a Test event. See [Automatic Picking](#automatic-picking).
|
||||
Notes | Free form text field for notes regarding this event. Text will be saved in the notes.txt file in the event folder.
|
||||
| Column | Description |
|
||||
|------------|--------------------------------------------------------------------------------------------------------------------|
|
||||
| Event | Full path to the events folder. |
|
||||
| Time | Time of event. |
|
||||
| Lat | Latitude in degrees of event location. |
|
||||
| Lon | Longitude in degrees of event location. |
|
||||
| Depth | Depth in km of event. |
|
||||
| Mag | Magnitude of event. |
|
||||
| [N] MP | Number of manual picks. |
|
||||
| [N] AP | Number of automatic picks. |
|
||||
| Tuning Set | Select whether this event is a Tuning event. See [Automatic Picking](#automatic-picking). |
|
||||
| Test Set | Select whether this event is a Test event. See [Automatic Picking](#automatic-picking). |
|
||||
| Notes | Free form text field for notes regarding this event. Text will be saved in the notes.txt file in the event folder. |
|
||||
|
||||
## Usage
|
||||
|
||||
@ -169,17 +168,16 @@ One example header and data line is shown below.
|
||||
|
||||
The meaning of the header entries is:
|
||||
|
||||
Header | description
|
||||
--- | ---
|
||||
event | Event id, has to be the same as the folder name in which waveform data for this event is kept.
|
||||
Data | Origin date of the event, format DD/MM/YY or DD/MM/YYYY.
|
||||
Time | Origin time of the event. Format HH:MM:SS.
|
||||
Lat, Long | Origin latitude and longitude in decimal degrees.
|
||||
Region | Flinn-Engdahl region name.
|
||||
Basis Lat, Basis Lon | Latitude and longitude of the basis of the station network in decimal degrees.
|
||||
Distance [km] | Distance from origin coordinates to basis coordinates in km.
|
||||
Distance [rad] | Distance from origin coordinates to basis coordinates in rad.
|
||||
|
||||
| Header | description |
|
||||
|----------------------|------------------------------------------------------------------------------------------------|
|
||||
| event | Event id, has to be the same as the folder name in which waveform data for this event is kept. |
|
||||
| Data | Origin date of the event, format DD/MM/YY or DD/MM/YYYY. |
|
||||
| Time | Origin time of the event. Format HH:MM:SS. |
|
||||
| Lat, Long | Origin latitude and longitude in decimal degrees. |
|
||||
| Region | Flinn-Engdahl region name. |
|
||||
| Basis Lat, Basis Lon | Latitude and longitude of the basis of the station network in decimal degrees. |
|
||||
| Distance [km] | Distance from origin coordinates to basis coordinates in km. |
|
||||
| Distance [rad] | Distance from origin coordinates to basis coordinates in rad. |
|
||||
|
||||
### Adding events to project
|
||||
|
||||
@ -220,24 +218,24 @@ Open the picking window of a station by leftclicking on any trace in the wavefor
|
||||
|
||||
#### Picking Window Settings
|
||||
|
||||
Icon | Shortcut | Menu Alternative | Description
|
||||
---|---|---|---
|
||||
<img src="../icons/filter_p.png" alt="Filter P" width="64" height="64"> | p | Filter->Apply P Filter | Filter all channels according to the options specified in Filter parameter, P Filter section.
|
||||
<img src="../icons/filter_s.png" alt="Filter S" width="64" height="64"> | s | Filter->Apply S Filter | Filter all channels according to the options specified in Filter parameter, S Filter section.
|
||||
<img src="../icons/key_A.png" alt="Filter Automatically" width="64" height="64"> | Ctrl + a | Filter->Automatic Filtering | If enabled, automatically select the correct filter option (P, S) depending on the selected phase to be picked.
|
||||
 | 1 (P) or 5 (S) | Picks->P or S | Select phase to pick. If Automatic Filtering is enabled, this will apply the appropriate filter depending on the phase.
|
||||
 | - | - | Zoom into waveform.
|
||||
 | - | - | Reset zoom to default view.
|
||||
 | - | - | Delete all manual picks on this station.
|
||||
 | - | - | Click this button and then the picked phase to rename it.
|
||||
 | - | - | If checked, after accepting the manual picks for this station with 'OK', the picking window for the next station will be opened. This option is useful for fast manual picking of a complete event.
|
||||
Estimated onsets | - | - | Show the theoretical onsets for this station. Needs metadata and origin information.
|
||||
Compare to channel | - | - | Select a data channel to compare against. The selected channel will be displayed in the picking window behind every channel allowing the analyst to visually compare signal correlation between different channels.
|
||||
Scaling | - | - | Individual means every channel is scaled to its own maximum. If a channel is selected here, all channels will be scaled relatively to this channel.
|
||||
| Icon | Shortcut | Menu Alternative | Description |
|
||||
|----------------------------------------------------------------------------------|----------------|-----------------------------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| <img src="../icons/filter_p.png" alt="Filter P" width="64" height="64"> | p | Filter->Apply P Filter | Filter all channels according to the options specified in Filter parameter, P Filter section. |
|
||||
| <img src="../icons/filter_s.png" alt="Filter S" width="64" height="64"> | s | Filter->Apply S Filter | Filter all channels according to the options specified in Filter parameter, S Filter section. |
|
||||
| <img src="../icons/key_A.png" alt="Filter Automatically" width="64" height="64"> | Ctrl + a | Filter->Automatic Filtering | If enabled, automatically select the correct filter option (P, S) depending on the selected phase to be picked. |
|
||||
|  | 1 (P) or 5 (S) | Picks->P or S | Select phase to pick. If Automatic Filtering is enabled, this will apply the appropriate filter depending on the phase. |
|
||||
|  | - | - | Zoom into waveform. |
|
||||
|  | - | - | Reset zoom to default view. |
|
||||
|  | - | - | Delete all manual picks on this station. |
|
||||
|  | - | - | Click this button and then the picked phase to rename it. |
|
||||
|  | - | - | If checked, after accepting the manual picks for this station with 'OK', the picking window for the next station will be opened. This option is useful for fast manual picking of a complete event. |
|
||||
| Estimated onsets | - | - | Show the theoretical onsets for this station. Needs metadata and origin information. |
|
||||
| Compare to channel | - | - | Select a data channel to compare against. The selected channel will be displayed in the picking window behind every channel allowing the analyst to visually compare signal correlation between different channels. |
|
||||
| Scaling | - | - | Individual means every channel is scaled to its own maximum. If a channel is selected here, all channels will be scaled relatively to this channel. |
|
||||
|
||||
Menu Command | Shortcut | Description
|
||||
---|---|---
|
||||
P Channels and S Channels | - | Select which channels should be treated as P or S channels during picking. When picking a phase, only the corresponding channels will be shown during the precise pick. Normally, the Z channel should be selected for the P phase and the N and E channel for the S phase.
|
||||
| Menu Command | Shortcut | Description |
|
||||
|---------------------------|----------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| P Channels and S Channels | - | Select which channels should be treated as P or S channels during picking. When picking a phase, only the corresponding channels will be shown during the precise pick. Normally, the Z channel should be selected for the P phase and the N and E channel for the S phase. |
|
||||
|
||||
### Filtering
|
||||
|
||||
|
||||
121
docs/tuning.md
121
docs/tuning.md
@ -6,95 +6,94 @@ A description of the parameters used for determining automatic picks.
|
||||
|
||||
Parameters applied to the traces before picking algorithm starts.
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*P Start*, *P Stop* | Define time interval relative to trace start time for CF calculation on vertical trace. Value is relative to theoretical onset time if 'Use TauPy' option is enabled in main settings of 'Tune Autopicker' dialogue.
|
||||
*S Start*, *S Stop* | Define time interval relative to trace start time for CF calculation on horizontal traces. Value is relative to theoretical onset time if 'Use TauPy' option is enabled in main settings of 'Tune Autopicker' dialogue.
|
||||
*Bandpass Z1* | Filter settings for Butterworth bandpass applied to vertical trace for calculation of initial P pick.
|
||||
*Bandpass Z2* | Filter settings for Butterworth bandpass applied to vertical trace for calculation of precise P pick.
|
||||
*Bandpass H1* | Filter settings for Butterworth bandpass applied to horizontal traces for calculation of initial S pick.
|
||||
*Bandpass H2* | Filter settings for Butterworth bandpass applied to horizontal traces for calculation of precise S pick.
|
||||
| Name | Description |
|
||||
|---------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *P Start*, *P Stop* | Define time interval relative to trace start time for CF calculation on vertical trace. Value is relative to theoretical onset time if 'Use TauPy' option is enabled in main settings of 'Tune Autopicker' dialogue. |
|
||||
| *S Start*, *S Stop* | Define time interval relative to trace start time for CF calculation on horizontal traces. Value is relative to theoretical onset time if 'Use TauPy' option is enabled in main settings of 'Tune Autopicker' dialogue. |
|
||||
| *Bandpass Z1* | Filter settings for Butterworth bandpass applied to vertical trace for calculation of initial P pick. |
|
||||
| *Bandpass Z2* | Filter settings for Butterworth bandpass applied to vertical trace for calculation of precise P pick. |
|
||||
| *Bandpass H1* | Filter settings for Butterworth bandpass applied to horizontal traces for calculation of initial S pick. |
|
||||
| *Bandpass H2* | Filter settings for Butterworth bandpass applied to horizontal traces for calculation of precise S pick. |
|
||||
|
||||
## Inital P pick
|
||||
|
||||
Parameters used for determination of initial P pick.
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*tLTA* | Size of gliding LTA window in seconds used for calculation of HOS-CF.
|
||||
*pickwin P* | Size of time window in seconds in which the minimum of the AIC-CF in front of the maximum of the HOS-CF is determined.
|
||||
*AICtsmooth* | Average of samples in this time window will be used for smoothing of the AIC-CF.
|
||||
*checkwinP* | Time in front of the global maximum of the HOS-CF in which to search for a second local extrema.
|
||||
*minfactorP* | Used with *checkwinP*. If a second local maximum is found, it has to be at least as big as the first maximum * *minfactorP*.
|
||||
*tsignal* | Time window in seconds after the initial P pick used for determining signal amplitude.
|
||||
*tnoise* | Time window in seconds in front of initial P pick used for determining noise amplitude.
|
||||
*tsafetey* | Time in seconds between *tsignal* and *tnoise*.
|
||||
*tslope* | Time window in seconds after initial P pick in which the slope of the onset is calculated.
|
||||
| Name | Description |
|
||||
|--------------|------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *tLTA* | Size of gliding LTA window in seconds used for calculation of HOS-CF. |
|
||||
| *pickwin P* | Size of time window in seconds in which the minimum of the AIC-CF in front of the maximum of the HOS-CF is determined. |
|
||||
| *AICtsmooth* | Average of samples in this time window will be used for smoothing of the AIC-CF. |
|
||||
| *checkwinP* | Time in front of the global maximum of the HOS-CF in which to search for a second local extrema. |
|
||||
| *minfactorP* | Used with *checkwinP*. If a second local maximum is found, it has to be at least as big as the first maximum * *minfactorP*. |
|
||||
| *tsignal* | Time window in seconds after the initial P pick used for determining signal amplitude. |
|
||||
| *tnoise* | Time window in seconds in front of initial P pick used for determining noise amplitude. |
|
||||
| *tsafetey* | Time in seconds between *tsignal* and *tnoise*. |
|
||||
| *tslope* | Time window in seconds after initial P pick in which the slope of the onset is calculated. |
|
||||
|
||||
## Inital S pick
|
||||
|
||||
Parameters used for determination of initial S pick
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*tdet1h* | Length of time window in seconds in which AR params of the waveform are determined.
|
||||
*tpred1h* | Length of time window in seconds in which the waveform is predicted using the AR model.
|
||||
*AICtsmoothS* | Average of samples in this time window is used for smoothing the AIC-CF.
|
||||
*pickwinS* | Time window in which the minimum in the AIC-CF in front of the maximum in the ARH-CF is determined.
|
||||
*checkwinS* | Time in front of the global maximum of the ARH-CF in which to search for a second local extrema.
|
||||
*minfactorP* | Used with *checkwinS*. If a second local maximum is found, it has to be at least as big as the first maximum * *minfactorS*.
|
||||
*tsignal* | Time window in seconds after the initial P pick used for determining signal amplitude.
|
||||
*tnoise* | Time window in seconds in front of initial P pick used for determining noise amplitude.
|
||||
*tsafetey* | Time in seconds between *tsignal* and *tnoise*.
|
||||
*tslope* | Time window in seconds after initial P pick in which the slope of the onset is calculated.
|
||||
| Name | Description |
|
||||
|---------------|------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *tdet1h* | Length of time window in seconds in which AR params of the waveform are determined. |
|
||||
| *tpred1h* | Length of time window in seconds in which the waveform is predicted using the AR model. |
|
||||
| *AICtsmoothS* | Average of samples in this time window is used for smoothing the AIC-CF. |
|
||||
| *pickwinS* | Time window in which the minimum in the AIC-CF in front of the maximum in the ARH-CF is determined. |
|
||||
| *checkwinS* | Time in front of the global maximum of the ARH-CF in which to search for a second local extrema. |
|
||||
| *minfactorP* | Used with *checkwinS*. If a second local maximum is found, it has to be at least as big as the first maximum * *minfactorS*. |
|
||||
| *tsignal* | Time window in seconds after the initial P pick used for determining signal amplitude. |
|
||||
| *tnoise* | Time window in seconds in front of initial P pick used for determining noise amplitude. |
|
||||
| *tsafetey* | Time in seconds between *tsignal* and *tnoise*. |
|
||||
| *tslope* | Time window in seconds after initial P pick in which the slope of the onset is calculated. |
|
||||
|
||||
## Precise P pick
|
||||
|
||||
Parameters used for determination of precise P pick.
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*Precalcwin* | Time window in seconds for recalculation of the HOS-CF. The new CF will be two times the size of *Precalcwin*, since it will be calculated from the initial pick to +/- *Precalcwin*.
|
||||
*tsmoothP* | Average of samples in this time window will be used for smoothing the second HOS-CF.
|
||||
*ausP* | Controls artificial uplift of samples during precise picking. A common local minimum of the smoothed and unsmoothed HOS-CF is found when the previous sample is larger or equal to the current sample times (1+*ausP*).
|
||||
| Name | Description |
|
||||
|--------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *Precalcwin* | Time window in seconds for recalculation of the HOS-CF. The new CF will be two times the size of *Precalcwin*, since it will be calculated from the initial pick to +/- *Precalcwin*. |
|
||||
| *tsmoothP* | Average of samples in this time window will be used for smoothing the second HOS-CF. |
|
||||
| *ausP* | Controls artificial uplift of samples during precise picking. A common local minimum of the smoothed and unsmoothed HOS-CF is found when the previous sample is larger or equal to the current sample times (1+*ausP*). |
|
||||
|
||||
## Precise S pick
|
||||
|
||||
Parameters used for determination of precise S pick.
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*tdet2h* | Time window for determination of AR coefficients.
|
||||
*tpred2h* | Time window in which the waveform is predicted using the determined AR parameters.
|
||||
*Srecalcwin* | Time window for recalculation of ARH-CF. New CF will be calculated from initial pick +/- *Srecalcwin*.
|
||||
*tsmoothS* | Average of samples in this time window will be used for smoothing the second ARH-CF.
|
||||
*ausS* | Controls artificial uplift of samples during precise picking. A common local minimum of the smoothed and unsmoothed ARH-CF is found when the previous sample is larger or equal to the current sample times (1+*ausS*).
|
||||
*pickwinS* | Time window around initial pick in which to look for a precise pick.
|
||||
|
||||
| Name | Description |
|
||||
|--------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *tdet2h* | Time window for determination of AR coefficients. |
|
||||
| *tpred2h* | Time window in which the waveform is predicted using the determined AR parameters. |
|
||||
| *Srecalcwin* | Time window for recalculation of ARH-CF. New CF will be calculated from initial pick +/- *Srecalcwin*. |
|
||||
| *tsmoothS* | Average of samples in this time window will be used for smoothing the second ARH-CF. |
|
||||
| *ausS* | Controls artificial uplift of samples during precise picking. A common local minimum of the smoothed and unsmoothed ARH-CF is found when the previous sample is larger or equal to the current sample times (1+*ausS*). |
|
||||
| *pickwinS* | Time window around initial pick in which to look for a precise pick. |
|
||||
|
||||
## Pick quality control
|
||||
|
||||
Parameters used for checking quality and integrity of automatic picks.
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*minAICPslope* | Initial P picks with a slope lower than this value will be discared.
|
||||
*minAICPSNR* | Initial P picks with a SNR below this value will be discarded.
|
||||
*minAICSslope* | Initial S picks with a slope lower than this value will be discarded.
|
||||
*minAICSSNR* | Initial S picks with a SNR below this value will be discarded.
|
||||
*minsiglength*, *noisefacor*. *minpercent* | Parameters for checking signal length. In the time window of size *minsiglength* after the initial P pick *minpercent* of samples have to be larger than the RMS value.
|
||||
*zfac* | To recognize misattributed S picks, the RMS amplitude of vertical and horizontal traces are compared. The RMS amplitude of the vertical traces has to be at least *zfac* higher than the RMS amplitude on the horizontal traces for the pick to be accepted as a valid P pick.
|
||||
*jackfactor* | A P pick is removed if the jackknife pseudo value of the variance of his subgroup is larger than the variance of all picks multiplied with the *jackfactor*.
|
||||
*mdttolerance* | Maximum allowed deviation of P onset times from the median. Value in seconds.
|
||||
*wdttolerance* | Maximum allowed deviation of S onset times from the line during the Wadati test. Value in seconds.
|
||||
| Name | Description |
|
||||
|--------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *minAICPslope* | Initial P picks with a slope lower than this value will be discared. |
|
||||
| *minAICPSNR* | Initial P picks with a SNR below this value will be discarded. |
|
||||
| *minAICSslope* | Initial S picks with a slope lower than this value will be discarded. |
|
||||
| *minAICSSNR* | Initial S picks with a SNR below this value will be discarded. |
|
||||
| *minsiglength*, *noisefacor*. *minpercent* | Parameters for checking signal length. In the time window of size *minsiglength* after the initial P pick *minpercent* of samples have to be larger than the RMS value. |
|
||||
| *zfac* | To recognize misattributed S picks, the RMS amplitude of vertical and horizontal traces are compared. The RMS amplitude of the vertical traces has to be at least *zfac* higher than the RMS amplitude on the horizontal traces for the pick to be accepted as a valid P pick. |
|
||||
| *jackfactor* | A P pick is removed if the jackknife pseudo value of the variance of his subgroup is larger than the variance of all picks multiplied with the *jackfactor*. |
|
||||
| *mdttolerance* | Maximum allowed deviation of P onset times from the median. Value in seconds. |
|
||||
| *wdttolerance* | Maximum allowed deviation of S onset times from the line during the Wadati test. Value in seconds. |
|
||||
|
||||
## Pick quality determination
|
||||
|
||||
Parameters for discrete quality classes.
|
||||
|
||||
Name | Description
|
||||
--- | ---
|
||||
*timeerrorsP* | Width of the time windows in seconds between earliest and latest possible pick which represent the quality classes 0, 1, 2, 3 for P onsets.
|
||||
*timeerrorsS* | Width of the time windows in seconds between earliest and latest possible pick which represent the quality classes 0, 1, 2, 3 for S onsets.
|
||||
*nfacP*, *nfacS* | For determination of latest possible onset time. The time when the signal reaches an amplitude of *nfac* * mean value of the RMS amplitude in the time window *tnoise* corresponds to the latest possible onset time.
|
||||
| Name | Description |
|
||||
|------------------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
|
||||
| *timeerrorsP* | Width of the time windows in seconds between earliest and latest possible pick which represent the quality classes 0, 1, 2, 3 for P onsets. |
|
||||
| *timeerrorsS* | Width of the time windows in seconds between earliest and latest possible pick which represent the quality classes 0, 1, 2, 3 for S onsets. |
|
||||
| *nfacP*, *nfacS* | For determination of latest possible onset time. The time when the signal reaches an amplitude of *nfac* * mean value of the RMS amplitude in the time window *tnoise* corresponds to the latest possible onset time. |
|
||||
|
||||
|
||||
@ -1,6 +1,5 @@
|
||||
#!/usr/bin/env python
|
||||
# -*- coding: utf-8 -*-
|
||||
import pdb
|
||||
import glob
|
||||
import matplotlib.pyplot as plt
|
||||
import numpy as np
|
||||
@ -393,9 +392,6 @@ def reassess_pilot_event(root_dir, db_dir, event_id, out_dir=None, fn_param=None
|
||||
from pylot.core.pick.utils import earllatepicker
|
||||
# TODO: change root to datapath
|
||||
|
||||
if fn_param is None:
|
||||
fn_param = defaults.AUTOMATIC_DEFAULTS
|
||||
|
||||
default = PylotParameter(fn_param, verbosity)
|
||||
|
||||
search_base = os.path.join(root_dir, db_dir, event_id)
|
||||
@ -518,7 +514,7 @@ def writephases(arrivals, fformat, filename, parameter=None, eventinfo=None):
|
||||
# write header
|
||||
fid.write('# EQEVENT: %s Label: EQ%s Loc: X 0.00 Y 0.00 Z 10.00 OT 0.00 \n' %
|
||||
(parameter.get('database'), parameter.get('eventID')))
|
||||
arrivals = chooseArrivals(arrivals)
|
||||
arrivals = chooseArrivals(arrivals) # MP MP what is chooseArrivals? It is not defined anywhere
|
||||
for key in arrivals:
|
||||
# P onsets
|
||||
if arrivals[key].has_key('P'):
|
||||
@ -594,7 +590,7 @@ def writephases(arrivals, fformat, filename, parameter=None, eventinfo=None):
|
||||
# write header
|
||||
fid.write(' %s\n' %
|
||||
parameter.get('eventID'))
|
||||
arrivals = chooseArrivals(arrivals)
|
||||
arrivals = chooseArrivals(arrivals) # MP MP what is chooseArrivals? It is not defined anywhere
|
||||
for key in arrivals:
|
||||
if arrivals[key]['P']['weight'] < 4:
|
||||
stat = key
|
||||
@ -671,7 +667,7 @@ def writephases(arrivals, fformat, filename, parameter=None, eventinfo=None):
|
||||
fid = open("%s" % filename, 'w')
|
||||
# write header
|
||||
fid.write('%s, event %s \n' % (parameter.get('database'), parameter.get('eventID')))
|
||||
arrivals = chooseArrivals(arrivals)
|
||||
arrivals = chooseArrivals(arrivals) # MP MP what is chooseArrivals? It is not defined anywhere
|
||||
for key in arrivals:
|
||||
# P onsets
|
||||
if arrivals[key].has_key('P') and arrivals[key]['P']['mpp'] is not None:
|
||||
@ -962,7 +958,7 @@ def writephases(arrivals, fformat, filename, parameter=None, eventinfo=None):
|
||||
erh, erz, eventinfo.magnitudes[0]['mag'],
|
||||
hashID))
|
||||
# Prefer Manual Picks over automatic ones if possible
|
||||
arrivals = chooseArrivals(arrivals)
|
||||
arrivals = chooseArrivals(arrivals) # MP MP what is chooseArrivals? It is not defined anywhere
|
||||
# write phase lines
|
||||
for key in arrivals:
|
||||
if arrivals[key].has_key('P'):
|
||||
|
||||
@ -509,7 +509,7 @@ class PragPicker(AutoPicker):
|
||||
self.Pick = pick_l
|
||||
pickflag = 1
|
||||
elif flagpick_l > 0 and flagpick_r > 0 and cfpick_l >= cfpick_r:
|
||||
self.Pick = pick_r
|
||||
self.Pick = pick_r # MP MP there is no pick_r defined, commented out after commit of LK on 29.07.2020 (see above)
|
||||
pickflag = 1
|
||||
elif flagpick_l == 0 and flagpick_r > 0 and cfpick_l >= cfpick_r:
|
||||
self.Pick = pick_l
|
||||
|
||||
@ -2,12 +2,13 @@
|
||||
# -*- coding: utf-8 -*-
|
||||
|
||||
try:
|
||||
# noinspection PyUnresolvedReferences
|
||||
from urllib2 import urlopen
|
||||
except:
|
||||
from urllib.request import urlopen
|
||||
|
||||
|
||||
def checkurl(url='https://ariadne.geophysik.ruhr-uni-bochum.de/trac/PyLoT/'):
|
||||
def checkurl(url='https://git.geophysik.ruhr-uni-bochum.de/marcel/pylot/'):
|
||||
"""
|
||||
check if URL is available
|
||||
:param url: url
|
||||
|
||||
@ -351,25 +351,25 @@ def check_time(datetime):
|
||||
:type datetime: list
|
||||
:return: returns True if Values are in supposed range, returns False otherwise
|
||||
|
||||
>>> check_time([1999, 01, 01, 23, 59, 59, 999000])
|
||||
>>> check_time([1999, 1, 1, 23, 59, 59, 999000])
|
||||
True
|
||||
>>> check_time([1999, 01, 01, 23, 59, 60, 999000])
|
||||
>>> check_time([1999, 1, 1, 23, 59, 60, 999000])
|
||||
False
|
||||
>>> check_time([1999, 01, 01, 23, 59, 59, 1000000])
|
||||
>>> check_time([1999, 1, 1, 23, 59, 59, 1000000])
|
||||
False
|
||||
>>> check_time([1999, 01, 01, 23, 60, 59, 999000])
|
||||
>>> check_time([1999, 1, 1, 23, 60, 59, 999000])
|
||||
False
|
||||
>>> check_time([1999, 01, 01, 23, 60, 59, 999000])
|
||||
>>> check_time([1999, 1, 1, 23, 60, 59, 999000])
|
||||
False
|
||||
>>> check_time([1999, 01, 01, 24, 59, 59, 999000])
|
||||
>>> check_time([1999, 1, 1, 24, 59, 59, 999000])
|
||||
False
|
||||
>>> check_time([1999, 01, 31, 23, 59, 59, 999000])
|
||||
>>> check_time([1999, 1, 31, 23, 59, 59, 999000])
|
||||
True
|
||||
>>> check_time([1999, 02, 30, 23, 59, 59, 999000])
|
||||
>>> check_time([1999, 2, 30, 23, 59, 59, 999000])
|
||||
False
|
||||
>>> check_time([1999, 02, 29, 23, 59, 59, 999000])
|
||||
>>> check_time([1999, 2, 29, 23, 59, 59, 999000])
|
||||
False
|
||||
>>> check_time([2000, 02, 29, 23, 59, 59, 999000])
|
||||
>>> check_time([2000, 2, 29, 23, 59, 59, 999000])
|
||||
True
|
||||
>>> check_time([2000, 13, 29, 23, 59, 59, 999000])
|
||||
False
|
||||
@ -444,7 +444,7 @@ def evt_head_check(root_dir, out_dir=None):
|
||||
"""
|
||||
if not out_dir:
|
||||
print('WARNING files are going to be overwritten!')
|
||||
inp = str(raw_input('Continue? [y/N]'))
|
||||
inp = str(input('Continue? [y/N]'))
|
||||
if not inp == 'y':
|
||||
sys.exit()
|
||||
filelist = get_file_list(root_dir)
|
||||
|
||||
@ -1,35 +0,0 @@
|
||||
import unittest
|
||||
from pylot.core.pick.autopick import PickingParameters
|
||||
|
||||
class TestPickingParameters(unittest.TestCase):
|
||||
|
||||
def setUp(self):
|
||||
self.simple_dict = {'a': 3, 'b': 14}
|
||||
self.nested_dict = {'a': self.simple_dict, 'b': self.simple_dict}
|
||||
|
||||
def assertParameterEquality(self, dic, instance):
|
||||
"""Test wether all parameters given in dic are found in instance"""
|
||||
for key, value in dic.items():
|
||||
self.assertEqual(value, getattr(instance, key))
|
||||
|
||||
def test_add_params_from_dict_simple(self):
|
||||
pickparam = PickingParameters()
|
||||
pickparam.add_params_from_dict(self.simple_dict)
|
||||
self.assertParameterEquality(self.simple_dict, pickparam)
|
||||
|
||||
def test_add_params_from_dict_nested(self):
|
||||
pickparam = PickingParameters()
|
||||
pickparam.add_params_from_dict(self.nested_dict)
|
||||
self.assertParameterEquality(self.nested_dict, pickparam)
|
||||
|
||||
def test_init(self):
|
||||
pickparam = PickingParameters(self.simple_dict)
|
||||
self.assertParameterEquality(self.simple_dict, pickparam)
|
||||
|
||||
def test_dot_access(self):
|
||||
pickparam = PickingParameters(self.simple_dict)
|
||||
self.assertEqual(pickparam.a, self.simple_dict['a'])
|
||||
|
||||
|
||||
if __name__ == '__main__':
|
||||
unittest.main()
|
||||
Loading…
Reference in New Issue
Block a user