rt_filter.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef RT_FILTER_RT_H
#define RT_FILTER_RT_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "../dsp_global.h"
 
#include "../filterkernel.h"
 
#include <fiff/fiff_info.h>
#include <fiff/fiff_evoked.h>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QSharedPointer>
#include <QtConcurrent/QtConcurrent>
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
#include <unsupported/Eigen/FFT>
 
//=============================================================================================================
// FORWARD DECLARATIONS
//=============================================================================================================
 
namespace FIFFLIB {
    class FiffRawData;
}
 
//=============================================================================================================
// DEFINE NAMESPACE RTPROCESSINGLIB
//=============================================================================================================
 
namespace RTPROCESSINGLIB
{

struct FilterObject {
    UTILSLIB::FilterKernel filterKernel;
    int iRow;
    Eigen::RowVectorXd vecData;
};
 
//=========================================================================================================
DSPSHARED_EXPORT bool filterFile(QIODevice& pIODevice,
                                         QSharedPointer<FIFFLIB::FiffRawData> pFiffRawData,
                                         int type,
                                         double dCenterfreq,
                                         double dBandwidth,
                                         double dTransition,
                                         double dSFreq,
                                         int iOrder = 4096,
                                         int designMethod = UTILSLIB::FilterKernel::m_designMethods.indexOf(UTILSLIB::FilterParameter("Cosine")),
                                         const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                                         bool bUseThreads = true);
 
//=========================================================================================================
DSPSHARED_EXPORT bool filterFile(QIODevice& pIODevice,
                                         QSharedPointer<FIFFLIB::FiffRawData> pFiffRawData,
                                         const UTILSLIB::FilterKernel& filterKernel,
                                         const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                                         bool bUseThreads = false);
 
//=========================================================================================================
DSPSHARED_EXPORT Eigen::MatrixXd filterData(const Eigen::MatrixXd& matData,
                                                    int type,
                                                    double dCenterfreq,
                                                    double dBandwidth,
                                                    double dTransition,
                                                    double dSFreq,
                                                    int iOrder = 1024,
                                                    int designMethod = UTILSLIB::FilterKernel::m_designMethods.indexOf(UTILSLIB::FilterParameter("Cosine")),
                                                    const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                                                    bool bUseThreads = true,
                                                    bool bKeepOverhead = false);
 
//=========================================================================================================
DSPSHARED_EXPORT Eigen::MatrixXd filterData(const Eigen::MatrixXd& matData,
                                                    const UTILSLIB::FilterKernel& filterKernel,
                                                    const Eigen::RowVectorXi& vecPicks = Eigen::RowVectorXi(),
                                                    bool bUseThreads = true,
                                                    bool bKeepOverhead = false);
 
//=========================================================================================================
DSPSHARED_EXPORT Eigen::MatrixXd filterDataBlock(const Eigen::MatrixXd& matData,
                                                         const Eigen::RowVectorXi& vecPicks,
                                                         const UTILSLIB::FilterKernel& filterKernel,
                                                         bool bUseThreads = true);
 
//=========================================================================================================
DSPSHARED_EXPORT void filterChannel(FilterObject &channelDataTime);
 
//=========================================================================================================
DSPSHARED_EXPORT FIFFLIB::FiffEvoked computeFilteredAverage(const FIFFLIB::FiffRawData& raw,
                                                                    const Eigen::MatrixXi& matEvents,
                                                                    float fTMinS,
                                                                    float fTMaxS,
                                                                    qint32 eventType,
                                                                    bool bApplyBaseline,
                                                                    float fTBaselineFromS,
                                                                    float fTBaselineToS,
                                                                    const QMap<QString,double>& mapReject,
                                                                    const UTILSLIB::FilterKernel& filterKernel,
                                                                    const QStringList &lExcludeChs = QStringList(),
                                                                    const Eigen::RowVectorXi& vecPicks = Eigen::RowVectorXi());
 
//=============================================================================================================
class DSPSHARED_EXPORT FilterOverlapAdd
{
public:
    typedef QSharedPointer<FilterOverlapAdd> SPtr;             
    typedef QSharedPointer<const FilterOverlapAdd> ConstSPtr;  
 
    //=========================================================================================================
    Eigen::MatrixXd calculate(const Eigen::MatrixXd& matData,
                              int type,
                              double dCenterfreq,
                              double dBandwidth,
                              double dTransition,
                              double dSFreq,
                              int iOrder = 1024,
                              int designMethod = UTILSLIB::FilterKernel::m_designMethods.indexOf(UTILSLIB::FilterParameter("Cosine")),
                              const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                              bool bFilterEnd = true,
                              bool bUseThreads = true,
                              bool bKeepOverhead = false);
 
    //=========================================================================================================
    Eigen::MatrixXd calculate(const Eigen::MatrixXd& matData,
                              const UTILSLIB::FilterKernel& filterKernel,
                              const Eigen::RowVectorXi& vecPicks = Eigen::RowVectorXi(),
                              bool bFilterEnd = true,
                              bool bUseThreads = true,
                              bool bKeepOverhead = false);
 
    //=========================================================================================================
    void reset();
 
private:
    Eigen::MatrixXd                 m_matOverlapBack;                   
    Eigen::MatrixXd                 m_matOverlapFront;                  
};
 
//=============================================================================================================
// INLINE DEFINITIONS
//=============================================================================================================
 
} // NAMESPACE
 
#endif // RT_FILTER_RT_H