old/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_h
#define __mitkExtendedStatisticsImageFilter_h

#include "itkStatisticsImageFilter.h"
#include "itkScalarImageToHistogramGenerator.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;

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


        typedef itk::Statistics::ScalarImageToHistogramGenerator< TInputImage >
          HistogramGeneratorType;

    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 HistogramGeneratorType::HistogramType*
        GetHistogram()
      {
        return m_HistogramGenerator->GetOutput();
      }

    void SetBinSize(int size);

  protected:

    ExtendedStatisticsImageFilter();

    virtual ~ExtendedStatisticsImageFilter(){};

    void AfterThreadedGenerateData();
    void CalculateHistogram();

    void ComputeSkewnessKurtosisAndMPP();
    void ComputeEntropyUniformityMedianAndUPP();

    typename  HistogramGeneratorType::Pointer  m_HistogramGenerator;
    int m_BinSize;
    bool m_HistogramCalculated;


    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 );


  }; // end of class

} // end namespace itk

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

#endif