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
148
149
|
/*=========================================================================
*
* 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 itkHistogramToImageFilter_h
#define itkHistogramToImageFilter_h
#include "itkImageSource.h"
#include "itkConceptChecking.h"
#include "itkHistogram.h"
#include "itkImageRegionIteratorWithIndex.h"
#include "itkSimpleDataObjectDecorator.h"
namespace itk
{
/** \class HistogramToImageFilter
* \brief This class takes a histogram as an input and returns an image of
* type specified by the functor.
*
* The dimension of the image is equal to the size of each measurement
* vector of the histogram. The size in the image along each dimension will be
* equal to the number of bins along each dimension of the histogram.
*
* The filter may be used in registration methods to plot the joint histogram
* after every iteration. A functor is used since it is customary to plot
* p log p where p is the probability of each measurement vector
* p is given by Number of occurrences of the measurement vector / total number
* of occurrences of all measurement vectors.
*
* \sa HistogramToProbabilityImageFilter, HistogramToLogProbabilityImageFilter,
* HistogramToIntensityImageFilter, HistogramToEntropyImageFilter
*
* \ingroup ITKStatistics
*/
template <typename THistogram, typename TImage, typename TFunction>
class ITK_TEMPLATE_EXPORT HistogramToImageFilter : public ImageSource<TImage>
{
public:
ITK_DISALLOW_COPY_AND_MOVE(HistogramToImageFilter);
/** Standard class type aliases. */
using FunctorType = TFunction;
using Self = HistogramToImageFilter;
using Superclass = ImageSource<TImage>;
using Pointer = SmartPointer<Self>;
using ConstPointer = SmartPointer<const Self>;
using OutputImageType = TImage;
using OutputImagePointer = typename Superclass::Pointer;
using SpacingType = typename OutputImageType::SpacingType;
using PointType = typename OutputImageType::PointType;
using OutputPixelType = typename OutputImageType::PixelType;
// Define an iterator to iterate through the image
using ImageIteratorType = itk::ImageRegionIteratorWithIndex<OutputImageType>;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** \see LightObject::GetNameOfClass() */
itkOverrideGetNameOfClassMacro(HistogramToImageFilter);
/** Superclass type alias. */
using typename Superclass::OutputImageRegionType;
/** Some convenient type alias. */
using HistogramType = THistogram;
using MeasurementVectorType = typename HistogramType::MeasurementVectorType;
using HistogramSizeType = typename HistogramType::SizeType;
using SizeType = typename OutputImageType::SizeType;
/** Determine the image dimension. */
static constexpr unsigned int ImageDimension = OutputImageType::ImageDimension;
/** Set/Get the input histogram. */
using Superclass::SetInput;
virtual void
SetInput(const HistogramType * input);
const HistogramType *
GetInput();
/** Set the functor object. This replaces the current Functor with a
* copy of the specified Functor. This allows the user to specify a
* functor that has ivars set differently than the default functor.
* This method requires an operator!=() be defined on the functor
* (or the compiler's default implementation of operator!=() being
* appropriate). */
void
SetFunctor(const FunctorType & functor)
{
m_Functor = functor;
this->Modified();
}
/** Get the functor object. The functor is returned by reference.
* (Functors do not have to derive from itk::LightObject, so they do
* not necessarily have a reference count. So we cannot return a
* SmartPointer.) */
FunctorType &
GetFunctor()
{
return m_Functor;
}
const FunctorType &
GetFunctor() const
{
return m_Functor;
}
void
SetTotalFrequency(SizeValueType n);
protected:
HistogramToImageFilter();
~HistogramToImageFilter() override = default;
void
GenerateOutputInformation() override;
void
GenerateData() override;
FunctorType m_Functor{};
void
PrintSelf(std::ostream & os, Indent indent) const override;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
# include "itkHistogramToImageFilter.hxx"
#endif
#endif
|
