itkMIDASMeanIntensityWithinARegionFilter.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 itkMIDASMeanIntensityWithinARegionFilter_h
#define itkMIDASMeanIntensityWithinARegionFilter_h

#include <itkImageToImageFilter.h>

namespace itk
{

  template <class TInputImage1, class TInputImage2, class TOutputImage>
  class ITK_EXPORT MIDASMeanIntensityWithinARegionFilter : public ImageToImageFilter<TInputImage1, TOutputImage>
  {
  public:
    typedef MIDASMeanIntensityWithinARegionFilter           Self;
    typedef ImageToImageFilter<TInputImage1, TOutputImage>  SuperClass;
    typedef SmartPointer<Self>                              Pointer;
    typedef SmartPointer<const Self>                        ConstPointer;

    itkNewMacro(Self);

    itkTypeMacro(MIDASMeanIntensityWithinARegionFilter, ImageToImageFilter);

    typedef TInputImage1                              InputMainImageType;
    typedef typename InputMainImageType::ConstPointer InputMainImageConstPointer;
    typedef typename InputMainImageType::Pointer      InputMainImagePointer;
    typedef typename InputMainImageType::RegionType   InputMainImageRegionType;
    typedef typename InputMainImageType::PixelType    InputMainImagePixelType;
    typedef typename InputMainImageType::IndexType    InputMainImageIndexType;
    typedef typename InputMainImageType::SizeType     InputMainImageSizeType;

    typedef TInputImage2                              InputMaskImageType;
    typedef typename InputMaskImageType::ConstPointer InputMaskImageConstPointer;
    typedef typename InputMaskImageType::Pointer      InputMaskImagePointer;
    typedef typename InputMaskImageType::RegionType   InputMaskImageRegionType;
    typedef typename InputMaskImageType::PixelType    InputMaskImagePixelType;
    typedef typename InputMaskImageType::IndexType    InputMaskImageIndexType;
    typedef typename InputMaskImageType::SizeType     InputMaskImageSizeType;

    typedef TInputImage1                              OutputImageType;
    typedef typename OutputImageType::Pointer         OutputImagePointer;
    typedef typename OutputImageType::RegionType      OutputImageRegionType;
    typedef typename OutputImageType::SizeType        OutputImageSizeType;
    typedef typename OutputImageType::IndexType       OutputImageIndexType;

    void SetGreyScaleImageInput(const InputMainImageType* image);

    void SetBinaryImageInput(const InputMaskImageType* image);

    double GetMeanIntensityMainImage();

    unsigned long int GetCount();

    itkSetMacro(InValue, InputMaskImagePixelType);
    itkGetConstMacro(InValue, InputMaskImagePixelType);

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

    void AllocateOutputs();

    virtual void BeforeThreadedGenerateData();
    
    // The main method to implement in derived classes, note, its threaded.
    virtual void ThreadedGenerateData(const InputMainImageRegionType &outputRegionForThread, ThreadIdType ThreadID);

    virtual void AfterThreadedGenerateData();

  private:
    MIDASMeanIntensityWithinARegionFilter(const Self&); //purposely not implemented
    void operator=(const Self&); //purposely not implemented
   
    double                         m_MeanIntensityMainImage;
    unsigned long int              m_Counter;
    std::vector<double>            m_TotalIntensityVector;
    std::vector<unsigned long int> m_CountPixelsVector;
    InputMaskImagePixelType        m_InValue;

  };

}

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

#endif