Polaris_MT_File.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2003- 2010, Hammersmith Imanet Ltd
    For internal GE use only.
*/
#ifndef __stir_Polaris_MT_File__
#define __stir_Polaris_MT_File__

#include "stir_experimental/Quaternion.h"
#include "stir/CartesianCoordinate3D.h"
#include "stir/Succeeded.h"

#include <vector>
#include <string>
#include <ctime>

#ifdef BOOST_NO_STDC_NAMESPACE
namespace std
{
using ::time_t;
using ::tm;
using ::localtime;
} // namespace std
#endif

START_NAMESPACE_STIR

class Polaris_MT_File
{

public:
  struct Record
  {
    double sample_time;
    unsigned int frame_num;
    unsigned int rand_num;
    char total_num;
    Quaternion<float> quat;
    CartesianCoordinate3D<float> trans;
    float rms;
    unsigned int out_of_FOV;
  };

  typedef std::vector<Record>::const_iterator const_iterator;

  ~Polaris_MT_File(){};
  Polaris_MT_File(const std::string& filename);


  Record operator[](unsigned int n) const;

  const_iterator begin() const { return vector_of_records.begin(); }
  const_iterator end() const { return vector_of_records.end(); }
  unsigned long num_samples() const { return vector_of_records.size(); }

  const_iterator begin_all_tags() const { return vector_of_tags.begin(); }
  const_iterator end_all_tags() const { return vector_of_tags.end(); }
  unsigned long num_tags() const { return vector_of_tags.size(); }

  std::time_t get_start_time_in_secs_since_1970();

private:
  std::time_t start_time_in_secs_since_1970;

  std::vector<Record> vector_of_records;
  // this contains all tags and times (even those with 'missing data')
  std::vector<Record> vector_of_tags;

  void read_Peter_Bloomfield_mt_file(const std::string& mt_filename, std::istream& mt_stream, const char* const first_line);
  void read_NDI_Toolviewer_mt_file(const std::string& mt_filename, std::istream& mt_stream, const char* const first_line);
};

END_NAMESPACE_STIR

#endif