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
|
/*=========================================================================
*
* 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.
*
*=========================================================================*/
// Disable warning for long symbol names in this file only
/*
* This is a test file for the itkImageMaskSpatialObject class.
* The supported pixel types does not include itkRGBPixel, itkRGBAPixel, etc...
* So far it only allows to manage images of simple types like unsigned short,
* unsigned int, or itk::Vector<...>.
*/
#include "itkImageRegionIterator.h"
#include "itkImageMaskSpatialObject.h"
#include "itkTestingMacros.h"
int
itkImageMaskSpatialObjectTest5(int, char *[])
{
constexpr unsigned int VDimension = 3;
using ImageMaskSpatialObject = itk::ImageMaskSpatialObject<VDimension>;
using PixelType = ImageMaskSpatialObject::PixelType;
using ImageType = ImageMaskSpatialObject::ImageType;
using Iterator = itk::ImageRegionIterator<ImageType>;
auto image = ImageType::New();
ImageType::SizeType size = { { 50, 50, 50 } };
ImageType::IndexType index = { { 0, 0, 0 } };
ImageType::RegionType region;
region.SetSize(size);
region.SetIndex(index);
image->SetRegions(region);
image->AllocateInitialized();
ImageType::RegionType insideRegion;
ImageType::SizeType insideSize = { { 30, 30, 30 } };
ImageType::IndexType insideIndex = { { 10, 10, 10 } };
insideRegion.SetSize(insideSize);
insideRegion.SetIndex(insideIndex);
Iterator it(image, insideRegion);
it.GoToBegin();
while (!it.IsAtEnd())
{
it.Set(itk::NumericTraits<PixelType>::max());
++it;
}
auto maskSO = ImageMaskSpatialObject::New();
ITK_EXERCISE_BASIC_OBJECT_METHODS(maskSO, ImageMaskSpatialObject, ImageSpatialObject);
maskSO->SetImage(image);
maskSO->Update();
Iterator itr(image, region);
itr.GoToBegin();
while (!itr.IsAtEnd())
{
const ImageType::IndexType constIndex = itr.GetIndex();
const bool reference = insideRegion.IsInside(constIndex);
ImageType::PointType point;
image->TransformIndexToPhysicalPoint(constIndex, point);
const bool test = maskSO->IsInsideInWorldSpace(point);
if (test != reference)
{
std::cerr << "Error in the evaluation of IsInside() - Part 1 " << std::endl;
std::cerr << "Index failed = " << constIndex << std::endl;
std::cerr << " Mask " << test << " != " << reference << std::endl;
return EXIT_FAILURE;
}
++itr;
}
// Repeat test using a mask value
maskSO->SetMaskValue(itk::NumericTraits<PixelType>::max());
maskSO->SetUseMaskValue(true);
itr.GoToBegin();
while (!itr.IsAtEnd())
{
const ImageType::IndexType constIndex = itr.GetIndex();
const bool reference = insideRegion.IsInside(constIndex);
ImageType::PointType point;
image->TransformIndexToPhysicalPoint(constIndex, point);
const bool test = maskSO->IsInsideInWorldSpace(point);
if (test != reference)
{
std::cerr << "Error in the evaluation of IsInside() - Part 2" << std::endl;
std::cerr << "Index failed = " << constIndex << std::endl;
std::cerr << " Mask " << test << " != " << reference << std::endl;
return EXIT_FAILURE;
}
++itr;
}
// The following should result in all IsInsideInWorldSpace calls
// returning false
maskSO->SetMaskValue(itk::NumericTraits<PixelType>::OneValue());
maskSO->SetUseMaskValue(true);
itr.GoToBegin();
while (!itr.IsAtEnd())
{
const ImageType::IndexType constIndex = itr.GetIndex();
const bool reference = insideRegion.IsInside(constIndex);
ImageType::PointType point;
image->TransformIndexToPhysicalPoint(constIndex, point);
const bool test = maskSO->IsInsideInWorldSpace(point);
if (test == true)
{
std::cerr << "Error in the evaluation of IsInside() - Part 3" << std::endl;
std::cerr << "Index failed = " << constIndex << std::endl;
std::cerr << " Mask " << test << " != " << reference << std::endl;
return EXIT_FAILURE;
}
++itr;
}
std::cout << "Test finished" << std::endl;
return EXIT_SUCCESS;
}
|
