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
150
151
152
153
154
155
156
157
|
/*=========================================================================
*
* 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 itkApproximateSignedDistanceMapImageFilter_h
#define itkApproximateSignedDistanceMapImageFilter_h
#include "itkFastChamferDistanceImageFilter.h"
#include "itkIsoContourDistanceImageFilter.h"
namespace itk
{
/**
* \class ApproximateSignedDistanceMapImageFilter
* \brief Create a map of the approximate signed distance from the boundaries of
* a binary image.
*
* The ApproximateSignedDistanceMapImageFilter takes as input a binary image and
* produces a signed distance map. Each pixel value in the output contains the
* approximate distance from that pixel to the nearest "object" in the binary
* image. This filter differs from the DanielssonDistanceMapImageFilter in that
* it calculates the distance to the "object edge" for pixels within the object.
*
* Negative values in the output indicate that the pixel at that position is
* within an object in the input image. The absolute value of a negative pixel
* represents the approximate distance to the nearest object boundary pixel.
*
* WARNING: This filter requires that the output type be floating-point. Otherwise
* internal calculations will not be performed to the appropriate precision,
* resulting in completely incorrect (read: zero-valued) output.
*
* The distances computed by this filter are Chamfer distances, which are only
* an approximation to Euclidean distances, and are not as exact approximations
* as those calculated by the DanielssonDistanceMapImageFilter. On the other hand,
* this filter is faster.
*
* This filter requires that an "inside value" and "outside value" be set as
* parameters. The "inside value" is the intensity value of the binary image
* which corresponds to objects, and the "outside value" is the intensity of the
* background. (A typical binary image often represents objects as black (0) and
* background as white (usually 255), or vice-versa.) Note that this filter is
* slightly faster if the inside value is less than the outside value. Otherwise
* an extra iteration through the image is required.
*
* This filter uses the FastChamferDistanceImageFilter and the
* IsoContourDistanceImageFilter internally to perform the distance calculations.
*
* \sa DanielssonDistanceMapImageFilter
* \sa SignedDanielssonDistanceMapImageFilter
* \sa SignedMaurerDistanceMapImageFilter
* \sa FastChamferDistanceImageFilter
* \sa IsoContourDistanceImageFilter
*
* \author Zach Pincus
* \ingroup ITKDistanceMap
*
* \sphinx
* \sphinxexample{Filtering/DistanceMap/ApproxDistanceMapOfBinary, Approximate Distance Map Of Binary Image}
* \endsphinx
*/
template <typename TInputImage, typename TOutputImage>
class ITK_TEMPLATE_EXPORT ApproximateSignedDistanceMapImageFilter : public ImageToImageFilter<TInputImage, TOutputImage>
{
public:
ITK_DISALLOW_COPY_AND_MOVE(ApproximateSignedDistanceMapImageFilter);
/** Standard type alias */
using Self = ApproximateSignedDistanceMapImageFilter;
using Superclass = ImageToImageFilter<TInputImage, TOutputImage>;
using Pointer = SmartPointer<Self>;
using ConstPointer = SmartPointer<const Self>;
/** \see LightObject::GetNameOfClass() */
itkOverrideGetNameOfClassMacro(ApproximateSignedDistanceMapImageFilter);
/** standard New() method support */
itkNewMacro(Self);
/** Type for input image. */
using InputImageType = TInputImage;
/** Type for the output image. */
using OutputImageType = TOutputImage;
/** Type for the pixels of the input image. */
using InputPixelType = typename InputImageType::PixelType;
/** Type for the pixels of the output image. */
using OutputPixelType = typename OutputImageType::PixelType;
/** Type of input image size and size value */
using OutputSizeType = typename OutputImageType::SizeType;
using OutputSizeValueType = typename OutputSizeType::SizeValueType;
/** The dimension of the input image. */
static constexpr unsigned int InputImageDimension = InputImageType::ImageDimension;
/** Pointer Type for input image. */
using InputImagePointer = typename InputImageType::ConstPointer;
/** Pointer Type for the output image. */
using OutputImagePointer = typename OutputImageType::Pointer;
/** Set/Get intensity value representing the interior of objects in the mask.
*/
itkSetMacro(InsideValue, InputPixelType);
itkGetConstMacro(InsideValue, InputPixelType);
/** Set/Get intensity value representing non-objects in the mask. */
itkSetMacro(OutsideValue, InputPixelType);
itkGetConstMacro(OutsideValue, InputPixelType);
#ifdef ITK_USE_CONCEPT_CHECKING
// Begin concept checking
itkConceptMacro(InputEqualityComparableCheck, (Concept::EqualityComparable<typename InputImageType::PixelType>));
// End concept checking
#endif
protected:
ApproximateSignedDistanceMapImageFilter();
~ApproximateSignedDistanceMapImageFilter() override = default;
void
GenerateData() override;
void
PrintSelf(std::ostream & os, Indent indent) const override;
private:
using IsoContourType = IsoContourDistanceImageFilter<InputImageType, OutputImageType>;
using ChamferType = FastChamferDistanceImageFilter<OutputImageType, OutputImageType>;
typename IsoContourType::Pointer m_IsoContourFilter{};
typename ChamferType::Pointer m_ChamferFilter{};
InputPixelType m_InsideValue{};
InputPixelType m_OutsideValue{};
};
} // end of namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
# include "itkApproximateSignedDistanceMapImageFilter.hxx"
#endif
#endif
|
