ProjDataInfoGenericNoArcCorr.inl

Go to the documentation of this file.

/*
    Copyright (C) 2021, 2026 University College London
    Copyright (C) 2017, ETH Zurich, Institute of Particle Physics and Astrophysics
    Copyright (C) 2018, 2021, University of Leeds
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/

 
// for sqrt
#include <cmath>
#include "stir/Bin.h"
#include "stir/Succeeded.h"
#include "stir/LORCoordinates.h"
#include "stir/is_null_ptr.h"
#include <algorithm>
#include "stir/CartesianCoordinate3D.h"
#include "stir/error.h"
 
START_NAMESPACE_STIR

float
ProjDataInfoGenericNoArcCorr::get_phi(const Bin& bin) const
{
  LORInAxialAndNoArcCorrSinogramCoordinates<float> lor;
  get_LOR(lor, bin);
  return lor.phi();
}

float
ProjDataInfoGenericNoArcCorr::get_m(const Bin& bin) const
{
  LORInAxialAndNoArcCorrSinogramCoordinates<float> lor;
  get_LOR(lor, bin);
  return (static_cast<float>(lor.z1() + lor.z2())) / 2.F;
}
 
float
ProjDataInfoGenericNoArcCorr::get_t(const Bin& bin) const
{
  return get_m(bin) * get_costheta(bin);
}
 
/*
        theta is copolar angle of normal to projection plane with z axis, i.e. copolar angle of lor with z axis.
        tan (theta) = dz/sqrt(dx2+dy2)
        cylindrical geometry:
                delta_z = delta_ring * ring spacing
        generic geometry:
                delta_z is calculated from lor
*/
float
ProjDataInfoGenericNoArcCorr::get_tantheta(const Bin& bin) const
{
  LORInAxialAndNoArcCorrSinogramCoordinates<float> lor;
  get_LOR(lor, bin);
 
  return (lor.z2() - lor.z1()) / (2 * lor.radius());
}
 
float
ProjDataInfoGenericNoArcCorr::get_sampling_in_m(const Bin& bin) const
{
  // TODOBLOCK
  return get_scanner_ptr()
      ->get_ring_spacing(); // TODO currently restricted to span=1 get_num_axial_poss_per_ring_inc(segment_num);
  // return get_axial_sampling(bin.segment_num());
}
 
float
ProjDataInfoGenericNoArcCorr::get_sampling_in_t(const Bin& bin) const
{
  // TODOBLOCK
  return get_scanner_ptr()->get_ring_spacing()
         * get_costheta(bin); // TODO currently restricted to span=1 get_num_axial_poss_per_ring_inc(segment_num);
  // return get_axial_sampling(bin.segment_num())*get_costheta(bin);
}
 
float
ProjDataInfoGenericNoArcCorr::get_axial_sampling(int segment_num) const
{
  // TODOBLOCK should check sampling
  return get_ring_spacing(); // TODO currently restricted to span=1 get_num_axial_poss_per_ring_inc(segment_num);
}
 
bool
ProjDataInfoGenericNoArcCorr::axial_sampling_is_uniform() const
{
  // TODOBLOCK should check sampling
  return true;
}
 
float
ProjDataInfoGenericNoArcCorr::get_ring_spacing() const
{
  return get_scanner_ptr()->get_ring_spacing();
}

float
ProjDataInfoGenericNoArcCorr::get_s(const Bin& bin) const
{
  LORInAxialAndNoArcCorrSinogramCoordinates<float> lor;
  get_LOR(lor, bin);
  return lor.s();
}
 
END_NAMESPACE_STIR