rtkProjectionsRegionConstIteratorRayBased.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright RTK Consortium
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         https://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/

#ifndef rtkProjectionsRegionConstIteratorRayBased_h
#define rtkProjectionsRegionConstIteratorRayBased_h

#include <itkImageConstIteratorWithIndex.h>

#include "rtkThreeDCircularProjectionGeometry.h"

namespace rtk
{
template <typename TImage>
class ProjectionsRegionConstIteratorRayBasedWithFlatPanel;
template <typename TImage>
class ProjectionsRegionConstIteratorRayBasedWithCylindricalPanel;
template <typename TImage>
class ProjectionsRegionConstIteratorRayBasedParallel;

template <typename TImage>
class ITK_TEMPLATE_EXPORT ProjectionsRegionConstIteratorRayBased : public itk::ImageConstIteratorWithIndex<TImage>
{
public:
  using Self = ProjectionsRegionConstIteratorRayBased;
  using Superclass = itk::ImageConstIteratorWithIndex<TImage>;

  using OffsetValueType = typename Superclass::OffsetValueType;
  using RegionType = typename Superclass::RegionType;
  using PointType = typename itk::Vector<double, 3>;
  using IndexValueType = typename Superclass::IndexValueType;

  using MatrixType = itk::Matrix<double, 3, 4>;
  using HomogeneousMatrixType = itk::Matrix<double, 4, 4>;

  ProjectionsRegionConstIteratorRayBased(const TImage *                           ptr,
                                         const RegionType &                       region,
                                         const ThreeDCircularProjectionGeometry * geometry,
                                         const MatrixType &                       postMat);

  static Self *
  New(const TImage *                           ptr,
      const RegionType &                       region,
      const ThreeDCircularProjectionGeometry * geometry,
      const MatrixType &                       postMat);

  static Self *
  New(const TImage *                           ptr,
      const RegionType &                       region,
      const ThreeDCircularProjectionGeometry * geometry,
      const HomogeneousMatrixType &            postMat);

  static Self *
  New(const TImage * ptr, const RegionType & region, const ThreeDCircularProjectionGeometry * geometry);

  Self &
  operator++();

  void
  Next()
  {
    ++*this;
  }

  const PointType &
  GetSourcePosition()
  {
    return this->m_SourcePosition;
  }
  const PointType &
  GetPixelPosition()
  {
    return this->m_PixelPosition;
  }
  const PointType &
  GetSourceToPixel()
  {
    return this->m_SourceToPixel;
  }

  const PointType
  GetDirection()
  {
    return m_SourceToPixel / m_SourceToPixel.GetNorm();
  }

protected:
  virtual void
  NewProjection() = 0;

  virtual void
  NewPixel() = 0;

  ThreeDCircularProjectionGeometry::ConstPointer m_Geometry;
  MatrixType                                     m_PostMultiplyMatrix;
  PointType                                      m_SourcePosition;
  PointType                                      m_PixelPosition;
  PointType                                      m_SourceToPixel;
};
} // namespace rtk

#ifndef ITK_MANUAL_INSTANTIATION
#  include "rtkProjectionsRegionConstIteratorRayBased.hxx"
#endif

#endif