You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
83 lines
3.1 KiB
C
83 lines
3.1 KiB
C
1 year ago
|
/***************************************************************************
|
||
|
|
* Copyright (C) gempa GmbH *
|
||
|
|
* All rights reserved. *
|
||
|
|
* Contact: gempa GmbH (seiscomp-dev@gempa.de) *
|
||
|
|
* *
|
||
|
|
* GNU Affero General Public License Usage *
|
||
|
|
* This file may be used under the terms of the GNU Affero *
|
||
|
|
* Public License version 3.0 as published by the Free Software Foundation *
|
||
|
|
* and appearing in the file LICENSE included in the packaging of this *
|
||
|
|
* file. Please review the following information to ensure the GNU Affero *
|
||
|
|
* Public License version 3.0 requirements will be met: *
|
||
|
|
* https://www.gnu.org/licenses/agpl-3.0.html. *
|
||
|
|
* *
|
||
|
|
* Other Usage *
|
||
|
|
* Alternatively, this file may be used in accordance with the terms and *
|
||
|
|
* conditions contained in a signed written agreement between you and *
|
||
|
|
* gempa GmbH. *
|
||
|
|
***************************************************************************/
|
||
|
|
|
||
|
|
|
||
|
|
#ifndef SEISCOMP_PROCESSING_BKPICKER_H
|
||
|
|
#define SEISCOMP_PROCESSING_BKPICKER_H
|
||
|
|
|
||
|
|
|
||
|
|
#include <seiscomp/processing/picker.h>
|
||
|
|
|
||
|
|
|
||
|
|
namespace Seiscomp {
|
||
|
|
namespace Processing {
|
||
|
|
|
||
|
|
|
||
|
|
class SC_SYSTEM_CLIENT_API BKPicker : public Picker {
|
||
|
|
// ----------------------------------------------------------------------
|
||
|
|
// X'truction
|
||
|
|
// ----------------------------------------------------------------------
|
||
|
|
public:
|
||
|
|
//! C'tor
|
||
|
|
BKPicker();
|
||
|
|
BKPicker(const Core::Time& trigger);
|
||
|
|
//! D'tor
|
||
|
|
~BKPicker();
|
||
|
|
|
||
|
|
|
||
|
|
// ----------------------------------------------------------------------
|
||
|
|
// Public Interface
|
||
|
|
// ----------------------------------------------------------------------
|
||
|
|
public:
|
||
|
|
virtual bool setup(const Settings &settings) override;
|
||
|
|
const std::string &methodID() const override;
|
||
|
|
const std::string &filterID() const override;
|
||
|
|
|
||
|
|
|
||
|
|
// ----------------------------------------------------------------------
|
||
|
|
// Protected Interface
|
||
|
|
// ----------------------------------------------------------------------
|
||
|
|
protected:
|
||
|
|
// filter settings
|
||
|
|
std::string filterType;
|
||
|
|
int filterPoles; // number of poles
|
||
|
|
double f1; // bandpass lower cutoff freq. in Hz
|
||
|
|
double f2; // bandpass upper cutoff freq. in Hz
|
||
|
|
std::string usedFilter;
|
||
|
|
|
||
|
|
// picker parameters
|
||
|
|
double thrshl1; // threshold to trigger for pick (c.f. paper), default 10
|
||
|
|
double thrshl2; // threshold for updating sigma (c.f. paper), default 20
|
||
|
|
|
||
|
|
int debugOutput;
|
||
|
|
|
||
|
|
void bk_wrapper(int n, double *data, int &kmin, double &snr, double samplespersec=120);
|
||
|
|
bool calculatePick(int n, const double *data,
|
||
|
|
int signalStartIdx, int signalEndIdx,
|
||
|
|
int &triggerIdx, int &lowerUncertainty,
|
||
|
|
int &upperUncertainty, double &snr,
|
||
|
|
OPT(Polarity) &polarity) override;
|
||
|
|
};
|
||
|
|
|
||
|
|
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
#endif
|