Mercurial > hg > orthanc-stone
comparison OrthancStone/UnitTestsSources/TestStructureSet.cpp @ 1908:affde38b84de
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()); |