ProjMatrixByBinUsingRayTracing.h

Go to the documentation of this file.

/*
    Copyright (C) 2000 PARAPET partners
    Copyright (C) 2000-2009, Hammersmith Imanet Ltd
    Copyright (C) 2014, University College London
    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0 AND License-ref-PARAPET-license

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

#include "stir/RegisteredParsingObject.h"
#include "stir/recon_buildblock/ProjMatrixByBin.h"
#include "stir/CartesianCoordinate3D.h"
#include "stir/shared_ptr.h"

START_NAMESPACE_STIR

template <int num_dimensions, typename elemT>
class DiscretisedDensity;
class ProjDataInfo;

class ProjMatrixByBinUsingRayTracing
    : public RegisteredParsingObject<ProjMatrixByBinUsingRayTracing, ProjMatrixByBin, ProjMatrixByBin>
{
public:
  static const char* const registered_name;

  ProjMatrixByBinUsingRayTracing();


  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;

  ProjMatrixByBinUsingRayTracing* clone() const override;

  bool get_restrict_to_cylindrical_FOV() const;
  void set_restrict_to_cylindrical_FOV(bool);
  int get_num_tangential_LORs() const;
  void set_num_tangential_LORs(int);

  bool get_use_actual_detector_boundaries() const;
  void set_use_actual_detector_boundaries(bool);

  bool get_do_symmetry_90degrees_min_phi() const;
  void set_do_symmetry_90degrees_min_phi(bool);

  bool get_do_symmetry_180degrees_min_phi() const;
  void set_do_symmetry_180degrees_min_phi(bool);

  bool get_do_symmetry_swap_segment() const;
  void set_do_symmetry_swap_segment(bool);

  bool get_do_symmetry_swap_s() const;
  void set_do_symmetry_swap_s(bool);

  bool get_do_symmetry_shift_z() const;
  void set_do_symmetry_shift_z(bool);

private:

  bool already_setup;

  bool restrict_to_cylindrical_FOV;
  int num_tangential_LORs;
  bool use_actual_detector_boundaries;
  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;

  // explicitly list necessary members for image details (should use an Info object instead)
  CartesianCoordinate3D<float> voxel_size;
  CartesianCoordinate3D<float> origin;
  CartesianCoordinate3D<int> min_index;
  CartesianCoordinate3D<int> max_index;

  void calculate_proj_matrix_elems_for_one_bin(ProjMatrixElemsForOneBin&) const override;

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

END_NAMESPACE_STIR

#endif