itkMidsagittalPlaneExtractionImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkMidsagittalPlaneExtractionImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2003/09/10 14:28:51 $
  Version:   $Revision: 1.4 $

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

#include "itkImageToImageFilter.h"
#include "itkPoint.h"

namespace itk
{
  
  typedef struct
  {
    double shiftX;
    double *shiftsY;
    double stepY;
    double n;
    double center;
    double angleEstimate;
//    itk::Point<double, InputImageType::ImageDimension> flippedCenterOfMass;
//    itk::Point<double, 3> flippedCenterOfMass;
//    itk::Point flippedCenterOfMass;
    double *flippedCenterOfMass;
  } estimateType;

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

  itkNewMacro(Self);
  
  itkTypeMacro(MidsagittalPlaneExtractionImageFilter, ImageToImageFilter);

  typedef TInputImage InputImageType;
  typedef typename    InputImageType::Pointer    InputImagePointer;
  typedef typename    InputImageType::RegionType InputImageRegionType; 
  typedef typename    InputImageType::PixelType  InputImagePixelType; 
  typedef TOutputImage OutputImageType;
  typedef typename     OutputImageType::Pointer    OutputImagePointer;
  typedef typename     OutputImageType::RegionType OutputImageRegionType;
  typedef typename     OutputImageType::PixelType  OutputImagePixelType;


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

protected:
  MidsagittalPlaneExtractionImageFilter();
  virtual ~MidsagittalPlaneExtractionImageFilter() {};
  void PrintSelf(std::ostream& os, Indent indent) const;
  double FindAngle(typename Superclass::InputImageConstPointer, int angleChoice, int numberOfPoints, double step, estimateType* anEstimate);
  double FindShift(typename Superclass::InputImageConstPointer, double numberOfPoints, double step, double shiftEstimate);
  double findHighestPoint(typename Superclass::InputImageConstPointer, int dimensionChoice, double startingHeight);
//  itk::Point<double, 3> *findCentroid(typename Superclass::InputImageConstPointer);
//  itk::Point<double, TInputImage::ImageDimension> *findCentroid(typename Superclass::InputImageConstPointer);

//  itk::Point<double, 3> *findCentroid(typename Superclass::InputImageConstPointer volume, TInputImage::RegionType::IndexType dimensionMask);
//  itk::Point<double, 3> *findCentroid(typename Superclass::InputImageConstPointer volume);
    double *findCentroid(typename Superclass::InputImageConstPointer volume, double* dimensionMask);

//  itk::Point *findCentroid(typename Superclass::InputImageConstPointer, TInputImage::RegionType::IndexType dimensionMask);
  estimateType* estimateAngle(typename Superclass::InputImageConstPointer);
  double estimateShift(typename Superclass::InputImageConstPointer);
  void GenerateData(void);

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

  unsigned int m_Direction;

};

} // end namespace itk

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

#endif

---