Mercurial > hg > orthanc
annotate OrthancServer/OrthancRestApi/RadiotherapyRestApi.cpp @ 767:c19552f604d5 dicom-rt
fix
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 16 Apr 2014 16:15:58 +0200 |
parents | 2f7802e95dd0 |
children | e57e08ed510f |
rev | line source |
---|---|
526 | 1 /** |
2 * Orthanc - A Lightweight, RESTful DICOM Store | |
3 * Copyright (C) 2012-2013 Medical Physics Department, CHU of Liege, | |
4 * Belgium | |
5 * | |
6 * This program is free software: you can redistribute it and/or | |
7 * modify it under the terms of the GNU General Public License as | |
8 * published by the Free Software Foundation, either version 3 of the | |
9 * License, or (at your option) any later version. | |
10 * | |
11 * In addition, as a special exception, the copyright holders of this | |
12 * program give permission to link the code of its release with the | |
13 * OpenSSL project's "OpenSSL" library (or with modified versions of it | |
14 * that use the same license as the "OpenSSL" library), and distribute | |
15 * the linked executables. You must obey the GNU General Public License | |
16 * in all respects for all of the code used other than "OpenSSL". If you | |
17 * modify file(s) with this exception, you may extend this exception to | |
18 * your version of the file(s), but you are not obligated to do so. If | |
19 * you do not wish to do so, delete this exception statement from your | |
20 * version. If you delete this exception statement from all source files | |
21 * in the program, then also delete it here. | |
22 * | |
23 * This program is distributed in the hope that it will be useful, but | |
24 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
25 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
26 * General Public License for more details. | |
27 * | |
28 * You should have received a copy of the GNU General Public License | |
29 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
30 **/ | |
31 | |
32 | |
764 | 33 #include "OrthancRestApi.h" |
526 | 34 |
764 | 35 #include "../ServerToolbox.h" |
36 | |
526 | 37 |
38 | |
532
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
39 // DICOM tags for RT-STRUCT |
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
40 |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
41 /** |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
42 * REFERENCE: http://www.dabsoft.ch/dicom/3/C.8.8.6/ |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
43 * |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
44 * IMPORTANT: The points/vertices coordinates are reported in [mm]. |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
45 * |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
46 * TODO: Support "Contour Offset Vector" |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
47 **/ |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
48 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
49 |
526 | 50 #define REFERENCED_STUDY_SEQUENCE "0008,1110" |
51 #define REFERENCED_SOP_INSTANCE_UID "0008,1155" | |
52 #define FRAME_OF_REFERENCE_UID "0020,0052" | |
53 #define REFERENCED_FRAME_OF_REFERENCE_SEQUENCE "3006,0010" | |
529 | 54 #define STRUCTURE_SET_ROI_SEQUENCE "3006,0020" |
55 #define ROI_NUMBER "3006,0022" | |
56 #define ROI_NAME "3006,0026" | |
57 #define ROI_GENERATION_ALGORITHM "3006,0036" | |
58 #define ROI_CONTOUR_SEQUENCE "3006,0039" | |
59 #define REFERENCED_ROI_NUMBER "3006,0084" | |
60 #define ROI_DISPLAY_COLOR "3006,002a" | |
532
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
61 #define CONTOUR_SEQUENCE "3006,0040" |
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
62 #define CONTOUR_IMAGE_SEQUENCE "3006,0016" |
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
63 #define CONTOUR_GEOMETRIC_TYPE "3006,0042" |
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
64 #define NUMBER_OF_CONTOUR_POINTS "3006,0046" |
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
65 #define CONTOUR_DATA "3006,0050" |
535 | 66 #define CONTOUR_SLAB_THICKNESS "3006,0044" |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
67 #define SLICE_THICKNESS "0018,0050" |
532
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
68 |
526 | 69 |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
70 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
71 #include <boost/geometry.hpp> |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
72 #include <boost/geometry/geometries/point.hpp> |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
73 #include <boost/geometry/geometries/point_xy.hpp> |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
74 #include <boost/geometry/geometries/polygon.hpp> |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
75 #include <boost/geometry/geometries/linestring.hpp> |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
76 |
526 | 77 namespace Orthanc |
78 { | |
79 static bool CheckSeriesModality(Json::Value& study, | |
80 Json::Value& series, | |
81 Json::Value& content, | |
82 ServerContext& context, | |
83 const std::string& seriesId, | |
84 const std::string& modality) | |
85 { | |
86 if (!context.GetIndex().LookupResource(series, seriesId, ResourceType_Series)) | |
87 { | |
88 return false; | |
89 } | |
90 | |
91 // Retrieve the parent study | |
92 std::string studyId = series["ParentStudy"].asString(); | |
93 if (!context.GetIndex().LookupResource(study, studyId, ResourceType_Study)) | |
94 { | |
95 return false; | |
96 } | |
97 | |
98 // Check the modality and that there is a single instance inside the series | |
99 if (!series["MainDicomTags"].isMember("Modality") || | |
100 series["MainDicomTags"]["Modality"].asString() != modality || | |
101 series["Instances"].size() != 1) | |
102 { | |
103 return false; | |
104 } | |
105 | |
106 // Retrieve the instance data | |
107 std::string instanceId = series["Instances"][0].asString(); | |
108 | |
109 context.ReadJson(content, instanceId); | |
110 | |
111 return true; | |
112 } | |
113 | |
114 | |
535 | 115 static bool ContourToPoints(Json::Value& result, |
116 const Json::Value& source) | |
117 { | |
118 std::vector<std::string> points; | |
764 | 119 Toolbox::TokenizeString(points, source.asString(), '\\'); |
535 | 120 |
121 if (points.size() % 3 != 0) | |
122 { | |
123 return false; | |
124 } | |
125 | |
126 result = Json::arrayValue; | |
127 | |
128 for (size_t k = 0; k < points.size(); k += 3) | |
129 { | |
130 Json::Value p = Json::arrayValue; | |
131 | |
132 try | |
133 { | |
134 p.append(boost::lexical_cast<float>(points[k])); | |
135 p.append(boost::lexical_cast<float>(points[k + 1])); | |
136 p.append(boost::lexical_cast<float>(points[k + 2])); | |
137 } | |
138 catch (boost::bad_lexical_cast) | |
139 { | |
140 return false; | |
141 } | |
142 | |
143 result.append(p); | |
144 } | |
145 | |
146 return true; | |
147 } | |
148 | |
149 | |
526 | 150 static bool GetRtStructuresInfo(Json::Value& study, |
151 Json::Value& series, | |
152 Json::Value& content, | |
153 std::string& frameOfReference, | |
154 ServerContext& context, | |
155 const std::string& seriesId) | |
156 { | |
157 if (!CheckSeriesModality(study, series, content, context, seriesId, "RTSTRUCT")) | |
158 { | |
159 return false; | |
160 } | |
161 | |
544 | 162 // Check that the "ReferencedStudySequence" (if any) is the same as the parent study. |
163 if (content.isMember(REFERENCED_STUDY_SEQUENCE)) | |
526 | 164 { |
544 | 165 if (content[REFERENCED_STUDY_SEQUENCE]["Value"].size() != 1 || |
166 !content[REFERENCED_STUDY_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID) || | |
167 content[REFERENCED_STUDY_SEQUENCE]["Value"][0][REFERENCED_SOP_INSTANCE_UID]["Value"].asString() != | |
168 study["MainDicomTags"]["StudyInstanceUID"].asString()) | |
169 { | |
170 return false; | |
171 } | |
526 | 172 } |
173 | |
174 // Lookup for the frame of reference. Orthanc does not support | |
175 // RTSTRUCT with multiple frames of reference. | |
176 if (!content.isMember(REFERENCED_FRAME_OF_REFERENCE_SEQUENCE) || | |
177 content[REFERENCED_FRAME_OF_REFERENCE_SEQUENCE]["Value"].size() != 1 || | |
178 !content[REFERENCED_FRAME_OF_REFERENCE_SEQUENCE]["Value"][0].isMember(FRAME_OF_REFERENCE_UID)) | |
179 { | |
180 return false; | |
181 } | |
182 | |
183 frameOfReference = content[REFERENCED_FRAME_OF_REFERENCE_SEQUENCE]["Value"][0][FRAME_OF_REFERENCE_UID]["Value"].asString(); | |
184 | |
185 return true; | |
186 } | |
187 | |
188 | |
535 | 189 static bool GetRtStructuresRoi(Json::Value& result, |
190 Json::Value& contourSequence, | |
191 std::string& instanceId, | |
192 ServerContext& context, | |
193 const std::string& seriesId, | |
194 const std::string& roiNumber) | |
195 { | |
196 Json::Value study, series, content; | |
197 std::string frameOfReference; | |
198 | |
199 if (!GetRtStructuresInfo(study, series, content, frameOfReference, context, seriesId)) | |
200 { | |
201 return false; | |
202 } | |
203 | |
204 if (!content.isMember(STRUCTURE_SET_ROI_SEQUENCE) || | |
205 !content.isMember(ROI_CONTOUR_SEQUENCE)) | |
206 { | |
207 return false; | |
208 } | |
209 | |
210 instanceId = series["Instances"][0].asString(); | |
211 | |
212 bool found = false; | |
213 | |
214 for (Json::Value::ArrayIndex i = 0; i < content[STRUCTURE_SET_ROI_SEQUENCE]["Value"].size(); i++) | |
215 { | |
216 const Json::Value& roi = content[STRUCTURE_SET_ROI_SEQUENCE]["Value"][i]; | |
217 | |
218 if (roi.isMember(ROI_NUMBER) && | |
219 roi.isMember(ROI_NAME) && | |
220 roi[ROI_NUMBER]["Value"].asString() == roiNumber) | |
221 { | |
222 result["InternalIndex"] = i; | |
223 result["Number"] = boost::lexical_cast<unsigned int>(roiNumber); | |
224 result["Name"] = roi[ROI_NAME]["Value"].asString(); | |
225 result["GenerationAlgorithm"] = roi[ROI_GENERATION_ALGORITHM]["Value"].asString(); | |
226 found = true; | |
227 } | |
228 } | |
229 | |
230 if (!found) | |
231 { | |
232 return false; | |
233 } | |
234 | |
235 for (Json::Value::ArrayIndex i = 0; i < content[ROI_CONTOUR_SEQUENCE]["Value"].size(); i++) | |
236 { | |
237 const Json::Value& contour = content[ROI_CONTOUR_SEQUENCE]["Value"][i]; | |
238 | |
239 if (contour.isMember(REFERENCED_ROI_NUMBER) && | |
240 contour.isMember(ROI_DISPLAY_COLOR) && | |
241 contour.isMember(CONTOUR_SEQUENCE) && | |
242 contour[REFERENCED_ROI_NUMBER]["Value"].asString() == roiNumber) | |
243 { | |
244 std::vector<std::string> color; | |
764 | 245 Toolbox::TokenizeString(color, contour[ROI_DISPLAY_COLOR]["Value"].asString(), '\\'); |
535 | 246 |
247 result["DisplayColor"] = Json::arrayValue; | |
248 if (color.size() != 3) | |
249 { | |
250 return false; | |
251 } | |
252 | |
253 for (size_t k = 0; k < color.size(); k++) | |
254 { | |
255 try | |
256 { | |
257 result["DisplayColor"].append(boost::lexical_cast<int>(color[k])); | |
258 } | |
259 catch (boost::bad_lexical_cast) | |
260 { | |
261 return false; | |
262 } | |
263 } | |
264 | |
265 contourSequence = contour[CONTOUR_SEQUENCE]["Value"]; | |
266 | |
267 return true; | |
268 } | |
269 } | |
270 | |
271 return false; | |
272 } | |
273 | |
274 | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
275 static bool GetClosedPlanarPoints(Json::Value& result, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
276 ServerContext& context, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
277 const std::string& instanceId, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
278 const Json::Value& roi, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
279 unsigned int index) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
280 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
281 boost::mutex::scoped_lock lock(context.GetDicomFileMutex()); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
282 ParsedDicomFile& dicom = context.GetDicomFile(instanceId); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
283 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
284 ParsedDicomFile::SequencePath path; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
285 path.push_back(std::make_pair(DicomTag(0x3006, 0x0039 /* ROIContourSequence */), roi["InternalIndex"].asInt())); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
286 path.push_back(std::make_pair(DicomTag(0x3006, 0x0040 /* ContourSequence */), index)); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
287 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
288 std::string contourData; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
289 std::string numberOfPoints; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
290 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
291 if (dicom.GetTagValue(contourData, path, DicomTag(0x3006, 0x0050 /* ContourData */)) && |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
292 dicom.GetTagValue(numberOfPoints, path, DicomTag(0x3006, 0x0046 /* NumberOfContourPoints */)) && |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
293 ContourToPoints(result, contourData) && |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
294 result.size() == boost::lexical_cast<unsigned int>(numberOfPoints)) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
295 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
296 return true; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
297 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
298 else |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
299 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
300 return false; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
301 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
302 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
303 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
304 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
305 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
306 static bool LookupReferencedInstance(Json::Value& result, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
307 ServerContext& context, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
308 const Json::Value& contour) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
309 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
310 if (contour.isMember(CONTOUR_IMAGE_SEQUENCE) && |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
311 contour[CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
312 contour[CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID)) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
313 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
314 std::string uid = contour[CONTOUR_IMAGE_SEQUENCE]["Value"][0][REFERENCED_SOP_INSTANCE_UID]["Value"].asString(); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
315 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
316 std::list<std::string> instance; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
317 context.GetIndex().LookupTagValue(instance, DICOM_TAG_SOP_INSTANCE_UID, uid); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
318 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
319 if (instance.size() == 1 && |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
320 context.GetIndex().LookupResource(result, instance.front(), ResourceType_Instance)) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
321 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
322 return true; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
323 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
324 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
325 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
326 return false; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
327 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
328 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
329 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
330 static bool GetRtStructuresClosedPlanarThickness(float& result, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
331 ServerContext& context, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
332 const Json::Value& contour) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
333 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
334 if (contour.isMember(CONTOUR_SLAB_THICKNESS)) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
335 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
336 result = boost::lexical_cast<float>(contour[CONTOUR_SLAB_THICKNESS]["Value"].asString()); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
337 return true; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
338 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
339 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
340 // No slab thickness is explicitely specified: Fallback to |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
341 // the thickness of the referred instance |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
342 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
343 Json::Value instance; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
344 if (LookupReferencedInstance(instance, context, contour)) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
345 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
346 Json::Value info; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
347 context.ReadJson(info, instance["ID"].asString()); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
348 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
349 if (info.isMember(SLICE_THICKNESS)) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
350 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
351 result = boost::lexical_cast<float>(info[SLICE_THICKNESS]["Value"].asString()); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
352 return true; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
353 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
354 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
355 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
356 return false; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
357 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
358 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
359 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
360 |
526 | 361 static void GetRtStructuresInfo(RestApi::GetCall& call) |
362 { | |
764 | 363 ServerContext& context = OrthancRestApi::GetContext(call); |
526 | 364 |
365 Json::Value study, series, content; | |
366 std::string frameOfReference; | |
367 if (GetRtStructuresInfo(study, series, content, frameOfReference, context, call.GetUriComponent("id", ""))) | |
368 { | |
369 Json::Value result; | |
370 | |
371 result["Study"] = study["ID"]; | |
372 | |
373 | |
374 // Lookup the series with the same frame of reference inside this study | |
375 result["RelatedSeries"] = Json::arrayValue; | |
376 | |
377 for (Json::Value::ArrayIndex i = 0; i < study["Series"].size(); i++) | |
378 { | |
379 Json::Value otherSeries; | |
380 if (context.GetIndex().LookupResource(otherSeries, study["Series"][i].asString(), ResourceType_Series) && | |
381 otherSeries["Instances"].size() > 0) | |
382 { | |
383 Json::Value info; | |
384 context.ReadJson(info, otherSeries["Instances"][0].asString()); | |
385 | |
386 if (info.isMember(FRAME_OF_REFERENCE_UID)) | |
387 { | |
388 result["RelatedSeries"].append(study["Series"][i].asString()); | |
389 } | |
390 } | |
391 } | |
392 | |
529 | 393 |
526 | 394 call.GetOutput().AnswerJson(result); |
395 } | |
396 } | |
397 | |
398 | |
529 | 399 static void GetRtStructuresListOfROIs(RestApi::GetCall& call) |
400 { | |
764 | 401 ServerContext& context = OrthancRestApi::GetContext(call); |
529 | 402 |
403 Json::Value study, series, content; | |
404 std::string frameOfReference; | |
405 if (GetRtStructuresInfo(study, series, content, frameOfReference, context, call.GetUriComponent("id", ""))) | |
406 { | |
407 Json::Value result(Json::arrayValue); | |
408 | |
409 if (content.isMember(STRUCTURE_SET_ROI_SEQUENCE)) | |
410 { | |
411 for (Json::Value::ArrayIndex i = 0; i < content[STRUCTURE_SET_ROI_SEQUENCE]["Value"].size(); i++) | |
412 { | |
413 if (content[STRUCTURE_SET_ROI_SEQUENCE]["Value"][i].isMember(ROI_NUMBER)) | |
414 { | |
535 | 415 result.append(boost::lexical_cast<int>(content[STRUCTURE_SET_ROI_SEQUENCE]["Value"][i][ROI_NUMBER]["Value"].asString())); |
529 | 416 } |
417 } | |
418 } | |
419 | |
420 call.GetOutput().AnswerJson(result); | |
421 } | |
422 } | |
423 | |
424 | |
425 static void GetRtStructuresROI(RestApi::GetCall& call) | |
426 { | |
764 | 427 ServerContext& context = OrthancRestApi::GetContext(call); |
535 | 428 |
764 | 429 Json::Value roi, contour; |
430 std::string instanceId; | |
535 | 431 |
764 | 432 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
433 call.GetUriComponent("id", ""), | |
434 call.GetUriComponent("roi", ""))) | |
435 { | |
436 roi.removeMember("InternalIndex"); | |
437 call.GetOutput().AnswerJson(roi); | |
438 } | |
535 | 439 } |
440 | |
441 | |
442 static void GetRtStructuresROIPoints(RestApi::GetCall& call) | |
443 { | |
764 | 444 ServerContext& context = OrthancRestApi::GetContext(call); |
535 | 445 |
764 | 446 Json::Value roi, contour; |
447 std::string instanceId; | |
535 | 448 |
764 | 449 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
450 call.GetUriComponent("id", ""), | |
451 call.GetUriComponent("roi", ""))) | |
452 { | |
453 Json::Value result = Json::arrayValue; | |
535 | 454 |
764 | 455 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
456 { | |
457 if (contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "POINT") | |
458 { | |
459 Json::Value p; | |
460 if (ContourToPoints(p, contour[i][CONTOUR_DATA]["Value"].asString()) && | |
461 p.size() == 1) | |
462 { | |
463 result.append(p[0]); | |
464 } | |
465 } | |
466 } | |
535 | 467 |
764 | 468 call.GetOutput().AnswerJson(result); |
469 } | |
535 | 470 } |
471 | |
472 | |
473 static void GetRtStructuresListOfClosedPlanars(RestApi::GetCall& call) | |
474 { | |
764 | 475 ServerContext& context = OrthancRestApi::GetContext(call); |
535 | 476 |
764 | 477 Json::Value roi, contour; |
478 std::string instanceId; | |
535 | 479 |
764 | 480 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
481 call.GetUriComponent("id", ""), | |
482 call.GetUriComponent("roi", ""))) | |
483 { | |
484 Json::Value result = Json::arrayValue; | |
535 | 485 |
764 | 486 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
487 { | |
488 if (contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
489 contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && | |
490 contour[i].isMember(NUMBER_OF_CONTOUR_POINTS) && | |
491 contour[i].isMember(CONTOUR_DATA) && | |
492 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
493 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID) && | |
494 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "CLOSED_PLANAR") | |
495 { | |
496 result.append(i); | |
497 } | |
498 } | |
535 | 499 |
764 | 500 call.GetOutput().AnswerJson(result); |
501 } | |
535 | 502 } |
503 | |
504 | |
505 static void GetRtStructuresSingleClosedPlanar(RestApi::GetCall& call) | |
506 { | |
764 | 507 ServerContext& context = OrthancRestApi::GetContext(call); |
535 | 508 |
764 | 509 Json::Value roi, contour, result; |
510 std::string instanceId; | |
535 | 511 |
764 | 512 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
513 call.GetUriComponent("id", ""), | |
514 call.GetUriComponent("roi", ""))) | |
515 { | |
516 unsigned int index = boost::lexical_cast<unsigned int>(call.GetUriComponent("polygon", "")); | |
535 | 517 |
764 | 518 if (GetClosedPlanarPoints(result, context, instanceId, roi, index)) |
519 { | |
520 call.GetOutput().AnswerJson(result); | |
521 } | |
522 } | |
535 | 523 } |
524 | |
525 | |
526 static void GetRtStructuresClosedPlanarThickness(RestApi::GetCall& call) | |
527 { | |
764 | 528 ServerContext& context = OrthancRestApi::GetContext(call); |
535 | 529 |
764 | 530 Json::Value roi, contour; |
531 std::string instanceId; | |
535 | 532 |
764 | 533 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
534 call.GetUriComponent("id", ""), | |
535 call.GetUriComponent("roi", ""))) | |
536 { | |
537 unsigned int index = boost::lexical_cast<unsigned int>(call.GetUriComponent("polygon", "")); | |
535 | 538 |
764 | 539 float thickness; |
540 if (GetRtStructuresClosedPlanarThickness(thickness, context, contour[index])) | |
541 { | |
542 call.GetOutput().AnswerBuffer(boost::lexical_cast<std::string>(thickness), "text/plain"); | |
543 } | |
544 } | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
545 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
546 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
547 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
548 static bool ComputeClosedPlanarArea(double& area, |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
549 const Json::Value& vertices) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
550 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
551 if (vertices.size() <= 1) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
552 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
553 area = 0; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
554 return true; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
555 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
556 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
557 // Check that all the points share the same z coordinates |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
558 for (Json::Value::ArrayIndex i = 1; i < vertices.size(); i++) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
559 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
560 static float THRESHOLD = 10.0f * std::numeric_limits<float>::epsilon(); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
561 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
562 assert(vertices[i].size() == 3); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
563 if (fabs(vertices[i][2].asFloat() - vertices[0][2].asFloat()) > THRESHOLD) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
564 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
565 // This point has not the same z coordinate |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
566 return false; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
567 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
568 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
569 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
570 // Calculate the area of a cartesian polygon |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
571 // TODO - What happens if self-crossing polygon? |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
572 typedef boost::geometry::model::d2::point_xy<float> point_type; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
573 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
574 boost::geometry::model::linestring<point_type> points; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
575 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
576 points.resize(vertices.size()); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
577 for (Json::Value::ArrayIndex i = 0; i < vertices.size(); i++) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
578 { |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
579 float x = vertices[i][0].asFloat(); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
580 float y = vertices[i][1].asFloat(); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
581 points[i] = boost::geometry::make<point_type>(x, y); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
582 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
583 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
584 boost::geometry::model::polygon<point_type> poly; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
585 boost::geometry::append(poly, points); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
586 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
587 area = abs(boost::geometry::area(poly)); |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
588 return true; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
589 } |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
590 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
591 |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
592 static void GetRtStructuresClosedPlanarArea(RestApi::GetCall& call) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
593 { |
764 | 594 ServerContext& context = OrthancRestApi::GetContext(call); |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
595 |
764 | 596 Json::Value roi, contour, vertices; |
597 std::string instanceId; | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
598 |
764 | 599 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
600 call.GetUriComponent("id", ""), | |
601 call.GetUriComponent("roi", ""))) | |
602 { | |
603 unsigned int index = boost::lexical_cast<unsigned int>(call.GetUriComponent("polygon", "")); | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
604 |
764 | 605 double area; |
606 if (GetClosedPlanarPoints(vertices, context, instanceId, roi, index) && | |
607 ComputeClosedPlanarArea(area, vertices)) | |
608 { | |
609 call.GetOutput().AnswerBuffer(boost::lexical_cast<std::string>(area), "text/plain"); | |
610 } | |
611 } | |
535 | 612 } |
613 | |
614 | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
615 static void GetRtStructuresInstanceOfClosedPlanar(RestApi::GetCall& call) |
535 | 616 { |
764 | 617 ServerContext& context = OrthancRestApi::GetContext(call); |
529 | 618 |
535 | 619 Json::Value roi, contour; |
620 std::string instanceId; | |
532
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
621 |
535 | 622 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
623 call.GetUriComponent("id", ""), | |
624 call.GetUriComponent("roi", ""))) | |
625 { | |
626 unsigned int index = boost::lexical_cast<unsigned int>(call.GetUriComponent("polygon", "")); | |
532
b22312081388
extract roi geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
529
diff
changeset
|
627 |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
628 Json::Value result; |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
629 if (LookupReferencedInstance(result, context, contour[index])) |
535 | 630 { |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
631 call.GetOutput().AnswerJson(result); |
529 | 632 } |
633 } | |
634 } | |
635 | |
636 | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
637 static void GetRtStructuresListOfInstances(RestApi::GetCall& call) |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
638 { |
764 | 639 ServerContext& context = OrthancRestApi::GetContext(call); |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
640 |
764 | 641 Json::Value roi, contour; |
642 std::string instanceId; | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
643 |
764 | 644 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
645 call.GetUriComponent("id", ""), | |
646 call.GetUriComponent("roi", ""))) | |
647 { | |
648 Json::Value result = Json::arrayValue; | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
649 |
764 | 650 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
651 { | |
652 if (contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
653 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
654 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID)) | |
655 { | |
656 std::string uid = contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0][REFERENCED_SOP_INSTANCE_UID]["Value"].asString(); | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
657 |
764 | 658 std::list<std::string> instance; |
659 context.GetIndex().LookupTagValue(instance, DICOM_TAG_SOP_INSTANCE_UID, uid); | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
660 |
764 | 661 if (instance.size() == 1) |
662 { | |
663 result.append(instance.front()); | |
664 } | |
665 } | |
666 } | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
667 |
764 | 668 call.GetOutput().AnswerJson(result); |
669 } | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
670 } |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
671 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
672 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
673 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
674 static void GetRtStructuresClosedPlanarsOfInstance(RestApi::GetCall& call) |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
675 { |
764 | 676 ServerContext& context = OrthancRestApi::GetContext(call); |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
677 |
764 | 678 Json::Value roi, contour, instance; |
679 std::string instanceId; | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
680 |
764 | 681 if (context.GetIndex().LookupResource(instance, call.GetUriComponent("instance", ""), ResourceType_Instance) && |
682 GetRtStructuresRoi(roi, contour, instanceId, context, | |
683 call.GetUriComponent("id", ""), | |
684 call.GetUriComponent("roi", ""))) | |
685 { | |
686 Json::Value result = Json::arrayValue; | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
687 |
764 | 688 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
689 { | |
690 if (contour[i].isMember(CONTOUR_DATA) && | |
691 contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
692 contour[i].isMember(NUMBER_OF_CONTOUR_POINTS) && | |
693 contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && | |
694 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
695 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID) && | |
696 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "CLOSED_PLANAR") | |
697 { | |
698 std::string uid = contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0][REFERENCED_SOP_INSTANCE_UID]["Value"].asString(); | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
699 |
764 | 700 Json::Value points; |
701 if (uid == instance["MainDicomTags"]["SOPInstanceUID"].asString() && | |
702 GetClosedPlanarPoints(points, context, instanceId, roi, i)) | |
703 { | |
704 result.append(points); | |
705 } | |
706 } | |
707 } | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
708 |
764 | 709 call.GetOutput().AnswerJson(result); |
710 } | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
711 } |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
712 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
713 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
714 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
715 static void GetRtStructuresPointsOfInstance(RestApi::GetCall& call) |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
716 { |
764 | 717 ServerContext& context = OrthancRestApi::GetContext(call); |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
718 |
764 | 719 Json::Value roi, contour, instance; |
720 std::string instanceId; | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
721 |
764 | 722 if (context.GetIndex().LookupResource(instance, call.GetUriComponent("instance", ""), ResourceType_Instance) && |
723 GetRtStructuresRoi(roi, contour, instanceId, context, | |
724 call.GetUriComponent("id", ""), | |
725 call.GetUriComponent("roi", ""))) | |
726 { | |
727 Json::Value result = Json::arrayValue; | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
728 |
764 | 729 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
730 { | |
731 if (contour[i].isMember(CONTOUR_DATA) && | |
732 contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
733 contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && | |
734 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
735 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID) && | |
736 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "POINT") | |
737 { | |
738 std::string uid = contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0][REFERENCED_SOP_INSTANCE_UID]["Value"].asString(); | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
739 |
764 | 740 if (uid == instance["MainDicomTags"]["SOPInstanceUID"].asString()) |
741 { | |
742 Json::Value p; | |
743 if (ContourToPoints(p, contour[i][CONTOUR_DATA]["Value"].asString()) && | |
744 p.size() == 1) | |
745 { | |
746 result.append(p[0]); | |
747 } | |
748 } | |
749 } | |
750 } | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
751 |
764 | 752 call.GetOutput().AnswerJson(result); |
753 } | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
754 } |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
755 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
756 |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
757 static void GetRtStructuresVolume(RestApi::GetCall& call) |
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
758 { |
764 | 759 ServerContext& context = OrthancRestApi::GetContext(call); |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
760 |
764 | 761 Json::Value roi, contour, vertices; |
762 std::string instanceId; | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
763 |
764 | 764 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
765 call.GetUriComponent("id", ""), | |
766 call.GetUriComponent("roi", ""))) | |
767 { | |
768 double volume = 0; | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
769 |
764 | 770 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
771 { | |
772 double area; | |
773 float thickness; | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
774 |
764 | 775 if (contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && |
776 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "CLOSED_PLANAR" && | |
777 GetClosedPlanarPoints(vertices, context, instanceId, roi, i) && | |
778 ComputeClosedPlanarArea(area, vertices) && | |
779 GetRtStructuresClosedPlanarThickness(thickness, context, contour[i])) | |
780 { | |
781 volume += area * static_cast<double>(thickness); | |
782 } | |
783 } | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
784 |
764 | 785 call.GetOutput().AnswerBuffer(boost::lexical_cast<std::string>(volume), "text/plain"); |
786 } | |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
787 } |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
788 |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
789 |
764 | 790 void OrthancRestApi::RegisterRadiotherapy() |
526 | 791 { |
792 Register("/series/{id}/rt-structures", GetRtStructuresInfo); | |
529 | 793 Register("/series/{id}/rt-structures/roi", GetRtStructuresListOfROIs); |
535 | 794 Register("/series/{id}/rt-structures/roi/{roi}/info", GetRtStructuresROI); |
795 Register("/series/{id}/rt-structures/roi/{roi}/points", GetRtStructuresROIPoints); | |
796 Register("/series/{id}/rt-structures/roi/{roi}/closed-planar", GetRtStructuresListOfClosedPlanars); | |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
797 Register("/series/{id}/rt-structures/roi/{roi}/closed-planar/{polygon}/vertices", GetRtStructuresSingleClosedPlanar); |
535 | 798 Register("/series/{id}/rt-structures/roi/{roi}/closed-planar/{polygon}/thickness", GetRtStructuresClosedPlanarThickness); |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
799 Register("/series/{id}/rt-structures/roi/{roi}/closed-planar/{polygon}/instance", GetRtStructuresInstanceOfClosedPlanar); |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
800 Register("/series/{id}/rt-structures/roi/{roi}/closed-planar/{polygon}/area", GetRtStructuresClosedPlanarArea); |
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
801 Register("/series/{id}/rt-structures/roi/{roi}/instances", GetRtStructuresListOfInstances); |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
802 Register("/series/{id}/rt-structures/roi/{roi}/instances/{instance}/closed-planar", GetRtStructuresClosedPlanarsOfInstance); |
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
803 Register("/series/{id}/rt-structures/roi/{roi}/instances/{instance}/points", GetRtStructuresPointsOfInstance); |
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
804 Register("/series/{id}/rt-structures/roi/{roi}/volume", GetRtStructuresVolume); |
526 | 805 } |
529 | 806 |
526 | 807 } |
808 | |
809 | |
767 | 810 /** |
811 storescu localhost 4242 ~/DICOM/Akos/data1/*.dcm | |
812 curl http://localhost:8042/series/10668f4a-fcc8fd1c-832e409c-5c7c018f-7ac8d3d9/rt-structures | |
813 **/ |