filter.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef FILTER_RTPROCESSING_H
#define FILTER_RTPROCESSING_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "rtprocessing_global.h"
 
#include "helpers/filterkernel.h"
 
#include <fiff/fiff_info.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
{

typedef struct {
    FilterKernel filterKernel;
    int iRow;
    Eigen::RowVectorXd vecData;
} FilterObject;
 
//=========================================================================================================
RTPROCESINGSHARED_EXPORT bool filterFile(QIODevice& pIODevice,
                                         QSharedPointer<FIFFLIB::FiffRawData> pFiffRawData,
                                         int type,
                                         double dCenterfreq,
                                         double dBandwidth,
                                         double dTransition,
                                         double dSFreq,
                                         int iOrder = 4096,
                                         int designMethod = FilterKernel::m_designMethods.indexOf(FilterParameter("Cosine")),
                                         const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                                         bool bUseThreads = true);
 
//=========================================================================================================
RTPROCESINGSHARED_EXPORT bool filterFile(QIODevice& pIODevice,
                                         QSharedPointer<FIFFLIB::FiffRawData> pFiffRawData,
                                         const RTPROCESSINGLIB::FilterKernel& filterKernel,
                                         const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                                         bool bUseThreads = false);
 
//=========================================================================================================
RTPROCESINGSHARED_EXPORT Eigen::MatrixXd filterData(const Eigen::MatrixXd& matData,
                                                    int type,
                                                    double dCenterfreq,
                                                    double dBandwidth,
                                                    double dTransition,
                                                    double dSFreq,
                                                    int iOrder = 1024,
                                                    int designMethod = FilterKernel::m_designMethods.indexOf(FilterParameter("Cosine")),
                                                    const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                                                    bool bUseThreads = true,
                                                    bool bKeepOverhead = false);
 
//=========================================================================================================
RTPROCESINGSHARED_EXPORT Eigen::MatrixXd filterData(const Eigen::MatrixXd& mataData,
                                                    const RTPROCESSINGLIB::FilterKernel& filterKernel,
                                                    const Eigen::RowVectorXi& vecPicks = Eigen::RowVectorXi(),
                                                    bool bUseThreads = true,
                                                    bool bKeepOverhead = false);
 
//=========================================================================================================
RTPROCESINGSHARED_EXPORT Eigen::MatrixXd filterDataBlock(const Eigen::MatrixXd& mataData,
                                                         const Eigen::RowVectorXi& vecPicks,
                                                         const RTPROCESSINGLIB::FilterKernel& filterKernel,
                                                         bool bUseThreads = true);
 
//=========================================================================================================
RTPROCESINGSHARED_EXPORT void filterChannel(FilterObject &channelDataTime);
 
//=============================================================================================================
class RTPROCESINGSHARED_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 = FilterKernel::m_designMethods.indexOf(FilterParameter("Cosine")),
                              const Eigen::RowVectorXi &vecPicks = Eigen::RowVectorXi(),
                              bool bFilterEnd = true,
                              bool bUseThreads = true,
                              bool bKeepOverhead = false);
 
    //=========================================================================================================
    Eigen::MatrixXd calculate(const Eigen::MatrixXd& mataData,
                              const RTPROCESSINGLIB::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 // FILTER_RTPROCESSING_H