packages/itk/src/bbitkImageProperties.cxx

   1 /*=========================================================================
   2   Program:   bbtk
   3   Module:    $RCSfile: bbitkImageProperties.cxx,v $
   4   Language:  C++
   5   Date:      $Date: 2011/03/04 08:33:59 $
   6   Version:   $Revision: 1.15 $
   7 =========================================================================*/
   8
   9 /* ---------------------------------------------------------------------
  10
  11 * Copyright (c) CREATIS-LRMN (Centre de Recherche en Imagerie Medicale)
  12 * Authors : Eduardo Davila, Laurent Guigues, Jean-Pierre Roux
  13 *
  14 *  This software is governed by the CeCILL-B license under French law and
  15 *  abiding by the rules of distribution of free software. You can  use,
  16 *  modify and/ or redistribute the software under the terms of the CeCILL-B
  17 *  license as circulated by CEA, CNRS and INRIA at the following URL
  18 *  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  19 *  or in the file LICENSE.txt.
  20 *
  21 *  As a counterpart to the access to the source code and  rights to copy,
  22 *  modify and redistribute granted by the license, users are provided only
  23 *  with a limited warranty  and the software's author,  the holder of the
  24 *  economic rights,  and the successive licensors  have only  limited
  25 *  liability.
  26 *
  27 *  The fact that you are presently reading this means that you have had
  28 *  knowledge of the CeCILL-B license and that you accept its terms.
  29 * ------------------------------------------------------------------------ */
  30
  31 /**
  32  *  \file
  33  *  \brief
  34  */
  35 #ifdef _USE_ITK_
  36
  37 #include "bbitkImageProperties.h"
  38 #include "bbitkPackage.h"
  39 #include "itkMinimumMaximumImageFilter.h"
  40
  41 namespace bbitk
  42 {
  43   BBTK_BLACK_BOX_IMPLEMENTATION(ImageProperties,bbtk::AtomicBlackBox);
  44
  45   BBTK_ADD_BLACK_BOX_TO_PACKAGE(itk,ImageProperties);
  46
  47         //-----------------------------------------------------------------
  48         void ImageProperties::bbUserSetDefaultValues()
  49         {
  50         }
  51
  52         //-----------------------------------------------------------------
  53         void ImageProperties::bbUserInitializeProcessing()
  54         {
  55         }
  56
  57         //-----------------------------------------------------------------
  58         void ImageProperties::bbUserFinalizeProcessing()
  59         {
  60         }
  61
  62   void ImageProperties::DoIt()
  63   {
  64     BBTK_TEMPLATE_ITK_IMAGE_SWITCH(bbGetInputIn().type(),DoIt);
  65   }
  66
  67   /**
  68       Template Processing
  69   */
  70   template<class itkImageType>
  71   void ImageProperties::DoIt()
  72   {
  73     bbtkDebugMessageInc("Core",9,"bbitk::ImageProperties::DoIt<"
  74                         <<bbtk::TypeName<itkImageType>()
  75                         <<">()"<<std::endl);
  76
  77     itkImageType* im = bbGetInputIn().get<itkImageType*>();
  78     unsigned int dim = im->GetImageDimension();
  79
  80     bbSetOutputTypeName(bbtk::TypeName<typename itkImageType::PixelType>());
  81     bbSetOutputDimension(dim);
  82     const typename itkImageType::RegionType& r =im->GetLargestPossibleRegion();
  83     bbSetOutputLargestPossibleRegion(r);
  84
  85     const typename itkImageType::IndexType& ind = r.GetIndex();
  86     std::vector<int> vind;
  87     for (unsigned int i=0;i<dim;++i)
  88        vind.push_back(ind[i]);
  89     bbSetOutputIndex(vind);
  90
  91     const typename itkImageType::SizeType& sz = r.GetSize();
  92     std::vector<int> vsz;
  93     for (unsigned int i=0;i<dim;++i)
  94        vsz.push_back(sz[i]);
  95
  96   // brute hack to avoid failure of most black boxes that expects 3D images. // JPR
  97    if (dim==2)
  98         vsz.push_back(1);
  99     bbSetOutputSize(vsz);
 100
 101     typename itkImageType::SpacingType sp = im->GetSpacing();
 102     std::vector<float> vsp;
 103     for (unsigned int i=0;i<dim;++i)
 104        vsp.push_back(sp[i]);
 105
 106  // brute hack to avoid failure of most black boxes that expects 3D images. // JPR
 107     if (dim==2)
 108        vsp.push_back(1.0);
 109
 110     bbSetOutputSpacing(vsp);
 111
 112
 113   // Min Max Value                     // EED
 114           std::vector<float> vminmax;
 115           typedef itk::MinimumMaximumImageFilter< itkImageType > MinMaxFilterType;
 116           typename MinMaxFilterType::Pointer minMaxCalculator = MinMaxFilterType::New();
 117           minMaxCalculator->SetInput(im);
 118           minMaxCalculator->Update();
 119           vminmax.push_back( (float) (minMaxCalculator->GetMinimum()) );
 120           vminmax.push_back( (float) (minMaxCalculator->GetMaximum()) );
 121           bbSetOutputMinMax(vminmax);
 122
 123     bbtkDebugDecTab("Core",9);
 124   }
 125 }
 126 // eo namespace bbtk
 127
 128 #endif