segmentation/clitkExtractPatientGenericFilter.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
  19 #include "clitkImageCommon.h"
  20
  21 //--------------------------------------------------------------------
  22 template<class ArgsInfoType>
  23 clitk::ExtractPatientGenericFilter<ArgsInfoType>::ExtractPatientGenericFilter():
  24   ImageToImageGenericFilter<Self>("ExtractPatient")
  25 {
  26   // Default values
  27   cmdline_parser_clitkExtractPatient_init(&mArgsInfo);
  28   InitializeImageType<3>();
  29 }
  30 //--------------------------------------------------------------------
  31
  32
  33 //--------------------------------------------------------------------
  34 template<class ArgsInfoType>
  35 template<unsigned int Dim>
  36 void clitk::ExtractPatientGenericFilter<ArgsInfoType>::InitializeImageType()
  37 {
  38   ADD_IMAGE_TYPE(Dim, short);
  39   // ADD_IMAGE_TYPE(Dim, int);
  40   // ADD_IMAGE_TYPE(Dim, float);
  41 }
  42 //--------------------------------------------------------------------
  43
  44
  45 //--------------------------------------------------------------------
  46 template<class ArgsInfoType>
  47 void clitk::ExtractPatientGenericFilter<ArgsInfoType>::SetArgsInfo(const ArgsInfoType & a)
  48 {
  49   mArgsInfo=a;
  50   SetIOVerbose(mArgsInfo.verbose_flag);
  51   if (mArgsInfo.imagetypes_flag) this->PrintAvailableImageTypes();
  52   if (mArgsInfo.input_given)  SetInputFilename(mArgsInfo.input_arg);
  53   if (mArgsInfo.output_given)  AddOutputFilename(mArgsInfo.output_arg);
  54 }
  55 //--------------------------------------------------------------------
  56
  57
  58 //--------------------------------------------------------------------
  59 template<class ArgsInfoType>
  60 template<class FilterType>
  61 void
  62 clitk::ExtractPatientGenericFilter<ArgsInfoType>::
  63 SetOptionsFromArgsInfoToFilter(FilterType * f)
  64 {
  65   f->SetVerboseOptionFlag(mArgsInfo.verbose_flag);
  66   f->SetVerboseStepFlag(mArgsInfo.verboseStep_flag);
  67   f->SetWriteStepFlag(mArgsInfo.writeStep_flag);
  68   f->SetVerboseWarningFlag(!mArgsInfo.verboseWarningOff_flag);
  69   f->SetVerboseMemoryFlag(mArgsInfo.verboseMemory_flag);
  70
  71   if (mArgsInfo.afdb_given)
  72     f->SetAFDBFilename(mArgsInfo.afdb_arg);
  73
  74   f->SetOutputPatientFilename(mArgsInfo.output_arg);
  75
  76   f->SetUpperThreshold(mArgsInfo.upper_arg);
  77   f->SetLowerThreshold(mArgsInfo.lower_arg);
  78
  79   f->SetDecomposeAndReconstructDuringFirstStep(mArgsInfo.erode1_flag);
  80
  81   // Convert to SizeType
  82   typename FilterType::InputImageSizeType s;
  83   if (mArgsInfo.radius1_given) {
  84     ConvertOptionMacro(mArgsInfo.radius1, s, 3, false);
  85     f->SetRadius1(s);
  86   }
  87
  88   f->SetMaximumNumberOfLabels1(mArgsInfo.max1_arg);
  89   f->SetNumberOfNewLabels1(mArgsInfo.new1_arg);
  90
  91   f->SetDecomposeAndReconstructDuringSecondStep(mArgsInfo.erode2_flag);
  92
  93   // Convert to SizeType
  94   if (mArgsInfo.radius2_given) {
  95     ConvertOptionMacro(mArgsInfo.radius2, s, 3, false);
  96     f->SetRadius2(s);
  97   }
  98
  99   f->SetMaximumNumberOfLabels2(mArgsInfo.max2_arg);
 100   f->SetNumberOfNewLabels2(mArgsInfo.new2_arg);
 101
 102   f->SetFirstKeep(mArgsInfo.firstKeep_arg);
 103   f->SetLastKeep(mArgsInfo.lastKeep_arg);
 104
 105   f->SetFinalOpenClose(mArgsInfo.openClose_flag);
 106   f->SetAutoCrop(!mArgsInfo.noAutoCrop_flag);
 107 }
 108 //--------------------------------------------------------------------
 109
 110
 111 //--------------------------------------------------------------------
 112 // Update with the number of dimensions and the pixeltype
 113 //--------------------------------------------------------------------
 114 template<class ArgsInfoType>
 115 template<class ImageType>
 116 void clitk::ExtractPatientGenericFilter<ArgsInfoType>::UpdateWithInputImageType()
 117 {
 118   // Output image type
 119   typedef itk::Image<uchar, ImageType::ImageDimension> OutputImageType;
 120
 121   // Reading input
 122   typename ImageType::Pointer input = this->template GetInput<ImageType>(0);
 123
 124   // Create filter
 125   typedef clitk::ExtractPatientFilter<ImageType> FilterType;
 126   typename FilterType::Pointer filter = FilterType::New();
 127
 128   // Set the filter (needed for example for threaded monitoring)
 129   this->SetFilterBase(filter);
 130
 131   // Set global Options
 132   filter->SetInput(input);
 133   SetOptionsFromArgsInfoToFilter<FilterType>(filter);
 134
 135   // Go !
 136   filter->Update();
 137
 138   // Write/Save results
 139   typename OutputImageType::Pointer output = filter->GetOutput();
 140   this->template SetNextOutput<OutputImageType>(output);
 141 }
 142 //--------------------------------------------------------------------
 143
 144