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 **/ |