tools/clitkUnsharpMaskGenericFilter.txx

   1 /*=========================================================================
   2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
   3
   4   Authors belong to:
   5   - University of LYON              http://www.universite-lyon.fr/
   6   - Léon Bérard cancer center       http://www.centreleonberard.fr
   7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
   8
   9   This software is distributed WITHOUT ANY WARRANTY; without even
  10   the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  11   PURPOSE.  See the copyright notices for more information.
  12
  13   It is distributed under dual licence
  14
  15   - BSD        See included LICENSE.txt file
  16   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  17 ===========================================================================**/
  18 #ifndef clitkUnsharpMaskGenericFilter_txx
  19 #define clitkUnsharpMaskGenericFilter_txx
  20
  21 /* =================================================
  22  * @file   clitkUnsharpMaskGenericFilter.txx
  23  * @author Jef Vandemeulebroucke <jef@creatis.insa-lyon.fr>
  24  * @date   29 june 2009
  25  *
  26  * @brief
  27  *
  28  ===================================================*/
  29
  30 // itk include
  31 #include "itkRecursiveGaussianImageFilter.h"
  32 #include "itkSubtractImageFilter.h"
  33 #include <clitkCommon.h>
  34
  35 namespace clitk
  36 {
  37
  38 //--------------------------------------------------------------------
  39 template<class args_info_type>
  40 UnsharpMaskGenericFilter<args_info_type>::UnsharpMaskGenericFilter():
  41   ImageToImageGenericFilter<Self>("UnsharpMask")
  42 {
  43   InitializeImageType<2>();
  44   InitializeImageType<3>();
  45   //InitializeImageType<4>();
  46 }
  47 //--------------------------------------------------------------------
  48
  49
  50 //--------------------------------------------------------------------
  51 template<class args_info_type>
  52 template<unsigned int Dim>
  53 void UnsharpMaskGenericFilter<args_info_type>::InitializeImageType()
  54 {
  55   ADD_DEFAULT_IMAGE_TYPES(Dim);
  56 }
  57 //--------------------------------------------------------------------
  58
  59
  60 //--------------------------------------------------------------------
  61 template<class args_info_type>
  62 void UnsharpMaskGenericFilter<args_info_type>::SetArgsInfo(const args_info_type & a)
  63 {
  64   mArgsInfo=a;
  65   this->SetIOVerbose(mArgsInfo.verbose_flag);
  66   if (mArgsInfo.imagetypes_flag) this->PrintAvailableImageTypes();
  67
  68   if (mArgsInfo.input_given) {
  69     this->SetInputFilename(mArgsInfo.input_arg);
  70   }
  71   if (mArgsInfo.output_given) {
  72     this->SetOutputFilename(mArgsInfo.output_arg);
  73   }
  74 }
  75 //--------------------------------------------------------------------
  76
  77 //--------------------------------------------------------------------
  78 // Update with the number of dimensions and the pixeltype
  79 //--------------------------------------------------------------------
  80 template<class args_info_type>
  81 template<class InputImageType>
  82 void
  83 UnsharpMaskGenericFilter<args_info_type>::UpdateWithInputImageType()
  84 {
  85
  86   // Reading input
  87   typename InputImageType::Pointer input = this->template GetInput<InputImageType>(0);
  88
  89   // Main filter
  90   typedef typename InputImageType::PixelType PixelType;
  91   typedef itk::Image<float, InputImageType::ImageDimension> OutputImageType;
  92
  93   // Filter
  94   typedef itk::RecursiveGaussianImageFilter<InputImageType, OutputImageType> RecursiveGaussianImageFilterType;
  95   typename RecursiveGaussianImageFilterType::Pointer gaussianFilter=RecursiveGaussianImageFilterType::New();
  96   gaussianFilter->SetInput(input);
  97   gaussianFilter->SetSigma(mArgsInfo.sigma_arg);
  98
  99   typedef itk::SubtractImageFilter<InputImageType, OutputImageType, OutputImageType> SubtractFilterType;
 100   typename SubtractFilterType::Pointer subtractFilter = SubtractFilterType::New();
 101   subtractFilter->SetInput1(input);
 102   subtractFilter->SetInput2(gaussianFilter->GetOutput());
 103   subtractFilter->Update();
 104
 105   this->template SetNextOutput<OutputImageType>(subtractFilter->GetOutput());
 106 }
 107 //--------------------------------------------------------------------
 108
 109
 110 }//end clitk
 111
 112 #endif //#define clitkUnsharpMaskGenericFilter_txx