itkShapeLabelMapFilter.h

Go to the documentation of this file.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkShapeLabelMapFilter.h,v $
  Language:  C++
  Date:      $Date: 2009-08-11 14:24:44 $
  Version:   $Revision: 1.6 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef __itkShapeLabelMapFilter_h
#define __itkShapeLabelMapFilter_h
#include "itkLabelPerimeterEstimationCalculator.h"

#include "itkInPlaceLabelMapFilter.h"

namespace itk {

template<class TImage, class TLabelImage=Image< ITK_TYPENAME TImage::PixelType, ::itk::GetImageDimension<TImage>::ImageDimension> > 
class ITK_EXPORT ShapeLabelMapFilter : 
public
InPlaceLabelMapFilter<TImage>
{
public:
  typedef ShapeLabelMapFilter           Self;
  typedef InPlaceLabelMapFilter<TImage> Superclass;
  typedef SmartPointer<Self>            Pointer;
  typedef SmartPointer<const Self>      ConstPointer;

  typedef TImage                               ImageType;
  typedef typename ImageType::Pointer          ImagePointer;
  typedef typename ImageType::ConstPointer     ImageConstPointer;
  typedef typename ImageType::PixelType        PixelType;
  typedef typename ImageType::IndexType        IndexType;
  typedef typename ImageType::SizeType         SizeType;
  typedef typename ImageType::LabelObjectType  LabelObjectType;
  typedef typename LabelObjectType::MatrixType MatrixType;
  typedef typename LabelObjectType::VectorType VectorType;

  typedef TLabelImage                              LabelImageType;
  typedef typename LabelImageType::Pointer         LabelImagePointer;
  typedef typename LabelImageType::ConstPointer    LabelImageConstPointer;
  typedef typename LabelImageType::PixelType       LabelPixelType;
  
  itkStaticConstMacro(ImageDimension, unsigned int, TImage::ImageDimension );

  typedef LabelPerimeterEstimationCalculator< LabelImageType > PerimeterCalculatorType;

  itkNewMacro(Self);  

  itkTypeMacro(ShapeLabelMapFilter, InPlaceLabelMapFilter );

#ifdef ITK_USE_CONCEPT_CHECKING

/*  itkConceptMacro(InputEqualityComparableCheck,
    (Concept::EqualityComparable<InputImagePixelType>));
  itkConceptMacro(IntConvertibleToInputCheck,
    (Concept::Convertible<int, InputImagePixelType>));
  itkConceptMacro(InputOStreamWritableCheck,
    (Concept::OStreamWritable<InputImagePixelType>));*/

#endif

  itkSetMacro( ComputeFeretDiameter, bool );
  itkGetConstReferenceMacro( ComputeFeretDiameter, bool );
  itkBooleanMacro( ComputeFeretDiameter );

  itkSetMacro( ComputePerimeter, bool );
  itkGetConstReferenceMacro( ComputePerimeter, bool );
  itkBooleanMacro( ComputePerimeter );

  void SetLabelImage( const TLabelImage *input )
    {
    m_LabelImage = input;
    }

protected:
  ShapeLabelMapFilter();
  ~ShapeLabelMapFilter() {};

  virtual void ThreadedProcessLabelObject( LabelObjectType * labelObject );
  virtual void BeforeThreadedGenerateData();
  virtual void AfterThreadedGenerateData();

  void PrintSelf(std::ostream& os, Indent indent) const;

private:
  ShapeLabelMapFilter(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

  bool                                      m_ComputeFeretDiameter;
  bool                                      m_ComputePerimeter;
  LabelImageConstPointer                    m_LabelImage;
  typename PerimeterCalculatorType::Pointer m_PerimeterCalculator;

  void ComputeFeretDiameter( LabelObjectType * labelObject );

  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 );

}; // end of class

} // end namespace itk
  
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkShapeLabelMapFilter.txx"
#endif

#endif