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