1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
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
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.
13 It is distributed under dual licence
15 - BSD See included LICENSE.txt file
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ===========================================================================**/
19 #ifndef CLITKREGIONGROWINGGENERICFILTER_TXX
20 #define CLITKREGIONGROWINGGENERICFILTER_TXX
22 #include <itkBinaryBallStructuringElement.h>
23 #include <itkConstShapedNeighborhoodIterator.h>
26 //--------------------------------------------------------------------
27 template<class ArgsInfoType>
28 clitk::RegionGrowingGenericFilter<ArgsInfoType>::RegionGrowingGenericFilter():
29 ImageToImageGenericFilter<Self>("RegionGrowing")
31 InitializeImageType<2>();
32 InitializeImageType<3>();
33 //InitializeImageType<4>();
35 //--------------------------------------------------------------------
38 //--------------------------------------------------------------------
39 template<class ArgsInfoType>
40 template<unsigned int Dim>
41 void clitk::RegionGrowingGenericFilter<ArgsInfoType>::InitializeImageType()
43 ADD_IMAGE_TYPE(Dim, uchar);
44 ADD_IMAGE_TYPE(Dim, short);
45 // ADD_IMAGE_TYPE(Dim, int);
46 ADD_IMAGE_TYPE(Dim, float);
48 //--------------------------------------------------------------------
51 //--------------------------------------------------------------------
52 template<class ArgsInfoType>
53 void clitk::RegionGrowingGenericFilter<ArgsInfoType>::SetArgsInfo(const ArgsInfoType & a)
56 this->SetIOVerbose(mArgsInfo.verbose_flag);
57 if (mArgsInfo.input_given) this->AddInputFilename(mArgsInfo.input_arg);
58 if (mArgsInfo.output_given) this->SetOutputFilename(mArgsInfo.output_arg);
60 //--------------------------------------------------------------------
63 //--------------------------------------------------------------------
64 // Update with the number of dimensions and the pixeltype
65 //--------------------------------------------------------------------
66 template<class ArgsInfoType>
67 template<class ImageType>
68 void clitk::RegionGrowingGenericFilter<ArgsInfoType>::UpdateWithInputImageType()
70 DD("UpdateWithInputImageType");
71 const int Dimension = ImageType::ImageDimension;
74 typedef ImageType InputImageType;
75 typedef ImageType OutputImageType;
76 typedef typename ImageType::PixelType PixelType;
79 typename ImageType::Pointer input = this->template GetInput<ImageType>(0);
82 typedef typename std::vector<typename InputImageType::IndexType> SeedsType;
84 if (mArgsInfo.seed_given==Dimension)
85 for (unsigned int i=0; i<Dimension;i++)
86 seeds[0][i]=mArgsInfo.seed_arg[i];
88 else if ( mArgsInfo.seed_given==1)
89 seeds[0].Fill(mArgsInfo.seed_arg[0]);
91 else seeds[0].Fill(mArgsInfo.seed_arg[0]);
92 if (mArgsInfo.verbose_flag) std::cout<<"Setting seed seeds to "<<seeds[0]<<"..."<<std::endl;
94 if (mArgsInfo.seedRadius_given)
96 typedef itk::BinaryBallStructuringElement<PixelType, Dimension> BallType;
97 typename BallType::RadiusType r;
99 if (mArgsInfo.seedRadius_given == Dimension)
100 for (unsigned i = 0; i < Dimension; i++)
101 r[i] = mArgsInfo.seedRadius_arg[i];
103 r.Fill(mArgsInfo.seedRadius_arg[0]);
107 ball.CreateStructuringElement();
109 typedef itk::ConstShapedNeighborhoodIterator<InputImageType> IteratorType;
110 IteratorType it(ball.GetRadius(),
112 input->GetLargestPossibleRegion());
113 #if ITK_VERSION_MAJOR < 4
114 typename BallType::ConstIterator nit;
116 for (nit = ball.Begin(); nit != ball.End(); ++nit, ++idx)
120 it.ActivateOffset(it.GetOffset(idx));
124 it.DeactivateOffset(it.GetOffset(idx));
128 it.CreateActiveListFromNeighborhood(ball);
129 it.NeedToUseBoundaryConditionOff();
132 it.SetLocation(seeds[0]);
133 for (typename IteratorType::ConstIterator i = it.Begin(); !i.IsAtEnd(); ++i)
135 typename InputImageType::IndexType id = seeds[0] + i.GetNeighborhoodOffset();
136 if (id != seeds[0] && input->GetLargestPossibleRegion().IsInside(id))
142 typedef itk::ImageToImageFilter<InputImageType, OutputImageType> ImageToImageFilterType;
143 typename ImageToImageFilterType::Pointer filter;
145 switch (mArgsInfo.type_arg)
149 typedef itk::ConnectedThresholdImageFilter<InputImageType, OutputImageType> ImageFilterType;
150 typename ImageFilterType::Pointer f= ImageFilterType::New();
152 f->SetLower(mArgsInfo.lower_arg);
153 f->SetUpper(mArgsInfo.upper_arg);
154 f->SetReplaceValue(static_cast<PixelType>(mArgsInfo.pad_arg));
155 for (typename SeedsType::const_iterator it = seeds.begin(); it != seeds.end(); ++it)
158 if (mArgsInfo.verbose_flag) std::cout<<"Using the connected threshold image filter..."<<std::endl;
165 typedef itk::NeighborhoodConnectedImageFilter<InputImageType, OutputImageType> ImageFilterType;
166 typename ImageFilterType::Pointer f= ImageFilterType::New();
169 typename InputImageType::SizeType size;
170 if (mArgsInfo.radius_given==Dimension)
171 for (unsigned int i=0; i<Dimension;i++)
172 size[i]=mArgsInfo.radius_arg[i];
174 else if ( mArgsInfo.radius_given==1)
175 size.Fill(mArgsInfo.radius_arg[0]);
177 else size.Fill(mArgsInfo.radius_arg[0]);
178 if (mArgsInfo.verbose_flag) std::cout<<"Setting neighborhood radius to "<<size<<"..."<<std::endl;
180 f->SetLower(mArgsInfo.lower_arg);
181 f->SetUpper(mArgsInfo.upper_arg);
182 f->SetReplaceValue(static_cast<PixelType>(mArgsInfo.pad_arg));
183 for (typename SeedsType::const_iterator it = seeds.begin(); it != seeds.end(); ++it)
187 if (mArgsInfo.verbose_flag) std::cout<<"Using the neighborhood threshold connected image filter..."<<std::endl;
194 typedef itk::ConfidenceConnectedImageFilter<InputImageType, OutputImageType> ImageFilterType;
195 typename ImageFilterType::Pointer f= ImageFilterType::New();
198 typename InputImageType::SizeType size;
199 if (mArgsInfo.radius_given==Dimension)
200 for (unsigned int i=0; i<Dimension;i++)
201 size[i]=mArgsInfo.radius_arg[i];
203 else if ( mArgsInfo.radius_given==1)
204 size.Fill(mArgsInfo.radius_arg[0]);
206 else size.Fill(mArgsInfo.radius_arg[0]);
207 if (mArgsInfo.verbose_flag) std::cout<<"Setting neighborhood radius to "<<size<<"..."<<std::endl;
209 f->SetMultiplier( mArgsInfo.multiplier_arg );
210 f->SetNumberOfIterations( mArgsInfo.multiplier_arg );
211 for (typename SeedsType::const_iterator it = seeds.begin(); it != seeds.end(); ++it)
213 f->SetNumberOfIterations( mArgsInfo.iter_arg);
214 f->SetReplaceValue(static_cast<PixelType>(mArgsInfo.pad_arg));
215 f->SetInitialNeighborhoodRadius(size[0]);
217 if (mArgsInfo.verbose_flag) std::cout<<"Using the confidence threshold connected image filter..."<<std::endl;
224 typedef clitk::LocallyAdaptiveThresholdConnectedImageFilter<InputImageType, OutputImageType> ImageFilterType;
225 typename ImageFilterType::Pointer f= ImageFilterType::New();
228 typename InputImageType::SizeType size;
229 if (mArgsInfo.radius_given==Dimension)
230 for (unsigned int i=0; i<Dimension;i++)
231 size[i]=mArgsInfo.radius_arg[i];
232 else size.Fill(mArgsInfo.radius_arg[0]);
233 if (mArgsInfo.verbose_flag) std::cout<<"Setting neighborhood radius to "<<size<<"..."<<std::endl;
236 f->SetLower(mArgsInfo.lower_arg);
237 f->SetUpper(mArgsInfo.upper_arg);
238 f->SetLowerBorderIsGiven(mArgsInfo.adaptLower_flag);
239 f->SetLowerBorderIsGiven(mArgsInfo.adaptUpper_flag);
240 f->SetReplaceValue(static_cast<PixelType>(mArgsInfo.pad_arg));
241 f->SetMultiplier(mArgsInfo.multiplier_arg);
242 f->SetMaximumSDIsGiven(mArgsInfo.maxSD_given);
243 if (mArgsInfo.maxSD_given) f->SetMaximumSD(mArgsInfo.maxSD_arg);
244 for (typename SeedsType::const_iterator it = seeds.begin(); it != seeds.end(); ++it)
248 if (mArgsInfo.verbose_flag) std::cout<<"Using the locally adaptive threshold connected image filter..."<<std::endl;
255 typedef clitk::ExplosionControlledThresholdConnectedImageFilter<InputImageType, OutputImageType> ImageFilterType;
256 typename ImageFilterType::Pointer f= ImageFilterType::New();
259 // typename InputImageType::SizeType size;
260 // if (mArgsInfo.radius_given==Dimension)
261 // for (unsigned int i=0; i<Dimension;i++)
262 // size[i]=mArgsInfo.radius_arg[i];
263 // else size.Fill(mArgsInfo.radius_arg[0]);
264 // if (mArgsInfo.verbose_flag) std::cout<<"Setting neighborhood radius to "<<size<<"..."<<std::endl;
267 f->SetVerbose(mArgsInfo.verbose_flag);
268 f->SetLower(mArgsInfo.lower_arg);
269 f->SetUpper(mArgsInfo.upper_arg);
270 f->SetMinimumLowerThreshold(mArgsInfo.minLower_arg);
271 f->SetMaximumUpperThreshold(mArgsInfo.maxUpper_arg);
272 f->SetAdaptLowerBorder(mArgsInfo.adaptLower_flag);
273 f->SetAdaptUpperBorder(mArgsInfo.adaptUpper_flag);
274 f->SetReplaceValue(static_cast<PixelType>(mArgsInfo.pad_arg));
275 f->SetMultiplier(mArgsInfo.multiplier_arg);
276 f->SetThresholdStepSize(mArgsInfo.step_arg);
277 f->SetMinimumThresholdStepSize(mArgsInfo.minStep_arg);
278 f->SetFullyConnected(mArgsInfo.full_flag);
279 for (typename SeedsType::const_iterator it = seeds.begin(); it != seeds.end(); ++it)
282 if (mArgsInfo.verbose_flag) std::cout<<"Using the explosion controlled threshold connected image filter..."<<std::endl;
290 filter->SetInput(input);
292 typename OutputImageType::Pointer output=filter->GetOutput();
294 // Write/Save results
295 this->template SetNextOutput<OutputImageType>(output);
297 //--------------------------------------------------------------------
300 #endif //#define CLITKREGIONGROWINGGENERICFILTER_TXX