itkSparseFieldLevelSetImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkSparseFieldLevelSetImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2002/09/24 16:48:17 $
  Version:   $Revision: 1.11 $

  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 __itkSparseFieldLevelSetImageFilter_h_
#define __itkSparseFieldLevelSetImageFilter_h_

#include "itkFiniteDifferenceImageFilter.h"
#include "itkMultiThreader.h"
#include "itkSparseFieldLayer.h"
#include "itkObjectStore.h"
#include <vector>
#include "itkSmartNeighborhoodIterator.h"

namespace itk {

template <class TValueType>
class SparseFieldLevelSetNode
{
public:
  TValueType m_Value;
  SparseFieldLevelSetNode *Next;
  SparseFieldLevelSetNode *Previous;
};

template <class TNeighborhoodType>
class SparseFieldCityBlockNeighborList
{
public:
  typedef TNeighborhoodType NeighborhoodType;
  typedef typename NeighborhoodType::OffsetType OffsetType;
  typedef typename NeighborhoodType::RadiusType RadiusType;
  itkStaticConstMacro(Dimension, unsigned int,
                      NeighborhoodType::Dimension );

  const RadiusType &GetRadius() const
  { return m_Radius; }
  
  const unsigned int &GetArrayIndex(unsigned int i) const
  { return m_ArrayIndex[i]; }

  const OffsetType &GetNeighborhoodOffset(unsigned int i) const
  { return m_NeighborhoodOffset[i]; }

  const unsigned int &GetSize() const
  { return m_Size; }
  
  SparseFieldCityBlockNeighborList();
  ~SparseFieldCityBlockNeighborList() {}

  void Print(std::ostream &os) const;
  
private:
  unsigned int m_Size;
  RadiusType m_Radius;
  std::vector<unsigned int> m_ArrayIndex;
  std::vector<OffsetType>   m_NeighborhoodOffset;
};

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

  typedef typename Superclass::TimeStepType TimeStepType;
  
  itkNewMacro(Self);

  itkTypeMacro(SparseFieldLevelSetImageFilter, FiniteDifferenceImageFilter);

  typedef TInputImage  InputImageType;
  typedef TOutputImage OutputImageType;
  typedef typename OutputImageType::IndexType IndexType;
  itkStaticConstMacro(ImageDimension, unsigned int,
                      TOutputImage::ImageDimension);

  typedef typename OutputImageType::ValueType ValueType;

  typedef SparseFieldLevelSetNode<IndexType> LayerNodeType;
  
  typedef SparseFieldLayer<LayerNodeType> LayerType;
  typedef typename LayerType::Pointer LayerPointerType;

  typedef std::vector<LayerPointerType> LayerListType;
  
  typedef char StatusType;
  
  typedef Image<StatusType, itkGetStaticConstMacro(ImageDimension)>
    StatusImageType;

  typedef ObjectStore<LayerNodeType> LayerNodeStorageType;

  typedef std::vector<ValueType> UpdateBufferType;

  itkSetMacro(NumberOfLayers, unsigned int);
  itkGetMacro(NumberOfLayers, unsigned int);

  itkSetMacro(IsoSurfaceValue, ValueType);
  itkGetMacro(IsoSurfaceValue, ValueType);

  itkGetMacro(RMSChange, ValueType);
    
protected:
  SparseFieldLevelSetImageFilter();
  ~SparseFieldLevelSetImageFilter();
  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)
  { return (value + dt * change); }

  virtual void PostProcessOutput();

  virtual void InitializeBackgroundPixels();
 
  void Initialize();

  void CopyInputToOutput(); 

  void AllocateUpdateBuffer();

  void ApplyUpdate(TimeStepType dt);

  TimeStepType CalculateChange();

  void ConstructLayer(StatusType from, StatusType to);

  void ConstructActiveLayer();

  void InitializeActiveLayerValues();
  
  void PropagateLayerValues(StatusType from, StatusType to,
                            StatusType promote, int InOrOut);

  void PropagateAllLayerValues();

  void UpdateActiveLayerValues(TimeStepType dt, LayerType *StatusUpList,
                               LayerType *StatusDownList);
  void ProcessStatusList(LayerType *InputList, LayerType *OutputList,
                         StatusType ChangeToStatus, StatusType SearchForStatus);

  void ProcessOutsideList(LayerType *OutsideList, StatusType ChangeToStatus);
  
  SparseFieldCityBlockNeighborList<SmartNeighborhoodIterator<OutputImageType> >
    m_NeighborList;
  
  static double m_ConstantGradientValue;

  static ValueType m_ValueOne;

  static ValueType m_ValueZero;

  static StatusType m_StatusChanging;

  static StatusType m_StatusActiveChangingUp;

  static StatusType m_StatusActiveChangingDown;

  static StatusType m_StatusNull;
    
  typename OutputImageType::Pointer m_ShiftedImage;

  LayerListType m_Layers;

  unsigned int m_NumberOfLayers;

  typename StatusImageType::Pointer m_StatusImage;

  typename LayerNodeStorageType::Pointer m_LayerNodeStore;
  
  ValueType m_IsoSurfaceValue;

  UpdateBufferType m_UpdateBuffer;

  ValueType m_RMSChange;  

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

};
  
  
} // end namespace itk

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

#endif

---