mitkExtendedStatisticsImageFilter.h

Go to the documentation of this file.

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

The Medical Imaging Interaction Toolkit (MITK)

Copyright (c) German Cancer Research Center,
Division of Medical and Biological Informatics.
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 or http://www.mitk.org for details.

===================================================================*/
#ifndef __mitkExtendedStatisticsImageFilter
#define __mitkExtendedStatisticsImageFilter

#include "itkStatisticsImageFilter.h"
#include <mbilog.h>
#include <mitkLogMacros.h>

namespace itk
{
  template< class TInputImage >
  class ExtendedStatisticsImageFilter : public StatisticsImageFilter< TInputImage >
  {
  public:
    typedef ExtendedStatisticsImageFilter                    Self;
    typedef StatisticsImageFilter< TInputImage >             Superclass;
    typedef SmartPointer< Self >                             Pointer;
    typedef SmartPointer< const Self >                       ConstPointer;
    typedef typename Superclass::RealType                    RealType;
    typedef typename Superclass::RealObjectType              RealObjectType;
    typedef typename Superclass::PixelType                   PixelType;

    typedef itk::Statistics::Histogram< RealType > HistogramType;
    typedef typename HistogramType::Pointer        HistogramPointer;

    itkFactorylessNewMacro( Self );
    itkCloneMacro( Self );
    itkTypeMacro( ExtendedStatisticsImageFilter, StatisticsImageFilter );

    double GetSkewness() const
    {
      return this->GetSkewnessOutput()->Get();
    }

    double GetMedian() const
    {
        return this->GetMedianOutput()->Get();
    }

    double GetKurtosis() const
    {
      return this->GetKurtosisOutput()->Get();
    }

    /* \brief Return the computed MPP.
      */
      double GetMPP() const
    {
      return this->GetMPPOutput()->Get();
    }

    double GetUniformity() const
    {
      return this->GetUniformityOutput()->Get();
    }

      double GetEntropy() const
    {
      return this->GetEntropyOutput()->Get();
    }

      double GetUPP() const
    {
      return this->GetUPPOutput()->Get();
    }


      const typename HistogramType::Pointer
        GetHistogram()
      {
          if (m_HistogramCalculated)
          {
              return m_Histogram;
          }
          else
          {
              return ITK_NULLPTR;
          }
      }


    void SetHistogramParameters(const int numBins, RealType lowerBound,
                                RealType upperBound);

  protected:

    ExtendedStatisticsImageFilter();

    virtual ~ExtendedStatisticsImageFilter(){};

    void BeforeThreadedGenerateData();

    void  ThreadedGenerateData(const typename StatisticsImageFilter<TInputImage>::RegionType &
                               outputRegionForThread,
                               ThreadIdType threadId);

    void AfterThreadedGenerateData();


    RealObjectType* GetSkewnessOutput();

    const RealObjectType* GetSkewnessOutput() const;

    RealObjectType* GetKurtosisOutput();

    const RealObjectType* GetKurtosisOutput() const;

    RealObjectType* GetMPPOutput();

    const RealObjectType* GetMPPOutput() const;

    RealObjectType* GetEntropyOutput();

    const RealObjectType* GetEntropyOutput() const;

    RealObjectType* GetUniformityOutput();

    const RealObjectType* GetUniformityOutput() const;

    RealObjectType* GetUPPOutput();

    const RealObjectType* GetUPPOutput() const;

    RealObjectType* GetMedianOutput();

    const RealObjectType* GetMedianOutput() const;

    virtual DataObject::Pointer MakeOutput( ProcessObject::DataObjectPointerArraySizeType idx );

private:
    Array< RealType >       m_ThreadSum;
    Array< RealType >       m_SumOfSquares;
    Array< RealType >       m_SumOfCubes;
    Array< RealType >       m_SumOfQuadruples;
    Array< SizeValueType >  m_Count;
    Array< SizeValueType >  m_PositivePixelCount;
    Array< RealType >       m_ThreadSumOfPositivePixels;
    Array< PixelType >      m_ThreadMin;
    Array< PixelType >      m_ThreadMax;
    std::vector< HistogramPointer > m_HistogramPerThread;
    HistogramPointer        m_Histogram;
    bool                    m_UseHistogram;
    bool                    m_HistogramCalculated;
    RealType                m_LowerBound, m_UpperBound;
    int                     m_NumBins;


  }; // end of class

} // end namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "mitkExtendedStatisticsImageFilter.hxx"
#endif

#endif