itkMammogramFatEstimationFitMetricForArray.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 __itkMammogramFatEstimationFitMetricForArray_h
#define __itkMammogramFatEstimationFitMetricForArray_h

#include <itkMammogramFatEstimationFitMetric.h>
#include <itkImageRegionIterator.h>
#include <itkImageRegionIteratorWithIndex.h>
#include <itkImageRegionConstIterator.h>
#include <itkImageRegionConstIteratorWithIndex.h>

namespace itk {
  
template <class TInputImage>
class  ITK_EXPORT MammogramFatEstimationFitMetricForArray :
  public MammogramFatEstimationFitMetric
{
//  Software Guide : EndCodeSnippet
public:
  typedef MammogramFatEstimationFitMetricForArray  Self;
  typedef MammogramFatEstimationFitMetric   Superclass;
  typedef SmartPointer< Self >       Pointer;
  typedef SmartPointer< const Self > ConstPointer;
  
  itkTypeMacro( MammogramFatEstimationFitMetricForArray, MammogramFatEstimationFitMetric );

  itkNewMacro(Self);

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

  itkStaticConstMacro( ParametricSpaceDimension, unsigned int, 4 );

  virtual unsigned int GetNumberOfParameters(void) const override  
  {
    return ParametricSpaceDimension;
  }

  typedef float DistancePixelType;

  void SetNumberOfDistances( unsigned int nDistances ) {
    m_NumberOfDistances = nDistances;
  }

  void SetInputArray( float *minIntensityVsEdgeDistance ) {
    m_MinIntensityVsEdgeDistance = minIntensityVsEdgeDistance;
  }

  DistancePixelType GetMaxDistance( void ) { return m_MaxDistance; }
  void SetMaxDistance( DistancePixelType distance ) { m_MaxDistance = distance; }


  MeasureType GetValue( const ParametersType &parameters ) const override;

  void WriteIntensityVsEdgeDistToFile( std::string fileOutputIntensityVsEdgeDist ) override;
  void WriteFitToFile( std::string fileOutputFit,
                       const ParametersType &parameters ) override;

  void GenerateFatArray( unsigned int nDistances, float *fatEstimate, 
                         const ParametersType &parameters );

protected:

  MammogramFatEstimationFitMetricForArray();
  virtual ~MammogramFatEstimationFitMetricForArray();
  MammogramFatEstimationFitMetricForArray(const Self &) {}
  void operator=(const Self &) {}
  void PrintSelf(std::ostream & os, Indent indent) const override;

  unsigned int m_NumberOfDistances;
  float *m_MinIntensityVsEdgeDistance;

  DistancePixelType      m_MaxDistance;

  double CalculateFit( double d, const ParametersType &parameters );

};

} // end namespace itk

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

#endif