DataSymmetriesForDensels_PET_CartesianGrid.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2001- 2011,  Hammersmith Imanet Ltd
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/
#ifndef __stir_recon_buildblock_DataSymmetriesForDensels_PET_CartesianGrid_H__
#define __stir_recon_buildblock_DataSymmetriesForDensels_PET_CartesianGrid_H__
 
#include "stir/recon_buildblock/DataSymmetriesForDensels.h"
#include "stir/ProjDataInfo.h"
//#include "stir/SymmetryOperations_PET_CartesianGrid.h"
//#include "stir/ViewSegmentNumbers.h"
//#include "stir/VoxelsOnCartesianGrid.h"
#include "stir/Densel.h"
#include "stir/shared_ptr.h"
#include <vector>
 
START_NAMESPACE_STIR
 
template <int num_dimensions, typename elemT>
class DiscretisedDensity;
template <int num_dimensions, typename elemT>
class DiscretisedDensityOnCartesianGrid;

class DataSymmetriesForDensels_PET_CartesianGrid : public DataSymmetriesForDensels
{
private:
  typedef DataSymmetriesForDensels base_type;
  typedef DataSymmetriesForDensels_PET_CartesianGrid self_type;
 
public:
  DataSymmetriesForDensels_PET_CartesianGrid(const shared_ptr<const ProjDataInfo>& proj_data_info_ptr,
                                             const shared_ptr<const DiscretisedDensity<3, float>>& image_info_ptr);
 
  DataSymmetriesForDensels_PET_CartesianGrid* clone() const override;
 
  bool operator==(const DataSymmetriesForDensels_PET_CartesianGrid&) const;
 
#if 0
  TODO!
  virtual DenselIndexRange
    get_basic_densel_index_range() const;
#endif
 
  inline void get_related_densels(std::vector<Densel>&, const Densel& b) const override;
 
  inline int num_related_densels(const Densel& b) const override;
 
  inline unique_ptr<SymmetryOperation> find_symmetry_operation_from_basic_densel(Densel&) const override;
 
  inline bool find_basic_densel(Densel& b) 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;
 
private:
  const shared_ptr<const ProjDataInfo>& proj_data_info_ptr;
  int num_planes;
  int num_independent_planes;
  int num_views;
  int num_planes_per_scanner_ring;
  VectorWithOffset<int> num_planes_per_axial_pos;
  VectorWithOffset<float> axial_pos_to_z_offset;
 
#if 0
  // at the moment, we don't need the following 2 members
 
  // TODO somehow store only the info
  shared_ptr<const DiscretisedDensity<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 SymmetryOperation* find_sym_op_general_densel(const int z, const int y, const int x) const;
};
 
END_NAMESPACE_STIR
 
#include "stir_experimental/recon_buildblock/DataSymmetriesForDensels_PET_CartesianGrid.inl"
 
#endif