]> Creatis software - bbtk.git/blob - packages/itk/src/bbitkImageProperties.cxx
no message
[bbtk.git] / packages / itk / src / bbitkImageProperties.cxx
1 /*=========================================================================                                                                               
2   Program:   bbtk
3   Module:    $RCSfile: bbitkImageProperties.cxx,v $
4   Language:  C++
5   Date:      $Date: 2011/02/28 06:56:57 $
6   Version:   $Revision: 1.14 $
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           //itk::MinimumMaximumImageFilter< TInputImage >
115           printf("EED ImageProperties::DoIt 1\n");
116           itk::MinimumMaximumImageFilter<itkImageType>   *minmaxfilter = itk::MinimumMaximumImageFilter< itkImageType >::New();
117           printf("EED ImageProperties::DoIt 2\n");
118           std::vector<float> vminmax;
119 /*        
120           minmaxfilter->SetInput(im);
121           printf("EED ImageProperties::DoIt 3\n");
122           minmaxfilter->Update();
123           printf("EED ImageProperties::DoIt 4\n");
124 //        vminmax.push_back( minmaxfilter->GetMinimumOutput()->Get() );
125 //        vminmax.push_back( minmaxfilter->GetMaximumOutput()->Get() );
126           vminmax.push_back( minmaxfilter->GetMinimum() );
127           printf("EED ImageProperties::DoIt 5\n");
128           vminmax.push_back( minmaxfilter->GetMaximum() );
129           printf("EED ImageProperties::DoIt 6\n");
130 */
131           vminmax.push_back( 22 );
132           vminmax.push_back( 33 );
133           
134           bbSetOutputMinMax(vminmax);
135           printf("EED ImageProperties::DoIt 7\n");
136           
137     bbtkDebugDecTab("Core",9);
138   }
139 }
140 // eo namespace bbtk
141
142 #endif