FourierRebinning.h

Go to the documentation of this file.

//
//

/*
    Copyright (C) 2000 PARAPET partners
    Copyright (C) 2003 - 2005, Hammersmith Imanet Ltd
    Copyright (C) 2004 - 2005 DKFZ Heidelberg, Germany
    Copyright (C) 2011-07-01 - 2012, Kris Thielemans

    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_FourierRebinning_H__
#define __stir_recon_buildblock_FourierRebinning_H__

#include "stir/recon_buildblock/ProjDataRebinning.h"
#include "stir/RegisteredParsingObject.h"
#include "stir/ArrayFwd.h"
#include <complex>

START_NAMESPACE_STIR
template <typename elemT>
class SegmentByView;
template <typename elemT>
class SegmentBySinogram;
// template <typename elemT> class Sinogram;
class Succeeded;

/*
  \class PETCount_rebinned
  \brief Class for rebinned elements counter
  \ingroup recon_buildblock
*/
class PETCount_rebinned
{

public:
  int total;
  int miss;
  int ssrb;

#ifdef PARALLEL
  friend PMessage& operator<<(PMessage&, PETCount_rebinned&);
  friend PMessage& operator>>(PMessage&, PETCount_rebinned&);

  PETCount_rebinned& operator+=(const PETCount_rebinned& rebin)
  {
    total += rebin.total;
    miss += rebin.miss;
    ssrb += rebin.ssrb;
    return *this;
  }
#endif
  // Default constructor by initialising all the elements conter to null
  explicit PETCount_rebinned(int total_v = 0, int miss_v = 0, int ssrb_v = 0)
      : total(total_v),
        miss(miss_v),
        ssrb(ssrb_v){}

  ;
};

class FourierRebinning : public RegisteredParsingObject<FourierRebinning, ProjDataRebinning, ProjDataRebinning>
{
private:
  typedef ProjDataRebinning base_type;

public:
  static const char* const registered_name;
  void set_defaults() override;

protected:
  int kmin;
  int wmin;
  int deltamin;
  int kc;
  int fore_debug_level;

public:
  FourierRebinning();

  std::string method_info() const override { return ("FORE"); }

  Succeeded rebin() override;

  inline void set_kmin(int km) { kmin = km; }
  inline void set_wmin(int wm) { wmin = wm; }
  inline void set_deltamin(int dm) { deltamin = dm; }
  inline void set_kc(int kcc) { kc = kcc; }
  inline void set_fore_debug_level(int fdebug) { fore_debug_level = fdebug; }

  inline int get_kmin() { return kmin; }
  inline int get_wmin() { return wmin; }
  inline int get_deltamin() { return deltamin; }
  inline int get_kc() { return kc; }
  inline int get_fore_debug_level() { return fore_debug_level; }

private:
  void rebinning(ArrayType<3, std::complex<float>>& FT_rebinned_data,
                 ArrayType<3, float>& Weights_for_FT_rebinned_data,
                 PETCount_rebinned& num_rebinned,
                 const ArrayType<2, std::complex<float>>& FT_current_sinogram,
                 const float z,
                 const float average_ring_difference_in_segment,
                 const int num_views_pow2,
                 const int num_tang_poss_pow2,
                 const float half_distance_between_rings,
                 const float sampling_distance_in_s,
                 const float radial_sampling_freq_w,
                 const float R_field_of_view_mm,
                 const float ratio_ring_spacing_to_ring_radius);

  void do_rebinning(ArrayType<3, std::complex<float>>& FT_rebinned_data,
                    ArrayType<3, float>& Weights_for_FT_rebinned_data,
                    PETCount_rebinned& count_rebinned,
                    const SegmentBySinogram<float>& segment,
                    const int num_tang_poss_pow2,
                    const int num_views_pow2,
                    const int num_planes,
                    const float average_ring_difference_in_segment,
                    const float half_distance_between_rings,
                    const float sampling_distance_in_s,
                    const float radial_sampling_freq_w,
                    const float R_field_of_view_mm,
                    const float ratio_ring_spacing_to_ring_radius);

  void do_log_file();

  void do_display_count(PETCount_rebinned& num_rebinned_total);

  void do_adjust_nb_views_to_pow2(SegmentBySinogram<float>& segment);

  Succeeded fore_check_parameters(int num_tang_poss_pow2, int num_views_pow2, int max_segment_num_to_process);

protected:
  bool post_processing() override;
  void initialise_keymap() override;
};

END_NAMESPACE_STIR
#endif