bad_channel_detect.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef BAD_CHANNEL_DETECT_DSP_H
#define BAD_CHANNEL_DETECT_DSP_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "dsp_global.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QVector>
 
//=============================================================================================================
// DEFINE NAMESPACE UTILSLIB
//=============================================================================================================
 
namespace UTILSLIB
{
 
//=============================================================================================================
struct DSPSHARED_EXPORT BadChannelDetectParams
{
    // Flat signal
    double dFlatThreshold = 1e-13; 
 
    // High variance
    double dVarZThresh = 4.0; 
 
    // Low correlation
    double dCorrThresh = 0.4; 
    int iNeighbours = 5;      
};
 
//=============================================================================================================
class DSPSHARED_EXPORT BadChannelDetect
{
public:
    using Params = BadChannelDetectParams; 
 
    //=========================================================================================================
    static QVector<int> detect(const Eigen::MatrixXd& matData,
                               const Params& params = Params());
 
    //=========================================================================================================
    static QVector<int> detectFlat(const Eigen::MatrixXd& matData,
                                   double dThreshold = 1e-13);
 
    //=========================================================================================================
    static QVector<int> detectHighVariance(const Eigen::MatrixXd& matData,
                                           double dZThresh = 4.0);
 
    //=========================================================================================================
    static QVector<int> detectLowCorrelation(const Eigen::MatrixXd& matData,
                                             double dCorrThresh = 0.4,
                                             int iNeighbours = 5);
 
private:
    //=========================================================================================================
    static double pearsonCorr(const Eigen::RowVectorXd& a, const Eigen::RowVectorXd& b);
 
    //=========================================================================================================
    static double median(QVector<double> values);
};
 
} // namespace UTILSLIB
 
#endif // BAD_CHANNEL_DETECT_DSP_H