itkSignedDanielssonDistanceMapImageFilter.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 itkSignedDanielssonDistanceMapImageFilter_h
#define itkSignedDanielssonDistanceMapImageFilter_h
 
#include "itkDanielssonDistanceMapImageFilter.h"
#include "itkSubtractImageFilter.h"
 
// Simple functor to invert an image for Outside Danielsson distance map
namespace itk
{
namespace Functor
{
template <typename InputPixelType>
class ITK_TEMPLATE_EXPORT InvertIntensityFunctor{ public: InputPixelType operator()(InputPixelType input)
                                                    const { if (input){ return InputPixelType {};
}
 
return NumericTraits<InputPixelType>::OneValue();
} // namespace itk
}
;
} // namespace Functor
} // namespace itk
 
namespace itk
{
template <typename TInputImage, typename TOutputImage, typename TVoronoiImage = TInputImage>
class ITK_TEMPLATE_EXPORT SignedDanielssonDistanceMapImageFilter : public ImageToImageFilter<TInputImage, TOutputImage>
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(SignedDanielssonDistanceMapImageFilter);
 
  using Self = SignedDanielssonDistanceMapImageFilter;
  using Superclass = ImageToImageFilter<TInputImage, TOutputImage>;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;
 
  itkNewMacro(Self);
 
  itkOverrideGetNameOfClassMacro(SignedDanielssonDistanceMapImageFilter);
 
  using InputImageType = TInputImage;
 
  using OutputImageType = TOutputImage;
 
  using RegionType = typename InputImageType::RegionType;
 
  using IndexType = typename RegionType::IndexType;
 
  using OffsetType = typename InputImageType::OffsetType;
 
  using PixelType = typename InputImageType::PixelType;
 
  using SizeType = typename RegionType::SizeType;
 
  static constexpr unsigned int InputImageDimension = InputImageType::ImageDimension;
 
  using VectorImageType = Image<OffsetType, Self::InputImageDimension>;
 
  using InputImagePointer = typename InputImageType::ConstPointer;
 
  using OutputImagePointer = typename OutputImageType::Pointer;
 
  using VectorImagePointer = typename VectorImageType::Pointer;
 
  using VoronoiImageType = TVoronoiImage;
  using VoronoiImagePointer = typename VoronoiImageType::Pointer;
  using VoronoiPixelType = typename VoronoiImageType::PixelType;
 
  using typename Superclass::DataObjectPointer;
 
  itkSetMacro(SquaredDistance, bool);
 
  itkGetConstReferenceMacro(SquaredDistance, bool);
 
  itkBooleanMacro(SquaredDistance);
 
  itkSetMacro(UseImageSpacing, bool);
 
  itkGetConstReferenceMacro(UseImageSpacing, bool);
 
  itkBooleanMacro(UseImageSpacing);
 
  itkSetMacro(InsideIsPositive, bool);
 
  itkGetConstReferenceMacro(InsideIsPositive, bool);
 
  itkBooleanMacro(InsideIsPositive);
 
  VoronoiImageType *
  GetVoronoiMap();
 
  OutputImageType *
  GetDistanceMap();
 
  VectorImageType *
  GetVectorDistanceMap();
 
  using DataObjectPointerArraySizeType = ProcessObject::DataObjectPointerArraySizeType;
  using Superclass::MakeOutput;
  DataObjectPointer
  MakeOutput(DataObjectPointerArraySizeType idx) override;
 
#ifdef ITK_USE_CONCEPT_CHECKING
  // Begin concept checking
  itkConceptMacro(IntConvertibleToInputCheck, (Concept::Convertible<int, PixelType>));
  itkConceptMacro(InputHasNumericTraitsCheck, (Concept::HasNumericTraits<PixelType>));
  itkConceptMacro(OutputImagePixelTypeIsFloatingPointCheck,
                  (Concept::IsFloatingPoint<typename OutputImageType::PixelType>));
  // End concept checking
#endif
 
protected:
  SignedDanielssonDistanceMapImageFilter();
  ~SignedDanielssonDistanceMapImageFilter() override = default;
  void
  PrintSelf(std::ostream & os, Indent indent) const override;
 
  void
  GenerateData() override;
 
private:
  bool m_SquaredDistance{};
  bool m_UseImageSpacing{ true };
  bool m_InsideIsPositive{}; // ON is treated as inside pixels
}; // end of SignedDanielssonDistanceMapImageFilter
   // class
} // end namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkSignedDanielssonDistanceMapImageFilter.hxx"
#endif
 
#endif