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