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