freImageMerger.txx

Go to the documentation of this file.

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

  Program:   F.R.E.E. - flexible registration evaluation engine
  Version:   v.1.0.0
  Date:      $Date: 2006/09/01 12:00:00 $
  Module:    $RCSfile: freImageMerger.txx,v $
  Language:  C++



  Copyright (c) 2007 Ralf o Floca (Department of Medical Informatics,
  Institute for Medical Biometry and Informatics, University of Heidelberg,
  Germany). All rights reserved.
  See FREECopyright.txt or http://www.mi.med.uni-hd.de/free/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 __freImageMerger_txx
#define __freImageMerger_txx

#include "freImageMerger.h"
#include "freExceptions.h"

#include "freDifferenceToRGBImageFilter.h"

#include "itkConstrainedValueAdditionImageFilter.h" 
#include "itkWeightedAddImageFilter.h" 
#include "itkCheckerBoardImageFilter.h"
#include "itkSquaredDifferenceImageFilter.h" 
#include "itkSqrtImageFilter.h"

namespace FREE
{

template <class TImage>
typename ImageMerger<TImage>::RGBImagePointer
ImageMerger<TImage>::
GenerateAlternationImage(ImageType* pInput1, ImageType* pInput2, ImageType* pInput3)
{
  if (!pInput1) throwStaticExceptionMacro("Cannot compute alternation image: Input1 is null");
  if (!pInput2) throwStaticExceptionMacro("Cannot compute alternation image: Input2 is null");
  if (!pInput3) throwStaticExceptionMacro("Cannot compute alternation image: Input3 is null");

  typedef FREE::DifferenceToRGBImageFilter<ImageType, ImageType, RGBImageType> DifRGBFilterType;

        ImagePointer smpDiff1 = GenerateDifferenceImage(pInput1, pInput2);
        ImagePointer smpDiff2 = GenerateDifferenceImage(pInput1, pInput3);

        typename DifRGBFilterType::Pointer difFilter = DifRGBFilterType::New();

        difFilter->SetInput1( smpDiff1 );
        difFilter->SetInput2( smpDiff2 );
        difFilter->Update();

        typename RGBImageType::Pointer smpResult = difFilter->GetOutput();

  return smpResult;
};

template <class TImage>
typename ImageMerger<TImage>::RGBImagePointer
ImageMerger<TImage>::
GenerateAlternationImage(ImageType* pInput3) const
{
  return GenerateAlternationImage(m_Input1, m_Input2, pInput3);
};

template <class TImage>
typename ImageMerger<TImage>::ImagePointer
ImageMerger<TImage>::
GenerateBlendImage(ImageType* pInput1, ImageType* pInput2, unsigned int iMergeValue)
{
  if (!pInput1) throwStaticExceptionMacro("Cannot compute blend image: Input1 is null");
  if (!pInput2) throwStaticExceptionMacro("Cannot compute blend image: Input2 is null");

  typedef itk::WeightedAddImageFilter<ImageType, ImageType, ImageType> FilterType;

        typename FilterType::Pointer filter = FilterType::New();
        filter->SetInput1( pInput1 );
        filter->SetInput2( pInput2 );
        filter->SetAlpha(iMergeValue/100.0);
        filter->Update();
        ImagePointer smpResult = filter->GetOutput();

  return smpResult;
};

template <class TImage>
typename ImageMerger<TImage>::ImagePointer
ImageMerger<TImage>::
GenerateBlendImage(unsigned int iMergeValue) const
{
  return GenerateBlendImage(m_Input1, m_Input2, iMergeValue);
};

template <class TImage>
typename ImageMerger<TImage>::ImagePointer
ImageMerger<TImage>::
GenerateCheckerImage(ImageType* pInput1, ImageType* pInput2, unsigned int iPatternSize)
{
  if (!pInput1) throwStaticExceptionMacro("Cannot compute checker: Input1 is null");
  if (!pInput2) throwStaticExceptionMacro("Cannot compute checker: Input2 is null");
        
  typedef itk::CheckerBoardImageFilter<ImageType> FilterType;

        typename FilterType::PatternArrayType pattern;
        pattern.Fill(iPatternSize);

        typename FilterType::Pointer filter = FilterType::New();
        filter->SetInput1( pInput1 );
        filter->SetInput2( pInput2 );
        filter->SetCheckerPattern(pattern);
        filter->Update();
        ImagePointer result = filter->GetOutput();
  return result;
};

template <class TImage>
typename ImageMerger<TImage>::ImagePointer
ImageMerger<TImage>::
GenerateCheckerImage(unsigned int iPatternSize) const
{
  return GenerateCheckerImage(m_Input1, m_Input2, iPatternSize);
};

template <class TImage>
typename ImageMerger<TImage>::ImagePointer
ImageMerger<TImage>::
GenerateDifferenceImage(ImageType* pInput1, ImageType* pInput2)
{
        if (!pInput1) throwStaticExceptionMacro("Cannot compute difference Input1 is null");
        if (!pInput2) throwStaticExceptionMacro("Cannot compute difference Input2 is null");

  typedef itk::SquaredDifferenceImageFilter<ImageType, ImageType, ImageType> DifFilterType;
        typedef itk::SqrtImageFilter<ImageType, ImageType> SqrtFilterType;

        typename DifFilterType::Pointer difFilter = DifFilterType::New();
        typename SqrtFilterType::Pointer sqrtFilter = SqrtFilterType::New();

        difFilter->SetInput1( pInput1 );
        difFilter->SetInput2( pInput2 );
        sqrtFilter->SetInput(difFilter->GetOutput());
        sqrtFilter->Update();
        return sqrtFilter->GetOutput();
};

template <class TImage>
typename ImageMerger<TImage>::ImagePointer
ImageMerger<TImage>::
GenerateDifferenceImage() const
{
  return GenerateDifferenceImage(m_Input1, m_Input2);
};

template <class TImage>
typename ImageMerger<TImage>::RGBImagePointer
ImageMerger<TImage>::
GenerateColoredDifferenceImage(ImageType* pInput1, ImageType* pInput2)
{
        if (!pInput1) throwStaticExceptionMacro("Cannot compute colored difference Input1 is null");
        if (!pInput2) throwStaticExceptionMacro("Cannot compute colored difference Input2 is null");

  typedef FREE::DifferenceToRGBImageFilter<ImageType, ImageType, RGBImageType> DifRGBFilterType;

        typename DifRGBFilterType::Pointer difFilter = DifRGBFilterType::New();

        difFilter->SetInput1( pInput1 );
        difFilter->SetInput2( pInput2 );
        difFilter->Update();

        RGBImagePointer smpResult = difFilter->GetOutput();
  return smpResult;
};

template <class TImage>
typename ImageMerger<TImage>::RGBImagePointer
ImageMerger<TImage>::
GenerateColoredDifferenceImage() const
{
  return GenerateColoredDifferenceImage(m_Input1, m_Input2);
};

template <class TImage>
ImageMerger<TImage>::
ImageMerger(ImageType* pInput1, ImageType* pInput2)
{
  m_Input1 = pInput1;
  m_Input2 = pInput2;
};

template <class TImage>
ImageMerger<TImage>::
~ImageMerger() {};

} //end of namespace free

#endif

---