ECAT7ImageInputFileFormat.h

Go to the documentation of this file.

//
//
#ifndef __stir_IO_ECAT7ImageInputFileFormat_h__
#define __stir_IO_ECAT7ImageInputFileFormat_h__
/*
    Copyright (C) 2006- 2013, Hammersmith Imanet Ltd
    This file is part of STIR.
    SPDX-License-Identifier: Apache-2.0

    See STIR/LICENSE.txt for details
*/
#include "stir/IO/InputFileFormat.h"
#include "stir/utilities.h"
#include "stir/VoxelsOnCartesianGrid.h"
#include "stir/error.h"
#include <fstream>
#include <string>

#ifndef HAVE_LLN_MATRIX
#  error HAVE_LLN_MATRIX not define: you need the lln ecat library.
#endif

#include "stir/IO/stir_ecat7.h"
START_NAMESPACE_STIR
START_NAMESPACE_ECAT
START_NAMESPACE_ECAT7


class ECAT7ImageInputFileFormat : public InputFileFormat<DiscretisedDensity<3, float>>
{
public:
  virtual const std::string get_name() const { return "ECAT7"; }

protected:
  virtual bool actual_can_read(const FileSignature& signature, std::istream& input) const
  {
    if (strncmp(signature.get_signature(), "MATRIX", 6) != 0)
      return false;

    // TODO
    // return (is_ECAT7_image_file(filename))
    return true;
  }

  virtual unique_ptr<data_type> read_from_file(std::istream& input) const
  {
    // TODO
    error("read_from_file for ECAT7 with istream not implemented %s:%d. Sorry", __FILE__, __LINE__);
    return unique_ptr<data_type>();
  }
  virtual unique_ptr<data_type> read_from_file(const std::string& filename) const
  {

    if (is_ECAT7_image_file(filename))
      {
        warning("\nReading frame 1, gate 1, data 0, bed 0 from file %s\n", filename.c_str());
        return unique_ptr<data_type>(ECAT7_to_VoxelsOnCartesianGrid(filename,
                                                                    /*frame_num, gate_num, data_num, bed_num*/ 1,
                                                                    1,
                                                                    0,
                                                                    0));
      }
    else
      {
        error("ECAT7 file %s is not an image file", filename.c_str());
        return unique_ptr<data_type>();
      }
  }
};

END_NAMESPACE_ECAT
END_NAMESPACE_ECAT7
END_NAMESPACE_STIR

#endif