fiff_digitizer_data.cpp

Go to the documentation of this file.

//=============================================================================================================
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "fiff_digitizer_data.h"
#include "fiff_coord_trans_old.h"
#include "../fiff_dig_point.h"
#include <iostream>
 
#ifndef TRUE
#define TRUE 1
#endif
 
#ifndef FALSE
#define FALSE 0
#endif
 
#define FREE_43(x) if ((char *)(x) != Q_NULLPTR) free((char *)(x))
 
#define MALLOC_43(x,t) (t *)malloc((x)*sizeof(t))
 
#define FREE_CMATRIX_43(m) mne_free_cmatrix_43((m))
 
void mne_free_cmatrix_43(float **m)
{
    if (m) {
        FREE_43(*m);
        FREE_43(m);
    }
}
 
//=============================================================================================================
// USED NAMESPACES
//=============================================================================================================
 
using namespace FIFFLIB;
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
FiffDigitizerData::FiffDigitizerData()
: head_mri_t(Q_NULLPTR)
, head_mri_t_adj(Q_NULLPTR)
, coord_frame(FIFFV_COORD_UNKNOWN)
, npoint(0)
, mri_fids(Q_NULLPTR)
, nfids(0)
, show(FALSE)
, show_minimal(FALSE)
, dist(Q_NULLPTR)
, closest(Q_NULLPTR)
, closest_point(Q_NULLPTR)
, dist_valid(FALSE)
{
}
 
//=============================================================================================================
 
FiffDigitizerData::FiffDigitizerData(const FiffDigitizerData& p_FiffDigitizerData)
: head_mri_t(p_FiffDigitizerData.head_mri_t)
, head_mri_t_adj(p_FiffDigitizerData.head_mri_t_adj)
, points(p_FiffDigitizerData.points)
, coord_frame(p_FiffDigitizerData.coord_frame)
, active(p_FiffDigitizerData.active)
, discard(p_FiffDigitizerData.discard)
, npoint(p_FiffDigitizerData.npoint)
, mri_fids(p_FiffDigitizerData.mri_fids)
, nfids(p_FiffDigitizerData.nfids)
, show(p_FiffDigitizerData.show)
, show_minimal(p_FiffDigitizerData.show_minimal)
, dist(p_FiffDigitizerData.dist)
, closest(p_FiffDigitizerData.closest)
, closest_point(p_FiffDigitizerData.closest_point)
, dist_valid(p_FiffDigitizerData.dist_valid)
{
}
 
//=============================================================================================================
 
FiffDigitizerData::FiffDigitizerData(QIODevice &p_IODevice)
: head_mri_t(Q_NULLPTR)
, head_mri_t_adj(Q_NULLPTR)
, coord_frame(FIFFV_COORD_UNKNOWN)
, npoint(0)
, mri_fids(Q_NULLPTR)
, nfids(0)
, show(FALSE)
, show_minimal(FALSE)
, dist(Q_NULLPTR)
, closest(Q_NULLPTR)
, closest_point(Q_NULLPTR)
, dist_valid(FALSE)
{
    // Open the io device
    FiffStream::SPtr t_pStream(new FiffStream(&p_IODevice));
    bool open_here = false;
 
    //Open if the device and stream have not been openend already
    if (!t_pStream->device()->isOpen()) {
        if(!t_pStream->open()) {
            qWarning() << "Warning in FiffDigitizerData::FiffDigitizerData - Could not open the didigitzer data file"; // ToDo throw error
            return;
        }
 
        open_here = true;
    }
 
    // If device is open read the data
    if(!t_pStream->read_digitizer_data(t_pStream->dirtree(), *this)) {
        qWarning() << "Warning in FiffDigitizerData::FiffDigitizerData - Could not read the FiffDigitizerData"; // ToDo throw error
    }
 
    // If stream has been opened in this function also close here again
    if(open_here) {
        t_pStream->close();
    }
}
 
//=============================================================================================================
 
FiffDigitizerData::~FiffDigitizerData()
{
    FREE_43(head_mri_t);
    FREE_43(head_mri_t_adj);
    FREE_43(dist);
    FREE_43(closest);
    FREE_CMATRIX_43(closest_point);
    FREE_43(mri_fids);
}
 
//=============================================================================================================
 
void FiffDigitizerData::print() const
{
    std::cout << "Number of digitizer points: " << points.size() << "\n";
 
    switch(coord_frame){
    case FIFFV_COORD_MRI:
        std::cout << "Coord. Frame: FIFFV_COORD_MRI \n";
         break;
    case FIFFV_COORD_HEAD:
        std::cout << "Coord. Frame: FIFFV_COORD_HEAD \n";
        break;
    }
 
    for (auto& point : points){
        if (point.kind == FIFFV_POINT_HPI){
            std::cout << "HPI Point " << point.ident << " - " << point.r[0] << ", " << point.r[1] << ", " << point.r[2] << "\n";
        }
    }
 
    std::cout << "Number of MRI fiducials: " << nfids << "\n";
 
    if (head_mri_t){
 
    }
}