ITK  4.3.0
Insight Segmentation and Registration Toolkit
itkDivideImageFilter.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright Insight Software Consortium
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * http://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 #ifndef __itkDivideImageFilter_h
19 #define __itkDivideImageFilter_h
20 
22 #include "itkNumericTraits.h"
23 
24 namespace itk
25 {
26 namespace Functor
27 {
33 template< class TInput1, class TInput2, class TOutput >
34 class Div
35 {
36 public:
37  Div() {}
38  ~Div() {}
39  bool operator!=(const Div &) const
40  {
41  return false;
42  }
44 
45  bool operator==(const Div & other) const
46  {
47  return !( *this != other );
48  }
49 
50  inline TOutput operator()(const TInput1 & A, const TInput2 & B) const
51  {
52  if ( B != (TInput2)0 )
53  {
54  return (TOutput)( A / B );
55  }
56  else
57  {
58  return NumericTraits< TOutput >::max( static_cast<TOutput>(A) );
59  }
60  }
61 };
62 }
80 template< class TInputImage1, class TInputImage2, class TOutputImage >
81 class ITK_EXPORT DivideImageFilter:
82  public
83  BinaryFunctorImageFilter< TInputImage1, TInputImage2, TOutputImage,
84  Functor::Div<
85  typename TInputImage1::PixelType,
86  typename TInputImage2::PixelType,
87  typename TOutputImage::PixelType > >
88 {
89 public:
94 
98  typedef BinaryFunctorImageFilter< TInputImage1, TInputImage2, TOutputImage,
100  typename TInputImage1::PixelType,
101  typename TInputImage2::PixelType,
102  typename TOutputImage::PixelType >
104 
110 
114  itkNewMacro(Self);
115 
117  itkTypeMacro(DivideImageFilter,
119 
120 #ifdef ITK_USE_CONCEPT_CHECKING
121 
122  itkConceptMacro( IntConvertibleToInput2Check,
124  itkConceptMacro( Input1Input2OutputDivisionOperatorsCheck,
125  ( Concept::DivisionOperators< typename TInputImage1::PixelType,
126  typename TInputImage2::PixelType,
127  typename TOutputImage::PixelType > ) );
128 
130 #endif
131 
132 protected:
134  virtual ~DivideImageFilter() {}
135 
136  void GenerateData()
137  {
138  const typename Superclass::DecoratedInput2ImagePixelType *input
139  = dynamic_cast< const typename Superclass::DecoratedInput2ImagePixelType * >(
140  this->ProcessObject::GetInput(1) );
142  {
143  itkGenericExceptionMacro(<<"The constant value used as denominator should not be set to zero");
144  }
145  else
146  {
147  Superclass::GenerateData();
148  }
149  }
150 
151 private:
152  DivideImageFilter(const Self &); //purposely not implemented
153  void operator=(const Self &); //purposely not implemented
154 
155 };
156 } // end namespace itk
157 
158 #endif
159