itkMeanVoxelwiseIntensityOfMultipleImages.h

Go to the documentation of this file.

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

  NifTK: A software platform for medical image computing.

  Copyright (c) University College London (UCL). All rights reserved.

  This software is distributed WITHOUT ANY WARRANTY; without even
  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  PURPOSE.

  See LICENSE.txt in the top level directory for details.

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

#ifndef itkMeanVoxelwiseIntensityOfMultipleImages_h
#define itkMeanVoxelwiseIntensityOfMultipleImages_h

#include <itkImageToImageFilter.h>
#include <itkEulerAffineTransform.h>


namespace itk {
  
template<class TInputImage, class TOutputImage>
class ITK_EXPORT MeanVoxelwiseIntensityOfMultipleImages:
    public ImageToImageFilter< TInputImage, TOutputImage >
{
public:
  typedef MeanVoxelwiseIntensityOfMultipleImages         Self;
  typedef ImageToImageFilter< TInputImage,TOutputImage > Superclass;
  typedef SmartPointer< Self >                           Pointer;
  typedef SmartPointer< const Self >                     ConstPointer;
  
  itkTypeMacro( MeanVoxelwiseIntensityOfMultipleImages, ImageToImageFilter );

  itkNewMacro(Self);

  itkStaticConstMacro(ImageDimension, unsigned int,
                      TInputImage::ImageDimension);

  typedef TInputImage                           InputImageType;
  typedef typename InputImageType::Pointer      InputImagePointer;
  typedef typename InputImageType::ConstPointer InputImageConstPointer;
  typedef typename InputImageType::RegionType   InputImageRegionType;
  typedef typename InputImageType::PixelType    InputImagePixelType;
  typedef typename InputImageType::SpacingType  InputImageSpacingType;
  typedef typename InputImageType::IndexType    InputImageIndexType;
  typedef typename InputImageType::PointType    InputImagePointType;

  typedef typename InputImageRegionType::SizeType    InputImageSizeType;

  typedef TOutputImage                          OutputImageType;
  typedef typename OutputImageType::Pointer     OutputImagePointer;
  typedef typename OutputImageType::RegionType  OutputImageRegionType;
  typedef typename OutputImageType::PixelType   OutputImagePixelType;
  typedef typename OutputImageType::SpacingType OutputImageSpacingType;
  typedef typename OutputImageType::IndexType   OutputImageIndexType;
  typedef typename OutputImageType::PointType   OutputImagePointType;

  typedef typename OutputImageRegionType::SizeType    OutputImageSizeType;

  typedef typename itk::EulerAffineTransform< double, ImageDimension, ImageDimension > TransformType;

  typedef typename TransformType::TranslationType TranslationType;
  typedef typename TransformType::InputPointType CenterType;
  typedef typename TransformType::ScaleType ScaleType;

  itkSetMacro( SubtractMinima, bool );
  itkGetMacro( SubtractMinima, bool );

  itkSetMacro( ExpandOutputRegion, double );
  itkGetMacro( ExpandOutputRegion, double );

  void SetTranslations( std::vector< TranslationType > &translations ) {
    m_TranslationVectors = translations;
    this->Modified();
  }
  void SetCenters( std::vector< CenterType > &centers ) {
    m_CenterVectors = centers;
    this->Modified();
  }
  void SetScales( std::vector< ScaleType > &scales ) {
    m_ScaleVectors = scales;
    this->Modified();
  }
  void ClearTransformation( void ) {
    m_TranslationVectors.clear();
    m_CenterVectors.clear();
    m_ScaleVectors.clear();
    this->Modified();
  }

#ifdef ITK_USE_CONCEPT_CHECKING

  itkConceptMacro(InputHasNumericTraitsCheck,
                  (Concept::HasNumericTraits<InputImagePixelType>));
  itkConceptMacro(OutputHasPixelTraitsCheck,
                  (Concept::HasPixelTraits<OutputImagePixelType>));
#endif

protected:
  MeanVoxelwiseIntensityOfMultipleImages();
  virtual ~MeanVoxelwiseIntensityOfMultipleImages() {};
  void PrintSelf(std::ostream& os, Indent indent) const;
  
  void GenerateOutputInformation();
  virtual void GenerateInputRequestedRegion();

  void GenerateData();

  std::vector< TranslationType > m_TranslationVectors;
  std::vector< CenterType > m_CenterVectors;
  std::vector< ScaleType > m_ScaleVectors;

  // Create a border around the mean image
  double m_ExpandOutputRegion;

  bool m_SubtractMinima;

  OutputImageRegionType  m_OutRegion;
  OutputImageSizeType    m_OutSize;
  OutputImageSpacingType m_OutSpacing;
  OutputImagePointType   m_OutOrigin;


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

} // end namespace itk

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

#endif