ProjMatrixByBinFromFile.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2004- 2008, Hammersmith Imanet Ltd
    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0

    See STIR/LICENSE.txt for details
*/
#ifndef __stir_recon_buildblock_ProjMatrixByBinFromFile__
#define __stir_recon_buildblock_ProjMatrixByBinFromFile__

#include "stir/RegisteredParsingObject.h"
#include "stir/recon_buildblock/ProjMatrixByBin.h"
#include "stir/ProjDataInfo.h"
#include "stir/CartesianCoordinate3D.h"
#include "stir/IndexRange.h"
#include "stir/shared_ptr.h"
#include <iostream>

START_NAMESPACE_STIR

template <int num_dimensions, typename elemT>
class DiscretisedDensity;
class Bin;
class ProjMatrixByBinFromFile : public RegisteredParsingObject<ProjMatrixByBinFromFile, ProjMatrixByBin, ProjMatrixByBin>
{
public:
  static const char* const registered_name;


  static Succeeded write_to_file(const std::string& output_filename_prefix,
                                 const ProjMatrixByBin& proj_matrix,
                                 const shared_ptr<const ProjDataInfo>& proj_data_info_sptr,
                                 const DiscretisedDensity<3, float>& template_density);

  ProjMatrixByBinFromFile();

  void set_up(const shared_ptr<const ProjDataInfo>& proj_data_info_ptr,
              const shared_ptr<const DiscretisedDensity<3, float>>& density_info_ptr // TODO should be Info only
              ) override;

  ProjMatrixByBinFromFile* clone() const override;

private:
  std::string parsed_version;
  std::string template_density_filename;
  std::string template_proj_data_filename;
  std::string data_filename;

  std::string symmetries_type;
  // should be in symmetries
  bool do_symmetry_90degrees_min_phi;
  bool do_symmetry_180degrees_min_phi;
  bool do_symmetry_swap_segment;
  bool do_symmetry_swap_s;
  bool do_symmetry_shift_z;

  // TODO this only works as long as we only have VoxelsOnCartesianGrid
  // explicitly list necessary members for image details (should use an Info object instead)
  CartesianCoordinate3D<float> voxel_size;
  CartesianCoordinate3D<float> origin;
  IndexRange<3> densel_range;

  shared_ptr<const ProjDataInfo> proj_data_info_ptr;

  void calculate_proj_matrix_elems_for_one_bin(ProjMatrixElemsForOneBin&) const override;

  void set_defaults() override;
  void initialise_keymap() override;
  bool post_processing() override;

  Succeeded read_data();
};

END_NAMESPACE_STIR

#endif