|
1
|
1 /**
|
|
|
2 * Orthanc - A Lightweight, RESTful DICOM Store
|
|
|
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 "../PrecompiledHeaders.h"
|
|
|
34 #include "DicomMap.h"
|
|
|
35
|
|
|
36 #include <stdio.h>
|
|
|
37 #include <memory>
|
|
|
38 #include "../Endianness.h"
|
|
|
39 #include "../OrthancException.h"
|
|
|
40
|
|
|
41
|
|
|
42 namespace Orthanc
|
|
|
43 {
|
|
|
44 static DicomTag patientTags[] =
|
|
|
45 {
|
|
|
46 //DicomTag(0x0010, 0x1010), // PatientAge
|
|
|
47 //DicomTag(0x0010, 0x1040) // PatientAddress
|
|
|
48 DicomTag(0x0010, 0x0010), // PatientName
|
|
|
49 DicomTag(0x0010, 0x0030), // PatientBirthDate
|
|
|
50 DicomTag(0x0010, 0x0040), // PatientSex
|
|
|
51 DicomTag(0x0010, 0x1000), // OtherPatientIDs
|
|
|
52 DICOM_TAG_PATIENT_ID
|
|
|
53 };
|
|
|
54
|
|
|
55 static DicomTag studyTags[] =
|
|
|
56 {
|
|
|
57 //DicomTag(0x0010, 0x1020), // PatientSize
|
|
|
58 //DicomTag(0x0010, 0x1030) // PatientWeight
|
|
|
59 DICOM_TAG_STUDY_DATE,
|
|
|
60 DicomTag(0x0008, 0x0030), // StudyTime
|
|
|
61 DicomTag(0x0020, 0x0010), // StudyID
|
|
|
62 DICOM_TAG_STUDY_DESCRIPTION,
|
|
|
63 DICOM_TAG_ACCESSION_NUMBER,
|
|
|
64 DICOM_TAG_STUDY_INSTANCE_UID,
|
|
|
65 DICOM_TAG_REQUESTED_PROCEDURE_DESCRIPTION, // New in db v6
|
|
|
66 DICOM_TAG_INSTITUTION_NAME, // New in db v6
|
|
|
67 DICOM_TAG_REQUESTING_PHYSICIAN, // New in db v6
|
|
|
68 DICOM_TAG_REFERRING_PHYSICIAN_NAME // New in db v6
|
|
|
69 };
|
|
|
70
|
|
|
71 static DicomTag seriesTags[] =
|
|
|
72 {
|
|
|
73 //DicomTag(0x0010, 0x1080), // MilitaryRank
|
|
|
74 DicomTag(0x0008, 0x0021), // SeriesDate
|
|
|
75 DicomTag(0x0008, 0x0031), // SeriesTime
|
|
|
76 DICOM_TAG_MODALITY,
|
|
|
77 DicomTag(0x0008, 0x0070), // Manufacturer
|
|
|
78 DicomTag(0x0008, 0x1010), // StationName
|
|
|
79 DICOM_TAG_SERIES_DESCRIPTION,
|
|
|
80 DicomTag(0x0018, 0x0015), // BodyPartExamined
|
|
|
81 DicomTag(0x0018, 0x0024), // SequenceName
|
|
|
82 DicomTag(0x0018, 0x1030), // ProtocolName
|
|
|
83 DicomTag(0x0020, 0x0011), // SeriesNumber
|
|
|
84 DICOM_TAG_CARDIAC_NUMBER_OF_IMAGES,
|
|
|
85 DICOM_TAG_IMAGES_IN_ACQUISITION,
|
|
|
86 DICOM_TAG_NUMBER_OF_TEMPORAL_POSITIONS,
|
|
|
87 DICOM_TAG_NUMBER_OF_SLICES,
|
|
|
88 DICOM_TAG_NUMBER_OF_TIME_SLICES,
|
|
|
89 DICOM_TAG_SERIES_INSTANCE_UID,
|
|
|
90 DICOM_TAG_IMAGE_ORIENTATION_PATIENT, // New in db v6
|
|
|
91 DICOM_TAG_SERIES_TYPE, // New in db v6
|
|
|
92 DICOM_TAG_OPERATOR_NAME, // New in db v6
|
|
|
93 DICOM_TAG_PERFORMED_PROCEDURE_STEP_DESCRIPTION, // New in db v6
|
|
|
94 DICOM_TAG_ACQUISITION_DEVICE_PROCESSING_DESCRIPTION, // New in db v6
|
|
|
95 DICOM_TAG_CONTRAST_BOLUS_AGENT // New in db v6
|
|
|
96 };
|
|
|
97
|
|
|
98 static DicomTag instanceTags[] =
|
|
|
99 {
|
|
|
100 DicomTag(0x0008, 0x0012), // InstanceCreationDate
|
|
|
101 DicomTag(0x0008, 0x0013), // InstanceCreationTime
|
|
|
102 DicomTag(0x0020, 0x0012), // AcquisitionNumber
|
|
|
103 DICOM_TAG_IMAGE_INDEX,
|
|
|
104 DICOM_TAG_INSTANCE_NUMBER,
|
|
|
105 DICOM_TAG_NUMBER_OF_FRAMES,
|
|
|
106 DICOM_TAG_TEMPORAL_POSITION_IDENTIFIER,
|
|
|
107 DICOM_TAG_SOP_INSTANCE_UID,
|
|
|
108 DICOM_TAG_IMAGE_POSITION_PATIENT, // New in db v6
|
|
|
109 DICOM_TAG_IMAGE_COMMENTS // New in db v6
|
|
|
110 };
|
|
|
111
|
|
|
112
|
|
|
113 void DicomMap::LoadMainDicomTags(const DicomTag*& tags,
|
|
|
114 size_t& size,
|
|
|
115 ResourceType level)
|
|
|
116 {
|
|
|
117 switch (level)
|
|
|
118 {
|
|
|
119 case ResourceType_Patient:
|
|
|
120 tags = patientTags;
|
|
|
121 size = sizeof(patientTags) / sizeof(DicomTag);
|
|
|
122 break;
|
|
|
123
|
|
|
124 case ResourceType_Study:
|
|
|
125 tags = studyTags;
|
|
|
126 size = sizeof(studyTags) / sizeof(DicomTag);
|
|
|
127 break;
|
|
|
128
|
|
|
129 case ResourceType_Series:
|
|
|
130 tags = seriesTags;
|
|
|
131 size = sizeof(seriesTags) / sizeof(DicomTag);
|
|
|
132 break;
|
|
|
133
|
|
|
134 case ResourceType_Instance:
|
|
|
135 tags = instanceTags;
|
|
|
136 size = sizeof(instanceTags) / sizeof(DicomTag);
|
|
|
137 break;
|
|
|
138
|
|
|
139 default:
|
|
|
140 throw OrthancException(ErrorCode_ParameterOutOfRange);
|
|
|
141 }
|
|
|
142 }
|
|
|
143
|
|
|
144
|
|
|
145 void DicomMap::SetValue(uint16_t group,
|
|
|
146 uint16_t element,
|
|
|
147 DicomValue* value)
|
|
|
148 {
|
|
|
149 DicomTag tag(group, element);
|
|
|
150 Map::iterator it = map_.find(tag);
|
|
|
151
|
|
|
152 if (it != map_.end())
|
|
|
153 {
|
|
|
154 delete it->second;
|
|
|
155 it->second = value;
|
|
|
156 }
|
|
|
157 else
|
|
|
158 {
|
|
|
159 map_.insert(std::make_pair(tag, value));
|
|
|
160 }
|
|
|
161 }
|
|
|
162
|
|
|
163 void DicomMap::SetValue(DicomTag tag,
|
|
|
164 DicomValue* value)
|
|
|
165 {
|
|
|
166 SetValue(tag.GetGroup(), tag.GetElement(), value);
|
|
|
167 }
|
|
|
168
|
|
|
169
|
|
|
170
|
|
|
171
|
|
|
172 void DicomMap::Clear()
|
|
|
173 {
|
|
|
174 for (Map::iterator it = map_.begin(); it != map_.end(); ++it)
|
|
|
175 {
|
|
|
176 delete it->second;
|
|
|
177 }
|
|
|
178
|
|
|
179 map_.clear();
|
|
|
180 }
|
|
|
181
|
|
|
182
|
|
|
183 void DicomMap::ExtractTags(DicomMap& result,
|
|
|
184 const DicomTag* tags,
|
|
|
185 size_t count) const
|
|
|
186 {
|
|
|
187 result.Clear();
|
|
|
188
|
|
|
189 for (unsigned int i = 0; i < count; i++)
|
|
|
190 {
|
|
|
191 Map::const_iterator it = map_.find(tags[i]);
|
|
|
192 if (it != map_.end())
|
|
|
193 {
|
|
|
194 result.SetValue(it->first, it->second->Clone());
|
|
|
195 }
|
|
|
196 }
|
|
|
197 }
|
|
|
198
|
|
|
199
|
|
|
200 void DicomMap::ExtractPatientInformation(DicomMap& result) const
|
|
|
201 {
|
|
|
202 ExtractTags(result, patientTags, sizeof(patientTags) / sizeof(DicomTag));
|
|
|
203 }
|
|
|
204
|
|
|
205 void DicomMap::ExtractStudyInformation(DicomMap& result) const
|
|
|
206 {
|
|
|
207 ExtractTags(result, studyTags, sizeof(studyTags) / sizeof(DicomTag));
|
|
|
208 }
|
|
|
209
|
|
|
210 void DicomMap::ExtractSeriesInformation(DicomMap& result) const
|
|
|
211 {
|
|
|
212 ExtractTags(result, seriesTags, sizeof(seriesTags) / sizeof(DicomTag));
|
|
|
213 }
|
|
|
214
|
|
|
215 void DicomMap::ExtractInstanceInformation(DicomMap& result) const
|
|
|
216 {
|
|
|
217 ExtractTags(result, instanceTags, sizeof(instanceTags) / sizeof(DicomTag));
|
|
|
218 }
|
|
|
219
|
|
|
220
|
|
|
221
|
|
|
222 DicomMap* DicomMap::Clone() const
|
|
|
223 {
|
|
|
224 std::auto_ptr<DicomMap> result(new DicomMap);
|
|
|
225
|
|
|
226 for (Map::const_iterator it = map_.begin(); it != map_.end(); ++it)
|
|
|
227 {
|
|
|
228 result->map_.insert(std::make_pair(it->first, it->second->Clone()));
|
|
|
229 }
|
|
|
230
|
|
|
231 return result.release();
|
|
|
232 }
|
|
|
233
|
|
|
234
|
|
|
235 void DicomMap::Assign(const DicomMap& other)
|
|
|
236 {
|
|
|
237 Clear();
|
|
|
238
|
|
|
239 for (Map::const_iterator it = other.map_.begin(); it != other.map_.end(); ++it)
|
|
|
240 {
|
|
|
241 map_.insert(std::make_pair(it->first, it->second->Clone()));
|
|
|
242 }
|
|
|
243 }
|
|
|
244
|
|
|
245
|
|
|
246 const DicomValue& DicomMap::GetValue(const DicomTag& tag) const
|
|
|
247 {
|
|
|
248 const DicomValue* value = TestAndGetValue(tag);
|
|
|
249
|
|
|
250 if (value)
|
|
|
251 {
|
|
|
252 return *value;
|
|
|
253 }
|
|
|
254 else
|
|
|
255 {
|
|
|
256 throw OrthancException(ErrorCode_InexistentTag);
|
|
|
257 }
|
|
|
258 }
|
|
|
259
|
|
|
260
|
|
|
261 const DicomValue* DicomMap::TestAndGetValue(const DicomTag& tag) const
|
|
|
262 {
|
|
|
263 Map::const_iterator it = map_.find(tag);
|
|
|
264
|
|
|
265 if (it == map_.end())
|
|
|
266 {
|
|
|
267 return NULL;
|
|
|
268 }
|
|
|
269 else
|
|
|
270 {
|
|
|
271 return it->second;
|
|
|
272 }
|
|
|
273 }
|
|
|
274
|
|
|
275
|
|
|
276 void DicomMap::Remove(const DicomTag& tag)
|
|
|
277 {
|
|
|
278 Map::iterator it = map_.find(tag);
|
|
|
279 if (it != map_.end())
|
|
|
280 {
|
|
|
281 delete it->second;
|
|
|
282 map_.erase(it);
|
|
|
283 }
|
|
|
284 }
|
|
|
285
|
|
|
286
|
|
|
287 static void SetupFindTemplate(DicomMap& result,
|
|
|
288 const DicomTag* tags,
|
|
|
289 size_t count)
|
|
|
290 {
|
|
|
291 result.Clear();
|
|
|
292
|
|
|
293 for (size_t i = 0; i < count; i++)
|
|
|
294 {
|
|
|
295 result.SetValue(tags[i], "", false);
|
|
|
296 }
|
|
|
297 }
|
|
|
298
|
|
|
299 void DicomMap::SetupFindPatientTemplate(DicomMap& result)
|
|
|
300 {
|
|
|
301 SetupFindTemplate(result, patientTags, sizeof(patientTags) / sizeof(DicomTag));
|
|
|
302 }
|
|
|
303
|
|
|
304 void DicomMap::SetupFindStudyTemplate(DicomMap& result)
|
|
|
305 {
|
|
|
306 SetupFindTemplate(result, studyTags, sizeof(studyTags) / sizeof(DicomTag));
|
|
|
307 result.SetValue(DICOM_TAG_ACCESSION_NUMBER, "", false);
|
|
|
308 result.SetValue(DICOM_TAG_PATIENT_ID, "", false);
|
|
|
309
|
|
|
310 // These main DICOM tags are only indirectly related to the
|
|
|
311 // General Study Module, remove them
|
|
|
312 result.Remove(DICOM_TAG_INSTITUTION_NAME);
|
|
|
313 result.Remove(DICOM_TAG_REQUESTING_PHYSICIAN);
|
|
|
314 result.Remove(DICOM_TAG_REQUESTED_PROCEDURE_DESCRIPTION);
|
|
|
315 }
|
|
|
316
|
|
|
317 void DicomMap::SetupFindSeriesTemplate(DicomMap& result)
|
|
|
318 {
|
|
|
319 SetupFindTemplate(result, seriesTags, sizeof(seriesTags) / sizeof(DicomTag));
|
|
|
320 result.SetValue(DICOM_TAG_ACCESSION_NUMBER, "", false);
|
|
|
321 result.SetValue(DICOM_TAG_PATIENT_ID, "", false);
|
|
|
322 result.SetValue(DICOM_TAG_STUDY_INSTANCE_UID, "", false);
|
|
|
323
|
|
|
324 // These tags are considered as "main" by Orthanc, but are not in the Series module
|
|
|
325 result.Remove(DicomTag(0x0008, 0x0070)); // Manufacturer
|
|
|
326 result.Remove(DicomTag(0x0008, 0x1010)); // Station name
|
|
|
327 result.Remove(DicomTag(0x0018, 0x0024)); // Sequence name
|
|
|
328 result.Remove(DICOM_TAG_CARDIAC_NUMBER_OF_IMAGES);
|
|
|
329 result.Remove(DICOM_TAG_IMAGES_IN_ACQUISITION);
|
|
|
330 result.Remove(DICOM_TAG_NUMBER_OF_SLICES);
|
|
|
331 result.Remove(DICOM_TAG_NUMBER_OF_TEMPORAL_POSITIONS);
|
|
|
332 result.Remove(DICOM_TAG_NUMBER_OF_TIME_SLICES);
|
|
|
333 result.Remove(DICOM_TAG_IMAGE_ORIENTATION_PATIENT);
|
|
|
334 result.Remove(DICOM_TAG_SERIES_TYPE);
|
|
|
335 result.Remove(DICOM_TAG_ACQUISITION_DEVICE_PROCESSING_DESCRIPTION);
|
|
|
336 result.Remove(DICOM_TAG_CONTRAST_BOLUS_AGENT);
|
|
|
337 }
|
|
|
338
|
|
|
339 void DicomMap::SetupFindInstanceTemplate(DicomMap& result)
|
|
|
340 {
|
|
|
341 SetupFindTemplate(result, instanceTags, sizeof(instanceTags) / sizeof(DicomTag));
|
|
|
342 result.SetValue(DICOM_TAG_ACCESSION_NUMBER, "", false);
|
|
|
343 result.SetValue(DICOM_TAG_PATIENT_ID, "", false);
|
|
|
344 result.SetValue(DICOM_TAG_STUDY_INSTANCE_UID, "", false);
|
|
|
345 result.SetValue(DICOM_TAG_SERIES_INSTANCE_UID, "", false);
|
|
|
346 }
|
|
|
347
|
|
|
348
|
|
|
349 void DicomMap::CopyTagIfExists(const DicomMap& source,
|
|
|
350 const DicomTag& tag)
|
|
|
351 {
|
|
|
352 if (source.HasTag(tag))
|
|
|
353 {
|
|
|
354 SetValue(tag, source.GetValue(tag));
|
|
|
355 }
|
|
|
356 }
|
|
|
357
|
|
|
358
|
|
|
359 bool DicomMap::IsMainDicomTag(const DicomTag& tag, ResourceType level)
|
|
|
360 {
|
|
|
361 DicomTag *tags = NULL;
|
|
|
362 size_t size;
|
|
|
363
|
|
|
364 switch (level)
|
|
|
365 {
|
|
|
366 case ResourceType_Patient:
|
|
|
367 tags = patientTags;
|
|
|
368 size = sizeof(patientTags) / sizeof(DicomTag);
|
|
|
369 break;
|
|
|
370
|
|
|
371 case ResourceType_Study:
|
|
|
372 tags = studyTags;
|
|
|
373 size = sizeof(studyTags) / sizeof(DicomTag);
|
|
|
374 break;
|
|
|
375
|
|
|
376 case ResourceType_Series:
|
|
|
377 tags = seriesTags;
|
|
|
378 size = sizeof(seriesTags) / sizeof(DicomTag);
|
|
|
379 break;
|
|
|
380
|
|
|
381 case ResourceType_Instance:
|
|
|
382 tags = instanceTags;
|
|
|
383 size = sizeof(instanceTags) / sizeof(DicomTag);
|
|
|
384 break;
|
|
|
385
|
|
|
386 default:
|
|
|
387 throw OrthancException(ErrorCode_ParameterOutOfRange);
|
|
|
388 }
|
|
|
389
|
|
|
390 for (size_t i = 0; i < size; i++)
|
|
|
391 {
|
|
|
392 if (tags[i] == tag)
|
|
|
393 {
|
|
|
394 return true;
|
|
|
395 }
|
|
|
396 }
|
|
|
397
|
|
|
398 return false;
|
|
|
399 }
|
|
|
400
|
|
|
401 bool DicomMap::IsMainDicomTag(const DicomTag& tag)
|
|
|
402 {
|
|
|
403 return (IsMainDicomTag(tag, ResourceType_Patient) ||
|
|
|
404 IsMainDicomTag(tag, ResourceType_Study) ||
|
|
|
405 IsMainDicomTag(tag, ResourceType_Series) ||
|
|
|
406 IsMainDicomTag(tag, ResourceType_Instance));
|
|
|
407 }
|
|
|
408
|
|
|
409
|
|
|
410 void DicomMap::GetMainDicomTagsInternal(std::set<DicomTag>& result, ResourceType level)
|
|
|
411 {
|
|
|
412 DicomTag *tags = NULL;
|
|
|
413 size_t size;
|
|
|
414
|
|
|
415 switch (level)
|
|
|
416 {
|
|
|
417 case ResourceType_Patient:
|
|
|
418 tags = patientTags;
|
|
|
419 size = sizeof(patientTags) / sizeof(DicomTag);
|
|
|
420 break;
|
|
|
421
|
|
|
422 case ResourceType_Study:
|
|
|
423 tags = studyTags;
|
|
|
424 size = sizeof(studyTags) / sizeof(DicomTag);
|
|
|
425 break;
|
|
|
426
|
|
|
427 case ResourceType_Series:
|
|
|
428 tags = seriesTags;
|
|
|
429 size = sizeof(seriesTags) / sizeof(DicomTag);
|
|
|
430 break;
|
|
|
431
|
|
|
432 case ResourceType_Instance:
|
|
|
433 tags = instanceTags;
|
|
|
434 size = sizeof(instanceTags) / sizeof(DicomTag);
|
|
|
435 break;
|
|
|
436
|
|
|
437 default:
|
|
|
438 throw OrthancException(ErrorCode_ParameterOutOfRange);
|
|
|
439 }
|
|
|
440
|
|
|
441 for (size_t i = 0; i < size; i++)
|
|
|
442 {
|
|
|
443 result.insert(tags[i]);
|
|
|
444 }
|
|
|
445 }
|
|
|
446
|
|
|
447
|
|
|
448 void DicomMap::GetMainDicomTags(std::set<DicomTag>& result, ResourceType level)
|
|
|
449 {
|
|
|
450 result.clear();
|
|
|
451 GetMainDicomTagsInternal(result, level);
|
|
|
452 }
|
|
|
453
|
|
|
454
|
|
|
455 void DicomMap::GetMainDicomTags(std::set<DicomTag>& result)
|
|
|
456 {
|
|
|
457 result.clear();
|
|
|
458 GetMainDicomTagsInternal(result, ResourceType_Patient);
|
|
|
459 GetMainDicomTagsInternal(result, ResourceType_Study);
|
|
|
460 GetMainDicomTagsInternal(result, ResourceType_Series);
|
|
|
461 GetMainDicomTagsInternal(result, ResourceType_Instance);
|
|
|
462 }
|
|
|
463
|
|
|
464
|
|
|
465 void DicomMap::GetTags(std::set<DicomTag>& tags) const
|
|
|
466 {
|
|
|
467 tags.clear();
|
|
|
468
|
|
|
469 for (Map::const_iterator it = map_.begin();
|
|
|
470 it != map_.end(); ++it)
|
|
|
471 {
|
|
|
472 tags.insert(it->first);
|
|
|
473 }
|
|
|
474 }
|
|
|
475
|
|
|
476
|
|
|
477 static uint16_t ReadUnsignedInteger16(const char* dicom)
|
|
|
478 {
|
|
|
479 return le16toh(*reinterpret_cast<const uint16_t*>(dicom));
|
|
|
480 }
|
|
|
481
|
|
|
482
|
|
|
483 static uint32_t ReadUnsignedInteger32(const char* dicom)
|
|
|
484 {
|
|
|
485 return le32toh(*reinterpret_cast<const uint32_t*>(dicom));
|
|
|
486 }
|
|
|
487
|
|
|
488
|
|
|
489 static bool ValidateTag(const ValueRepresentation& vr,
|
|
|
490 const std::string& value)
|
|
|
491 {
|
|
|
492 switch (vr)
|
|
|
493 {
|
|
|
494 case ValueRepresentation_ApplicationEntity:
|
|
|
495 return value.size() <= 16;
|
|
|
496
|
|
|
497 case ValueRepresentation_AgeString:
|
|
|
498 return (value.size() == 4 &&
|
|
|
499 isdigit(value[0]) &&
|
|
|
500 isdigit(value[1]) &&
|
|
|
501 isdigit(value[2]) &&
|
|
|
502 (value[3] == 'D' || value[3] == 'W' || value[3] == 'M' || value[3] == 'Y'));
|
|
|
503
|
|
|
504 case ValueRepresentation_AttributeTag:
|
|
|
505 return value.size() == 4;
|
|
|
506
|
|
|
507 case ValueRepresentation_CodeString:
|
|
|
508 return value.size() <= 16;
|
|
|
509
|
|
|
510 case ValueRepresentation_Date:
|
|
|
511 return value.size() <= 18;
|
|
|
512
|
|
|
513 case ValueRepresentation_DecimalString:
|
|
|
514 return value.size() <= 16;
|
|
|
515
|
|
|
516 case ValueRepresentation_DateTime:
|
|
|
517 return value.size() <= 54;
|
|
|
518
|
|
|
519 case ValueRepresentation_FloatingPointSingle:
|
|
|
520 return value.size() == 4;
|
|
|
521
|
|
|
522 case ValueRepresentation_FloatingPointDouble:
|
|
|
523 return value.size() == 8;
|
|
|
524
|
|
|
525 case ValueRepresentation_IntegerString:
|
|
|
526 return value.size() <= 12;
|
|
|
527
|
|
|
528 case ValueRepresentation_LongString:
|
|
|
529 return value.size() <= 64;
|
|
|
530
|
|
|
531 case ValueRepresentation_LongText:
|
|
|
532 return value.size() <= 10240;
|
|
|
533
|
|
|
534 case ValueRepresentation_OtherByte:
|
|
|
535 return true;
|
|
|
536
|
|
|
537 case ValueRepresentation_OtherDouble:
|
|
|
538 return value.size() <= (static_cast<uint64_t>(1) << 32) - 8;
|
|
|
539
|
|
|
540 case ValueRepresentation_OtherFloat:
|
|
|
541 return value.size() <= (static_cast<uint64_t>(1) << 32) - 4;
|
|
|
542
|
|
|
543 case ValueRepresentation_OtherLong:
|
|
|
544 return true;
|
|
|
545
|
|
|
546 case ValueRepresentation_OtherWord:
|
|
|
547 return true;
|
|
|
548
|
|
|
549 case ValueRepresentation_PersonName:
|
|
|
550 return true;
|
|
|
551
|
|
|
552 case ValueRepresentation_ShortString:
|
|
|
553 return value.size() <= 16;
|
|
|
554
|
|
|
555 case ValueRepresentation_SignedLong:
|
|
|
556 return value.size() == 4;
|
|
|
557
|
|
|
558 case ValueRepresentation_Sequence:
|
|
|
559 return true;
|
|
|
560
|
|
|
561 case ValueRepresentation_SignedShort:
|
|
|
562 return value.size() == 2;
|
|
|
563
|
|
|
564 case ValueRepresentation_ShortText:
|
|
|
565 return value.size() <= 1024;
|
|
|
566
|
|
|
567 case ValueRepresentation_Time:
|
|
|
568 return value.size() <= 28;
|
|
|
569
|
|
|
570 case ValueRepresentation_UnlimitedCharacters:
|
|
|
571 return value.size() <= (static_cast<uint64_t>(1) << 32) - 2;
|
|
|
572
|
|
|
573 case ValueRepresentation_UniqueIdentifier:
|
|
|
574 return value.size() <= 64;
|
|
|
575
|
|
|
576 case ValueRepresentation_UnsignedLong:
|
|
|
577 return value.size() == 4;
|
|
|
578
|
|
|
579 case ValueRepresentation_Unknown:
|
|
|
580 return true;
|
|
|
581
|
|
|
582 case ValueRepresentation_UniversalResource:
|
|
|
583 return value.size() <= (static_cast<uint64_t>(1) << 32) - 2;
|
|
|
584
|
|
|
585 case ValueRepresentation_UnsignedShort:
|
|
|
586 return value.size() == 2;
|
|
|
587
|
|
|
588 case ValueRepresentation_UnlimitedText:
|
|
|
589 return value.size() <= (static_cast<uint64_t>(1) << 32) - 2;
|
|
|
590
|
|
|
591 default:
|
|
|
592 // Assume unsupported tags are OK
|
|
|
593 return true;
|
|
|
594 }
|
|
|
595 }
|
|
|
596
|
|
|
597
|
|
|
598 static void RemoveTagPadding(std::string& value,
|
|
|
599 const ValueRepresentation& vr)
|
|
|
600 {
|
|
|
601 /**
|
|
|
602 * Remove padding from character strings, if need be. For the time
|
|
|
603 * being, only the UI VR is supported.
|
|
|
604 * http://dicom.nema.org/medical/dicom/current/output/chtml/part05/sect_6.2.html
|
|
|
605 **/
|
|
|
606
|
|
|
607 switch (vr)
|
|
|
608 {
|
|
|
609 case ValueRepresentation_UniqueIdentifier:
|
|
|
610 {
|
|
|
611 /**
|
|
|
612 * "Values with a VR of UI shall be padded with a single
|
|
|
613 * trailing NULL (00H) character when necessary to achieve even
|
|
|
614 * length."
|
|
|
615 **/
|
|
|
616
|
|
|
617 if (!value.empty() &&
|
|
|
618 value[value.size() - 1] == '\0')
|
|
|
619 {
|
|
|
620 value.resize(value.size() - 1);
|
|
|
621 }
|
|
|
622
|
|
|
623 break;
|
|
|
624 }
|
|
|
625
|
|
|
626 /**
|
|
|
627 * TODO implement other VR
|
|
|
628 **/
|
|
|
629
|
|
|
630 default:
|
|
|
631 // No padding is applicable to this VR
|
|
|
632 break;
|
|
|
633 }
|
|
|
634 }
|
|
|
635
|
|
|
636
|
|
|
637 static bool ReadNextTag(DicomTag& tag,
|
|
|
638 ValueRepresentation& vr,
|
|
|
639 std::string& value,
|
|
|
640 const char* dicom,
|
|
|
641 size_t size,
|
|
|
642 size_t& position)
|
|
|
643 {
|
|
|
644 /**
|
|
|
645 * http://dicom.nema.org/medical/dicom/current/output/chtml/part05/chapter_7.html#sect_7.1.2
|
|
|
646 * This function reads a data element with Explicit VR encoded using Little-Endian.
|
|
|
647 **/
|
|
|
648
|
|
|
649 if (position + 6 > size)
|
|
|
650 {
|
|
|
651 return false;
|
|
|
652 }
|
|
|
653
|
|
|
654 tag = DicomTag(ReadUnsignedInteger16(dicom + position),
|
|
|
655 ReadUnsignedInteger16(dicom + position + 2));
|
|
|
656
|
|
|
657 vr = StringToValueRepresentation(std::string(dicom + position + 4, 2), true);
|
|
|
658 if (vr == ValueRepresentation_NotSupported)
|
|
|
659 {
|
|
|
660 return false;
|
|
|
661 }
|
|
|
662
|
|
|
663 if (vr == ValueRepresentation_OtherByte ||
|
|
|
664 vr == ValueRepresentation_OtherDouble ||
|
|
|
665 vr == ValueRepresentation_OtherFloat ||
|
|
|
666 vr == ValueRepresentation_OtherLong ||
|
|
|
667 vr == ValueRepresentation_OtherWord ||
|
|
|
668 vr == ValueRepresentation_Sequence ||
|
|
|
669 vr == ValueRepresentation_UnlimitedCharacters ||
|
|
|
670 vr == ValueRepresentation_UniversalResource ||
|
|
|
671 vr == ValueRepresentation_UnlimitedText ||
|
|
|
672 vr == ValueRepresentation_Unknown) // Note that "UN" should never appear in the Meta Information
|
|
|
673 {
|
|
|
674 if (position + 12 > size)
|
|
|
675 {
|
|
|
676 return false;
|
|
|
677 }
|
|
|
678
|
|
|
679 uint32_t length = ReadUnsignedInteger32(dicom + position + 8);
|
|
|
680
|
|
|
681 if (position + 12 + length > size)
|
|
|
682 {
|
|
|
683 return false;
|
|
|
684 }
|
|
|
685
|
|
|
686 value.assign(dicom + position + 12, length);
|
|
|
687 position += (12 + length);
|
|
|
688 }
|
|
|
689 else
|
|
|
690 {
|
|
|
691 if (position + 8 > size)
|
|
|
692 {
|
|
|
693 return false;
|
|
|
694 }
|
|
|
695
|
|
|
696 uint16_t length = ReadUnsignedInteger16(dicom + position + 6);
|
|
|
697
|
|
|
698 if (position + 8 + length > size)
|
|
|
699 {
|
|
|
700 return false;
|
|
|
701 }
|
|
|
702
|
|
|
703 value.assign(dicom + position + 8, length);
|
|
|
704 position += (8 + length);
|
|
|
705 }
|
|
|
706
|
|
|
707 if (!ValidateTag(vr, value))
|
|
|
708 {
|
|
|
709 return false;
|
|
|
710 }
|
|
|
711
|
|
|
712 RemoveTagPadding(value, vr);
|
|
|
713
|
|
|
714 return true;
|
|
|
715 }
|
|
|
716
|
|
|
717
|
|
|
718 bool DicomMap::ParseDicomMetaInformation(DicomMap& result,
|
|
|
719 const char* dicom,
|
|
|
720 size_t size)
|
|
|
721 {
|
|
|
722 /**
|
|
|
723 * http://dicom.nema.org/medical/dicom/current/output/chtml/part10/chapter_7.html
|
|
|
724 * According to Table 7.1-1, besides the "DICM" DICOM prefix, the
|
|
|
725 * file preamble (i.e. dicom[0..127]) should not be taken into
|
|
|
726 * account to determine whether the file is or is not a DICOM file.
|
|
|
727 **/
|
|
|
728
|
|
|
729 if (size < 132 ||
|
|
|
730 dicom[128] != 'D' ||
|
|
|
731 dicom[129] != 'I' ||
|
|
|
732 dicom[130] != 'C' ||
|
|
|
733 dicom[131] != 'M')
|
|
|
734 {
|
|
|
735 return false;
|
|
|
736 }
|
|
|
737
|
|
|
738
|
|
|
739 /**
|
|
|
740 * The DICOM File Meta Information must be encoded using the
|
|
|
741 * Explicit VR Little Endian Transfer Syntax
|
|
|
742 * (UID=1.2.840.10008.1.2.1).
|
|
|
743 **/
|
|
|
744
|
|
|
745 result.Clear();
|
|
|
746
|
|
|
747 // First, we read the "File Meta Information Group Length" tag
|
|
|
748 // (0002,0000) to know where to stop reading the meta header
|
|
|
749 size_t position = 132;
|
|
|
750
|
|
|
751 DicomTag tag(0x0000, 0x0000); // Dummy initialization
|
|
|
752 ValueRepresentation vr;
|
|
|
753 std::string value;
|
|
|
754 if (!ReadNextTag(tag, vr, value, dicom, size, position) ||
|
|
|
755 tag.GetGroup() != 0x0002 ||
|
|
|
756 tag.GetElement() != 0x0000 ||
|
|
|
757 vr != ValueRepresentation_UnsignedLong ||
|
|
|
758 value.size() != 4)
|
|
|
759 {
|
|
|
760 return false;
|
|
|
761 }
|
|
|
762
|
|
|
763 size_t stopPosition = position + ReadUnsignedInteger32(value.c_str());
|
|
|
764 if (stopPosition > size)
|
|
|
765 {
|
|
|
766 return false;
|
|
|
767 }
|
|
|
768
|
|
|
769 while (position < stopPosition)
|
|
|
770 {
|
|
|
771 if (ReadNextTag(tag, vr, value, dicom, size, position))
|
|
|
772 {
|
|
|
773 result.SetValue(tag, value, IsBinaryValueRepresentation(vr));
|
|
|
774 }
|
|
|
775 else
|
|
|
776 {
|
|
|
777 return false;
|
|
|
778 }
|
|
|
779 }
|
|
|
780
|
|
|
781 return true;
|
|
|
782 }
|
|
|
783 }
|
