itkGPUBinaryThresholdImageFilter.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright NumFOCUS
 *
 *  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 itkGPUBinaryThresholdImageFilter_h
#define itkGPUBinaryThresholdImageFilter_h
 
#include "itkOpenCLUtil.h"
#include "itkGPUFunctorBase.h"
#include "itkGPUKernelManager.h"
#include "itkBinaryThresholdImageFilter.h"
#include "itkGPUUnaryFunctorImageFilter.h"
 
namespace itk
{
 
namespace Functor
{
template <typename TInput, typename TOutput>
class ITK_TEMPLATE_EXPORT GPUBinaryThreshold : public GPUFunctorBase
{
public:
  GPUBinaryThreshold()
  {
    m_LowerThreshold = NumericTraits<TInput>::NonpositiveMin();
    m_UpperThreshold = NumericTraits<TInput>::max();
    m_OutsideValue = TOutput{};
    m_InsideValue = NumericTraits<TOutput>::max();
  }
 
  ~GPUBinaryThreshold() override = default;
 
  void
  SetLowerThreshold(const TInput & thresh)
  {
    m_LowerThreshold = thresh;
  }
  void
  SetUpperThreshold(const TInput & thresh)
  {
    m_UpperThreshold = thresh;
  }
  void
  SetInsideValue(const TOutput & value)
  {
    m_InsideValue = value;
  }
  void
  SetOutsideValue(const TOutput & value)
  {
    m_OutsideValue = value;
  }
 
  int
  SetGPUKernelArguments(GPUKernelManager::Pointer KernelManager, int KernelHandle) override
  {
    KernelManager->SetKernelArg(KernelHandle, 0, sizeof(TInput), &(m_LowerThreshold));
    KernelManager->SetKernelArg(KernelHandle, 1, sizeof(TInput), &(m_UpperThreshold));
    KernelManager->SetKernelArg(KernelHandle, 2, sizeof(TOutput), &(m_InsideValue));
    KernelManager->SetKernelArg(KernelHandle, 3, sizeof(TOutput), &(m_OutsideValue));
    return 4;
  }
private:
  TInput  m_LowerThreshold{};
  TInput  m_UpperThreshold{};
  TOutput m_InsideValue{};
  TOutput m_OutsideValue{};
};
} // end of namespace Functor
 
itkGPUKernelClassMacro(GPUBinaryThresholdImageFilterKernel);
 
template <typename TInputImage, typename TOutputImage>
class ITK_TEMPLATE_EXPORT GPUBinaryThresholdImageFilter
  : public GPUUnaryFunctorImageFilter<
      TInputImage,
      TOutputImage,
      Functor::GPUBinaryThreshold<typename TInputImage::PixelType, typename TOutputImage::PixelType>,
      BinaryThresholdImageFilter<TInputImage, TOutputImage>>
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(GPUBinaryThresholdImageFilter);
 
  using Self = GPUBinaryThresholdImageFilter;
  using GPUSuperclass = GPUUnaryFunctorImageFilter<
    TInputImage,
    TOutputImage,
    Functor::GPUBinaryThreshold<typename TInputImage::PixelType, typename TOutputImage::PixelType>,
    BinaryThresholdImageFilter<TInputImage, TOutputImage>>;
  using CPUSuperclass = BinaryThresholdImageFilter<TInputImage, TOutputImage>;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;
 
  itkNewMacro(Self);
 
  itkOverrideGetNameOfClassMacro(GPUBinaryThresholdImageFilter);
 
  using InputPixelType = typename TInputImage::PixelType;
  using OutputPixelType = typename TOutputImage::PixelType;
 
  using InputPixelObjectType = SimpleDataObjectDecorator<InputPixelType>;
 
  itkGetOpenCLSourceFromKernelMacro(GPUBinaryThresholdImageFilterKernel);
 
protected:
  GPUBinaryThresholdImageFilter();
  ~GPUBinaryThresholdImageFilter() override = default;
 
  // virtual void BeforeThreadedGenerateData();
 
  void
  GPUGenerateData() override;
};
 
class GPUBinaryThresholdImageFilterFactory : public ObjectFactoryBase
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(GPUBinaryThresholdImageFilterFactory);
 
  using Self = GPUBinaryThresholdImageFilterFactory;
  using Superclass = ObjectFactoryBase;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;
 
  const char *
  GetITKSourceVersion() const override
  {
    return ITK_SOURCE_VERSION;
  }
  const char *
  GetDescription() const override
  {
    return "A Factory for GPUBinaryThresholdImageFilter";
  }
  itkFactorylessNewMacro(Self);
 
  itkOverrideGetNameOfClassMacro(GPUBinaryThresholdImageFilterFactory);
 
  static void
  RegisterOneFactory()
  {
    auto factory = GPUBinaryThresholdImageFilterFactory::New();
 
    itk::ObjectFactoryBase::RegisterFactory(factory);
  }
 
private:
#define OverrideThresholdFilterTypeMacro(ipt, opt, dm)                                                   \
  {                                                                                                      \
    using InputImageType = itk::Image<ipt, dm>;                                                          \
    using OutputImageType = itk::Image<opt, dm>;                                                         \
    this->RegisterOverride(                                                                              \
      typeid(itk::BinaryThresholdImageFilter<InputImageType, OutputImageType>).name(),                   \
      typeid(itk::GPUBinaryThresholdImageFilter<InputImageType, OutputImageType>).name(),                \
      "GPU Binary Threshold Image Filter Override",                                                      \
      true,                                                                                              \
      itk::CreateObjectFunction<GPUBinaryThresholdImageFilter<InputImageType, OutputImageType>>::New()); \
  }                                                                                                      \
  ITK_MACROEND_NOOP_STATEMENT
 
  GPUBinaryThresholdImageFilterFactory()
  {
    if (IsGPUAvailable())
    {
      OverrideThresholdFilterTypeMacro(unsigned char, unsigned char, 1);
      OverrideThresholdFilterTypeMacro(char, char, 1);
      OverrideThresholdFilterTypeMacro(float, float, 1);
      OverrideThresholdFilterTypeMacro(int, int, 1);
      OverrideThresholdFilterTypeMacro(unsigned int, unsigned int, 1);
      OverrideThresholdFilterTypeMacro(double, double, 1);
 
      OverrideThresholdFilterTypeMacro(unsigned char, unsigned char, 2);
      OverrideThresholdFilterTypeMacro(char, char, 2);
      OverrideThresholdFilterTypeMacro(float, float, 2);
      OverrideThresholdFilterTypeMacro(int, int, 2);
      OverrideThresholdFilterTypeMacro(unsigned int, unsigned int, 2);
      OverrideThresholdFilterTypeMacro(double, double, 2);
 
      OverrideThresholdFilterTypeMacro(unsigned char, unsigned char, 3);
      OverrideThresholdFilterTypeMacro(unsigned short, unsigned short, 3);
      OverrideThresholdFilterTypeMacro(char, char, 3);
      OverrideThresholdFilterTypeMacro(float, float, 3);
      OverrideThresholdFilterTypeMacro(int, int, 3);
      OverrideThresholdFilterTypeMacro(unsigned int, unsigned int, 3);
      OverrideThresholdFilterTypeMacro(double, double, 3);
    }
  }
};
 
} // end of namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkGPUBinaryThresholdImageFilter.hxx"
#endif
 
#endif