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://oncora1.lyon.fnclcc.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 -------------------------------------------------------------------
20 * @file clitkGuerreroVentilationGenericFilter.cxx
21 * @author Joël Schaerer
25 -------------------------------------------------------------------*/
27 #include "clitkGuerreroVentilationGenericFilter.h"
28 #include <itkBinaryGuerreroFilter.h>
29 #include <itkImageDuplicator.h>
30 #include <itkExtractImageFilter.h>
32 //--------------------------------------------------------------------
33 clitk::GuerreroVentilationGenericFilter::GuerreroVentilationGenericFilter()
34 :ImageToImageGenericFilter<Self>("GuerreroVentilationGenericFilter") {
36 InitializeImageType<2>();
37 InitializeImageType<3>();
39 //--------------------------------------------------------------------
41 //--------------------------------------------------------------------
42 template<unsigned int Dim>
43 void clitk::GuerreroVentilationGenericFilter::InitializeImageType() {
44 ADD_IMAGE_TYPE(Dim, short);
46 //--------------------------------------------------------------------
49 //--------------------------------------------------------------------
50 template<class ImageType>
51 void clitk::GuerreroVentilationGenericFilter::UpdateWithInputImageType() {
54 assert(mInputFilenames.size() == 2);
57 typedef ImageType InputImageType;
58 typename InputImageType::Pointer input = this->template GetInput<InputImageType>(0);
59 typename InputImageType::Pointer ref = this->template GetInput<InputImageType>(1);
61 typedef itk::Image<float,InputImageType::ImageDimension> OutputImageType;
62 // typename ImageType::Pointer input = clitk::readImage<ImageType>(mInputFilenames[0], mIOVerbose);
63 //typename ImageType::Pointer ref = clitk::readImage<ImageType>(mInputFilenames[1], mIOVerbose);
65 typedef itk::BinaryGuerreroFilter<ImageType,ImageType,OutputImageType> GFilterType;
66 typename GFilterType::Pointer filter = GFilterType::New();
67 filter->SetInput1(ref);
68 filter->SetInput2(input);
69 filter->SetBloodCorrectionFactor(blood_mass_factor);
70 filter->SetUseCorrectFormula(use_correct_formula);
73 this->SetNextOutput<OutputImageType>(filter->GetOutput());
74 //clitk::writeImage<OutputImageType>(filter->GetOutput(), mOutputFilename, mIOVerbose);
75 //std::cout << "Warning: removed " << filter->GetFunctor().aberant_voxels << " aberant voxels from the ventilation image"