itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2003/02/10 16:58:51 $
  Version:   $Revision: 1.8 $

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

#include "itkBloxBoundaryPointImage.h"
#include "itkBloxBoundaryProfileImage.h"
#include "itkImageToImageFilter.h"
#include "itkSize.h"

namespace itk
{

template< typename TSourceImage >
class ITK_EXPORT BloxBoundaryPointImageToBloxBoundaryProfileImageFilter :
public ImageToImageFilter<TSourceImage,BloxBoundaryProfileImage< ::itk::GetImageDimension<TSourceImage>::ImageDimension> >
{
public:
  itkStaticConstMacro(NDimensions, unsigned int, TSourceImage::ImageDimension);

  typedef BloxBoundaryPointImageToBloxBoundaryProfileImageFilter  Self;
  typedef ImageToImageFilter<TSourceImage,BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)> >  Superclass;
  typedef SmartPointer<Self>   Pointer;
  typedef SmartPointer<const Self>  ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(itkBloxBoundaryPointImageToBloxBoundaryProfileImageFilter, ImageToImageFilter);

  typedef BloxBoundaryPointImage<itkGetStaticConstMacro(NDimensions)> BoundaryPointImageType;
  typedef typename BoundaryPointImageType::Pointer BoundaryPointImagePointer;
  typedef typename BoundaryPointImageType::RegionType BoundaryPointImageRegionType; 
  typedef typename BoundaryPointImageType::PixelType BoundaryPointImagePixelType; 
  typedef typename BoundaryPointImageType::ConstPointer BoundaryPointImageConstPointer;

  typedef TSourceImage SourceImageType;
  typedef typename SourceImageType::Pointer SourceImagePointer;
  typedef typename SourceImageType::RegionType SourceImageRegionType; 
  typedef typename SourceImageType::PixelType SourceImagePixelType; 
  typedef typename SourceImageType::ConstPointer SourceImageConstPointer;

  typedef BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)> OutputImageType;
  typedef typename OutputImageType::Pointer OutputImagePointer;
  typedef typename OutputImageType::RegionType OutputImageRegionType;
  typedef typename OutputImageType::PixelType OutputImagePixelType;

  typedef typename BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)>::IndexType IndexType;

  typedef typename BloxBoundaryProfileImage<itkGetStaticConstMacro(NDimensions)>::PixelType PixelType;

  typedef Point<double, itkGetStaticConstMacro(NDimensions)> PositionType;

  typedef typename PositionType::VectorType VectorType;

  void SetInput1( const SourceImageType * image1 );

  void SetInput2( const BoundaryPointImageType * image2 );

  double FindAccumulatorMaximum();
  
  double FindAccumulatorMinimum();

  void FindBoundaryProfiles();

  bool AddSplatToAccumulatorAndNormalizer(int binNumber, double weight, double sourcePixelValue);

  void NormalizeSplatAccumulator();

  int FitProfile();
  
  void Initialize(double setUniqueAxis, double setSymmetricAxes, unsigned int numberOfBins,
                  unsigned int splatMethod, unsigned int spaceDimension)
    {
    m_UniqueAxis = setUniqueAxis;
    m_SymmetricAxes = setSymmetricAxes;
    m_NumberOfBins = numberOfBins;
    m_SplatMethod = splatMethod;
    m_SpaceDimension = spaceDimension;

    m_Accumulator = new double[m_NumberOfBins];
    m_Normalizer = new double[m_NumberOfBins];
    }

protected:
  BloxBoundaryPointImageToBloxBoundaryProfileImageFilter();
  ~BloxBoundaryPointImageToBloxBoundaryProfileImageFilter()
  {
    delete [] m_Accumulator;
    delete [] m_Normalizer;
    delete [] m_NormalizedAccumulator;
    delete [] m_FinalParameters;
  };
  void PrintSelf(std::ostream& os, Indent indent) const;

  void GenerateData();

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

  BoundaryPointImagePointer m_BoundaryPointImage;
  const double* m_BPImageOrigin;
  const double* m_BPImageSpacing;
    
  typename TSourceImage::Pointer m_OriginalImage;
  const double* m_OriginalImageOrigin;
  const double* m_OriginalImageSpacing;
   
  double m_UniqueAxis;

  double m_SymmetricAxes;

  unsigned int m_NumberOfBins;

  unsigned int m_SplatMethod;

  unsigned long int m_NumBoundaryProfiles;

  double * m_Accumulator;
  
  double * m_Normalizer;

  double * m_NormalizedAccumulator;

  double * m_FinalParameters;

  unsigned int m_SpaceDimension;
};

} // end namespace itk

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

#endif

---