itkAntiAliasBinaryImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkAntiAliasBinaryImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2002/05/16 20:07:11 $
  Version:   $Revision: 1.4 $

  Copyright (c) 2002 Insight 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 __itkAntiAliasBinaryImageFilter_h_
#define __itkAntiAliasBinaryImageFilter_h_
#include "itkSparseFieldLevelSetImageFilter.h"
#include "itkCurvatureFlowFunction.h"

namespace itk {

template <class TInputImage, class TOutputImage>
class AntiAliasBinaryImageFilter
  : public SparseFieldLevelSetImageFilter<TInputImage, TOutputImage>
{
public:
  typedef AntiAliasBinaryImageFilter Self;
  typedef SparseFieldLevelSetImageFilter<TInputImage, TOutputImage> Superclass;
  typedef SmartPointer<Self>  Pointer;
  typedef SmartPointer<const Self>  ConstPointer;

  typedef typename Superclass::ValueType ValueType;
  typedef typename Superclass::IndexType IndexType;
  typedef typename Superclass::TimeStepType TimeStepType;
  typedef typename Superclass::OutputImageType OutputImageType;
  typedef typename Superclass::InputImageType  InputImageType;

  
  typedef CurvatureFlowFunction<OutputImageType> CurvatureFunctionType;
  
  typedef typename TInputImage::ValueType  BinaryValueType;
  
  itkNewMacro(Self);

  itkTypeMacro(AntiAliasBinaryImageFilter, SparseFieldLevelSetImageFilter);

  itkSetMacro(MaximumRMSError, ValueType);
  itkGetMacro(MaximumRMSError, ValueType);

  itkSetMacro(MaximumIterations, unsigned int);
  itkGetMacro(MaximumIterations, unsigned int);
  
  itkGetMacro(UpperBinaryValue, BinaryValueType);
  itkGetMacro(LowerBinaryValue, BinaryValueType);
  
protected:
  AntiAliasBinaryImageFilter();
  ~AntiAliasBinaryImageFilter() {}
  virtual void PrintSelf(std::ostream& os, Indent indent) const;

  inline virtual ValueType CalculateUpdateValue(const IndexType &idx,
                                                const TimeStepType &dt,
                                                const ValueType &value,
                                                const ValueType &change);
  
  bool Halt()
  {
    if (this->GetElapsedIterations() >= m_MaximumIterations)
      {
        itkWarningMacro("This filter has passed the maximum number of iterations allowed and will be halted.  This means that the solution did not converge to the MaximumRMSError you specified.  Try setting m_MaximumRMSError to a higher value, or set m_MaxmimumIterations to a higher value.");
        return true;
      }
    else if ( this->GetElapsedIterations() == 0)
      { return false; }
    else if ( this->GetRMSChange() <= m_MaximumRMSError )
      { return true; }
    else
      { return false; }
  }

  void GenerateData();
  
private:
  ValueType m_MaximumRMSError;
  BinaryValueType m_UpperBinaryValue;
  BinaryValueType m_LowerBinaryValue;
  typename CurvatureFunctionType::Pointer m_CurvatureFunction;
  
  unsigned int m_MaximumIterations;
  
};

} // end namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkAntiAliasBinaryImageFilter.txx"
#endif

#endif

---