doxy/html/ForwardProjectorByBinUsingRayTracing_8h_source.html

//

//

/*

    Copyright (C) 2000 PARAPET partners

    Copyright (C) 2000- 2011, Hammersmith Imanet Ltd

    This file is part of STIR.


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


    See STIR/LICENSE.txt for details

*/


#ifndef __stir_recon_buildblock_ForwardProjectorByBinUsingRayTracing__H__

#define __stir_recon_buildblock_ForwardProjectorByBinUsingRayTracing__H__


#include "stir/recon_buildblock/ForwardProjectorByBin.h"

#include "stir/RegisteredParsingObject.h"

#include "stir/shared_ptr.h"

#include "stir/recon_buildblock/DataSymmetriesForBins_PET_CartesianGrid.h"

#include "stir/ArrayFwd.h"


START_NAMESPACE_STIR


template <typename elemT>

class Viewgram;

template <typename elemT>

class RelatedViewgrams;

template <typename elemT>

class VoxelsOnCartesianGrid;

class ProjDataInfo;

class ProjDataInfoCylindrical;


class ForwardProjectorByBinUsingRayTracing

    : public RegisteredParsingObject<ForwardProjectorByBinUsingRayTracing, ForwardProjectorByBin>

{

public:

  static const char* const registered_name;


  ForwardProjectorByBinUsingRayTracing();


  ForwardProjectorByBinUsingRayTracing(const shared_ptr<const ProjDataInfo>&,

                                       const shared_ptr<const DiscretisedDensity<3, float>>&);


  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 DataSymmetriesForViewSegmentNumbers* get_symmetries_used() const override;


protected:

  bool restrict_to_cylindrical_FOV;


private:

  void actual_forward_project(RelatedViewgrams<float>&,

                              const DiscretisedDensity<3, float>&,

                              const int min_axial_pos_num,

                              const int max_axial_pos_num,

                              const int min_tangential_pos_num,

                              const int max_tangential_pos_num) override;

#if 0 // disabled as currently not used. needs to be written in the new style anyway

  void actual_forward_project(Bin&,

                              const DiscretisedDensity<3,float>&);

#endif


  // KT 20/06/2001 changed type from 'const DataSymmetriesForViewSegmentNumbers *'

  shared_ptr<DataSymmetriesForBins_PET_CartesianGrid> symmetries_ptr;

  /*

    The version which uses all possible symmetries.

    Here 0<=view < num_views/4 (= 45 degrees)

    */


  void forward_project_all_symmetries(Viewgram<float>& pos_view,

                                      Viewgram<float>& neg_view,

                                      Viewgram<float>& pos_plus90,

                                      Viewgram<float>& neg_plus90,

                                      Viewgram<float>& pos_min180,

                                      Viewgram<float>& neg_min180,

                                      Viewgram<float>& pos_min90,

                                      Viewgram<float>& neg_min90,

                                      const VoxelsOnCartesianGrid<float>& image,

                                      const int min_axial_pos_num,

                                      const int max_axial_pos_num,

                                      const int min_tangential_pos_num,

                                      const int max_tangential_pos_num) const;


  /*

    This function projects 4 viewgrams related by symmetry.

    It will be used for view=0 or 45 degrees

    (or others if the number of views is not a multiple of 4)

    Here 0<=view < num_views/2 (= 90 degrees)

    */

  void forward_project_view_plus_90_and_delta(Viewgram<float>& pos_view,

                                              Viewgram<float>& neg_view,

                                              Viewgram<float>& pos_plus90,

                                              Viewgram<float>& neg_plus90,

                                              const VoxelsOnCartesianGrid<float>& image,

                                              const int min_axial_pos_num,

                                              const int max_axial_pos_num,

                                              const int min_tangential_pos_num,

                                              const int max_tangential_pos_num) const;

  /*

    This function projects 4 viewgrams related by symmetry.

    It will be used for view=0 or 45 degrees

    (or others if the number of views is not a multiple of 4)

    Here 0<=view < num_views/2 (= 90 degrees)

    */

  void forward_project_view_min_180_and_delta(Viewgram<float>& pos_view,

                                              Viewgram<float>& neg_view,

                                              Viewgram<float>& pos_min180,

                                              Viewgram<float>& neg_min180,

                                              const VoxelsOnCartesianGrid<float>& image,

                                              const int min_axial_pos_num,

                                              const int max_axial_pos_num,

                                              const int min_tangential_pos_num,

                                              const int max_tangential_pos_num) const;


  /*

    This function projects 4 viewgrams related by symmetry.

    It will be used for view=0 or 45 degrees

    (or others if the number of views is not a multiple of 4)

    Here 0<=view < num_views/2 (= 90 degrees)

    */

  void forward_project_delta(Viewgram<float>& pos_view,

                             Viewgram<float>& neg_view,

                             const VoxelsOnCartesianGrid<float>& image,

                             const int min_axial_pos_num,

                             const int max_axial_pos_num,

                             const int min_tangential_pos_num,

                             const int max_tangential_pos_num) const;


  void forward_project_all_symmetries_2D(Viewgram<float>& pos_view,

                                         Viewgram<float>& pos_plus90,

                                         Viewgram<float>& pos_min180,

                                         Viewgram<float>& pos_min90,

                                         const VoxelsOnCartesianGrid<float>& image,

                                         const int min_axial_pos_num,

                                         const int max_axial_pos_num,

                                         const int min_tangential_pos_num,

                                         const int max_tangential_pos_num) const;

  void forward_project_view_plus_90_2D(Viewgram<float>& pos_view,

                                       Viewgram<float>& pos_plus90,

                                       const VoxelsOnCartesianGrid<float>& image,

                                       const int min_axial_pos_num,

                                       const int max_axial_pos_num,

                                       const int min_tangential_pos_num,

                                       const int max_tangential_pos_num) const;

  void forward_project_view_min_180_2D(Viewgram<float>& pos_view,

                                       Viewgram<float>& pos_min180,

                                       const VoxelsOnCartesianGrid<float>& image,

                                       const int min_axial_pos_num,

                                       const int max_axial_pos_num,

                                       const int min_tangential_pos_num,

                                       const int max_tangential_pos_num) const;

  // no symmetries

  void forward_project_view_2D(Viewgram<float>& pos_view,

                               const VoxelsOnCartesianGrid<float>& image,

                               const int min_axial_pos_num,

                               const int max_axial_pos_num,

                               const int min_tangential_pos_num,

                               const int max_tangential_pos_num) const;


  template <int symmetry_type>

  static bool proj_Siddon(Array<4, float>& Projptr,

                          const VoxelsOnCartesianGrid<float>&,

                          const shared_ptr<const ProjDataInfoCylindrical> proj_data_info_sptr,

                          const float cphi,

                          const float sphi,

                          const float delta,

                          const float s_in_mm,

                          const float R,

                          const int min_ax_pos_num,

                          const int max_ax_pos_num,

                          const float offset,

                          const int num_planes_per_axial_pos,

                          const float axial_pos_to_z_offset,

                          const float norm_factor,

                          const bool restrict_to_cylindrical_FOV);


  void set_defaults() override;

  void initialise_keymap() override;

};


END_NAMESPACE_STIR

#endif

ArrayFwd.h
forward declaration of stir::Array class for multi-dimensional (numeric) arrays

DataSymmetriesForBins_PET_CartesianGrid.h
Declaration of class stir::DataSymmetriesForBins_PET_CartesianGrid.

ForwardProjectorByBin.h
Base class for forward projectors which work on 'large' collections of bins: given the whole image,...

RegisteredParsingObject.h
Declaration of class stir::RegisteredParsingObject.

stir::Array
This class defines multi-dimensional (numeric) arrays.
Definition Array.h:78

stir::Bin
A class for storing coordinates and value of a single projection bin.
Definition Bin.h:49

stir::DataSymmetriesForViewSegmentNumbers
A class for encoding/finding symmetries. Works only on ViewSegmentNumbers (instead of Bin).
Definition DataSymmetriesForViewSegmentNumbers.h:50

stir::DiscretisedDensity
This abstract class is the basis for all image representations.
Definition DiscretisedDensity.h:99

stir::ForwardProjectorByBinUsingRayTracing::registered_name
static const char *const registered_name
Name which will be used when parsing a ForwardProjectorByBin object.
Definition ForwardProjectorByBinUsingRayTracing.h:78

stir::ForwardProjectorByBinUsingRayTracing::restrict_to_cylindrical_FOV
bool restrict_to_cylindrical_FOV
variable that determines if a cylindrical FOV or the whole image will be handled
Definition ForwardProjectorByBinUsingRayTracing.h:97

stir::ForwardProjectorByBinUsingRayTracing::set_up
void set_up(const shared_ptr< const ProjDataInfo > &proj_data_info_ptr, const shared_ptr< const DiscretisedDensity< 3, float > > &density_info_ptr) override
Stores all necessary geometric info.
Definition ForwardProjectorByBinUsingRayTracing.cxx:103

stir::ForwardProjectorByBinUsingRayTracing::get_symmetries_used
const DataSymmetriesForViewSegmentNumbers * get_symmetries_used() const override
Informs on which symmetries the projector handles.
Definition ForwardProjectorByBinUsingRayTracing.cxx:164

stir::RegisteredParsingObject
Parent class for all leaves in a RegisteredObject hierarchy that do parsing of parameter files.
Definition RegisteredParsingObject.h:78

stir::RelatedViewgrams
A class for storing viewgrams which are related by symmetry.
Definition RelatedViewgrams.h:41

stir::Viewgram
A class for 2d projection data.
Definition Viewgram.h:49

stir::VoxelsOnCartesianGrid
This class is used to represent voxelised densities on a cuboid grid (3D).
Definition VoxelsOnCartesianGrid.h:46

shared_ptr.h
Import of std::shared_ptr, std::dynamic_pointer_cast and std::static_pointer_cast into the stir names...