1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
/*=========================================================================
*
* Copyright NumFOCUS
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* https://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#ifndef itkInvertIntensityImageFilter_h
#define itkInvertIntensityImageFilter_h
#include "itkUnaryFunctorImageFilter.h"
namespace itk
{
namespace Functor
{
/**
* \class InvertIntensityTransform
* \brief
* \ingroup ITKImageIntensity
*/
template <typename TInput, typename TOutput>
class ITK_TEMPLATE_EXPORT InvertIntensityTransform
{
public:
using RealType = typename NumericTraits<TInput>::RealType;
InvertIntensityTransform() { m_Maximum = NumericTraits<TInput>::max(); }
~InvertIntensityTransform() = default;
void
SetMaximum(TOutput max)
{
m_Maximum = max;
}
bool
operator==(const InvertIntensityTransform & other) const
{
return m_Maximum == other.m_Maximum;
}
ITK_UNEQUAL_OPERATOR_MEMBER_FUNCTION(InvertIntensityTransform);
inline TOutput
operator()(const TInput & x) const
{
auto result = static_cast<TOutput>(m_Maximum - x);
return result;
}
private:
TInput m_Maximum;
};
} // end namespace Functor
/** \class InvertIntensityImageFilter
* \brief Invert the intensity of an image.
*
* InvertIntensityImageFilter inverts intensity of pixels by
* subtracting pixel value to a maximum value. The maximum value can
* be set with SetMaximum and defaults the maximum of input pixel
* type. This filter can be used to invert, for example, a binary
* image, a distance map, etc.
*
* \author Gaetan Lehmann. Biologie du Developpement et de la Reproduction, INRA de Jouy-en-Josas, France.
*
* \sa IntensityWindowingImageFilter ShiftScaleImageFilter
* \ingroup IntensityImageFilters MultiThreaded
*
* \ingroup ITKImageIntensity
*
* \sphinx
* \sphinxexample{Filtering/ImageIntensity/InvertImage,Invert Image}
* \endsphinx
*/
template <typename TInputImage, typename TOutputImage = TInputImage>
class ITK_TEMPLATE_EXPORT InvertIntensityImageFilter
: public UnaryFunctorImageFilter<
TInputImage,
TOutputImage,
Functor::InvertIntensityTransform<typename TInputImage::PixelType, typename TOutputImage::PixelType>>
{
public:
ITK_DISALLOW_COPY_AND_MOVE(InvertIntensityImageFilter);
/** Standard class type aliases. */
using Self = InvertIntensityImageFilter;
using Superclass = UnaryFunctorImageFilter<
TInputImage,
TOutputImage,
Functor::InvertIntensityTransform<typename TInputImage::PixelType, typename TOutputImage::PixelType>>;
using Pointer = SmartPointer<Self>;
using ConstPointer = SmartPointer<const Self>;
using OutputPixelType = typename TOutputImage::PixelType;
using InputPixelType = typename TInputImage::PixelType;
using RealType = typename NumericTraits<InputPixelType>::RealType;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** \see LightObject::GetNameOfClass() */
itkOverrideGetNameOfClassMacro(InvertIntensityImageFilter);
/** Set/Get the maximum intensity value for the inversion. */
itkSetMacro(Maximum, InputPixelType);
itkGetConstReferenceMacro(Maximum, InputPixelType);
void
PrintSelf(std::ostream & os, Indent indent) const override;
/** Process to execute before entering the multithreaded section */
void
BeforeThreadedGenerateData() override;
#ifdef ITK_USE_CONCEPT_CHECKING
// Begin concept checking
itkConceptMacro(InputHasNumericTraitsCheck, (Concept::HasNumericTraits<InputPixelType>));
// End concept checking
#endif
protected:
InvertIntensityImageFilter();
~InvertIntensityImageFilter() override = default;
private:
InputPixelType m_Maximum{};
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
# include "itkInvertIntensityImageFilter.hxx"
#endif
#endif
|
