Mercurial > hg > orthanc-stone

comparison OrthancStone/UnitTestsSources/TestStructureSet.cpp @ 1908:affde38b84de

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
moved tentative bgo reimplementation of rt-struct into graveyard
author Sebastien Jodogne <s.jodogne@gmail.com>
date Tue, 01 Feb 2022 08:38:32 +0100
parents 14c8f339d480
children 07964689cb0b
comparison
equal deleted inserted replaced
1907:0208f99b8bde 1908:affde38b84de
19 * License along with this program. If not, see 19 * License along with this program. If not, see
20 * <http://www.gnu.org/licenses/>. 20 * <http://www.gnu.org/licenses/>.
21 **/ 21 **/
22 22
23 23
24 // working around a bug where the Visual C++ compiler would get
25 // stuck trying to compile this cpp file in release mode
26 // (versions: https://en.wikipedia.org/wiki/Microsoft_Visual_C%2B%2B)
27 #ifdef _MSC_VER
28 # pragma optimize("", off)
29 // warning C4748: /GS can not protect parameters and local variables from
30 // local buffer overrun because optimizations are disabled in function
31 # pragma warning(disable: 4748)
32 #endif
33
34 #include "../Sources/Toolbox/DicomStructureSet.h" 24 #include "../Sources/Toolbox/DicomStructureSet.h"
35 #include "../Sources/Toolbox/DicomStructureSet2.h"
36 #include "../Sources/Toolbox/DicomStructureSetUtils.h"
37 #include "../Sources/Toolbox/DisjointDataSet.h"
38 25
39 #include <EmbeddedResources.h> 26 #include <EmbeddedResources.h>
40 27
41 #include <boost/date_time/posix_time/posix_time.hpp>
42 #include <boost/make_shared.hpp>
43
44 #include <gtest/gtest.h> 28 #include <gtest/gtest.h>
45 #include <string>
46
47 using namespace OrthancStone;
48
49 static const double DELTA_MAX = 10.0 * std::numeric_limits<float>::epsilon();
50
51
52
53 #define STONE_ARRAY_SIZE(arr) (sizeof(arr) / sizeof((arr)[0]))
54
55 #ifdef BGO_ENABLE_DICOMSTRUCTURESETLOADER2
56
57 static void CheckGroundTruth(
58 const std::vector<DicomStructure2>& structures,
59 const size_t structureIndex,
60 const size_t sliceIndex,
61 std::vector<double> groundTruth)
62 {
63 const std::vector<DicomStructurePolygon2>& polygonsForThisStruct = structures.at(structureIndex).GetPolygons();
64 const DicomStructurePolygon2& polygon = polygonsForThisStruct.at(sliceIndex);
65
66 //double groundTruth[] = { 7.657838, 108.2725, 304.01, 6.826687, 107.4413, 304.01, 6.152492, 106.4785, 304.01, 5.655735, 105.4132, 304.01, 5.351513, 104.2778, 304.01, 5.249068, 103.1069, 304.01, 5.351513, 101.9359, 304.01, 5.655735, 100.8005, 304.01, 6.152492, 99.73524, 304.01, 6.826687, 98.77239, 304.01, 7.657838, 97.94124, 304.01, 8.620689, 97.26704, 304.01, 9.685987, 96.77029, 304.01, 10.82136, 96.46606, 304.01, 11.99231, 96.36362, 304.01, 13.16326, 96.46606, 304.01, 14.29864, 96.77029, 304.01, 15.36393, 97.26704, 304.01, 16.32678, 97.94124, 304.01, 17.15794, 98.77239, 304.01, 17.83213, 99.73524, 304.01, 18.32889, 100.8005, 304.01, 18.63311, 101.9359, 304.01, 18.73555, 103.1069, 304.01, 18.63311, 104.2778, 304.01, 18.32889, 105.4132, 304.01, 17.83213, 106.4785, 304.01, 17.15794, 107.4413, 304.01, 16.32678, 108.2725, 304.01, 15.36393, 108.9467, 304.01, 14.29864, 109.4434, 304.01, 13.16326, 109.7477, 304.01, 11.99231, 109.8501, 304.01, 10.82136, 109.7477, 304.01, 9.685987, 109.4434, 304.01, 8.620689, 108.9467, 304.01 };
67 size_t groundTruthItems = groundTruth.size();
68
69 size_t pointCount = 3 * polygon.GetPointCount();
70
71 EXPECT_EQ(groundTruthItems, pointCount);
72
73 for (size_t i = 0; i < polygon.GetPointCount(); ++i)
74 {
75 const Vector& point = polygon.GetPoint(i);
76
77 // loop over X, Y then Z.
78 for (size_t j = 0; j < 3; ++j)
79 {
80 size_t index = 3 * i + j;
81 ASSERT_LT(index, groundTruthItems);
82 bool isNear = LinearAlgebra::IsNear(groundTruth[index], point[j]);
83 EXPECT_TRUE(isNear);
84 }
85 }
86 }
87
88
89 #include <Toolbox.h>
90
91 TEST(StructureSet2, ReadFromJsonThatsAll)
92 {
93 /*
94 The "RT_STRUCT_00" string is the reply to the following Orthanc request:
95
96 http://localhost:8042/instances/1aa5f84b-c32a03b4-3c1857da-da2e69f3-3ef6e2b3/tags?ignore-length=3006-0050
97
98 The tag hierarchy can be found here: https://dicom.innolitics.com/ciods/rt-dose
99 */
100
101 DicomStructureSet2 structureSet;
102
103 FullOrthancDataset dicom(Orthanc::EmbeddedResources::GetFileResourceBuffer(Orthanc::EmbeddedResources::RT_STRUCT_00),
104 Orthanc::EmbeddedResources::GetFileResourceSize(Orthanc::EmbeddedResources::RT_STRUCT_00));
105 structureSet.Clear();
106
107 structureSet.FillStructuresFromDataset(dicom);
108 structureSet.ComputeDependentProperties();
109
110 const std::vector<DicomStructure2>& structures = structureSet.structures_;
111
112 /*
113
114 ██████╗ █████╗ ███████╗██╗ ██████╗ ██████╗██╗ ██╗███████╗ ██████╗██╗ ██╗███████╗
115 ██╔══██╗██╔══██╗██╔════╝██║██╔════╝ ██╔════╝██║ ██║██╔════╝██╔════╝██║ ██╔╝██╔════╝
116 ██████╔╝███████║███████╗██║██║ ██║ ███████║█████╗ ██║ █████╔╝ ███████╗
117 ██╔══██╗██╔══██║╚════██║██║██║ ██║ ██╔══██║██╔══╝ ██║ ██╔═██╗ ╚════██║
118 ██████╔╝██║ ██║███████║██║╚██████╗ ╚██████╗██║ ██║███████╗╚██████╗██║ ██╗███████║
119 ╚═════╝ ╚═╝ ╚═╝╚══════╝╚═╝ ╚═════╝ ╚═════╝╚═╝ ╚═╝╚══════╝ ╚═════╝╚═╝ ╚═╝╚══════╝
120 http://patorjk.com/software/taag/#p=display&f=ANSI%20Shadow&t=BASIC%20CHECKS
121 */
122
123 // (0x3006, 0x0080) seq. size
124 EXPECT_EQ(7u, structures.size());
125
126 // (0x3006, 0x0080)[i]/(0x3006, 0x00a4)
127 for (size_t i = 0; i < 5; ++i)
128 {
129 EXPECT_EQ(std::string("ORGAN"), structures[i].interpretation_);
130 }
131 EXPECT_EQ(std::string("EXTERNAL"), structures[5].interpretation_);
132 EXPECT_EQ(std::string("PTV"), structures[6].interpretation_);
133
134 // (0x3006, 0x0020)[i]/(0x3006, 0x0026)
135 EXPECT_EQ(std::string("LN300"), structures[0].name_);
136 EXPECT_EQ(std::string("Cortical Bone"), structures[1].name_);
137 EXPECT_EQ(std::string("Adipose"), structures[2].name_);
138 EXPECT_EQ(std::string("CB2-50%"), structures[3].name_);
139 EXPECT_EQ(std::string("Water"), structures[4].name_);
140 EXPECT_EQ(std::string("External"), structures[5].name_);
141 EXPECT_EQ(std::string("PTV"), structures[6].name_);
142
143 // (0x3006, 0x0039)[i]/(0x3006, 0x002a)
144 EXPECT_EQ(0xff, structures[0].red_);
145 EXPECT_EQ(0x00, structures[0].green_);
146 EXPECT_EQ(0x00, structures[0].blue_);
147
148 EXPECT_EQ(0x00, structures[1].red_);
149 EXPECT_EQ(0xff, structures[1].green_);
150 EXPECT_EQ(0xff, structures[1].blue_);
151
152 // ...
153
154 EXPECT_EQ(0x00, structures[5].red_);
155 EXPECT_EQ(0x80, structures[5].green_);
156 EXPECT_EQ(0x00, structures[5].blue_);
157
158 EXPECT_EQ(0xff, structures[6].red_);
159 EXPECT_EQ(0x00, structures[6].green_);
160 EXPECT_EQ(0xff, structures[6].blue_);
161
162 /*
163
164 ██████╗ ███████╗ ██████╗ ███╗ ███╗███████╗████████╗██████╗ ██╗ ██╗
165 ██╔════╝ ██╔════╝██╔═══██╗████╗ ████║██╔════╝╚══██╔══╝██╔══██╗╚██╗ ██╔╝
166 ██║ ███╗█████╗ ██║ ██║██╔████╔██║█████╗ ██║ ██████╔╝ ╚████╔╝
167 ██║ ██║██╔══╝ ██║ ██║██║╚██╔╝██║██╔══╝ ██║ ██╔══██╗ ╚██╔╝
168 ╚██████╔╝███████╗╚██████╔╝██║ ╚═╝ ██║███████╗ ██║ ██║ ██║ ██║
169 ╚═════╝ ╚══════╝ ╚═════╝ ╚═╝ ╚═╝╚══════╝ ╚═╝ ╚═╝ ╚═╝ ╚═╝
170 http://patorjk.com/software/taag/#p=display&f=ANSI%20Shadow&t=BASIC%20CHECKS
171 */
172
173
174 {
175 double groundTruthRaw[] = { 7.657838, 108.2725, 304.01, 6.826687, 107.4413, 304.01, 6.152492, 106.4785, 304.01, 5.655735, 105.4132, 304.01, 5.351513, 104.2778, 304.01, 5.249068, 103.1069, 304.01, 5.351513, 101.9359, 304.01, 5.655735, 100.8005, 304.01, 6.152492, 99.73524, 304.01, 6.826687, 98.77239, 304.01, 7.657838, 97.94124, 304.01, 8.620689, 97.26704, 304.01, 9.685987, 96.77029, 304.01, 10.82136, 96.46606, 304.01, 11.99231, 96.36362, 304.01, 13.16326, 96.46606, 304.01, 14.29864, 96.77029, 304.01, 15.36393, 97.26704, 304.01, 16.32678, 97.94124, 304.01, 17.15794, 98.77239, 304.01, 17.83213, 99.73524, 304.01, 18.32889, 100.8005, 304.01, 18.63311, 101.9359, 304.01, 18.73555, 103.1069, 304.01, 18.63311, 104.2778, 304.01, 18.32889, 105.4132, 304.01, 17.83213, 106.4785, 304.01, 17.15794, 107.4413, 304.01, 16.32678, 108.2725, 304.01, 15.36393, 108.9467, 304.01, 14.29864, 109.4434, 304.01, 13.16326, 109.7477, 304.01, 11.99231, 109.8501, 304.01, 10.82136, 109.7477, 304.01, 9.685987, 109.4434, 304.01, 8.620689, 108.9467, 304.01 };
176 size_t n = sizeof(groundTruthRaw) / sizeof(groundTruthRaw[0]);
177 std::vector<double> groundTruth(groundTruthRaw, groundTruthRaw+n);
178 CheckGroundTruth(structures, 0, 0, groundTruth);
179 }
180 {
181 double groundTruthRaw[] = { 7.657838, 108.2725, 310.01, 6.826687, 107.4413, 310.01, 6.152492, 106.4785, 310.01, 5.655735, 105.4132, 310.01, 5.351513, 104.2778, 310.01, 5.249068, 103.1069, 310.01, 5.351513, 101.9359, 310.01, 5.655735, 100.8005, 310.01, 6.152492, 99.73524, 310.01, 6.826687, 98.77239, 310.01, 7.657838, 97.94124, 310.01, 8.620689, 97.26704, 310.01, 9.685987, 96.77029, 310.01, 10.82136, 96.46606, 310.01, 11.99231, 96.36362, 310.01, 13.16326, 96.46606, 310.01, 14.29864, 96.77029, 310.01, 15.36393, 97.26704, 310.01, 16.32678, 97.94124, 310.01, 17.15794, 98.77239, 310.01, 17.83213, 99.73524, 310.01, 18.32889, 100.8005, 310.01, 18.63311, 101.9359, 310.01, 18.73555, 103.1069, 310.01, 18.63311, 104.2778, 310.01, 18.32889, 105.4132, 310.01, 17.83213, 106.4785, 310.01, 17.15794, 107.4413, 310.01, 16.32678, 108.2725, 310.01, 15.36393, 108.9467, 310.01, 14.29864, 109.4434, 310.01, 13.16326, 109.7477, 310.01, 11.99231, 109.8501, 310.01, 10.82136, 109.7477, 310.01, 9.685987, 109.4434, 310.01, 8.620689, 108.9467, 310.01 };
182 size_t n = sizeof(groundTruthRaw) / sizeof(groundTruthRaw[0]);
183 std::vector<double> groundTruth(groundTruthRaw, groundTruthRaw+n);
184 CheckGroundTruth(structures, 0, 2, groundTruth);
185 }
186 {
187 double groundTruthRaw[] = { -37.967, 161.9664, 304.01, -39.10237, 161.6622, 304.01, -40.16767, 161.1655, 304.01, -41.13052, 160.4913, 304.01, -41.96167, 159.6601, 304.01, -42.63587, 158.6973, 304.01, -43.13263, 157.632, 304.01, -43.43685, 156.4966, 304.01, -43.53929, 155.3257, 304.01, -43.43685, 154.1547, 304.01, -43.13263, 153.0193, 304.01, -42.63587, 151.954, 304.01, -41.96167, 150.9912, 304.01, -41.13052, 150.16, 304.01, -40.16767, 149.4858, 304.01, -39.10237, 148.9891, 304.01, -37.967, 148.6849, 304.01, -36.79605, 148.5824, 304.01, -35.6251, 148.6849, 304.01, -34.48972, 148.9891, 304.01, -33.42443, 149.4858, 304.01, -32.46157, 150.16, 304.01, -31.63042, 150.9912, 304.01, -30.95623, 151.954, 304.01, -30.45947, 153.0193, 304.01, -30.15525, 154.1547, 304.01, -30.0528, 155.3257, 304.01, -30.15525, 156.4966, 304.01, -30.45947, 157.632, 304.01, -30.95623, 158.6973, 304.01, -31.63042, 159.6601, 304.01, -32.46157, 160.4913, 304.01, -33.42443, 161.1655, 304.01, -34.48972, 161.6622, 304.01, -35.6251, 161.9664, 304.01, -36.79605, 162.0689, 304.01 };
188 size_t n = sizeof(groundTruthRaw) / sizeof(groundTruthRaw[0]);
189 std::vector<double> groundTruth(groundTruthRaw, groundTruthRaw+n);
190 CheckGroundTruth(structures, 1, 0, groundTruth);
191 }
192 {
193 double groundTruthRaw[] = { 69.4042, 150.7324, 307.01, 69.70842, 151.8678, 307.01, 69.81087, 153.0387, 307.01, 69.70842, 154.2097, 307.01, 69.4042, 155.345, 307.01, 68.90745, 156.4103, 307.01, 68.23325, 157.3732, 307.01, 67.4021, 158.2043, 307.01, 66.43925, 158.8785, 307.01, 65.37395, 159.3753, 307.01, 64.23858, 159.6795, 307.01, 63.06762, 159.7819, 307.01, 61.89667, 159.6795, 307.01, 60.7613, 159.3753, 307.01, 59.696, 158.8785, 307.01, 58.73315, 158.2043, 307.01, 57.902, 157.3732, 307.01, 57.22781, 156.4103, 307.01, 56.73105, 155.345, 307.01, 56.42683, 154.2097, 307.01, 56.32438, 153.0387, 307.01, 56.42683, 151.8678, 307.01, 56.73105, 150.7324, 307.01, 57.22781, 149.6671, 307.01, 57.902, 148.7042, 307.01, 58.73315, 147.8731, 307.01, 59.696, 147.1989, 307.01, 60.7613, 146.7021, 307.01, 61.89667, 146.3979, 307.01, 63.06762, 146.2955, 307.01, 64.23858, 146.3979, 307.01, 65.37395, 146.7021, 307.01, 66.43925, 147.1989, 307.01, 67.4021, 147.8731, 307.01, 68.23325, 148.7042, 307.01, 68.90745, 149.6671, 307.01 };
194 size_t n = sizeof(groundTruthRaw) / sizeof(groundTruthRaw[0]);
195 std::vector<double> groundTruth(groundTruthRaw, groundTruthRaw+n);
196 CheckGroundTruth(structures, 2, 1, groundTruth);
197 }
198
199 {
200 double groundTruthRaw[] = { 108.3984, 232.7406, 274.01, 106.0547, 231.7948, 274.01, 103.7109, 232.8407, 274.01, 96.67969, 232.8757, 274.01, 77.92969, 232.887, 274.01, 47.46094, 232.8902, 274.01, 38.08594, 232.7537, 274.01, 37.6668, 232.3734, 274.01, 38.08594, 231.9774, 274.01, 40.42969, 231.8475, 274.01, 41.76413, 230.0297, 274.01, 42.77344, 229.1388, 274.01, 45.11719, 228.5069, 274.01, 47.46094, 227.1533, 274.01, 49.80469, 226.3505, 274.01, 52.14844, 224.6564, 274.01, 54.49219, 223.923, 274.01, 56.83594, 222.0692, 274.01, 59.17969, 220.3438, 274.01, 61.52344, 219.3888, 274.01, 63.86719, 217.1287, 274.01, 65.83488, 215.9672, 274.01, 68.55469, 213.2383, 274.01, 70.89844, 211.2328, 274.01, 72.8125, 208.9359, 274.01, 75.58594, 206.3615, 274.01, 76.91445, 204.2484, 274.01, 78.89509, 201.9047, 274.01, 80.51276, 199.5609, 274.01, 81.51955, 197.2172, 274.01, 83.67448, 194.8734, 274.01, 84.60938, 192.5297, 274.01, 85.86986, 190.1859, 274.01, 86.57623, 187.8422, 274.01, 88.30051, 185.4984, 274.01, 88.94002, 183.1547, 274.01, 89.23261, 180.8109, 274.01, 89.64844, 180.3263, 274.01, 90.71885, 178.4672, 274.01, 90.97656, 176.1234, 274.01, 91.99219, 174.4794, 274.01, 92.56773, 173.7797, 274.01, 92.80016, 171.4359, 274.01, 93.23473, 169.0922, 274.01, 93.37606, 166.7484, 274.01, 93.60748, 157.3734, 274.01, 93.6341, 152.6859, 274.01, 93.35742, 140.9672, 274.01, 92.89317, 138.6234, 274.01, 92.7069, 136.2797, 274.01, 92.03726, 133.9359, 274.01, 90.84009, 131.5922, 274.01, 90.3769, 129.2484, 274.01, 89.09074, 126.9047, 274.01, 88.13225, 122.2172, 274.01, 86.17828, 119.8734, 274.01, 84.96094, 117.4163, 274.01, 83.99619, 115.1859, 274.01, 83.13079, 112.8422, 274.01, 82.61719, 112.2984, 274.01, 80.27344, 108.8454, 274.01, 79.64514, 108.1547, 274.01, 77.21497, 105.8109, 274.01, 76.47787, 103.4672, 274.01, 75.58594, 102.6177, 274.01, 73.24219, 100.0077, 274.01, 69.54492, 96.43594, 274.01, 67.34096, 94.09219, 274.01, 64.66306, 91.74844, 274.01, 63.86719, 90.92619, 274.01, 61.52344, 90.20454, 274.01, 59.17969, 87.78574, 274.01, 56.83594, 86.48566, 274.01, 54.49219, 84.31388, 274.01, 52.14844, 83.44438, 274.01, 49.80469, 82.75121, 274.01, 49.37617, 82.37344, 274.01, 47.46094, 81.26244, 274.01, 45.71391, 80.02969, 274.01, 45.11719, 79.45415, 274.01, 42.77344, 79.08185, 274.01, 40.42969, 78.51941, 274.01, 38.08594, 78.27534, 274.01, 37.36932, 77.68594, 274.01, 35.74219, 76.67624, 274.01, 33.39844, 76.49941, 274.01, 31.05469, 76.03495, 274.01, 28.71094, 74.83174, 274.01, 26.36719, 74.62859, 274.01, 24.02344, 74.55463, 274.01, 21.67969, 74.22861, 274.01, 19.33594, 74.05312, 274.01, 12.30469, 73.99397, 274.01, 5.273438, 74.0736, 274.01, 2.929688, 74.55463, 274.01, 0.5859375, 74.68513, 274.01, -1.757813, 74.914, 274.01, -2.319131, 75.34219, 274.01, -4.101563, 76.31516, 274.01, -8.789063, 76.74514, 274.01, -11.13281, 78.39038, 274.01, -13.47656, 78.6124, 274.01, -15.82031, 79.19784, 274.01, -18.16406, 81.11024, 274.01, -20.50781, 82.03296, 274.01, -22.85156, 83.13991, 274.01, -25.19531, 83.70732, 274.01, -27.53906, 85.85863, 274.01, -29.88281, 87.03368, 274.01, -32.22656, 88.3274, 274.01, -34.57031, 90.53674, 274.01, -36.91406, 92.5602, 274.01, -39.25781, 93.55952, 274.01, -41.60156, 95.74537, 274.01, -43.94531, 98.26609, 274.01, -46.28906, 100.3701, 274.01, -47.02621, 101.1234, 274.01, -47.86611, 103.4672, 274.01, -49.83594, 105.8109, 274.01, -51.98182, 108.1547, 274.01, -53.06448, 110.4984, 274.01, -53.32031, 110.7675, 274.01, -54.53804, 112.8422, 274.01, -55.66406, 114.273, 274.01, -56.55722, 115.1859, 274.01, -57.13953, 117.5297, 274.01, -58.29264, 119.8734, 274.01, -59.26869, 122.2172, 274.01, -60.35156, 124.0119, 274.01, -60.84229, 124.5609, 274.01, -61.54484, 126.9047, 274.01, -61.71691, 129.2484, 274.01, -63.62281, 131.5922, 274.01, -63.81256, 133.9359, 274.01, -64.12511, 136.2797, 274.01, -64.84515, 138.6234, 274.01, -65.13599, 140.9672, 274.01, -65.33604, 143.3109, 274.01, -65.87358, 145.6547, 274.01, -66.10577, 147.9984, 274.01, -66.17618, 155.0297, 274.01, -66.09933, 162.0609, 274.01, -65.40382, 164.4047, 274.01, -65.24833, 166.7484, 274.01, -64.71442, 171.4359, 274.01, -63.88171, 173.7797, 274.01, -63.69299, 176.1234, 274.01, -61.79081, 178.4672, 274.01, -61.59269, 180.8109, 274.01, -61.19405, 183.1547, 274.01, -60.35156, 185.2055, 274.01, -59.08288, 187.8422, 274.01, -58.00781, 189.3499, 274.01, -57.25858, 190.1859, 274.01, -56.64558, 192.5297, 274.01, -55.29191, 194.8734, 274.01, -54.28698, 197.2172, 274.01, -52.28595, 199.5609, 274.01, -51.47569, 201.9047, 274.01, -48.63281, 204.6417, 274.01, -47.10181, 206.5922, 274.01, -44.64154, 208.9359, 274.01, -42.38504, 211.2797, 274.01, -39.25781, 214.4025, 274.01, -37.42723, 215.9672, 274.01, -34.57031, 218.9107, 274.01, -32.22656, 219.7277, 274.01, -29.88281, 221.6934, 274.01, -27.53906, 222.852, 274.01, -25.19531, 224.5168, 274.01, -22.85156, 225.9419, 274.01, -20.50781, 226.7359, 274.01, -18.16406, 228.3332, 274.01, -15.82031, 229.065, 274.01, -13.47656, 229.267, 274.01, -12.63854, 230.0297, 274.01, -11.13281, 231.9201, 274.01, -10.65505, 232.3734, 274.01, -11.13281, 232.7794, 274.01, -15.82031, 232.792, 274.01, -18.16406, 232.8902, 274.01, -36.91406, 232.9015, 274.01, -39.25781, 232.8902, 274.01, -50.97656, 232.9236, 274.01, -60.35156, 232.9126, 274.01, -67.38281, 232.8407, 274.01, -72.07031, 232.8642, 274.01, -79.10156, 232.8555, 274.01, -83.78906, 232.8788, 274.01, -95.50781, 232.8902, 274.01, -97.85156, 233.4886, 274.01, -100.1953, 233.647, 274.01, -102.5391, 232.9858, 274.01, -104.8828, 233.6969, 274.01, -109.5703, 233.722, 274.01, -125.9766, 233.7086, 274.01, -128.3203, 233.2849, 274.01, -130.6641, 233.702, 274.01, -135.3516, 233.727, 274.01, -149.4141, 233.7135, 274.01, -156.4453, 233.727, 274.01, -163.4766, 233.7119, 274.01, -168.1641, 233.7643, 274.01, -191.6016, 233.7809, 274.01, -210.3516, 233.7716, 274.01, -224.4141, 233.7998, 274.01, -233.7891, 233.7647, 274.01, -243.1641, 233.7785, 274.01, -247.8516, 233.7378, 274.01, -254.8828, 233.8578, 274.01, -257.2266, 235.2519, 274.01, -259.5703, 236.0817, 274.01, -260.7617, 237.0609, 274.01, -261.9141, 238.2262, 274.01, -262.8989, 239.4047, 274.01, -262.9743, 241.7484, 274.01, -262.5977, 244.0922, 274.01, -260.6675, 246.4359, 274.01, -259.6161, 248.7797, 274.01, -257.2266, 251.0035, 274.01, -255.0361, 253.4672, 274.01, -252.5391, 256.0995, 274.01, -251.2277, 258.1547, 274.01, -246.7444, 262.8422, 274.01, -243.1641, 266.3515, 274.01, -239.7411, 269.8734, 274.01, -238.4766, 270.9495, 274.01, -237.2269, 272.2172, 274.01, -236.1328, 273.5215, 274.01, -235.0934, 274.5609, 274.01, -233.7891, 275.6655, 274.01, -232.5319, 276.9047, 274.01, -231.4453, 278.1693, 274.01, -227.917, 281.5922, 274.01, -224.4141, 285.1802, 274.01, -222.0703, 287.4025, 274.01, -218.6841, 290.9672, 274.01, -217.3828, 291.9709, 274.01, -215.0391, 293.1788, 274.01, -212.6953, 294.5138, 274.01, -210.3516, 295.2614, 274.01, -209.8994, 295.6547, 274.01, -208.0078, 296.7083, 274.01, -203.3203, 296.9372, 274.01, -196.2891, 296.9317, 274.01, -193.9453, 296.8988, 274.01, -172.8516, 296.8482, 274.01, -161.1328, 296.843, 274.01, -137.6953, 296.8542, 274.01, -130.6641, 296.8378, 274.01, -107.2266, 296.8379, 274.01, -93.16406, 296.8208, 274.01, -74.41406, 296.838, 274.01, -65.03906, 296.8609, 274.01, -50.97656, 296.8556, 274.01, -46.28906, 296.9051, 274.01, -41.60156, 298.5331, 274.01, -39.25781, 298.5624, 274.01, -36.91406, 297.1455, 274.01, -34.57031, 297.0498, 274.01, -32.22656, 298.5589, 274.01, -25.19531, 298.5624, 274.01, -22.85156, 297.2842, 274.01, -20.50781, 298.5624, 274.01, -1.757813, 298.5624, 274.01, 0.5859375, 297.2104, 274.01, 2.929688, 298.5624, 274.01, 5.273438, 297.6946, 274.01, 7.617188, 298.5168, 274.01, 9.960938, 298.5512, 274.01, 12.30469, 296.937, 274.01, 14.64844, 298.5478, 274.01, 16.99219, 298.5478, 274.01, 19.33594, 297.0782, 274.01, 21.67969, 296.844, 274.01, 23.54531, 297.9984, 274.01, 24.02344, 298.4023, 274.01, 24.50156, 297.9984, 274.01, 26.36719, 296.844, 274.01, 38.08594, 296.8381, 274.01, 52.14844, 296.8033, 274.01, 59.17969, 296.8033, 274.01, 73.24219, 296.7682, 274.01, 99.02344, 296.7566, 274.01, 117.7734, 296.7216, 274.01, 129.4922, 296.7152, 274.01, 131.8359, 295.9083, 274.01, 134.1797, 295.5245, 274.01, 138.8672, 295.4763, 274.01, 155.2734, 295.4763, 274.01, 176.3672, 295.3861, 274.01, 190.4297, 295.3718, 274.01, 197.4609, 295.4763, 274.01, 202.1484, 295.4454, 274.01, 204.4922, 295.3438, 274.01, 206.8359, 295.0757, 274.01, 209.1797, 294.4124, 274.01, 211.5234, 292.3133, 274.01, 213.8672, 291.0809, 274.01, 216.2109, 289.6743, 274.01, 217.3081, 288.6234, 274.01, 219.3558, 286.2797, 274.01, 221.8608, 283.9359, 274.01, 225.5859, 280.045, 274.01, 227.9297, 277.8885, 274.01, 230.2734, 275.2857, 274.01, 232.6172, 273.2225, 274.01, 233.6225, 272.2172, 274.01, 234.9609, 270.5822, 274.01, 238.2254, 267.5297, 274.01, 240.3691, 265.1859, 274.01, 244.3359, 261.3326, 274.01, 246.6797, 258.8034, 274.01, 249.0234, 256.7196, 274.01, 251.3672, 254.0746, 274.01, 254.5313, 251.1234, 274.01, 255.333, 248.7797, 274.01, 257.3723, 246.4359, 274.01, 259.7201, 244.0922, 274.01, 260.106, 241.7484, 274.01, 261.6423, 239.4047, 274.01, 261.0804, 237.0609, 274.01, 259.3552, 234.7172, 274.01, 258.3984, 233.7696, 274.01, 256.0547, 232.8757, 274.01, 253.7109, 232.792, 274.01, 251.3672, 232.8161, 274.01, 246.6797, 232.6981, 274.01, 244.3359, 232.725, 274.01, 239.6484, 232.9137, 274.01, 234.9609, 232.8525, 274.01, 225.5859, 232.8757, 274.01, 209.1797, 232.8757, 274.01, 204.4922, 232.7537, 274.01, 195.1172, 232.7794, 274.01, 171.6797, 232.792, 274.01, 164.6484, 232.7666, 274.01, 152.9297, 232.7666, 274.01, 148.2422, 232.792, 274.01, 138.8672, 232.7406, 274.01 };
201 size_t n = sizeof(groundTruthRaw) / sizeof(groundTruthRaw[0]);
202 std::vector<double> groundTruth(groundTruthRaw, groundTruthRaw+n);
203 EXPECT_EQ(340u * 3, groundTruth.size());
204 CheckGroundTruth(structures, 5, 0, groundTruth);
205 }
206
207 {
208 double groundTruthRaw[] = { -18.16406, 233.0632, 298.01, -27.53906, 233.1042, 298.01, -29.88281, 233.0819, 298.01, -34.57031, 233.131, 298.01, -43.94531, 233.1221, 298.01, -50.97656, 233.1736, 298.01, -62.69531, 233.1397, 298.01, -65.03906, 232.8376, 298.01, -69.72656, 232.9839, 298.01, -79.10156, 233.0245, 298.01, -90.82031, 233.0382, 298.01, -93.16406, 233.0859, 298.01, -109.5703, 233.1132, 298.01, -111.9141, 233.1791, 298.01, -114.2578, 233.7139, 298.01, -118.9453, 233.9793, 298.01, -128.3203, 234.0284, 298.01, -130.6641, 233.9793, 298.01, -135.3516, 234.0591, 298.01, -137.6953, 234.0284, 298.01, -142.3828, 234.0855, 298.01, -144.7266, 234.0284, 298.01, -151.7578, 234.002, 298.01, -158.7891, 234.0263, 298.01, -163.4766, 233.9784, 298.01, -165.8203, 234.0072, 298.01, -168.1641, 234.1756, 298.01, -170.5078, 234.2214, 298.01, -179.8828, 234.1934, 298.01, -186.9141, 234.2721, 298.01, -189.2578, 234.2289, 298.01, -193.9453, 234.2431, 298.01, -198.6328, 234.1692, 298.01, -200.9766, 234.2326, 298.01, -205.6641, 234.1271, 298.01, -212.6953, 234.2224, 298.01, -215.0391, 234.1992, 298.01, -222.0703, 234.3115, 298.01, -224.4141, 234.2224, 298.01, -226.7578, 234.2502, 298.01, -233.7891, 234.0906, 298.01, -238.4766, 234.0329, 298.01, -243.1641, 234.0283, 298.01, -247.8516, 233.7949, 298.01, -250.1953, 233.8681, 298.01, -252.5391, 234.7626, 298.01, -254.3469, 237.0609, 298.01, -255.6034, 239.4047, 298.01, -254.5181, 241.7484, 298.01, -254.2274, 244.0922, 298.01, -254.181, 248.7797, 298.01, -253.9355, 251.1234, 298.01, -253.5926, 253.4672, 298.01, -252.7483, 255.8109, 298.01, -250.8092, 258.1547, 298.01, -248.713, 260.4984, 298.01, -246.263, 262.8422, 298.01, -244.1406, 265.1859, 298.01, -241.6671, 267.5297, 298.01, -239.4754, 269.8734, 298.01, -237.0156, 272.2172, 298.01, -233.7891, 275.382, 298.01, -231.4453, 277.8249, 298.01, -229.1016, 279.9981, 298.01, -226.7578, 282.5281, 298.01, -224.4141, 284.6784, 298.01, -222.0703, 287.2355, 298.01, -220.5414, 288.6234, 298.01, -218.2745, 290.9672, 298.01, -217.3828, 291.6508, 298.01, -212.6953, 294.5949, 298.01, -210.3516, 295.3142, 298.01, -208.0078, 296.4674, 298.01, -205.6641, 296.8852, 298.01, -203.3203, 297.1563, 298.01, -196.2891, 297.1488, 298.01, -193.9453, 297.0597, 298.01, -182.2266, 296.9529, 298.01, -168.1641, 296.8576, 298.01, -154.1016, 296.9249, 298.01, -149.4141, 296.8921, 298.01, -128.3203, 296.9228, 298.01, -121.2891, 296.8623, 298.01, -111.9141, 296.8549, 298.01, -107.2266, 296.8266, 298.01, -102.5391, 296.8731, 298.01, -95.50781, 296.8453, 298.01, -88.47656, 296.9218, 298.01, -83.78906, 296.9016, 298.01, -69.72656, 296.979, 298.01, -67.38281, 296.9514, 298.01, -65.03906, 297.2199, 298.01, -62.69531, 296.9622, 298.01, -55.66406, 296.9926, 298.01, -50.97656, 296.9467, 298.01, -48.63281, 297.3652, 298.01, -46.28906, 297.0439, 298.01, -43.94531, 297.2875, 298.01, -39.25781, 297.0121, 298.01, -34.57031, 297.1564, 298.01, -32.22656, 297.3612, 298.01, -29.88281, 297.4229, 298.01, -27.53906, 297.1687, 298.01, -25.19531, 297.4334, 298.01, -18.16406, 297.3612, 298.01, -15.82031, 297.4441, 298.01, -13.47656, 297.4125, 298.01, -11.13281, 297.2468, 298.01, -8.789063, 297.4125, 298.01, -6.445313, 297.373, 298.01, -4.101563, 297.4195, 298.01, -1.757813, 297.077, 298.01, 0.5859375, 297.4229, 298.01, 2.929688, 297.4125, 298.01, 5.273438, 296.9489, 298.01, 7.617188, 297.3168, 298.01, 9.960938, 296.9377, 298.01, 12.30469, 296.8998, 298.01, 14.64844, 297.1975, 298.01, 16.99219, 296.8579, 298.01, 28.71094, 296.878, 298.01, 40.42969, 296.8163, 298.01, 42.77344, 296.8369, 298.01, 49.80469, 296.734, 298.01, 59.17969, 296.6906, 298.01, 61.52344, 296.6365, 298.01, 68.55469, 296.6278, 298.01, 73.24219, 296.5777, 298.01, 75.58594, 296.6191, 298.01, 84.96094, 296.5284, 298.01, 96.67969, 296.5538, 298.01, 103.7109, 296.479, 298.01, 115.4297, 296.4259, 298.01, 122.4609, 296.3434, 298.01, 129.4922, 296.3495, 298.01, 131.8359, 295.9141, 298.01, 136.5234, 296.2256, 298.01, 138.8672, 295.833, 298.01, 143.5547, 295.9857, 298.01, 145.8984, 295.8791, 298.01, 152.9297, 295.833, 298.01, 164.6484, 295.6819, 298.01, 171.6797, 295.6819, 298.01, 181.0547, 295.5401, 298.01, 185.7422, 295.5742, 298.01, 192.7734, 295.557, 298.01, 197.4609, 295.8012, 298.01, 202.1484, 295.6819, 298.01, 204.4922, 295.3698, 298.01, 206.8359, 294.803, 298.01, 209.1797, 294.3656, 298.01, 211.5234, 292.4764, 298.01, 213.8672, 291.1765, 298.01, 216.2109, 289.5873, 298.01, 217.229, 288.6234, 298.01, 218.5547, 287.0752, 298.01, 221.7097, 283.9359, 298.01, 225.5859, 279.8775, 298.01, 227.9297, 277.5633, 298.01, 230.2734, 275.0808, 298.01, 233.1989, 272.2172, 298.01, 234.9609, 270.2887, 298.01, 237.7384, 267.5297, 298.01, 241.9922, 263.0843, 298.01, 244.3359, 260.7643, 298.01, 246.788, 258.1547, 298.01, 249.0234, 255.451, 298.01, 250.3651, 253.4672, 298.01, 251.5297, 251.1234, 298.01, 252.1947, 248.7797, 298.01, 252.4915, 246.4359, 298.01, 252.5755, 241.7484, 298.01, 252.8592, 239.4047, 298.01, 252.9236, 237.0609, 298.01, 252.2924, 234.7172, 298.01, 251.3672, 233.4697, 298.01, 249.0234, 232.882, 298.01, 244.3359, 232.9048, 298.01, 241.9922, 233.0145, 298.01, 232.6172, 232.9048, 298.01, 227.9297, 233.0007, 298.01, 216.2109, 233.0632, 298.01, 211.5234, 233.0537, 298.01, 206.8359, 232.9699, 298.01, 204.4922, 232.7322, 298.01, 199.8047, 232.7186, 298.01, 190.4297, 232.7719, 298.01, 183.3984, 232.7719, 298.01, 181.0547, 232.7322, 298.01, 174.0234, 232.7048, 298.01, 171.6797, 232.7322, 298.01, 166.9922, 232.6908, 298.01, 157.6172, 232.7975, 298.01, 155.2734, 232.7588, 298.01, 148.2422, 232.7875, 298.01, 143.5547, 232.7614, 298.01, 138.8672, 232.6477, 298.01, 124.8047, 232.6179, 298.01, 122.4609, 232.6477, 298.01, 113.0859, 232.6027, 298.01, 110.7422, 232.4552, 298.01, 108.3984, 232.2192, 298.01, 106.0547, 231.6764, 298.01, 103.7109, 231.8559, 298.01, 102.8237, 232.3734, 298.01, 101.3672, 232.9839, 298.01, 99.02344, 233.0951, 298.01, 87.30469, 233.0819, 298.01, 84.96094, 233.1091, 298.01, 80.27344, 233.0726, 298.01, 77.92969, 233.1132, 298.01, 70.89844, 233.1397, 298.01, 68.55469, 233.1132, 298.01, 52.14844, 233.131, 298.01, 45.11719, 233.0859, 298.01, 44.16726, 232.3734, 298.01, 42.77344, 231.0206, 298.01, 42.04498, 230.0297, 298.01, 42.77344, 229.2462, 298.01, 45.11719, 228.5664, 298.01, 47.46094, 227.0695, 298.01, 49.80469, 226.0552, 298.01, 52.14844, 224.5723, 298.01, 54.49219, 223.6857, 298.01, 56.83594, 221.8519, 298.01, 59.17969, 220.2086, 298.01, 61.52344, 218.8854, 298.01, 64.94469, 215.9672, 298.01, 66.21094, 215.0191, 298.01, 67.72036, 213.6234, 298.01, 68.55469, 212.6986, 298.01, 70.89844, 210.5055, 298.01, 74.53191, 206.5922, 298.01, 76.54903, 204.2484, 298.01, 78.26105, 201.9047, 298.01, 80.27344, 198.9262, 298.01, 82.61719, 195.2822, 298.01, 82.98087, 194.8734, 298.01, 84.96094, 190.9255, 298.01, 85.43701, 190.1859, 298.01, 86.33423, 187.8422, 298.01, 87.78722, 185.4984, 298.01, 88.60233, 183.1547, 298.01, 89.10253, 180.8109, 298.01, 90.17504, 178.4672, 298.01, 90.88959, 176.1234, 298.01, 91.43783, 173.7797, 298.01, 92.39601, 171.4359, 298.01, 92.95762, 169.0922, 298.01, 93.55695, 159.7172, 298.01, 93.65527, 157.3734, 298.01, 93.67542, 152.6859, 298.01, 93.61213, 150.3422, 298.01, 93.22542, 143.3109, 298.01, 93.06345, 140.9672, 298.01, 92.77563, 138.6234, 298.01, 91.21714, 133.9359, 298.01, 90.67235, 131.5922, 298.01, 89.88776, 129.2484, 298.01, 88.8737, 126.9047, 298.01, 88.44087, 124.5609, 298.01, 86.09712, 119.8734, 298.01, 85.05786, 117.5297, 298.01, 83.87151, 115.1859, 298.01, 82.22388, 112.8422, 298.01, 81.09117, 110.4984, 298.01, 77.92969, 106.4052, 298.01, 77.3894, 105.8109, 298.01, 75.94332, 103.4672, 298.01, 71.71799, 98.77969, 298.01, 68.55469, 95.65721, 298.01, 63.86719, 91.54878, 298.01, 61.52344, 90.1121, 298.01, 59.17969, 88.15762, 298.01, 56.83594, 86.51503, 298.01, 54.49219, 85.42721, 298.01, 52.14844, 83.64907, 298.01, 49.80469, 82.89023, 298.01, 47.46094, 81.50237, 298.01, 45.11719, 80.62591, 298.01, 42.77344, 79.18153, 298.01, 40.42969, 78.7203, 298.01, 38.08594, 78.1349, 298.01, 35.74219, 77.11755, 298.01, 33.39844, 76.51949, 298.01, 31.05469, 76.07934, 298.01, 26.36719, 74.67744, 298.01, 24.02344, 74.42056, 298.01, 14.64844, 74.07317, 298.01, 9.960938, 74.11538, 298.01, 2.929688, 74.40105, 298.01, 0.5859375, 74.67952, 298.01, -1.757813, 75.31406, 298.01, -4.101563, 76.07065, 298.01, -6.445313, 76.49051, 298.01, -8.789063, 77.17276, 298.01, -11.13281, 78.20097, 298.01, -15.82031, 79.31967, 298.01, -18.16406, 80.76948, 298.01, -20.50781, 81.64266, 298.01, -22.85156, 83.0305, 298.01, -25.19531, 83.7937, 298.01, -27.53906, 85.63515, 298.01, -29.88281, 86.7363, 298.01, -32.22656, 88.36089, 298.01, -34.57031, 90.3302, 298.01, -36.56719, 91.74844, 298.01, -41.60156, 95.93605, 298.01, -46.58845, 101.1234, 298.01, -50.17995, 105.8109, 298.01, -52.10386, 108.1547, 298.01, -53.63992, 110.4984, 298.01, -54.95532, 112.8422, 298.01, -56.64794, 115.1859, 298.01, -57.4403, 117.5297, 298.01, -58.91927, 119.8734, 298.01, -59.78655, 122.2172, 298.01, -61.11754, 124.5609, 298.01, -61.58921, 126.9047, 298.01, -62.38012, 129.2484, 298.01, -63.49118, 131.5922, 298.01, -64.02599, 133.9359, 298.01, -64.3932, 136.2797, 298.01, -65.11897, 138.6234, 298.01, -65.64544, 140.9672, 298.01, -66.23938, 147.9984, 298.01, -66.46289, 152.6859, 298.01, -66.48911, 155.0297, 298.01, -66.34437, 159.7172, 298.01, -65.99894, 164.4047, 298.01, -65.49149, 169.0922, 298.01, -64.6875, 171.4359, 298.01, -63.7739, 176.1234, 298.01, -62.9398, 178.4672, 298.01, -61.86011, 180.8109, 298.01, -61.33423, 183.1547, 298.01, -60.43332, 185.4984, 298.01, -58.00781, 190.0632, 298.01, -56.85406, 192.5297, 298.01, -55.66406, 194.7283, 298.01, -54.11692, 197.2172, 298.01, -50.97656, 201.8369, 298.01, -47.36435, 206.5922, 298.01, -45.04395, 208.9359, 298.01, -42.83026, 211.2797, 298.01, -39.25781, 214.7435, 298.01, -34.57031, 218.4974, 298.01, -32.22656, 219.9595, 298.01, -28.02053, 222.9984, 298.01, -27.53906, 223.4238, 298.01, -25.19531, 224.4187, 298.01, -22.85156, 225.8252, 298.01, -20.50781, 226.9067, 298.01, -18.16406, 228.4286, 298.01, -15.82031, 229.1235, 298.01, -14.9447, 230.0297, 298.01, -15.82031, 231.3969, 298.01, -16.94484, 232.3734, 298.01 };
209 size_t n = sizeof(groundTruthRaw) / sizeof(groundTruthRaw[0]);
210 std::vector<double> groundTruth(groundTruthRaw, groundTruthRaw+n);
211 EXPECT_EQ(358u * 3, groundTruth.size());
212 CheckGroundTruth(structures, 5, 8, groundTruth);
213 }
214 }
215
216 #endif
217 // BGO_ENABLE_DICOMSTRUCTURESETLOADER2
218
219 #if 0
220
221 TEST(StructureSet2, ReadFromJsonAndCompute1)
222 {
223 DicomStructureSet2 structureSet;
224
225 OrthancPlugins::FullOrthancDataset dicom(GetTestJson());
226 structureSet.Clear();
227
228 structureSet.FillStructuresFromDataset(dicom);
229
230 structureSet.ComputeDependentProperties();
231 }
232
233 TEST(StructureSet2, ReadFromJsonAndCompute2)
234 {
235 DicomStructureSet2 structureSet;
236
237 OrthancPlugins::FullOrthancDataset dicom(GetTestJson());
238 structureSet.Clear();
239
240 structureSet.SetContents(dicom);
241 }
242 #endif
243
244 #ifdef BGO_ENABLE_DICOMSTRUCTURESETLOADER2
245
246 static bool CutStructureWithPlane(
247 std::vector< std::pair<ScenePoint2D, ScenePoint2D> >& segments,
248 const DicomStructure2& structure,
249 const double originX, const double originY, const double originZ,
250 const double axisX_X, const double axisX_Y, const double axisX_Z,
251 const double axisY_X, const double axisY_Y, const double axisY_Z
252 )
253 {
254 // create an AXIAL cutting plane, too far away from the volume
255 // (> sliceThickness/2)
256 Vector origin, axisX, axisY;
257 LinearAlgebra::AssignVector(origin, originX, originY, originZ);
258 LinearAlgebra::AssignVector(axisX, axisX_X, axisX_Y, axisX_Z);
259 LinearAlgebra::AssignVector(axisY, axisY_X, axisY_Y, axisY_Z);
260 CoordinateSystem3D cuttingPlane(origin, axisX, axisY);
261
262 // compute intersection
263 bool ok = structure.Project(segments, cuttingPlane);
264 return ok;
265 }
266
267 #endif
268 // BGO_ENABLE_DICOMSTRUCTURESETLOADER2
269
270 static double pointsCoord1[] = { 2, 2, 3, 3, 6, 8, 8, 7, 8, 8, 6 };
271 static double pointsCoord2[] = { 2, 6, 8, 10, 12, 10, 8, 6, 4, 2, 4 };
272 static const size_t pointsCoord1Count = STONE_ARRAY_SIZE(pointsCoord1);
273 static const size_t pointsCoord2Count = STONE_ARRAY_SIZE(pointsCoord2);
274 const size_t POLYGON_POINT_COUNT = pointsCoord1Count;
275
276 #ifdef BGO_ENABLE_DICOMSTRUCTURESETLOADER2
277
278 static void CreateBasicStructure(DicomStructure2& structure)
279 {
280 // see https://www.dropbox.com/s/1o1vg53hsbvx4cc/test-rtstruct-polygons.jpg?dl=0
281 EXPECT_EQ(pointsCoord1Count, pointsCoord2Count);
282 EXPECT_EQ(11u, pointsCoord2Count);
283
284 for (size_t slice = 0; slice < 3; ++slice)
285 {
286 DicomStructurePolygon2 polygon("Oblomptu", "CLOSED_PLANAR");
287 for (size_t ip = 0; ip < pointsCoord1Count; ++ip)
288 {
289 Vector pt;
290 double pt0 = pointsCoord1[ip];
291 double pt1 = pointsCoord2[ip];
292 double pt2 = 4 * (static_cast<double>(slice) - 1); // -4, 0, 4
293 LinearAlgebra::AssignVector(pt, pt0, pt1, pt2);
294 polygon.AddPoint(pt);
295 }
296 structure.AddPolygon(polygon);
297 }
298 structure.ComputeDependentProperties();
299 }
300
301
302 TEST(StructureSet2, CutAxialOutsideTop)
303 {
304 DicomStructure2 structure;
305 CreateBasicStructure(structure);
306 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
307
308 // create an AXIAL cutting plane, too far away from the volume
309 // (> sliceThickness/2)
310 bool ok = CutStructureWithPlane(segments, structure,
311 0, 0, 7,
312 1, 0, 0,
313 0, 1, 0);
314 EXPECT_FALSE(ok);
315 }
316
317
318 TEST(StructureSet2, CutAxialOutsideBottom)
319 {
320 DicomStructure2 structure;
321 CreateBasicStructure(structure);
322 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
323
324 // create an AXIAL cutting plane, too far away from the volume
325 // (> sliceThickness/2)
326 bool ok = CutStructureWithPlane(segments, structure,
327 0, 0, -6.66,
328 1, 0, 0,
329 0, 1, 0);
330 EXPECT_FALSE(ok);
331 }
332
333 TEST(StructureSet2, CutAxialInsideClose)
334 {
335 DicomStructure2 structure;
336 CreateBasicStructure(structure);
337 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
338
339 // create an AXIAL cutting plane in the volume
340 bool ok = CutStructureWithPlane(segments, structure,
341 0, 0, 1.1,
342 1, 0, 0,
343 0, 1, 0);
344 EXPECT_TRUE(ok);
345 EXPECT_EQ(POLYGON_POINT_COUNT, segments.size());
346
347 for (size_t i = 0; i < segments.size(); ++i)
348 {
349 EXPECT_LT(i, POLYGON_POINT_COUNT);
350 EXPECT_LT(i, POLYGON_POINT_COUNT);
351
352 const ScenePoint2D& pt = segments[i].first;
353
354 // ...should be at the same location as the 3D coords since the plane
355 // is rooted at 0,0,0 with normal 0,0,1
356 EXPECT_NEAR(pt.GetX(), pointsCoord1[i], DELTA_MAX);
357 EXPECT_NEAR(pt.GetY(), pointsCoord2[i], DELTA_MAX);
358 }
359 }
360
361
362 TEST(StructureSet2, CutAxialInsideFar)
363 {
364 DicomStructure2 structure;
365 CreateBasicStructure(structure);
366 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
367
368 // create an AXIAL cutting plane
369 Vector origin, axisX, axisY;
370 LinearAlgebra::AssignVector(origin, 0, 0, 0);
371 LinearAlgebra::AssignVector(axisX, 1, 0, 0);
372 LinearAlgebra::AssignVector(axisY, 0, 1, 0);
373 CoordinateSystem3D cuttingPlane(origin, axisX, axisY);
374
375 // compute intersection
376 bool ok = structure.Project(segments, cuttingPlane);
377 EXPECT_TRUE(ok);
378
379 EXPECT_EQ(11u, segments.size());
380 for (size_t i = 0; i < segments.size(); ++i)
381 {
382 EXPECT_LT(i, pointsCoord1Count);
383 EXPECT_LT(i, pointsCoord2Count);
384
385 // the 2D points of the projected polygon
386 const ScenePoint2D& pt = segments[i].first;
387
388 // ...should be at the same location as the 3D coords since the plane
389 // is rooted at 0,0,0 with normal 0,0,1
390 EXPECT_NEAR(pt.GetX(), pointsCoord1[i], DELTA_MAX);
391 EXPECT_NEAR(pt.GetY(), pointsCoord2[i], DELTA_MAX);
392 }
393 }
394
395 TEST(StructureSet2, CutCoronalOutsideClose)
396 {
397 DicomStructure2 structure;
398 CreateBasicStructure(structure);
399 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
400
401 // create an X,Z cutting plane, outside of the volume
402 Vector origin, axisX, axisY;
403 LinearAlgebra::AssignVector(origin, 0, 0, 0);
404 LinearAlgebra::AssignVector(axisX, 1, 0, 0);
405 LinearAlgebra::AssignVector(axisY, 0, 0, 1);
406 CoordinateSystem3D cuttingPlane(origin, axisX, axisY);
407
408 // compute intersection
409 bool ok = structure.Project(segments, cuttingPlane);
410 EXPECT_FALSE(ok);
411 }
412
413 TEST(StructureSet2, CutCoronalInsideClose1DTest)
414 {
415 DicomStructure2 structure;
416 CreateBasicStructure(structure);
417
418 // create an X,Z cutting plane, outside of the volume
419 Vector origin, axisX, axisY;
420 LinearAlgebra::AssignVector(origin, 0, 3, 0);
421 LinearAlgebra::AssignVector(axisX, 1, 0, 0);
422 LinearAlgebra::AssignVector(axisY, 0, 0, 1);
423 CoordinateSystem3D cuttingPlane(origin, axisX, axisY);
424
425 ASSERT_EQ(3u, structure.GetPolygons().size());
426
427 for (int i = 0; i < 3; ++i)
428 {
429 double polygonZ = static_cast<double>(i - 1) * 4.0;
430
431 const DicomStructurePolygon2& topSlab = structure.GetPolygons()[i];
432
433 // let's compute the intersection between the polygon and the plane
434 // intersections are in plane coords
435 std::vector<ScenePoint2D> intersects;
436 topSlab.ProjectOnConstantPlane(intersects, cuttingPlane);
437
438 ASSERT_EQ(4u, intersects.size());
439
440 EXPECT_NEAR(2, intersects[0].GetX(), DELTA_MAX);
441 EXPECT_NEAR(4, intersects[1].GetX(), DELTA_MAX);
442 EXPECT_NEAR(7, intersects[2].GetX(), DELTA_MAX);
443 EXPECT_NEAR(8, intersects[3].GetX(), DELTA_MAX);
444
445 for (size_t i = 0; i < 4u; ++i)
446 {
447 EXPECT_NEAR(polygonZ, intersects[i].GetY(), DELTA_MAX);
448 }
449 }
450 }
451
452 TEST(StructureSet2, CutCoronalInsideClose)
453 {
454 DicomStructure2 structure;
455 CreateBasicStructure(structure);
456 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
457
458 // create an X,Z cutting plane, outside of the volume
459 Vector origin, axisX, axisY;
460 LinearAlgebra::AssignVector(origin, 0, 3, 0);
461 LinearAlgebra::AssignVector(axisX, 1, 0, 0);
462 LinearAlgebra::AssignVector(axisY, 0, 0, 1);
463 CoordinateSystem3D cuttingPlane(origin, axisX, axisY);
464
465 // compute intersection
466 ASSERT_TRUE(structure.Project(segments, cuttingPlane));
467 EXPECT_EQ(24u, segments.size());
468
469 size_t numberOfVeryShortSegments = 0;
470 for (size_t iSegment = 0; iSegment < segments.size(); ++iSegment)
471 {
472 // count the NON vertical very short segments
473 if (LinearAlgebra::IsNear(segments[iSegment].first.GetX(), segments[iSegment].second.GetX()))
474 {
475 if (LinearAlgebra::IsNear(segments[iSegment].first.GetY(), segments[iSegment].second.GetY()))
476 {
477 numberOfVeryShortSegments++;
478 }
479 }
480 }
481 EXPECT_EQ(8u, numberOfVeryShortSegments);
482 }
483
484 #endif
485 // BGO_ENABLE_DICOMSTRUCTURESETLOADER2
486
487
488 TEST(DisjointDataSet, BasicTest)
489 {
490 const size_t ITEM_COUNT = 10;
491 DisjointDataSet ds(ITEM_COUNT);
492
493 for (size_t i = 0; i < ITEM_COUNT; ++i)
494 {
495 EXPECT_EQ(i, ds.Find(i));
496 }
497
498 ds.Union(0, 4);
499 EXPECT_EQ(0u, ds.Find(0));
500 EXPECT_EQ(0u, ds.Find(4));
501
502 ds.Union(4, 6);
503 ds.Union(8, 9);
504 ds.Union(0, 8);
505
506 for (size_t i = 0; i < ITEM_COUNT; ++i)
507 {
508 size_t parent = ds.Find(i);
509 EXPECT_TRUE(0 == parent || 1 == parent || 2 == parent || 3 == parent || 5 == parent || 7 == parent);
510 }
511
512 ds.Union(1, 2);
513 ds.Union(1, 7);
514 for (size_t i = 0; i < ITEM_COUNT; ++i)
515 {
516 size_t parent = ds.Find(i);
517 EXPECT_TRUE(0 == parent || 1 == parent || 3 == parent || 5 == parent);
518 }
519
520 ds.Union(3, 5);
521 for (size_t i = 0; i < ITEM_COUNT; ++i)
522 {
523 size_t parent = ds.Find(i);
524 EXPECT_TRUE(0 == parent || 1 == parent || 3 == parent);
525 }
526
527 EXPECT_EQ(ds.Find(0), ds.Find(0));
528 EXPECT_EQ(ds.Find(0), ds.Find(4));
529 EXPECT_EQ(ds.Find(0), ds.Find(6));
530 EXPECT_EQ(ds.Find(0), ds.Find(8));
531 EXPECT_EQ(ds.Find(0), ds.Find(8));
532
533 EXPECT_EQ(ds.Find(1), ds.Find(7));
534 EXPECT_EQ(ds.Find(2), ds.Find(1));
535 EXPECT_EQ(ds.Find(7), ds.Find(2));
536
537 EXPECT_EQ(ds.Find(3), ds.Find(5));
538 EXPECT_EQ(ds.Find(5), ds.Find(3));
539
540 ds.Union(0, 1);
541 ds.Union(3, 1);
542 for (size_t i = 0; i < ITEM_COUNT; ++i)
543 {
544 EXPECT_EQ(ds.Find(0), ds.Find(i));
545 }
546 }
547
548 #ifdef BGO_ENABLE_DICOMSTRUCTURESETLOADER2
549
550 TEST(StructureSet2, CutSagittalInsideClose)
551 {
552 DicomStructure2 structure;
553 CreateBasicStructure(structure);
554 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
555
556 // create an X,Z cutting plane, inside of the volume
557 Vector origin, axisX, axisY;
558 LinearAlgebra::AssignVector(origin, 0, 3, 0);
559 LinearAlgebra::AssignVector(axisX, 1, 0, 0);
560 LinearAlgebra::AssignVector(axisY, 0, 0, 1);
561 CoordinateSystem3D cuttingPlane(origin, axisX, axisY);
562
563 // compute intersection
564 bool ok = structure.Project(segments, cuttingPlane);
565 EXPECT_TRUE(ok);
566 }
567
568 #endif
569 // BGO_ENABLE_DICOMSTRUCTURESETLOADER2
570
571
572 static size_t ConvertListOfSlabsToSegments_Add(RtStructRectanglesInSlab& rectangles, int row, double xmin, double xmax)
573 {
574 double ymin = static_cast<double>(row) * 5.0;
575 double ymax = static_cast<double>(row + 1) * 5.0;
576
577 RtStructRectangleInSlab rectangle;
578 rectangle.xmin = xmin;
579 rectangle.xmax = xmax;
580 rectangle.ymin = ymin;
581 rectangle.ymax = ymax;
582
583 rectangles.push_back(rectangle);
584
585 return 1u;
586 }
587
588 static size_t FillTestRectangleList(std::vector< RtStructRectanglesInSlab >& rectanglesForEachSlab)
589 {
590 // ConvertListOfSlabsToSegments
591 size_t rectCount = 0;
592
593 rectanglesForEachSlab.push_back(RtStructRectanglesInSlab());
594 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 0, 5, 31);
595 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 0, 36, 50);
596
597 rectanglesForEachSlab.push_back(RtStructRectanglesInSlab());
598 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 1, 20, 45);
599 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 1, 52, 70);
600
601 rectanglesForEachSlab.push_back(RtStructRectanglesInSlab());
602 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 2, 0, 32);
603 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 2, 35, 44);
604 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 2, 60, 75);
605
606 rectanglesForEachSlab.push_back(RtStructRectanglesInSlab());
607 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 3, 10, 41);
608 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 3, 46, 80);
609
610 rectanglesForEachSlab.push_back(RtStructRectanglesInSlab());
611 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 4, 34, 42);
612 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 4, 90, 96);
613
614 rectanglesForEachSlab.push_back(RtStructRectanglesInSlab());
615 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 5, 1, 33);
616 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 5, 40, 43);
617 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 5, 51, 61);
618 rectCount += ConvertListOfSlabsToSegments_Add(rectanglesForEachSlab.back(), 5, 76, 95);
619
620 return rectCount;
621 }
622
623 /*
624 void AddSlabBoundaries(
625 std::vector<std::pair<double, RectangleBoundaryKind> >& boundaries,
626 const std::vector<RtStructRectanglesInSlab>& slabCuts, size_t iSlab)
627 */
628
629
630 /*
631 void ProcessBoundaryList(
632 std::vector< std::pair<ScenePoint2D, ScenePoint2D> >& segments,
633 const std::vector<std::pair<double, RectangleBoundaryKind> >& boundaries,
634 double y)
635 */
636
637
638 TEST(StructureSet2, ProcessBoundaryList_Empty)
639 {
640 std::vector< RtStructRectanglesInSlab > slabCuts;
641 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
642
643 boundaries.clear();
644 EXPECT_NO_THROW(AddSlabBoundaries(boundaries, slabCuts, 0));
645 ASSERT_EQ(0u, boundaries.size());
646 }
647
648 TEST(StructureSet2, ProcessBoundaryListTopRow)
649 {
650 std::vector< RtStructRectanglesInSlab > slabCuts;
651 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
652 FillTestRectangleList(slabCuts);
653
654 boundaries.clear();
655 AddSlabBoundaries(boundaries, slabCuts, 0);
656
657 {
658 size_t i = 0;
659 ASSERT_EQ(4u, boundaries.size());
660
661 ASSERT_EQ(RectangleBoundaryKind_Start, boundaries[i].second);
662 ASSERT_NEAR(5, boundaries[i].first, DELTA_MAX);
663 i++;
664
665 ASSERT_EQ(RectangleBoundaryKind_End, boundaries[i].second);
666 ASSERT_NEAR(31, boundaries[i].first, DELTA_MAX);
667 i++;
668
669 ASSERT_EQ(RectangleBoundaryKind_Start, boundaries[i].second);
670 ASSERT_NEAR(36, boundaries[i].first, DELTA_MAX);
671 i++;
672
673 ASSERT_EQ(RectangleBoundaryKind_End, boundaries[i].second);
674 ASSERT_NEAR(50, boundaries[i].first, DELTA_MAX);
675 i++;
676 }
677 }
678
679 TEST(StructureSet2, ProcessBoundaryListRows_0_and_1)
680 {
681 std::vector< RtStructRectanglesInSlab > slabCuts;
682 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
683 FillTestRectangleList(slabCuts);
684
685 boundaries.clear();
686 AddSlabBoundaries(boundaries, slabCuts, 0);
687 AddSlabBoundaries(boundaries, slabCuts, 1);
688
689 ASSERT_EQ(8u, boundaries.size());
690
691 {
692 size_t i = 0;
693
694 ASSERT_EQ(RectangleBoundaryKind_Start, boundaries[i].second);
695 ASSERT_NEAR(5, boundaries[i].first, DELTA_MAX);
696 i++;
697
698 ASSERT_EQ(RectangleBoundaryKind_Start, boundaries[i].second);
699 ASSERT_NEAR(20, boundaries[i].first, DELTA_MAX);
700 i++;
701
702 ASSERT_EQ(RectangleBoundaryKind_End, boundaries[i].second);
703 ASSERT_NEAR(31, boundaries[i].first, DELTA_MAX);
704 i++;
705
706 ASSERT_EQ(RectangleBoundaryKind_Start, boundaries[i].second);
707 ASSERT_NEAR(36, boundaries[i].first, DELTA_MAX);
708 i++;
709
710 ASSERT_EQ(RectangleBoundaryKind_End, boundaries[i].second);
711 ASSERT_NEAR(45, boundaries[i].first, DELTA_MAX);
712 i++;
713
714 ASSERT_EQ(RectangleBoundaryKind_End, boundaries[i].second);
715 ASSERT_NEAR(50, boundaries[i].first, DELTA_MAX);
716 i++;
717
718 ASSERT_EQ(RectangleBoundaryKind_Start, boundaries[i].second);
719 ASSERT_NEAR(52, boundaries[i].first, DELTA_MAX);
720 i++;
721
722 ASSERT_EQ(RectangleBoundaryKind_End, boundaries[i].second);
723 ASSERT_NEAR(70, boundaries[i].first, DELTA_MAX);
724 i++;
725 }
726 }
727
728 TEST(StructureSet2, ConvertListOfSlabsToSegments_EmptyBoundaries)
729 {
730 std::vector< RtStructRectanglesInSlab > slabCuts;
731 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
732 FillTestRectangleList(slabCuts);
733 boundaries.clear();
734 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
735 ASSERT_NO_THROW(ProcessBoundaryList(segments, boundaries, 42.0));
736 ASSERT_EQ(0u, segments.size());
737 }
738
739 TEST(StructureSet2, ConvertListOfSlabsToSegments_TopRow_Horizontal)
740 {
741 std::vector< RtStructRectanglesInSlab > slabCuts;
742 FillTestRectangleList(slabCuts);
743
744 // top row
745 {
746 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
747 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
748 AddSlabBoundaries(boundaries, slabCuts, 0);
749 ProcessBoundaryList(segments, boundaries, slabCuts[0][0].ymin);
750
751 ASSERT_EQ(2u, segments.size());
752
753 ASSERT_NEAR( 5.0, segments[0].first.GetX(), DELTA_MAX);
754 ASSERT_NEAR(31.0, segments[0].second.GetX(), DELTA_MAX);
755 ASSERT_NEAR( 0.0, segments[0].first.GetY(), DELTA_MAX);
756 ASSERT_NEAR( 0.0, segments[0].second.GetY(), DELTA_MAX);
757
758 ASSERT_NEAR(36.0, segments[1].first.GetX(), DELTA_MAX);
759 ASSERT_NEAR(50.0, segments[1].second.GetX(), DELTA_MAX);
760 ASSERT_NEAR( 0.0, segments[1].first.GetY(), DELTA_MAX);
761 ASSERT_NEAR( 0.0, segments[1].second.GetY(), DELTA_MAX);
762 }
763 }
764
765 TEST(StructureSet2, ConvertListOfSlabsToSegments_All_Horizontal)
766 {
767 std::vector< RtStructRectanglesInSlab > slabCuts;
768 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
769 FillTestRectangleList(slabCuts);
770
771 // top row
772 {
773 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
774 AddSlabBoundaries(boundaries, slabCuts, 0);
775 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
776 ProcessBoundaryList(segments, boundaries, slabCuts[0][0].ymin);
777 }
778
779 // mids
780 {
781 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
782 AddSlabBoundaries(boundaries, slabCuts, 0);
783 AddSlabBoundaries(boundaries, slabCuts, 1);
784 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
785 ProcessBoundaryList(segments, boundaries, slabCuts[0][0].ymax);
786
787 ASSERT_EQ(4u, segments.size());
788
789 ASSERT_NEAR(05.0, segments[0].first.GetX(), DELTA_MAX);
790 ASSERT_NEAR(20.0, segments[0].second.GetX(), DELTA_MAX);
791 ASSERT_NEAR(05.0, segments[0].first.GetY(), DELTA_MAX);
792 ASSERT_NEAR(05.0, segments[0].second.GetY(), DELTA_MAX);
793
794 ASSERT_NEAR(31.0, segments[1].first.GetX(), DELTA_MAX);
795 ASSERT_NEAR(36.0, segments[1].second.GetX(), DELTA_MAX);
796 ASSERT_NEAR(05.0, segments[1].first.GetY(), DELTA_MAX);
797 ASSERT_NEAR(05.0, segments[1].second.GetY(), DELTA_MAX);
798
799 ASSERT_NEAR(45.0, segments[2].first.GetX(), DELTA_MAX);
800 ASSERT_NEAR(50.0, segments[2].second.GetX(), DELTA_MAX);
801 ASSERT_NEAR(05.0, segments[2].first.GetY(), DELTA_MAX);
802 ASSERT_NEAR(05.0, segments[2].second.GetY(), DELTA_MAX);
803
804 ASSERT_NEAR(52.0, segments[3].first.GetX(), DELTA_MAX);
805 ASSERT_NEAR(70.0, segments[3].second.GetX(), DELTA_MAX);
806 ASSERT_NEAR(05.0, segments[3].first.GetY(), DELTA_MAX);
807 ASSERT_NEAR(05.0, segments[3].second.GetY(), DELTA_MAX);
808 }
809
810 // bottom row
811 {
812 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
813 AddSlabBoundaries(boundaries, slabCuts, 1);
814 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
815 ProcessBoundaryList(segments, boundaries, slabCuts[1][0].ymax);
816
817 ASSERT_EQ(2u, segments.size());
818
819 ASSERT_NEAR(20.0, segments[0].first.GetX(), DELTA_MAX);
820 ASSERT_NEAR(45.0, segments[0].second.GetX(), DELTA_MAX);
821 ASSERT_NEAR(10.0, segments[0].first.GetY(), DELTA_MAX);
822 ASSERT_NEAR(10.0, segments[0].second.GetY(), DELTA_MAX);
823
824 ASSERT_NEAR(52.0, segments[1].first.GetX(), DELTA_MAX);
825 ASSERT_NEAR(70.0, segments[1].second.GetX(), DELTA_MAX);
826 ASSERT_NEAR(10.0, segments[1].first.GetY(), DELTA_MAX);
827 ASSERT_NEAR(10.0, segments[1].second.GetY(), DELTA_MAX);
828 }
829
830 }
831
832 TEST(StructureSet2, ConvertListOfSlabsToSegments_Complete_Empty)
833 {
834 std::vector< RtStructRectanglesInSlab > slabCuts;
835 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
836
837 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
838
839 ASSERT_NO_THROW(ConvertListOfSlabsToSegments(segments, slabCuts, 0));
840 ASSERT_EQ(0u, segments.size());
841 }
842
843 TEST(StructureSet2, ConvertListOfSlabsToSegments_Complete_Regular)
844 {
845 std::vector< RtStructRectanglesInSlab > slabCuts;
846 std::vector<std::pair<double, RectangleBoundaryKind> > boundaries;
847 size_t totalRectCount = FillTestRectangleList(slabCuts);
848
849 std::vector< std::pair<ScenePoint2D, ScenePoint2D> > segments;
850
851 ASSERT_NO_THROW(ConvertListOfSlabsToSegments(segments, slabCuts, totalRectCount));
852 ASSERT_EQ(60u, segments.size());
853
854 size_t i = 0;
855
856 ASSERT_NEAR(segments[i].first.GetX(), 5.0000000000000000, DELTA_MAX);
857 ASSERT_NEAR(segments[i].first.GetY(), 0.0000000000000000, DELTA_MAX);
858 ASSERT_NEAR(segments[i].second.GetX(), 5.0000000000000000, DELTA_MAX);
859 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
860 i++;
861 ASSERT_NEAR(segments[i].first.GetX(), 31.000000000000000, DELTA_MAX);
862 ASSERT_NEAR(segments[i].first.GetY(), 0.0000000000000000, DELTA_MAX);
863 ASSERT_NEAR(segments[i].second.GetX(), 31.000000000000000, DELTA_MAX);
864 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
865 i++;
866 ASSERT_NEAR(segments[i].first.GetX(), 36.000000000000000, DELTA_MAX);
867 ASSERT_NEAR(segments[i].first.GetY(), 0.0000000000000000, DELTA_MAX);
868 ASSERT_NEAR(segments[i].second.GetX(), 36.000000000000000, DELTA_MAX);
869 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
870 i++;
871 ASSERT_NEAR(segments[i].first.GetX(), 50.000000000000000, DELTA_MAX);
872 ASSERT_NEAR(segments[i].first.GetY(), 0.0000000000000000, DELTA_MAX);
873 ASSERT_NEAR(segments[i].second.GetX(), 50.000000000000000, DELTA_MAX);
874 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
875 i++;
876 ASSERT_NEAR(segments[i].first.GetX(), 20.000000000000000, DELTA_MAX);
877 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
878 ASSERT_NEAR(segments[i].second.GetX(), 20.000000000000000, DELTA_MAX);
879 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
880 i++;
881 ASSERT_NEAR(segments[i].first.GetX(), 45.000000000000000, DELTA_MAX);
882 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
883 ASSERT_NEAR(segments[i].second.GetX(), 45.000000000000000, DELTA_MAX);
884 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
885 i++;
886 ASSERT_NEAR(segments[i].first.GetX(), 52.000000000000000, DELTA_MAX);
887 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
888 ASSERT_NEAR(segments[i].second.GetX(), 52.000000000000000, DELTA_MAX);
889 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
890 i++;
891 ASSERT_NEAR(segments[i].first.GetX(), 70.000000000000000, DELTA_MAX);
892 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
893 ASSERT_NEAR(segments[i].second.GetX(), 70.000000000000000, DELTA_MAX);
894 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
895 i++;
896 ASSERT_NEAR(segments[i].first.GetX(), 0.0000000000000000, DELTA_MAX);
897 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
898 ASSERT_NEAR(segments[i].second.GetX(), 0.0000000000000000, DELTA_MAX);
899 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
900 i++;
901 ASSERT_NEAR(segments[i].first.GetX(), 32.000000000000000, DELTA_MAX);
902 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
903 ASSERT_NEAR(segments[i].second.GetX(), 32.000000000000000, DELTA_MAX);
904 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
905 i++;
906 ASSERT_NEAR(segments[i].first.GetX(), 35.000000000000000, DELTA_MAX);
907 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
908 ASSERT_NEAR(segments[i].second.GetX(), 35.000000000000000, DELTA_MAX);
909 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
910 i++;
911 ASSERT_NEAR(segments[i].first.GetX(), 44.000000000000000, DELTA_MAX);
912 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
913 ASSERT_NEAR(segments[i].second.GetX(), 44.000000000000000, DELTA_MAX);
914 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
915 i++;
916 ASSERT_NEAR(segments[i].first.GetX(), 60.000000000000000, DELTA_MAX);
917 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
918 ASSERT_NEAR(segments[i].second.GetX(), 60.000000000000000, DELTA_MAX);
919 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
920 i++;
921 ASSERT_NEAR(segments[i].first.GetX(), 75.000000000000000, DELTA_MAX);
922 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
923 ASSERT_NEAR(segments[i].second.GetX(), 75.000000000000000, DELTA_MAX);
924 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
925 i++;
926 ASSERT_NEAR(segments[i].first.GetX(), 10.000000000000000, DELTA_MAX);
927 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
928 ASSERT_NEAR(segments[i].second.GetX(), 10.000000000000000, DELTA_MAX);
929 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
930 i++;
931 ASSERT_NEAR(segments[i].first.GetX(), 41.000000000000000, DELTA_MAX);
932 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
933 ASSERT_NEAR(segments[i].second.GetX(), 41.000000000000000, DELTA_MAX);
934 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
935 i++;
936 ASSERT_NEAR(segments[i].first.GetX(), 46.000000000000000, DELTA_MAX);
937 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
938 ASSERT_NEAR(segments[i].second.GetX(), 46.000000000000000, DELTA_MAX);
939 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
940 i++;
941 ASSERT_NEAR(segments[i].first.GetX(), 80.000000000000000, DELTA_MAX);
942 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
943 ASSERT_NEAR(segments[i].second.GetX(), 80.000000000000000, DELTA_MAX);
944 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
945 i++;
946 ASSERT_NEAR(segments[i].first.GetX(), 34.000000000000000, DELTA_MAX);
947 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
948 ASSERT_NEAR(segments[i].second.GetX(), 34.000000000000000, DELTA_MAX);
949 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
950 i++;
951 ASSERT_NEAR(segments[i].first.GetX(), 42.000000000000000, DELTA_MAX);
952 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
953 ASSERT_NEAR(segments[i].second.GetX(), 42.000000000000000, DELTA_MAX);
954 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
955 i++;
956 ASSERT_NEAR(segments[i].first.GetX(), 90.000000000000000, DELTA_MAX);
957 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
958 ASSERT_NEAR(segments[i].second.GetX(), 90.000000000000000, DELTA_MAX);
959 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
960 i++;
961 ASSERT_NEAR(segments[i].first.GetX(), 96.000000000000000, DELTA_MAX);
962 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
963 ASSERT_NEAR(segments[i].second.GetX(), 96.000000000000000, DELTA_MAX);
964 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
965 i++;
966 ASSERT_NEAR(segments[i].first.GetX(), 1.0000000000000000, DELTA_MAX);
967 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
968 ASSERT_NEAR(segments[i].second.GetX(), 1.0000000000000000, DELTA_MAX);
969 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
970 i++;
971 ASSERT_NEAR(segments[i].first.GetX(), 33.000000000000000, DELTA_MAX);
972 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
973 ASSERT_NEAR(segments[i].second.GetX(), 33.000000000000000, DELTA_MAX);
974 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
975 i++;
976 ASSERT_NEAR(segments[i].first.GetX(), 40.000000000000000, DELTA_MAX);
977 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
978 ASSERT_NEAR(segments[i].second.GetX(), 40.000000000000000, DELTA_MAX);
979 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
980 i++;
981 ASSERT_NEAR(segments[i].first.GetX(), 43.000000000000000, DELTA_MAX);
982 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
983 ASSERT_NEAR(segments[i].second.GetX(), 43.000000000000000, DELTA_MAX);
984 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
985 i++;
986 ASSERT_NEAR(segments[i].first.GetX(), 51.000000000000000, DELTA_MAX);
987 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
988 ASSERT_NEAR(segments[i].second.GetX(), 51.000000000000000, DELTA_MAX);
989 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
990 i++;
991 ASSERT_NEAR(segments[i].first.GetX(), 61.000000000000000, DELTA_MAX);
992 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
993 ASSERT_NEAR(segments[i].second.GetX(), 61.000000000000000, DELTA_MAX);
994 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
995 i++;
996 ASSERT_NEAR(segments[i].first.GetX(), 76.000000000000000, DELTA_MAX);
997 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
998 ASSERT_NEAR(segments[i].second.GetX(), 76.000000000000000, DELTA_MAX);
999 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
1000 i++;
1001 ASSERT_NEAR(segments[i].first.GetX(), 95.000000000000000, DELTA_MAX);
1002 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1003 ASSERT_NEAR(segments[i].second.GetX(), 95.000000000000000, DELTA_MAX);
1004 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
1005 i++;
1006 ASSERT_NEAR(segments[i].first.GetX(), 5.0000000000000000, DELTA_MAX);
1007 ASSERT_NEAR(segments[i].first.GetY(), 0.0000000000000000, DELTA_MAX);
1008 ASSERT_NEAR(segments[i].second.GetX(), 31.000000000000000, DELTA_MAX);
1009 ASSERT_NEAR(segments[i].second.GetY(), 0.0000000000000000, DELTA_MAX);
1010 i++;
1011 ASSERT_NEAR(segments[i].first.GetX(), 36.000000000000000, DELTA_MAX);
1012 ASSERT_NEAR(segments[i].first.GetY(), 0.0000000000000000, DELTA_MAX);
1013 ASSERT_NEAR(segments[i].second.GetX(), 50.000000000000000, DELTA_MAX);
1014 ASSERT_NEAR(segments[i].second.GetY(), 0.0000000000000000, DELTA_MAX);
1015 i++;
1016 ASSERT_NEAR(segments[i].first.GetX(), 5.0000000000000000, DELTA_MAX);
1017 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
1018 ASSERT_NEAR(segments[i].second.GetX(), 20.000000000000000, DELTA_MAX);
1019 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
1020 i++;
1021 ASSERT_NEAR(segments[i].first.GetX(), 31.000000000000000, DELTA_MAX);
1022 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
1023 ASSERT_NEAR(segments[i].second.GetX(), 36.000000000000000, DELTA_MAX);
1024 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
1025 i++;
1026 ASSERT_NEAR(segments[i].first.GetX(), 45.000000000000000, DELTA_MAX);
1027 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
1028 ASSERT_NEAR(segments[i].second.GetX(), 50.000000000000000, DELTA_MAX);
1029 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
1030 i++;
1031 ASSERT_NEAR(segments[i].first.GetX(), 52.000000000000000, DELTA_MAX);
1032 ASSERT_NEAR(segments[i].first.GetY(), 5.0000000000000000, DELTA_MAX);
1033 ASSERT_NEAR(segments[i].second.GetX(), 70.000000000000000, DELTA_MAX);
1034 ASSERT_NEAR(segments[i].second.GetY(), 5.0000000000000000, DELTA_MAX);
1035 i++;
1036 ASSERT_NEAR(segments[i].first.GetX(), 0.0000000000000000, DELTA_MAX);
1037 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
1038 ASSERT_NEAR(segments[i].second.GetX(), 20.000000000000000, DELTA_MAX);
1039 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
1040 i++;
1041 ASSERT_NEAR(segments[i].first.GetX(), 32.000000000000000, DELTA_MAX);
1042 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
1043 ASSERT_NEAR(segments[i].second.GetX(), 35.000000000000000, DELTA_MAX);
1044 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
1045 i++;
1046 ASSERT_NEAR(segments[i].first.GetX(), 44.000000000000000, DELTA_MAX);
1047 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
1048 ASSERT_NEAR(segments[i].second.GetX(), 45.000000000000000, DELTA_MAX);
1049 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
1050 i++;
1051 ASSERT_NEAR(segments[i].first.GetX(), 52.000000000000000, DELTA_MAX);
1052 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
1053 ASSERT_NEAR(segments[i].second.GetX(), 60.000000000000000, DELTA_MAX);
1054 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
1055 i++;
1056 ASSERT_NEAR(segments[i].first.GetX(), 70.000000000000000, DELTA_MAX);
1057 ASSERT_NEAR(segments[i].first.GetY(), 10.000000000000000, DELTA_MAX);
1058 ASSERT_NEAR(segments[i].second.GetX(), 75.000000000000000, DELTA_MAX);
1059 ASSERT_NEAR(segments[i].second.GetY(), 10.000000000000000, DELTA_MAX);
1060 i++;
1061 ASSERT_NEAR(segments[i].first.GetX(), 0.0000000000000000, DELTA_MAX);
1062 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
1063 ASSERT_NEAR(segments[i].second.GetX(), 10.000000000000000, DELTA_MAX);
1064 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
1065 i++;
1066 ASSERT_NEAR(segments[i].first.GetX(), 32.000000000000000, DELTA_MAX);
1067 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
1068 ASSERT_NEAR(segments[i].second.GetX(), 35.000000000000000, DELTA_MAX);
1069 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
1070 i++;
1071 ASSERT_NEAR(segments[i].first.GetX(), 41.000000000000000, DELTA_MAX);
1072 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
1073 ASSERT_NEAR(segments[i].second.GetX(), 44.000000000000000, DELTA_MAX);
1074 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
1075 i++;
1076 ASSERT_NEAR(segments[i].first.GetX(), 46.000000000000000, DELTA_MAX);
1077 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
1078 ASSERT_NEAR(segments[i].second.GetX(), 60.000000000000000, DELTA_MAX);
1079 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
1080 i++;
1081 ASSERT_NEAR(segments[i].first.GetX(), 75.000000000000000, DELTA_MAX);
1082 ASSERT_NEAR(segments[i].first.GetY(), 15.000000000000000, DELTA_MAX);
1083 ASSERT_NEAR(segments[i].second.GetX(), 80.000000000000000, DELTA_MAX);
1084 ASSERT_NEAR(segments[i].second.GetY(), 15.000000000000000, DELTA_MAX);
1085 i++;
1086 ASSERT_NEAR(segments[i].first.GetX(), 10.000000000000000, DELTA_MAX);
1087 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
1088 ASSERT_NEAR(segments[i].second.GetX(), 34.000000000000000, DELTA_MAX);
1089 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
1090 i++;
1091 ASSERT_NEAR(segments[i].first.GetX(), 41.000000000000000, DELTA_MAX);
1092 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
1093 ASSERT_NEAR(segments[i].second.GetX(), 42.000000000000000, DELTA_MAX);
1094 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
1095 i++;
1096 ASSERT_NEAR(segments[i].first.GetX(), 46.000000000000000, DELTA_MAX);
1097 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
1098 ASSERT_NEAR(segments[i].second.GetX(), 80.000000000000000, DELTA_MAX);
1099 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
1100 i++;
1101 ASSERT_NEAR(segments[i].first.GetX(), 90.000000000000000, DELTA_MAX);
1102 ASSERT_NEAR(segments[i].first.GetY(), 20.000000000000000, DELTA_MAX);
1103 ASSERT_NEAR(segments[i].second.GetX(), 96.000000000000000, DELTA_MAX);
1104 ASSERT_NEAR(segments[i].second.GetY(), 20.000000000000000, DELTA_MAX);
1105 i++;
1106 ASSERT_NEAR(segments[i].first.GetX(), 1.0000000000000000, DELTA_MAX);
1107 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1108 ASSERT_NEAR(segments[i].second.GetX(), 33.000000000000000, DELTA_MAX);
1109 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
1110 i++;
1111 ASSERT_NEAR(segments[i].first.GetX(), 34.000000000000000, DELTA_MAX);
1112 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1113 ASSERT_NEAR(segments[i].second.GetX(), 40.000000000000000, DELTA_MAX);
1114 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
1115 i++;
1116 ASSERT_NEAR(segments[i].first.GetX(), 42.000000000000000, DELTA_MAX);
1117 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1118 ASSERT_NEAR(segments[i].second.GetX(), 43.000000000000000, DELTA_MAX);
1119 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
1120 i++;
1121 ASSERT_NEAR(segments[i].first.GetX(), 51.000000000000000, DELTA_MAX);
1122 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1123 ASSERT_NEAR(segments[i].second.GetX(), 61.000000000000000, DELTA_MAX);
1124 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
1125 i++;
1126 ASSERT_NEAR(segments[i].first.GetX(), 76.000000000000000, DELTA_MAX);
1127 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1128 ASSERT_NEAR(segments[i].second.GetX(), 90.000000000000000, DELTA_MAX);
1129 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
1130 i++;
1131 ASSERT_NEAR(segments[i].first.GetX(), 95.000000000000000, DELTA_MAX);
1132 ASSERT_NEAR(segments[i].first.GetY(), 25.000000000000000, DELTA_MAX);
1133 ASSERT_NEAR(segments[i].second.GetX(), 96.000000000000000, DELTA_MAX);
1134 ASSERT_NEAR(segments[i].second.GetY(), 25.000000000000000, DELTA_MAX);
1135 i++;
1136 ASSERT_NEAR(segments[i].first.GetX(), 1.0000000000000000, DELTA_MAX);
1137 ASSERT_NEAR(segments[i].first.GetY(), 30.000000000000000, DELTA_MAX);
1138 ASSERT_NEAR(segments[i].second.GetX(), 33.000000000000000, DELTA_MAX);
1139 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
1140 i++;
1141 ASSERT_NEAR(segments[i].first.GetX(), 40.000000000000000, DELTA_MAX);
1142 ASSERT_NEAR(segments[i].first.GetY(), 30.000000000000000, DELTA_MAX);
1143 ASSERT_NEAR(segments[i].second.GetX(), 43.000000000000000, DELTA_MAX);
1144 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
1145 i++;
1146 ASSERT_NEAR(segments[i].first.GetX(), 51.000000000000000, DELTA_MAX);
1147 ASSERT_NEAR(segments[i].first.GetY(), 30.000000000000000, DELTA_MAX);
1148 ASSERT_NEAR(segments[i].second.GetX(), 61.000000000000000, DELTA_MAX);
1149 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
1150 i++;
1151 ASSERT_NEAR(segments[i].first.GetX(), 76.000000000000000, DELTA_MAX);
1152 ASSERT_NEAR(segments[i].first.GetY(), 30.000000000000000, DELTA_MAX);
1153 ASSERT_NEAR(segments[i].second.GetX(), 95.000000000000000, DELTA_MAX);
1154 ASSERT_NEAR(segments[i].second.GetY(), 30.000000000000000, DELTA_MAX);
1155 }
1156
1157 #if defined(BGO_ENABLE_DICOMSTRUCTURESETLOADER2) && (ORTHANC_SANDBOXED != 1)
1158
1159 #include <SystemToolbox.h>
1160
1161 TEST(StructureSet2, ReadFromJsonPart2)
1162 {
1163 DicomStructureSet2 structureSet;
1164 std::string jsonText;
1165
1166 Orthanc::SystemToolbox::ReadFile(jsonText, "72c773ac-5059f2c4-2e6a9120-4fd4bca1-45701661.json");
1167
1168 FullOrthancDataset dicom(jsonText);
1169 structureSet.Clear();
1170
1171 structureSet.FillStructuresFromDataset(dicom);
1172 structureSet.ComputeDependentProperties();
1173
1174 //const std::vector<DicomStructure2>& structures = structureSet.structures_;
1175 }
1176
1177 #endif
1178 // BGO_ENABLE_DICOMSTRUCTURESETLOADER2
1179 29
1180 30
1181 TEST(StructureSet, ReadFromJson) 31 TEST(StructureSet, ReadFromJson)
1182 { 32 {
1183 FullOrthancDataset dicom(Orthanc::EmbeddedResources::GetFileResourceBuffer(Orthanc::EmbeddedResources::RT_STRUCT_00), 33 OrthancStone::FullOrthancDataset dicom(
1184 Orthanc::EmbeddedResources::GetFileResourceSize(Orthanc::EmbeddedResources::RT_STRUCT_00)); 34 Orthanc::EmbeddedResources::GetFileResourceBuffer(Orthanc::EmbeddedResources::RT_STRUCT_00),
35 Orthanc::EmbeddedResources::GetFileResourceSize(Orthanc::EmbeddedResources::RT_STRUCT_00));
1185 36
1186 DicomStructureSet rtstruct(dicom); 37 OrthancStone::DicomStructureSet rtstruct(dicom);
1187 38
1188 ASSERT_DOUBLE_EQ(0.0, rtstruct.GetEstimatedNormal() [0]); 39 ASSERT_DOUBLE_EQ(0.0, rtstruct.GetEstimatedNormal() [0]);
1189 ASSERT_DOUBLE_EQ(0.0, rtstruct.GetEstimatedNormal() [1]); 40 ASSERT_DOUBLE_EQ(0.0, rtstruct.GetEstimatedNormal() [1]);
1190 ASSERT_DOUBLE_EQ(1.0, rtstruct.GetEstimatedNormal() [2]); 41 ASSERT_DOUBLE_EQ(1.0, rtstruct.GetEstimatedNormal() [2]);
1191 ASSERT_DOUBLE_EQ(3.0, rtstruct.GetEstimatedSliceThickness()); 42 ASSERT_DOUBLE_EQ(3.0, rtstruct.GetEstimatedSliceThickness());