Mercurial > hg > orthanc
annotate OrthancServer/OrthancRestApi/RadiotherapyRestApi.cpp @ 948:e57e08ed510f dicom-rt
integration mainline -> dicom-rt
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Wed, 25 Jun 2014 13:57:05 +0200 |
| parents | c19552f604d5 |
| children |
| 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 { |
|
948
e57e08ed510f
integration mainline -> dicom-rt
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
767
diff
changeset
|
281 std::string content; |
|
e57e08ed510f
integration mainline -> dicom-rt
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
767
diff
changeset
|
282 context.ReadFile(content, instanceId, FileContentType_Dicom); |
|
e57e08ed510f
integration mainline -> dicom-rt
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
767
diff
changeset
|
283 |
|
e57e08ed510f
integration mainline -> dicom-rt
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
767
diff
changeset
|
284 ParsedDicomFile dicom(content); |
|
541
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 { | |
| 764 | 365 ServerContext& context = OrthancRestApi::GetContext(call); |
| 526 | 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 { | |
| 764 | 403 ServerContext& context = OrthancRestApi::GetContext(call); |
| 529 | 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 { | |
| 764 | 429 ServerContext& context = OrthancRestApi::GetContext(call); |
| 535 | 430 |
| 764 | 431 Json::Value roi, contour; |
| 432 std::string instanceId; | |
| 535 | 433 |
| 764 | 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 } | |
| 535 | 441 } |
| 442 | |
| 443 | |
| 444 static void GetRtStructuresROIPoints(RestApi::GetCall& call) | |
| 445 { | |
| 764 | 446 ServerContext& context = OrthancRestApi::GetContext(call); |
| 535 | 447 |
| 764 | 448 Json::Value roi, contour; |
| 449 std::string instanceId; | |
| 535 | 450 |
| 764 | 451 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
| 452 call.GetUriComponent("id", ""), | |
| 453 call.GetUriComponent("roi", ""))) | |
| 454 { | |
| 455 Json::Value result = Json::arrayValue; | |
| 535 | 456 |
| 764 | 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 } | |
| 535 | 469 |
| 764 | 470 call.GetOutput().AnswerJson(result); |
| 471 } | |
| 535 | 472 } |
| 473 | |
| 474 | |
| 475 static void GetRtStructuresListOfClosedPlanars(RestApi::GetCall& call) | |
| 476 { | |
| 764 | 477 ServerContext& context = OrthancRestApi::GetContext(call); |
| 535 | 478 |
| 764 | 479 Json::Value roi, contour; |
| 480 std::string instanceId; | |
| 535 | 481 |
| 764 | 482 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
| 483 call.GetUriComponent("id", ""), | |
| 484 call.GetUriComponent("roi", ""))) | |
| 485 { | |
| 486 Json::Value result = Json::arrayValue; | |
| 535 | 487 |
| 764 | 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 } | |
| 535 | 501 |
| 764 | 502 call.GetOutput().AnswerJson(result); |
| 503 } | |
| 535 | 504 } |
| 505 | |
| 506 | |
| 507 static void GetRtStructuresSingleClosedPlanar(RestApi::GetCall& call) | |
| 508 { | |
| 764 | 509 ServerContext& context = OrthancRestApi::GetContext(call); |
| 535 | 510 |
| 764 | 511 Json::Value roi, contour, result; |
| 512 std::string instanceId; | |
| 535 | 513 |
| 764 | 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", "")); | |
| 535 | 519 |
| 764 | 520 if (GetClosedPlanarPoints(result, context, instanceId, roi, index)) |
| 521 { | |
| 522 call.GetOutput().AnswerJson(result); | |
| 523 } | |
| 524 } | |
| 535 | 525 } |
| 526 | |
| 527 | |
| 528 static void GetRtStructuresClosedPlanarThickness(RestApi::GetCall& call) | |
| 529 { | |
| 764 | 530 ServerContext& context = OrthancRestApi::GetContext(call); |
| 535 | 531 |
| 764 | 532 Json::Value roi, contour; |
| 533 std::string instanceId; | |
| 535 | 534 |
| 764 | 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", "")); | |
| 535 | 540 |
| 764 | 541 float thickness; |
| 542 if (GetRtStructuresClosedPlanarThickness(thickness, context, contour[index])) | |
| 543 { | |
| 544 call.GetOutput().AnswerBuffer(boost::lexical_cast<std::string>(thickness), "text/plain"); | |
| 545 } | |
| 546 } | |
|
541
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 { |
| 764 | 596 ServerContext& context = OrthancRestApi::GetContext(call); |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
597 |
| 764 | 598 Json::Value roi, contour, vertices; |
| 599 std::string instanceId; | |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
600 |
| 764 | 601 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
| 602 call.GetUriComponent("id", ""), | |
| 603 call.GetUriComponent("roi", ""))) | |
| 604 { | |
| 605 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
|
606 |
| 764 | 607 double area; |
| 608 if (GetClosedPlanarPoints(vertices, context, instanceId, roi, index) && | |
| 609 ComputeClosedPlanarArea(area, vertices)) | |
| 610 { | |
| 611 call.GetOutput().AnswerBuffer(boost::lexical_cast<std::string>(area), "text/plain"); | |
| 612 } | |
| 613 } | |
| 535 | 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 { |
| 764 | 619 ServerContext& context = OrthancRestApi::GetContext(call); |
| 529 | 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 { |
| 764 | 641 ServerContext& context = OrthancRestApi::GetContext(call); |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
642 |
| 764 | 643 Json::Value roi, contour; |
| 644 std::string instanceId; | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
645 |
| 764 | 646 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
| 647 call.GetUriComponent("id", ""), | |
| 648 call.GetUriComponent("roi", ""))) | |
| 649 { | |
| 650 Json::Value result = Json::arrayValue; | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
651 |
| 764 | 652 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
| 653 { | |
| 654 if (contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
| 655 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
| 656 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID)) | |
| 657 { | |
| 658 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
|
659 |
| 764 | 660 std::list<std::string> instance; |
| 661 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
|
662 |
| 764 | 663 if (instance.size() == 1) |
| 664 { | |
| 665 result.append(instance.front()); | |
| 666 } | |
| 667 } | |
| 668 } | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
669 |
| 764 | 670 call.GetOutput().AnswerJson(result); |
| 671 } | |
|
538
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 { |
| 764 | 678 ServerContext& context = OrthancRestApi::GetContext(call); |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
679 |
| 764 | 680 Json::Value roi, contour, instance; |
| 681 std::string instanceId; | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
682 |
| 764 | 683 if (context.GetIndex().LookupResource(instance, call.GetUriComponent("instance", ""), ResourceType_Instance) && |
| 684 GetRtStructuresRoi(roi, contour, instanceId, context, | |
| 685 call.GetUriComponent("id", ""), | |
| 686 call.GetUriComponent("roi", ""))) | |
| 687 { | |
| 688 Json::Value result = Json::arrayValue; | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
689 |
| 764 | 690 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
| 691 { | |
| 692 if (contour[i].isMember(CONTOUR_DATA) && | |
| 693 contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
| 694 contour[i].isMember(NUMBER_OF_CONTOUR_POINTS) && | |
| 695 contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && | |
| 696 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
| 697 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID) && | |
| 698 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "CLOSED_PLANAR") | |
| 699 { | |
| 700 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
|
701 |
| 764 | 702 Json::Value points; |
| 703 if (uid == instance["MainDicomTags"]["SOPInstanceUID"].asString() && | |
| 704 GetClosedPlanarPoints(points, context, instanceId, roi, i)) | |
| 705 { | |
| 706 result.append(points); | |
| 707 } | |
| 708 } | |
| 709 } | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
710 |
| 764 | 711 call.GetOutput().AnswerJson(result); |
| 712 } | |
|
538
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 { |
| 764 | 719 ServerContext& context = OrthancRestApi::GetContext(call); |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
720 |
| 764 | 721 Json::Value roi, contour, instance; |
| 722 std::string instanceId; | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
723 |
| 764 | 724 if (context.GetIndex().LookupResource(instance, call.GetUriComponent("instance", ""), ResourceType_Instance) && |
| 725 GetRtStructuresRoi(roi, contour, instanceId, context, | |
| 726 call.GetUriComponent("id", ""), | |
| 727 call.GetUriComponent("roi", ""))) | |
| 728 { | |
| 729 Json::Value result = Json::arrayValue; | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
730 |
| 764 | 731 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
| 732 { | |
| 733 if (contour[i].isMember(CONTOUR_DATA) && | |
| 734 contour[i].isMember(CONTOUR_IMAGE_SEQUENCE) && | |
| 735 contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && | |
| 736 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"].size() == 1 && | |
| 737 contour[i][CONTOUR_IMAGE_SEQUENCE]["Value"][0].isMember(REFERENCED_SOP_INSTANCE_UID) && | |
| 738 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "POINT") | |
| 739 { | |
| 740 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
|
741 |
| 764 | 742 if (uid == instance["MainDicomTags"]["SOPInstanceUID"].asString()) |
| 743 { | |
| 744 Json::Value p; | |
| 745 if (ContourToPoints(p, contour[i][CONTOUR_DATA]["Value"].asString()) && | |
| 746 p.size() == 1) | |
| 747 { | |
| 748 result.append(p[0]); | |
| 749 } | |
| 750 } | |
| 751 } | |
| 752 } | |
|
538
505d6deb9947
access contours by instances
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
535
diff
changeset
|
753 |
| 764 | 754 call.GetOutput().AnswerJson(result); |
| 755 } | |
|
538
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 { |
| 764 | 761 ServerContext& context = OrthancRestApi::GetContext(call); |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
762 |
| 764 | 763 Json::Value roi, contour, vertices; |
| 764 std::string instanceId; | |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
765 |
| 764 | 766 if (GetRtStructuresRoi(roi, contour, instanceId, context, |
| 767 call.GetUriComponent("id", ""), | |
| 768 call.GetUriComponent("roi", ""))) | |
| 769 { | |
| 770 double volume = 0; | |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
771 |
| 764 | 772 for (Json::Value::ArrayIndex i = 0; i < contour.size(); i++) |
| 773 { | |
| 774 double area; | |
| 775 float thickness; | |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
776 |
| 764 | 777 if (contour[i].isMember(CONTOUR_GEOMETRIC_TYPE) && |
| 778 contour[i][CONTOUR_GEOMETRIC_TYPE]["Value"].asString() == "CLOSED_PLANAR" && | |
| 779 GetClosedPlanarPoints(vertices, context, instanceId, roi, i) && | |
| 780 ComputeClosedPlanarArea(area, vertices) && | |
| 781 GetRtStructuresClosedPlanarThickness(thickness, context, contour[i])) | |
| 782 { | |
| 783 volume += area * static_cast<double>(thickness); | |
| 784 } | |
| 785 } | |
|
541
c8347eef225b
volume of structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
538
diff
changeset
|
786 |
| 764 | 787 call.GetOutput().AnswerBuffer(boost::lexical_cast<std::string>(volume), "text/plain"); |
| 788 } | |
|
541
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 |
| 764 | 792 void OrthancRestApi::RegisterRadiotherapy() |
| 526 | 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 | |
| 767 | 812 /** |
| 813 storescu localhost 4242 ~/DICOM/Akos/data1/*.dcm | |
| 814 curl http://localhost:8042/series/10668f4a-fcc8fd1c-832e409c-5c7c018f-7ac8d3d9/rt-structures | |
|
948
e57e08ed510f
integration mainline -> dicom-rt
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
767
diff
changeset
|
815 curl http://localhost:8042/series/10668f4a-fcc8fd1c-832e409c-5c7c018f-7ac8d3d9/rt-structures/roi/5/closed-planar/19/vertices |
| 767 | 816 **/ |