]> Creatis software - creaRigidRegistration.git/blob - lib/Convolution.cxx
#3271 creaRigidRegistration Bug New Normal - Convolution
[creaRigidRegistration.git] / lib / Convolution.cxx
1 /*
2 # ---------------------------------------------------------------------
3 #
4 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image 
5 #                        pour la Santé)
6 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
7 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
8 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
9 #
10 #  This software is governed by the CeCILL-B license under French law and 
11 #  abiding by the rules of distribution of free software. You can  use, 
12 #  modify and/ or redistribute the software under the terms of the CeCILL-B 
13 #  license as circulated by CEA, CNRS and INRIA at the following URL 
14 #  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html 
15 #  or in the file LICENSE.txt.
16 #
17 #  As a counterpart to the access to the source code and  rights to copy,
18 #  modify and redistribute granted by the license, users are provided only
19 #  with a limited warranty  and the software's author,  the holder of the
20 #  economic rights,  and the successive licensors  have only  limited
21 #  liability. 
22 #
23 #  The fact that you are presently reading this means that you have had
24 #  knowledge of the CeCILL-B license and that you accept its terms.
25 # ------------------------------------------------------------------------      */                                                                    
26
27 #include "Convolution.h"
28
29 /*
30 * Constructor
31 */
32 //------------------------------------------------------------
33 Convolution::Convolution()
34 {
35         _image=NULL;
36         _convolve = vtkImageConvolve::New();
37         _cast = vtkImageCast::New();
38 }
39
40 /*
41 * Destructor
42 */
43 //------------------------------------------------------------
44 Convolution::~Convolution()
45 {
46         if (_convolve   != NULL ) { _convolve->Delete();        }
47         if (_cast               != NULL ) { _cast->Delete();            }
48 }
49
50 vtkImageData *Convolution::getImage()
51 {
52         return _cast->GetOutput();
53 }
54
55 void Convolution::setImage(vtkImageData *image)
56 {
57         _image = image;
58 //EED 2017-01-01 Migration VTK7
59 #if VTK_MAJOR_VERSION <= 5
60         _convolve->SetInput(_image);
61 #else
62         _convolve->SetInputData(_image);
63 #endif
64 }
65
66 void Convolution::setFactor(double factor)
67 {
68
69 //      _factor = (factor/100.0)*5.0;   
70         _factor = factor;       
71 }
72
73 void Convolution::setOn(bool on)
74 {
75         _on = on;
76 }
77
78 void Convolution::Run()
79 {
80         if(_on==true)
81         {
82                 double kernel[] = {0.0,1.0,0.0,1.0,-_factor,1.0,0.0,1.0,0.0};
83                 _convolve->SetKernel3x3(kernel);
84                 _convolve->Update();
85 //EED 2017-01-01 Migration VTK7
86 #if VTK_MAJOR_VERSION <= 5
87                 _cast->SetInput(_convolve->GetOutput());
88 #else
89                 _cast->SetInputData(_convolve->GetOutput());
90 #endif
91                 //_cast->SetOutputScalarTypeToDouble();
92                 _cast->Update();
93         } else {
94 //EED 2017-01-01 Migration VTK7
95 #if VTK_MAJOR_VERSION <= 5
96                 _cast->SetInput(_image);
97 #else
98                 _cast->SetInputData(_image);
99 #endif
100                 _cast->SetOutputScalarType(_image->GetScalarType());
101                 _cast->Update();
102         }       
103 }