otbShapeAttributesLabelMapFilter.h

Go to the documentation of this file.

/*
 * Copyright (C) 1999-2011 Insight Software Consortium
 * 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 otbShapeAttributesLabelMapFilter_h
#define otbShapeAttributesLabelMapFilter_h
 
#include "itkConfigure.h"
 
#include "otbLabelMapFeaturesFunctorImageFilter.h"
#include "otbImage.h"
#include "otbPolygon.h"
#include "otbLabelObjectToPolygonFunctor.h"
#include "otbFlusserPathFunction.h"
#include "otbSimplifyPathFunctor.h"
 
namespace otb
{
 
namespace Functor
{
template< unsigned int VOffsetDimension >
class OffsetLexicographicCompare
{
public:
  bool operator()(itk::Offset< VOffsetDimension > const & l,
                  itk::Offset< VOffsetDimension > const & r) const
  {
    for ( unsigned int i = 0; i < VOffsetDimension; ++i )  {
      if ( l[i] < r[i] )  {
        return true;
      }
      else if ( l[i] > r[i] ) {
        return false;
      }
    }
    return false;
  }
};
 
template <class TLabelObject, class TLabelImage>
class ShapeAttributesLabelObjectFunctor
{
public:
  typedef ShapeAttributesLabelObjectFunctor Self;
 
  typedef TLabelObject LabelObjectType;
 
  typedef typename LabelObjectType::ConstLineIterator ConstLineIteratorType;
 
  typedef TLabelImage LabelImageType;
 
  typedef typename TLabelObject::PolygonType PolygonType;
  typedef Functor::LabelObjectToPolygonFunctor<LabelObjectType, PolygonType> PolygonFunctorType;
  typedef FlusserPathFunction<PolygonType> FlusserPathFunctionType;
  typedef SimplifyPathFunctor<PolygonType, PolygonType> SimplifyPolygonFunctorType;
 
  typedef unsigned int DimensionType;
 
  itkStaticConstMacro(ImageDimension, unsigned int, TLabelObject::ImageDimension);
  typedef itk::ImageRegion< TLabelObject::ImageDimension > RegionType;
  typedef itk::Offset< TLabelObject::ImageDimension > OffsetType;
  typedef otb::Functor::OffsetLexicographicCompare<TLabelObject::ImageDimension> OffsetLexicographicCompare;
 
  ShapeAttributesLabelObjectFunctor();
 
  virtual ~ShapeAttributesLabelObjectFunctor()
  {
  }
 
  bool operator!=(const Self& self);
  bool operator==(const Self& self);
 
  void SetComputePerimeter(bool flag);
 
  bool GetComputePerimeter() const;
 
  void SetComputeFlusser(bool flag);
 
  bool GetComputeFlusser() const;
 
  void SetComputePolygon(bool flag);
 
  bool GetComputePolygon() const;
 
  void SetComputeFeretDiameter(bool flag);
 
  bool GetComputeFeretDiameter() const;
 
  void SetReducedAttributeSet(bool flag);
 
  bool GetReducedAttributeSet() const;
 
  void SetLabelImage(const TLabelImage* image);
 
  const TLabelImage* GetLabelImage() const;
 
  inline void operator()(LabelObjectType* lo);
 
private:
  static long factorial(long n);
 
  static long doubleFactorial(long n);
 
  static double gammaN2p1(long n);
 
  static double hyperSphereVolume(double radius);
 
  static double hyperSpherePerimeter(double radius);
 
  static double hyperSphereRadiusFromVolume(double volume);
 
  double ComputePerimeter(LabelObjectType* labelObject, const RegionType& region);
 
  typedef itk::Offset<2>                                                          Offset2Type;
  typedef itk::Offset<3>                                                          Offset3Type;
  typedef itk::Vector<double, 2>                                                  Spacing2Type;
  typedef itk::Vector<double, 3>                                                  Spacing3Type;
  typedef otb::Functor::OffsetLexicographicCompare<2>                             OffsetLexicographicCompare2;
  typedef otb::Functor::OffsetLexicographicCompare<3>                             OffsetLexicographicCompare3;
  typedef std::map<Offset2Type, itk::SizeValueType, OffsetLexicographicCompare2>  MapIntercept2Type;
  typedef std::map<Offset3Type, itk::SizeValueType, OffsetLexicographicCompare3>  MapIntercept3Type;
  template<class TMapIntercept, class TSpacing> double PerimeterFromInterceptCount( TMapIntercept & intercepts, const TSpacing & spacing );
 
#if !defined(ITK_DO_NOT_USE_PERIMETER_SPECIALIZATION)
  double PerimeterFromInterceptCount(MapIntercept2Type& intercepts, const Spacing2Type spacing);
  double PerimeterFromInterceptCount(MapIntercept3Type& intercepts, const Spacing3Type spacing);
#endif
 
  bool m_ComputeFeretDiameter;
 
  bool m_ComputePerimeter;
 
  bool m_ComputeFlusser;
 
  bool m_ComputePolygon;
 
  bool m_ReducedAttributeSet;
 
  typename LabelImageType::ConstPointer m_LabelImage;
};
 
} // End namespace Functor
 
template <class TImage, class TLabelImage = Image<typename TImage::PixelType, TImage::ImageDimension>>
class ITK_EXPORT ShapeAttributesLabelMapFilter
    : public otb::LabelMapFeaturesFunctorImageFilter<TImage, typename Functor::ShapeAttributesLabelObjectFunctor<typename TImage::LabelObjectType, TLabelImage>>
{
public:
  typedef TImage                              ImageType;
  typedef typename ImageType::LabelObjectType LabelObjectType;
  typedef typename ImageType::RegionType      InputImageRegionType;
  typedef TLabelImage                         LabelImageType;
  typedef Functor::ShapeAttributesLabelObjectFunctor<LabelObjectType, LabelImageType> FunctorType;
 
  typedef ShapeAttributesLabelMapFilter Self;
  typedef LabelMapFeaturesFunctorImageFilter<ImageType, FunctorType> Superclass;
  typedef itk::SmartPointer<Self>       Pointer;
  typedef itk::SmartPointer<const Self> ConstPointer;
 
  typedef typename ImageType::Pointer ImagePointer;
 
  itkStaticConstMacro(ImageDimension, unsigned int, TImage::ImageDimension);
 
  itkNewMacro(Self);
 
  itkTypeMacro(ShapeAttributesLabelMapFilter, LabelMapFeaturesFunctorImageFilter);
 
  void SetComputeFeretDiameter(bool flag);
  bool GetComputeFeretDiameter() const;
  itkBooleanMacro(ComputeFeretDiameter);
 
  void SetComputePerimeter(bool flag);
  bool GetComputePerimeter() const;
  itkBooleanMacro(ComputePerimeter);
 
  void SetComputePolygon(bool flag);
  bool GetComputePolygon() const;
  itkBooleanMacro(ComputePolygon);
 
  void SetComputeFlusser(bool flag);
  bool GetComputeFlusser() const;
  itkBooleanMacro(ComputeFlusser);
 
 
  void SetReducedAttributeSet(bool flag);
  bool GetReducedAttributeSet() const;
  itkBooleanMacro(ReducedAttributeSet);
 
  void               SetLabelImage(const TLabelImage*);
  const TLabelImage* GetLabelImage() const;
 
protected:
  ShapeAttributesLabelMapFilter()
  {
  }
 
  ~ShapeAttributesLabelMapFilter() override
  {
  }
 
  void PrintSelf(std::ostream& os, itk::Indent indent) const override;
 
  void AllocateOutputs() override;
 
  void BeforeThreadedGenerateData() override;
 
  void GenerateInputRequestedRegion() override;
 
  void EnlargeOutputRequestedRegion(itk::DataObject*) override{};
 
private:
  ShapeAttributesLabelMapFilter(const Self&) = delete;
  void operator=(const Self&) = delete;
 
}; // end of class
 
} // end namespace otb
 
#ifndef OTB_MANUAL_INSTANTIATION
#include "otbShapeAttributesLabelMapFilter.hxx"
#endif
 
#endif