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