ecdview.cpp

Go to the documentation of this file.

//=============================================================================================================
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "ecdview.h"
 
#include "../engine/model/data3Dtreemodel.h"
 
#include <inverse/dipoleFit/dipole_fit_settings.h>
#include <inverse/dipoleFit/ecd_set.h>
#include <mne/mne_bem.h>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace DISP3DLIB;
using namespace INVERSELIB;
using namespace MNELIB;
using namespace Eigen;
using namespace FIFFLIB;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
ECDView::ECDView(const DipoleFitSettings& dipFitSettings,
                 const ECDSet& ecdSet,
                 QWidget* parent,
                 Qt::WindowFlags f)
: AbstractView(parent, f)
{
    //Read mri transform
    QFile file(dipFitSettings.mriname);
    ECDSet ecdSetTrans = ecdSet;
 
    if(file.exists()) {
        FiffCoordTrans coordTrans(file);
 
        for(int i = 0; i < ecdSet.size() ; ++i) {
            MatrixX3f dipoles(1, 3);
            //transform location
            dipoles(0,0) = ecdSet[i].rd(0);
            dipoles(0,1) = ecdSet[i].rd(1);
            dipoles(0,2) = ecdSet[i].rd(2);
 
            dipoles = coordTrans.apply_trans(dipoles);
 
            ecdSetTrans[i].rd(0) = dipoles(0,0);
            ecdSetTrans[i].rd(1) = dipoles(0,1);
            ecdSetTrans[i].rd(2) = dipoles(0,2);
 
            //transform orientation
            dipoles(0,0) = ecdSet[i].Q(0);
            dipoles(0,1) = ecdSet[i].Q(1);
            dipoles(0,2) = ecdSet[i].Q(2);
 
            dipoles = coordTrans.apply_trans(dipoles, false);
 
            ecdSetTrans[i].Q(0) = dipoles(0,0);
            ecdSetTrans[i].Q(1) = dipoles(0,1);
            ecdSetTrans[i].Q(2) = dipoles(0,2);
        }
    } else {
        qCritical("ECDView::ECDView - Cannot open FiffCoordTrans file");
    }
 
    //Add ECD data
    m_pData3DModel->addDipoleFitData("sample", QString("Set %1").arg(dipFitSettings.setno), ecdSetTrans);
 
    //Read and show BEM
    QFile t_fileBem(dipFitSettings.bemname);
 
    if(t_fileBem.exists()) {
        MNEBem t_Bem(t_fileBem);
        m_pData3DModel->addBemData("sample", "BEM", t_Bem);
    } else {
        qCritical("ECDView::ECDView - Cannot open MNEBem file");
    }
}
 
//=============================================================================================================
 
ECDView::~ECDView()
{
}