rtkTotalVariationDenoisingBPDQImageFilter.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 rtkTotalVariationDenoisingBPDQImageFilter_h
#define rtkTotalVariationDenoisingBPDQImageFilter_h

#include "rtkDenoisingBPDQImageFilter.h"
#include "rtkMagnitudeThresholdImageFilter.h"

#include <itkPeriodicBoundaryCondition.h>

namespace rtk
{
template <typename TOutputImage,
          typename TGradientImage =
            itk::Image<itk::CovariantVector<typename TOutputImage::ValueType, TOutputImage::ImageDimension>,
                       TOutputImage::ImageDimension>>
class ITK_TEMPLATE_EXPORT TotalVariationDenoisingBPDQImageFilter
  : public rtk::DenoisingBPDQImageFilter<TOutputImage, TGradientImage>
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(TotalVariationDenoisingBPDQImageFilter);

  using Self = TotalVariationDenoisingBPDQImageFilter;
  using Superclass = rtk::DenoisingBPDQImageFilter<TOutputImage, TGradientImage>;
  using Pointer = itk::SmartPointer<Self>;
  using ConstPointer = itk::SmartPointer<const Self>;

  itkNewMacro(Self);

#ifdef itkOverrideGetNameOfClassMacro
  itkOverrideGetNameOfClassMacro(TotalVariationDenoisingBPDQImageFilter);
#else
  itkTypeMacro(TotalVariationDenoisingBPDQImageFilter, DenoisingBPDQImageFilter);
#endif


  using MagnitudeThresholdFilterType =
    MagnitudeThresholdImageFilter<TGradientImage, typename TOutputImage::ValueType, TGradientImage>;

  void
  SetDimensionsProcessed(bool * arg);

  void
  SetBoundaryConditionToPeriodic();

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

  void
  GenerateOutputInformation() override;

  typename MagnitudeThresholdFilterType::Pointer m_ThresholdFilter;
  typename Superclass::ThresholdFilterType *
  GetThresholdFilter() override
  {
    return dynamic_cast<typename Superclass::ThresholdFilterType *>(this->m_ThresholdFilter.GetPointer());
  }
};

} // namespace rtk

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

#endif //__rtkTotalVariationDenoisingBPDQImageFilter__