DataSymmetriesForBins_PET_CartesianGrid.h

Go to the documentation of this file.

//
//
/*
    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_DataSymmetriesForBins_PET_CartesianGrid_H__
#define __stir_recon_buildblock_DataSymmetriesForBins_PET_CartesianGrid_H__

#include "stir/recon_buildblock/DataSymmetriesForBins.h"
//#include "stir/SymmetryOperations_PET_CartesianGrid.h"
//#include "stir/ViewSegmentNumbers.h"
//#include "stir/VoxelsOnCartesianGrid.h"
#include "stir/Bin.h"
#include "stir/shared_ptr.h"

START_NAMESPACE_STIR

template <int num_dimensions, typename elemT>
class DiscretisedDensity;
template <int num_dimensions, typename elemT>
class DiscretisedDensityOnCartesianGrid;
class ProjDataInfoCylindrical;

class DataSymmetriesForBins_PET_CartesianGrid : public DataSymmetriesForBins
{
private:
  typedef DataSymmetriesForBins base_type;
  typedef DataSymmetriesForBins_PET_CartesianGrid self_type;

public:

  DataSymmetriesForBins_PET_CartesianGrid(const shared_ptr<const ProjDataInfo>& proj_data_info_ptr,
                                          const shared_ptr<const DiscretisedDensity<3, float>>& image_info_ptr,
                                          const bool do_symmetry_90degrees_min_phi = true,
                                          const bool do_symmetry_180degrees_min_phi = true,
                                          const bool do_symmetry_swap_segment = true,
                                          const bool do_symmetry_swap_s = true,
                                          const bool do_symmetry_shift_z = true);

#ifndef STIR_NO_COVARIANT_RETURN_TYPES
  DataSymmetriesForBins_PET_CartesianGrid
#else
  DataSymmetriesForViewSegmentNumbers
#endif
      *
      clone() const override;

  virtual bool operator==(const DataSymmetriesForBins_PET_CartesianGrid&) const;

#if 0
  TODO!
  virtual BinIndexRange
    get_basic_bin_index_range() const;
#endif

  inline void get_related_bins_factorised(std::vector<AxTangPosNumbers>&,
                                          const Bin& b,
                                          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 override;

  inline int num_related_bins(const Bin& b) const override;

  inline unique_ptr<SymmetryOperation> find_symmetry_operation_from_basic_bin(Bin&) const override;

  inline bool find_basic_bin(Bin& b) const override;

  inline int num_related_view_segment_numbers(const ViewSegmentNumbers& vs) const override;

  inline void get_related_view_segment_numbers(std::vector<ViewSegmentNumbers>& rel_vs,
                                               const ViewSegmentNumbers& vs) const override;

  inline bool find_basic_view_segment_numbers(ViewSegmentNumbers& v_s) const override;

  inline float get_num_planes_per_scanner_ring() const;


  inline float get_num_planes_per_axial_pos(const int segment_num) const;
  inline float get_axial_pos_to_z_offset(const int segment_num) const;


  inline bool using_symmetry_90degrees_min_phi() const
  {
    return do_symmetry_90degrees_min_phi;
  }
  inline bool using_symmetry_180degrees_min_phi() const
  {
    return do_symmetry_180degrees_min_phi;
  }
  inline bool using_symmetry_swap_segment() const
  {
    return do_symmetry_swap_segment;
  }
  inline bool using_symmetry_swap_s() const
  {
    return do_symmetry_swap_s;
  }
  inline bool using_symmetry_shift_z() const
  {
    return do_symmetry_shift_z;
  }

private:
  bool do_symmetry_90degrees_min_phi;
  bool do_symmetry_180degrees_min_phi;
  bool do_symmetry_swap_segment;
  bool do_symmetry_swap_s;
  bool do_symmetry_shift_z;
  // const shared_ptr<ProjDataInfo>& proj_data_info_ptr;
  int num_views;
  int num_planes_per_scanner_ring;
  VectorWithOffset<int> num_planes_per_axial_pos;
  VectorWithOffset<float> axial_pos_to_z_offset;
  VectorWithOffset<float> deltas;

  void initialise_deltas(const ProjDataInfoCylindrical*);

#if 0
  // at the moment, we don't need the following 2 members

  // TODO somehow store only the info
  shared_ptr<constDiscretisedDensity<3,float> > image_info_ptr;

  // a convenience function that does the dynamic_cast from the above
  inline const DiscretisedDensityOnCartesianGrid<3,float> *
    cartesian_grid_info_ptr() const;
#endif

  bool blindly_equals(const root_type* const) const override;

  inline bool
  find_basic_bin(int& segment_num, int& view_num, int& axial_pos_num, int& tangential_pos_num, int& timing_pos_num) const;

  inline int find_transform_z(const int segment_num, const int axial_pos_num) const;

  inline SymmetryOperation* find_sym_op_general_bin(int s, int seg, int view_num, int axial_pos_num) const;

  inline SymmetryOperation* find_sym_op_bin0(int seg, int view_num, int axial_pos_num) const;
};

END_NAMESPACE_STIR

#include "stir/recon_buildblock/DataSymmetriesForBins_PET_CartesianGrid.inl"

#endif