inv_sensor_set.cpp

Go to the documentation of this file.

//=============================================================================================================
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "inv_sensor_set.h"
#include <iostream>
#include <fwd/fwd_coil_set.h>
#include <fiff/fiff_ch_info.h>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QCoreApplication>
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Dense>
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace INVLIB;
using namespace FIFFLIB;
using namespace FWDLIB;
using namespace Eigen;
 
//=============================================================================================================
// DEFINE GLOBAL METHODS
//=============================================================================================================
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
InvSensorSet::InvSensorSet(const FwdCoilSet::SPtr pFwdSensorSet)
{
    if (pFwdSensorSet != nullptr) {
        initFromFwdCoilSet(pFwdSensorSet);
    }
}
 
//=============================================================================================================
 
void InvSensorSet::initFromFwdCoilSet(const QSharedPointer<FWDLIB::FwdCoilSet> pFwdSensorSet)
{
    m_ncoils = pFwdSensorSet->ncoil();
    m_np = pFwdSensorSet->coils[0]->np;
    initMatrices(m_ncoils, m_np);
 
    // get data froms Fwd Coilset
    for (int i = 0; i < m_ncoils; i++) {
        const FwdCoil* coil = pFwdSensorSet->coils[i].get();
        MatrixXd matRmag = MatrixXd::Zero(m_np, 3);
        MatrixXd matCosmag = MatrixXd::Zero(m_np, 3);
        RowVectorXd vecW(m_np);
 
        m_r0(i, 0) = coil->r0[0];
        m_r0(i, 1) = coil->r0[1];
        m_r0(i, 2) = coil->r0[2];
 
        m_ez(i, 0) = coil->ez[0];
        m_ez(i, 1) = coil->ez[1];
        m_ez(i, 2) = coil->ez[2];
 
        for (int p = 0; p < m_np; p++) {
            m_w(i * m_np + p) = coil->w[p];
            for (int c = 0; c < 3; c++) {
                matRmag(p, c) = coil->rmag(p, c);
                matCosmag(p, c) = coil->cosmag(p, c);
            }
        }
 
        m_cosmag.block(i * m_np, 0, m_np, 3) = matCosmag;
        m_rmag.block(i * m_np, 0, m_np, 3) = matRmag;
    }
    m_tra = MatrixXd::Identity(m_ncoils, m_ncoils);
}
 
//=============================================================================================================
 
void InvSensorSet::initMatrices(int ncoils, int np)
{
    m_ez = MatrixXd(ncoils, 3);
    m_r0 = MatrixXd(ncoils, 3);
    m_rmag = MatrixXd(ncoils * np, 3);
    m_cosmag = MatrixXd(ncoils * np, 3);
    m_cosmag = MatrixXd(ncoils * np, 3);
    m_tra = MatrixXd(ncoils, ncoils);
    m_w = RowVectorXd(ncoils * np);
}
 
//=============================================================================================================
 
InvSensorSetCreator::InvSensorSetCreator()
{
    const QString qPath = QString(QCoreApplication::applicationDirPath() + "/../resources/general/coilDefinitions/coil_def.dat");
    m_pCoilDefinitions = FwdCoilSet::SPtr(FwdCoilSet::read_coil_defs(qPath).release());
}
 
//=============================================================================================================
 
InvSensorSet InvSensorSetCreator::updateSensorSet(const QList<FIFFLIB::FiffChInfo>& channelList,
                                                  const Accuracy& accuracy)
{
    if (channelList.isEmpty()) {
        return InvSensorSet();
    } else {
        auto pCoilMeg = FwdCoilSet::SPtr(m_pCoilDefinitions->create_meg_coils(channelList, channelList.size(), static_cast<int>(accuracy)).release());
        return InvSensorSet(pCoilMeg);
    }
}
 
//=============================================================================================================