Mercurial > hg > orthanc-stone
annotate Resources/Graveyard/Toolbox/DicomDataset.cpp @ 40:7207a407bcd8
shared copyright with osimis
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 04 Jan 2017 16:37:42 +0100 |
parents | 9ee7e2f5f1a3 |
children | 28956ed68280 |
rev | line source |
---|---|
39 | 1 /** |
2 * Stone of Orthanc | |
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
4 * Department, University Hospital of Liege, Belgium | |
40
7207a407bcd8
shared copyright with osimis
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
39
diff
changeset
|
5 * Copyright (C) 2017 Osimis, Belgium |
39 | 6 * |
7 * This program is free software: you can redistribute it and/or | |
8 * modify it under the terms of the GNU General Public License as | |
9 * published by the Free Software Foundation, either version 3 of the | |
10 * License, or (at your option) any later version. | |
11 * | |
12 * In addition, as a special exception, the copyright holders of this | |
13 * program give permission to link the code of its release with the | |
14 * OpenSSL project's "OpenSSL" library (or with modified versions of it | |
15 * that use the same license as the "OpenSSL" library), and distribute | |
16 * the linked executables. You must obey the GNU General Public License | |
17 * in all respects for all of the code used other than "OpenSSL". If you | |
18 * modify file(s) with this exception, you may extend this exception to | |
19 * your version of the file(s), but you are not obligated to do so. If | |
20 * you do not wish to do so, delete this exception statement from your | |
21 * version. If you delete this exception statement from all source files | |
22 * in the program, then also delete it here. | |
23 * | |
24 * This program is distributed in the hope that it will be useful, but | |
25 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
26 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
27 * General Public License for more details. | |
28 * | |
29 * You should have received a copy of the GNU General Public License | |
30 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
31 **/ | |
32 | |
33 | |
34 #include "DicomDataset.h" | |
35 | |
36 #include "../../Resources/Orthanc/Core/OrthancException.h" | |
37 #include "../../Resources/Orthanc/Core/Logging.h" | |
38 #include "../../Resources/Orthanc/Core/Toolbox.h" | |
39 | |
40 #include <boost/lexical_cast.hpp> | |
41 #include <json/value.h> | |
42 #include <json/reader.h> | |
43 | |
44 namespace OrthancStone | |
45 { | |
46 static uint16_t GetCharValue(char c) | |
47 { | |
48 if (c >= '0' && c <= '9') | |
49 return c - '0'; | |
50 else if (c >= 'a' && c <= 'f') | |
51 return c - 'a' + 10; | |
52 else if (c >= 'A' && c <= 'F') | |
53 return c - 'A' + 10; | |
54 else | |
55 return 0; | |
56 } | |
57 | |
58 | |
59 static uint16_t GetHexadecimalValue(const char* c) | |
60 { | |
61 return ((GetCharValue(c[0]) << 12) + | |
62 (GetCharValue(c[1]) << 8) + | |
63 (GetCharValue(c[2]) << 4) + | |
64 GetCharValue(c[3])); | |
65 } | |
66 | |
67 | |
68 static DicomDataset::Tag ParseTag(const std::string& tag) | |
69 { | |
70 if (tag.size() == 9 && | |
71 isxdigit(tag[0]) && | |
72 isxdigit(tag[1]) && | |
73 isxdigit(tag[2]) && | |
74 isxdigit(tag[3]) && | |
75 (tag[4] == '-' || tag[4] == ',') && | |
76 isxdigit(tag[5]) && | |
77 isxdigit(tag[6]) && | |
78 isxdigit(tag[7]) && | |
79 isxdigit(tag[8])) | |
80 { | |
81 uint16_t group = GetHexadecimalValue(tag.c_str()); | |
82 uint16_t element = GetHexadecimalValue(tag.c_str() + 5); | |
83 return std::make_pair(group, element); | |
84 } | |
85 else if (tag.size() == 8 && | |
86 isxdigit(tag[0]) && | |
87 isxdigit(tag[1]) && | |
88 isxdigit(tag[2]) && | |
89 isxdigit(tag[3]) && | |
90 isxdigit(tag[4]) && | |
91 isxdigit(tag[5]) && | |
92 isxdigit(tag[6]) && | |
93 isxdigit(tag[7])) | |
94 { | |
95 uint16_t group = GetHexadecimalValue(tag.c_str()); | |
96 uint16_t element = GetHexadecimalValue(tag.c_str() + 4); | |
97 return std::make_pair(group, element); | |
98 } | |
99 else | |
100 { | |
101 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
102 } | |
103 } | |
104 | |
105 void DicomDataset::Parse(const std::string& content) | |
106 { | |
107 Json::Value json; | |
108 Json::Reader reader; | |
109 if (!reader.parse(content, json)) | |
110 { | |
111 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
112 } | |
113 | |
114 Parse(json); | |
115 } | |
116 | |
117 | |
118 void DicomDataset::Parse(const Json::Value& content) | |
119 { | |
120 if (content.type() != Json::objectValue) | |
121 { | |
122 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
123 } | |
124 | |
125 Json::Value::Members members = content.getMemberNames(); | |
126 for (size_t i = 0; i < members.size(); i++) | |
127 { | |
128 Tag tag = ParseTag(members[i]); | |
129 | |
130 const Json::Value& item = content[members[i]]; | |
131 | |
132 if (item.type() != Json::objectValue || | |
133 !item.isMember("Type") || | |
134 !item.isMember("Value") || | |
135 !item.isMember("Name") || | |
136 item["Type"].type() != Json::stringValue || | |
137 item["Name"].type() != Json::stringValue) | |
138 { | |
139 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
140 } | |
141 | |
142 if (item["Type"].asString() == "String") | |
143 { | |
144 if (item["Value"].type() == Json::stringValue) | |
145 { | |
146 values_[tag] = item["Value"].asString(); | |
147 } | |
148 else | |
149 { | |
150 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
151 } | |
152 } | |
153 } | |
154 } | |
155 | |
156 | |
157 DicomDataset::DicomDataset(OrthancPlugins::IOrthancConnection& orthanc, | |
158 const std::string& instanceId) | |
159 { | |
160 std::string content; | |
161 orthanc.RestApiGet(content, "/instances/" + instanceId + "/tags"); | |
162 | |
163 Parse(content); | |
164 } | |
165 | |
166 | |
167 std::string DicomDataset::GetStringValue(const Tag& tag) const | |
168 { | |
169 Values::const_iterator it = values_.find(tag); | |
170 | |
171 if (it == values_.end()) | |
172 { | |
173 LOG(ERROR) << "Trying to access a DICOM tag that is not set in a DICOM dataset"; | |
174 throw Orthanc::OrthancException(Orthanc::ErrorCode_InexistentItem); | |
175 } | |
176 else | |
177 { | |
178 return it->second; | |
179 } | |
180 } | |
181 | |
182 | |
183 std::string DicomDataset::GetStringValue(const Tag& tag, | |
184 const std::string& defaultValue) const | |
185 { | |
186 Values::const_iterator it = values_.find(tag); | |
187 | |
188 if (it == values_.end()) | |
189 { | |
190 return defaultValue; | |
191 } | |
192 else | |
193 { | |
194 return it->second; | |
195 } | |
196 } | |
197 | |
198 | |
199 float DicomDataset::GetFloatValue(const Tag& tag) const | |
200 { | |
201 try | |
202 { | |
203 return boost::lexical_cast<float>(Orthanc::Toolbox::StripSpaces(GetStringValue(tag))); | |
204 } | |
205 catch (boost::bad_lexical_cast&) | |
206 { | |
207 LOG(ERROR) << "Trying to access a DICOM tag that is not a float"; | |
208 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
209 } | |
210 } | |
211 | |
212 | |
213 double DicomDataset::GetDoubleValue(const Tag& tag) const | |
214 { | |
215 try | |
216 { | |
217 return boost::lexical_cast<double>(Orthanc::Toolbox::StripSpaces(GetStringValue(tag))); | |
218 } | |
219 catch (boost::bad_lexical_cast&) | |
220 { | |
221 LOG(ERROR) << "Trying to access a DICOM tag that is not a float"; | |
222 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
223 } | |
224 } | |
225 | |
226 | |
227 int DicomDataset::GetIntegerValue(const Tag& tag) const | |
228 { | |
229 try | |
230 { | |
231 return boost::lexical_cast<int>(Orthanc::Toolbox::StripSpaces(GetStringValue(tag))); | |
232 } | |
233 catch (boost::bad_lexical_cast&) | |
234 { | |
235 LOG(ERROR) << "Trying to access a DICOM tag that is not an integer"; | |
236 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
237 } | |
238 } | |
239 | |
240 | |
241 unsigned int DicomDataset::GetUnsignedIntegerValue(const Tag& tag) const | |
242 { | |
243 int v = GetIntegerValue(tag); | |
244 | |
245 if (v < 0) | |
246 { | |
247 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
248 } | |
249 else | |
250 { | |
251 return static_cast<unsigned int>(v); | |
252 } | |
253 } | |
254 | |
255 | |
256 void DicomDataset::GetVectorValue(Vector& vector, | |
257 const Tag& tag) const | |
258 { | |
259 if (!GeometryToolbox::ParseVector(vector, Orthanc::Toolbox::StripSpaces(GetStringValue(tag)))) | |
260 { | |
261 LOG(ERROR) << "Trying to access a DICOM tag that is not a vector"; | |
262 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
263 } | |
264 } | |
265 | |
266 | |
267 void DicomDataset::GetVectorValue(Vector& vector, | |
268 const Tag& tag, | |
269 size_t expectedSize) const | |
270 { | |
271 GetVectorValue(vector, tag); | |
272 | |
273 if (vector.size() != expectedSize) | |
274 { | |
275 LOG(ERROR) << "A vector in a DICOM tag has a bad size"; | |
276 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
277 } | |
278 } | |
279 | |
280 | |
281 void DicomDataset::Print() const | |
282 { | |
283 for (Values::const_iterator it = values_.begin(); it != values_.end(); ++it) | |
284 { | |
285 printf("%04x,%04x = [%s]\n", it->first.first, it->first.second, it->second.c_str()); | |
286 } | |
287 printf("\n"); | |
288 } | |
289 | |
290 | |
291 bool DicomDataset::IsGrayscale() const | |
292 { | |
293 std::string photometric = Orthanc::Toolbox::StripSpaces(GetStringValue(DICOM_TAG_PHOTOMETRIC_INTERPRETATION)); | |
294 | |
295 return (photometric == "MONOCHROME1" || | |
296 photometric == "MONOCHROME2"); | |
297 } | |
298 | |
299 | |
300 void DicomDataset::GetPixelSpacing(double& spacingX, | |
301 double& spacingY) const | |
302 { | |
303 if (HasTag(DICOM_TAG_PIXEL_SPACING)) | |
304 { | |
305 Vector spacing; | |
306 GetVectorValue(spacing, DICOM_TAG_PIXEL_SPACING, 2); | |
307 spacingX = spacing[0]; | |
308 spacingY = spacing[1]; | |
309 } | |
310 else | |
311 { | |
312 spacingX = 1.0; | |
313 spacingY = 1.0; | |
314 } | |
315 } | |
316 } |