rtkFDKVarianceReconstructionFilter.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 rtkFDKVarianceReconstructionFilter_h
#define rtkFDKVarianceReconstructionFilter_h

#include "rtkFDKWeightProjectionFilter.h"
#include "rtkFFTVarianceRampImageFilter.h"
#include "rtkFDKBackProjectionImageFilter.h"
#include "rtkConfiguration.h"

#include <itkExtractImageFilter.h>

namespace rtk
{

template <class TInputImage, class TOutputImage = TInputImage, class TFFTPrecision = double>
class ITK_TEMPLATE_EXPORT FDKVarianceReconstructionFilter : public itk::InPlaceImageFilter<TInputImage, TOutputImage>
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(FDKVarianceReconstructionFilter);

  using Self = FDKVarianceReconstructionFilter;
  using Superclass = itk::InPlaceImageFilter<TInputImage, TOutputImage>;
  using Pointer = itk::SmartPointer<Self>;
  using ConstPointer = itk::SmartPointer<const Self>;

  using InputImageType = TInputImage;
  using OutputImageType = TOutputImage;

  using ExtractFilterType = itk::ExtractImageFilter<InputImageType, OutputImageType>;
  using WeightFilterType = rtk::FDKWeightProjectionFilter<InputImageType, OutputImageType>;
  using VarianceRampFilterType = rtk::FFTVarianceRampImageFilter<OutputImageType, OutputImageType, TFFTPrecision>;
  using BackProjectionFilterType = rtk::FDKBackProjectionImageFilter<OutputImageType, OutputImageType>;
  using BackProjectionFilterPointer = typename BackProjectionFilterType::Pointer;

  itkNewMacro(Self);

  itkOverrideGetNameOfClassMacro(FDKVarianceReconstructionFilter);

  itkGetModifiableObjectMacro(Geometry, ThreeDCircularProjectionGeometry);
  itkSetObjectMacro(Geometry, ThreeDCircularProjectionGeometry);

  typename VarianceRampFilterType::Pointer
  GetVarianceRampFilter()
  {
    return m_VarianceRampFilter;
  }

  itkGetMacro(ProjectionSubsetSize, unsigned int);
  itkSetMacro(ProjectionSubsetSize, unsigned int);

  itkGetMacro(BackProjectionFilter, BackProjectionFilterPointer);
  virtual void
  SetBackProjectionFilter(const BackProjectionFilterPointer _arg);

protected:
  FDKVarianceReconstructionFilter();
  ~FDKVarianceReconstructionFilter() override = default;

  void
  VerifyPreconditions() const override;

  void
  GenerateInputRequestedRegion() override;

  void
  GenerateOutputInformation() override;

  void
  GenerateData() override;

  void
  VerifyInputInformation() const override
  {}

  typename ExtractFilterType::Pointer      m_ExtractFilter;
  typename WeightFilterType::Pointer       m_WeightFilter1;
  typename WeightFilterType::Pointer       m_WeightFilter2;
  typename VarianceRampFilterType::Pointer m_VarianceRampFilter;
  BackProjectionFilterPointer              m_BackProjectionFilter;

private:
  unsigned int m_ProjectionSubsetSize{ 16 };

  ThreeDCircularProjectionGeometry::Pointer m_Geometry;
}; // end of class

} // end namespace rtk

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

#endif