otbNeighborhoodMajorityVotingImageFilter.h

Go to the documentation of this file.

/*
 * Copyright (C) 2005-2024 Centre National d'Etudes Spatiales (CNES)
 *
 * This file is part of Orfeo Toolbox
 *
 *     https://www.orfeo-toolbox.org/
 *
 * 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
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * 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 otbNeighborhoodMajorityVotingImageFilter_h
#define otbNeighborhoodMajorityVotingImageFilter_h
 
// First make sure that the configuration is available.
// This line can be removed once the optimized versions
// gets integrated into the main directories.
#include "itkConfigure.h"
 
#include "itkMorphologyImageFilter.h"
#include "itkBinaryBallStructuringElement.h"
 
namespace otb
{
 
template <class TInputImage, class TOutputImage = TInputImage,
          class TKernel = typename itk::BinaryBallStructuringElement<typename TInputImage::PixelType, TInputImage::ImageDimension>>
class ITK_EXPORT NeighborhoodMajorityVotingImageFilter : public itk::MorphologyImageFilter<TInputImage, TOutputImage, TKernel>
{
public:
 
  typedef NeighborhoodMajorityVotingImageFilter Self;
  typedef itk::MorphologyImageFilter<TInputImage, TOutputImage, TKernel> Superclass;
  typedef itk::SmartPointer<Self>       Pointer;
  typedef itk::SmartPointer<const Self> ConstPointer;
 
  itkNewMacro(Self);
 
  itkTypeMacro(NeighborhoodMajorityVotingImageFilter, itk::MorphologyImageFilter);
 
  typedef typename Superclass::PixelType PixelType;
 
 
  typedef typename Superclass::KernelIteratorType KernelIteratorType;
 
  typedef typename Superclass::NeighborhoodIteratorType NeighborhoodIteratorType;
 
  typedef typename Superclass::KernelType KernelType;
 
 
  typedef typename Superclass::DefaultBoundaryConditionType DefaultBoundaryConditionType;
 
  itkStaticConstMacro(InputImageDimension, unsigned int, TInputImage::ImageDimension);
  itkStaticConstMacro(OutputImageDimension, unsigned int, TOutputImage::ImageDimension);
  itkStaticConstMacro(KernelDimension, unsigned int, TKernel::NeighborhoodDimension);
 
 
  typedef typename TKernel::PixelType KernelPixelType;
 
#ifdef ITK_USE_CONCEPT_CHECKING
  itkConceptMacro(InputConvertibleToOutputCheck, (itk::Concept::Convertible<PixelType, typename TOutputImage::PixelType>));
  itkConceptMacro(SameDimensionCheck1, (itk::Concept::SameDimension<InputImageDimension, OutputImageDimension>));
  itkConceptMacro(SameDimensionCheck2, (itk::Concept::SameDimension<InputImageDimension, KernelDimension>));
  itkConceptMacro(InputGreaterThanComparableCheck, (itk::Concept::GreaterThanComparable<PixelType>));
  itkConceptMacro(KernelGreaterThanComparableCheck, (itk::Concept::GreaterThanComparable<KernelPixelType>));
 
#endif
 
 
  // Creates a SetLabelForNoDataPixels method
  virtual void SetLabelForNoDataPixels(const PixelType _arg)
  {
    itkDebugMacro("setting LabelForNoDataPixels to " << _arg);
    if (this->m_LabelForNoDataPixels != _arg)
    {
      this->m_LabelForNoDataPixels = _arg;
 
      m_MajorityVotingBoundaryCondition.SetConstant(m_LabelForNoDataPixels);
      this->OverrideBoundaryCondition(&m_MajorityVotingBoundaryCondition);
 
      this->Modified();
    }
  }
 
  // Creates a SetLabelForUndecidedPixels method
  itkSetMacro(LabelForUndecidedPixels, PixelType);
 
  // Creates a SetKeepOriginalLabelBool method
  itkSetMacro(KeepOriginalLabelBool, bool);
 
  // Process only isolated pixels
  itkSetMacro(OnlyIsolatedPixels, bool);
  itkSetMacro(IsolatedThreshold, unsigned int);
 
 
protected:
  NeighborhoodMajorityVotingImageFilter();
  ~NeighborhoodMajorityVotingImageFilter() override{};
 
  PixelType Evaluate(const NeighborhoodIteratorType& nit, const KernelIteratorType kernelBegin, const KernelIteratorType kernelEnd) override;
 
  void GenerateOutputInformation() override;
 
 
  // Type to store the useful information from the label histogram
  struct HistoSummary
  {
    unsigned int freqCenterLabel;
    PixelType    majorityLabel;
    bool         majorityUnique;
  };
 
  struct CompareHistoFequencies
  {
    typedef std::pair<PixelType, unsigned int> HistoValueType;
    bool operator()(const HistoValueType& a, const HistoValueType& b)
    {
      return a.second > b.second;
    }
  };
 
  // Get a histogram of frequencies of labels with the 2 highest
  // frequencies sorted in decreasing order and return the frequency
  // of the label of the center pixel
  const HistoSummary ComputeNeighborhoodHistogramSummary(const NeighborhoodIteratorType& nit, const KernelIteratorType kernelBegin,
                                                         const KernelIteratorType kernelEnd) const;
 
private:
  NeighborhoodMajorityVotingImageFilter(const Self&) = delete;
  void operator=(const Self&) = delete;
 
  // Default boundary condition for majority voting filter, defaults to
  DefaultBoundaryConditionType m_MajorityVotingBoundaryCondition;
 
  PixelType m_LabelForNoDataPixels;
  PixelType m_LabelForUndecidedPixels;
  bool      m_KeepOriginalLabelBool;
 
  // Choose to filter only isolated pixels
  bool m_OnlyIsolatedPixels;
  // The center pixel is isolated if there are fewer neighbours than
  // this threshold with the same label
  unsigned int m_IsolatedThreshold;
 
}; // end of class
 
} // end namespace otb
 
#ifndef ITK_MANUAL_INSTANTIATION
#include "otbNeighborhoodMajorityVotingImageFilter.hxx"
#endif
 
 
#endif