itkLMMSEVectorImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkLMMSEVectorImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2006/03/29 14:53:40 $
  Version:   $Revision: 1.5 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef __itkLMMSEVectorImageFilter_h
#define __itkLMMSEVectorImageFilter_h

#include "itkImageToImageFilter.h"
#include "itkImage.h"
#include "itkVectorImage.h"
#include "itkVectorImageCastFilter.h"
#include "itkLMMSEVectorImageFilterStep.h"
#include "itkExtractVolumeFilter.h"
#include "itkComputeStatisticsWherePositiveFilter.h"
#include "itkComputeRestrictedHistogram.h"
#include "itkMaskedMeanImageFilter.h"
namespace itk
{
template <class TInputImage, class TOutputImage>
class ITK_EXPORT LMMSEVectorImageFilter : public ImageToImageFilter<TInputImage, TOutputImage>
{
public:
  typedef TInputImage  InputImageType;
  typedef TOutputImage OutputImageType;

  typedef LMMSEVectorImageFilter                              Self;
  typedef ImageToImageFilter<InputImageType, OutputImageType> Superclass;
  typedef SmartPointer<Self>                                  Pointer;
  typedef SmartPointer<const Self>                            ConstPointer;

  typedef typename OutputImageType::Pointer     OutputImagePointer;
  typedef typename InputImageType::Pointer      InputImagePointer;
  typedef typename InputImageType::ConstPointer InputImageConstPointer;

  typedef itk::VectorImage<float, TInputImage::ImageDimension> InternalImageType;
  typedef typename InternalImageType::Pointer                  InternalImagePointer;
  typedef typename InternalImageType::ConstPointer             InternalImageConstPointer;

  itkNewMacro( Self );

  itkTypeMacro( LMMSEVectorImageFilter, ImageToImageFilter );

  typedef typename InputImageType::PixelType  InputPixelType;
  typedef typename OutputImageType::PixelType OutputPixelType;

  typedef typename InputImageType::RegionType  InputImageRegionType;
  typedef typename OutputImageType::RegionType OutputImageRegionType;
  typedef typename InputImageType::SizeType    InputSizeType;

  typedef LMMSEVectorImageFilterStep<InternalImageType, InternalImageType> LMMSEStepType;
  typedef typename LMMSEStepType::Pointer                                  LMMSEStepPointer;

  typedef itk::VectorImageCastFilter<InputImageType, InternalImageType>               InputCastType;
  typedef typename InputCastType::Pointer                                             InputCastPointer;
  typedef itk::VectorImageCastFilter<InternalImageType, OutputImageType>              OutputCastType;
  typedef typename OutputCastType::Pointer                                            OutputCastPointer;
  typedef float                                                                       ScalarPixelType;
  typedef itk::Image<ScalarPixelType, TInputImage::ImageDimension>                    ScalarImageType;
  typedef itk::ExtractVolumeFilter<InternalImageType, ScalarImageType>                ExtractType;
  typedef typename ExtractType::Pointer                                               ExtractPointer;
  typedef itk::MaskedMeanImageFilter<ScalarImageType, ScalarImageType>                LocalMeanType;
  typedef typename LocalMeanType::Pointer                                             LocalMeanPointer;
  typedef itk::ComputeStatisticsWherePositiveFilter<ScalarImageType, ScalarImageType> StatsType;
  typedef typename StatsType::Pointer                                                 StatsPointer;
  typedef itk::ComputeRestrictedHistogram<ScalarImageType, ScalarImageType>           HistogramType;
  typedef typename HistogramType::Pointer                                             HistogramPointer;

  itkSetMacro( RadiusFiltering, InputSizeType );
  itkSetMacro( RadiusEstimation, InputSizeType );
  itkGetConstReferenceMacro(RadiusFiltering, InputSizeType);
  itkGetConstReferenceMacro(RadiusEstimation, InputSizeType);

  itkSetMacro( Iterations, unsigned int );
  itkGetMacro( Iterations, unsigned int );

  itkSetMacro( FirstBaseline, unsigned int );
  itkGetMacro( FirstBaseline, unsigned int );

  itkSetMacro( MaximumNumberOfBins, unsigned long );
  itkGetMacro( MaximumNumberOfBins, unsigned long );

  itkSetMacro( UseAbsoluteValue, bool );
  itkBooleanMacro( UseAbsoluteValue );
  itkGetMacro( KeepValue, bool );
  itkSetMacro( KeepValue, bool );
  itkBooleanMacro( KeepValue );

  itkSetMacro( MinimumNoiseSTD, double );
  itkGetMacro( MinimumNoiseSTD, double );
  itkSetMacro( MaximumNoiseSTD, double );
  itkGetMacro( MaximumNoiseSTD, double );

  itkSetMacro(MinimumNumberOfUsedVoxelsEstimation, int );
  itkSetMacro(MinimumNumberOfUsedVoxelsFiltering, int );
  itkGetMacro(MinimumNumberOfUsedVoxelsEstimation, int );
  itkGetMacro(MinimumNumberOfUsedVoxelsFiltering, int );

  itkSetMacro(HistogramResolutionFactor, double );
  itkGetMacro(HistogramResolutionFactor, double );

  itkSetMacro( Channels, unsigned int );
  itkGetConstReferenceMacro( Channels, unsigned int );

  virtual void GenerateInputRequestedRegion()
  throw (InvalidRequestedRegionError);

protected:
  LMMSEVectorImageFilter();
  virtual ~LMMSEVectorImageFilter()
  {
  }
  void GenerateData();

  void PrintSelf( std::ostream& os, Indent indent) const;

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

  InputSizeType m_RadiusEstimation;
  InputSizeType m_RadiusFiltering;
  unsigned int  m_Iterations;
  bool          m_UseAbsoluteValue;
  bool          m_KeepValue;

  int m_MinimumNumberOfUsedVoxelsEstimation;
  int m_MinimumNumberOfUsedVoxelsFiltering;

  double m_MinimumNoiseSTD;
  double m_MaximumNoiseSTD;

  unsigned int m_FirstBaseline;
  double       m_HistogramResolutionFactor;

  unsigned int  m_Channels;
  unsigned long m_MaximumNumberOfBins;
};

} // end namespace itk

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

#endif