PostsmoothingForwardProjectorByBin.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2000- 2001, IRSL
    Copyright (C) 2019, University College London
    See STIR/LICENSE.txt for details
*/
#ifndef __stir_recon_buildblock_PostsmoothingForwardProjectorByBin__H__
#define __stir_recon_buildblock_PostsmoothingForwardProjectorByBin__H__

#include "stir/RegisteredParsingObject.h"
#include "stir/recon_buildblock/ForwardProjectorByBin.h"
#include "stir/shared_ptr.h"
#include "stir/VectorWithOffset.h"

START_NAMESPACE_STIR

template <typename elemT>
class Viewgram;

class PostsmoothingForwardProjectorByBin
    : public RegisteredParsingObject<PostsmoothingForwardProjectorByBin, ForwardProjectorByBin>
{
public:
  static const char* const registered_name;

  PostsmoothingForwardProjectorByBin();


  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;

  PostsmoothingForwardProjectorByBin(const shared_ptr<ForwardProjectorByBin>& original_forward_projector_ptr,
                                     const VectorWithOffset<float>& tangential_kernel,
                                     const VectorWithOffset<float>& axial_kernel,
                                     const bool smooth_segment_0_axially = false);

  // Informs on which symmetries the projector handles
  // It should get data related by at least those symmetries.
  // Otherwise, a run-time error will occur (unless the derived
  // class has other behaviour).
  const DataSymmetriesForViewSegmentNumbers* get_symmetries_used() const override;

private:
  shared_ptr<ForwardProjectorByBin> original_forward_projector_ptr;
  VectorWithOffset<float> tang_kernel;
  VectorWithOffset<float> ax_kernel;
  bool smooth_segment_0_axially;

  // next 2 necessary for parsing because of limitation in KeyParser
  std::vector<double> tang_kernel_double;
  std::vector<double> ax_kernel_double;

#ifdef STIR_PROJECTORS_AS_V3
  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);
#endif
  void actual_forward_project(RelatedViewgrams<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;
  void smooth(Viewgram<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) const;

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

END_NAMESPACE_STIR

#endif