ProjMatrixByDenselUsingRayTracing.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2000- 2002, IRSL
    See STIR/LICENSE.txt for details
*/
#ifndef __ProjMatrixByDenselUsingRayTracing__
#define __ProjMatrixByDenselUsingRayTracing__

#include "stir/RegisteredParsingObject.h"
#include "stir_experimental/recon_buildblock/ProjMatrixByDenselOnCartesianGridUsingElement.h"
#include "stir/CartesianCoordinate3D.h"

START_NAMESPACE_STIR

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

class ProjMatrixByDenselUsingRayTracing
    : public RegisteredParsingObject<ProjMatrixByDenselUsingRayTracing, ProjMatrixByDenselOnCartesianGridUsingElement>
{
  typedef ProjMatrixByDenselOnCartesianGridUsingElement base_type;

public:
  static const char* const registered_name;

  ProjMatrixByDenselUsingRayTracing();


  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;

  const DataSymmetriesForDensels* get_symmetries_ptr() const override;

  float get_element(const Bin&, const CartesianCoordinate3D<float>&) const override;

private:
  shared_ptr<DataSymmetriesForDensels_PET_CartesianGrid> symmetries_ptr;
#if 0
  bool restrict_to_cylindrical_FOV;
#endif
  int num_tangential_LORs;
  bool use_actual_detector_boundaries;

  // explicitly list necessary members for image details (should use an Info object instead)
  // ideally these should be const, but I have some trouble initialising them in that case
  CartesianCoordinate3D<float> voxel_size;
  CartesianCoordinate3D<float> origin;
  CartesianCoordinate3D<int> min_index;
  CartesianCoordinate3D<int> max_index;

  float xhalfsize;
  float yhalfsize;
  float zhalfsize;

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

END_NAMESPACE_STIR

#endif