0
|
1 /**
|
|
2 * Palantir - A Lightweight, RESTful DICOM Store
|
|
3 * Copyright (C) 2012 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 * This program is distributed in the hope that it will be useful, but
|
|
12 * WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
14 * General Public License for more details.
|
|
15 *
|
|
16 * You should have received a copy of the GNU General Public License
|
|
17 * along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
18 **/
|
|
19
|
|
20
|
|
21 #include "FromDcmtkBridge.h"
|
|
22
|
|
23 #include "ToDcmtkBridge.h"
|
|
24 #include "DicomIntegerPixelAccessor.h"
|
|
25 #include "../Core/PalantirException.h"
|
|
26 #include "../Core/PngWriter.h"
|
|
27 #include "../Core/DicomFormat/DicomString.h"
|
|
28 #include "../Core/DicomFormat/DicomNullValue.h"
|
|
29
|
|
30 #include <boost/locale.hpp>
|
|
31 #include <boost/lexical_cast.hpp>
|
|
32
|
|
33 #include <dcmtk/dcmdata/dcdicent.h>
|
|
34 #include <dcmtk/dcmdata/dcdict.h>
|
|
35 #include <dcmtk/dcmdata/dcelem.h>
|
|
36 #include <dcmtk/dcmdata/dcfilefo.h>
|
|
37 #include <dcmtk/dcmdata/dcistrmb.h>
|
|
38 #include <dcmtk/dcmdata/dcsequen.h>
|
|
39 #include <dcmtk/dcmdata/dcvrfd.h>
|
|
40 #include <dcmtk/dcmdata/dcvrfl.h>
|
|
41 #include <dcmtk/dcmdata/dcvrsl.h>
|
|
42 #include <dcmtk/dcmdata/dcvrss.h>
|
|
43 #include <dcmtk/dcmdata/dcvrul.h>
|
|
44 #include <dcmtk/dcmdata/dcvrus.h>
|
|
45
|
2
|
46 #include <boost/math/special_functions/round.hpp>
|
0
|
47
|
|
48 namespace Palantir
|
|
49 {
|
|
50 void FromDcmtkBridge::Convert(DicomMap& target, DcmDataset& dataset)
|
|
51 {
|
|
52 target.Clear();
|
|
53 for (unsigned long i = 0; i < dataset.card(); i++)
|
|
54 {
|
|
55 DcmElement* element = dataset.getElement(i);
|
|
56 if (element && element->isLeaf())
|
|
57 {
|
|
58 target.SetValue(element->getTag().getGTag(),
|
|
59 element->getTag().getETag(),
|
|
60 ConvertLeafElement(*element));
|
|
61 }
|
|
62 }
|
|
63 }
|
|
64
|
|
65
|
|
66 DicomTag FromDcmtkBridge::GetTag(const DcmElement& element)
|
|
67 {
|
|
68 return DicomTag(element.getGTag(), element.getETag());
|
|
69 }
|
|
70
|
|
71
|
|
72 DicomValue* FromDcmtkBridge::ConvertLeafElement(DcmElement& element)
|
|
73 {
|
|
74 if (!element.isLeaf())
|
|
75 {
|
|
76 throw PalantirException("Only applicable to leaf elements");
|
|
77 }
|
|
78
|
|
79 if (element.isaString())
|
|
80 {
|
|
81 char *c;
|
|
82 if (element.getString(c).good() &&
|
|
83 c != NULL)
|
|
84 {
|
|
85 std::string s(c);
|
|
86 std::string utf8;
|
|
87 try
|
|
88 {
|
|
89 utf8 = boost::locale::conv::to_utf<char>(s, "ISO-8859-1"); // TODO Parameter?
|
|
90 }
|
2
|
91 catch (std::runtime_error&)
|
0
|
92 {
|
|
93 // Bad input string or bad encoding
|
|
94 utf8 = s;
|
|
95 }
|
|
96
|
|
97 return new DicomString(utf8);
|
|
98 }
|
|
99 else
|
|
100 {
|
|
101 return new DicomNullValue;
|
|
102 }
|
|
103 }
|
|
104
|
|
105 try
|
|
106 {
|
|
107 // http://support.dcmtk.org/docs/dcvr_8h-source.html
|
|
108 switch (element.getVR())
|
|
109 {
|
|
110
|
|
111 /**
|
|
112 * TODO.
|
|
113 **/
|
|
114
|
|
115 case EVR_DS: // decimal string
|
|
116 case EVR_IS: // integer string
|
|
117 case EVR_OB: // other byte
|
|
118 case EVR_OF: // other float
|
|
119 case EVR_OW: // other word
|
|
120 case EVR_AS: // age string
|
|
121 case EVR_AT: // attribute tag
|
|
122 case EVR_DA: // date string
|
|
123 case EVR_DT: // date time string
|
|
124 case EVR_TM: // time string
|
|
125 case EVR_UN: // unknown value representation
|
|
126 return new DicomNullValue();
|
|
127
|
|
128
|
|
129 /**
|
|
130 * String types, should never happen at this point because of
|
|
131 * "element.isaString()".
|
|
132 **/
|
|
133
|
|
134 case EVR_AE: // application entity title
|
|
135 case EVR_CS: // code string
|
|
136 case EVR_SH: // short string
|
|
137 case EVR_LO: // long string
|
|
138 case EVR_ST: // short text
|
|
139 case EVR_LT: // long text
|
|
140 case EVR_UT: // unlimited text
|
|
141 case EVR_PN: // person name
|
|
142 case EVR_UI: // unique identifier
|
|
143 return new DicomNullValue();
|
|
144
|
|
145
|
|
146 /**
|
|
147 * Numerical types
|
|
148 **/
|
|
149
|
|
150 case EVR_SL: // signed long
|
|
151 {
|
|
152 Sint32 f;
|
|
153 if (dynamic_cast<DcmSignedLong&>(element).getSint32(f).good())
|
|
154 {
|
|
155 return new DicomString(boost::lexical_cast<std::string>(f));
|
|
156 }
|
|
157 else
|
|
158 {
|
|
159 return new DicomNullValue();
|
|
160 }
|
|
161 }
|
|
162
|
|
163 case EVR_SS: // signed short
|
|
164 {
|
|
165 Sint16 f;
|
|
166 if (dynamic_cast<DcmSignedShort&>(element).getSint16(f).good())
|
|
167 {
|
|
168 return new DicomString(boost::lexical_cast<std::string>(f));
|
|
169 }
|
|
170 else
|
|
171 {
|
|
172 return new DicomNullValue();
|
|
173 }
|
|
174 }
|
|
175
|
|
176 case EVR_UL: // unsigned long
|
|
177 {
|
|
178 Uint32 f;
|
|
179 if (dynamic_cast<DcmUnsignedLong&>(element).getUint32(f).good())
|
|
180 {
|
|
181 return new DicomString(boost::lexical_cast<std::string>(f));
|
|
182 }
|
|
183 else
|
|
184 {
|
|
185 return new DicomNullValue();
|
|
186 }
|
|
187 }
|
|
188
|
|
189 case EVR_US: // unsigned short
|
|
190 {
|
|
191 Uint16 f;
|
|
192 if (dynamic_cast<DcmUnsignedShort&>(element).getUint16(f).good())
|
|
193 {
|
|
194 return new DicomString(boost::lexical_cast<std::string>(f));
|
|
195 }
|
|
196 else
|
|
197 {
|
|
198 return new DicomNullValue();
|
|
199 }
|
|
200 }
|
|
201
|
|
202 case EVR_FL: // float single-precision
|
|
203 {
|
|
204 Float32 f;
|
|
205 if (dynamic_cast<DcmFloatingPointSingle&>(element).getFloat32(f).good())
|
|
206 {
|
|
207 return new DicomString(boost::lexical_cast<std::string>(f));
|
|
208 }
|
|
209 else
|
|
210 {
|
|
211 return new DicomNullValue();
|
|
212 }
|
|
213 }
|
|
214
|
|
215 case EVR_FD: // float double-precision
|
|
216 {
|
|
217 Float64 f;
|
|
218 if (dynamic_cast<DcmFloatingPointDouble&>(element).getFloat64(f).good())
|
|
219 {
|
|
220 return new DicomString(boost::lexical_cast<std::string>(f));
|
|
221 }
|
|
222 else
|
|
223 {
|
|
224 return new DicomNullValue();
|
|
225 }
|
|
226 }
|
|
227
|
|
228
|
|
229 /**
|
|
230 * Sequence types, should never occur at this point because of
|
|
231 * "element.isLeaf()".
|
|
232 **/
|
|
233
|
|
234 case EVR_SQ: // sequence of items
|
|
235 return new DicomNullValue;
|
|
236
|
|
237
|
|
238 /**
|
|
239 * Internal to DCMTK.
|
|
240 **/
|
|
241
|
|
242 case EVR_ox: // OB or OW depending on context
|
|
243 case EVR_xs: // SS or US depending on context
|
|
244 case EVR_lt: // US, SS or OW depending on context, used for LUT Data (thus the name)
|
|
245 case EVR_na: // na="not applicable", for data which has no VR
|
|
246 case EVR_up: // up="unsigned pointer", used internally for DICOMDIR suppor
|
|
247 case EVR_item: // used internally for items
|
|
248 case EVR_metainfo: // used internally for meta info datasets
|
|
249 case EVR_dataset: // used internally for datasets
|
|
250 case EVR_fileFormat: // used internally for DICOM files
|
|
251 case EVR_dicomDir: // used internally for DICOMDIR objects
|
|
252 case EVR_dirRecord: // used internally for DICOMDIR records
|
|
253 case EVR_pixelSQ: // used internally for pixel sequences in a compressed image
|
|
254 case EVR_pixelItem: // used internally for pixel items in a compressed image
|
|
255 case EVR_UNKNOWN: // used internally for elements with unknown VR (encoded with 4-byte length field in explicit VR)
|
|
256 case EVR_PixelData: // used internally for uncompressed pixeld data
|
|
257 case EVR_OverlayData: // used internally for overlay data
|
|
258 case EVR_UNKNOWN2B: // used internally for elements with unknown VR with 2-byte length field in explicit VR
|
|
259 return new DicomNullValue;
|
|
260
|
|
261
|
|
262 /**
|
|
263 * Default case.
|
|
264 **/
|
|
265
|
|
266 default:
|
|
267 return new DicomNullValue;
|
|
268 }
|
|
269 }
|
|
270 catch (boost::bad_lexical_cast)
|
|
271 {
|
|
272 return new DicomNullValue;
|
|
273 }
|
|
274 }
|
|
275
|
|
276
|
|
277 static void StoreElement(Json::Value& target,
|
|
278 DcmElement& element,
|
|
279 unsigned int maxStringLength);
|
|
280
|
|
281 static void StoreItem(Json::Value& target,
|
|
282 DcmItem& item,
|
|
283 unsigned int maxStringLength)
|
|
284 {
|
|
285 target = Json::Value(Json::objectValue);
|
|
286
|
|
287 for (unsigned long i = 0; i < item.card(); i++)
|
|
288 {
|
|
289 DcmElement* element = item.getElement(i);
|
|
290 StoreElement(target, *element, maxStringLength);
|
|
291 }
|
|
292 }
|
|
293
|
|
294
|
|
295 static void StoreElement(Json::Value& target,
|
|
296 DcmElement& element,
|
|
297 unsigned int maxStringLength)
|
|
298 {
|
|
299 assert(target.type() == Json::objectValue);
|
|
300
|
|
301 const std::string tagName = FromDcmtkBridge::GetName(FromDcmtkBridge::GetTag(element));
|
|
302
|
|
303 if (element.isLeaf())
|
|
304 {
|
|
305 std::auto_ptr<DicomValue> v(FromDcmtkBridge::ConvertLeafElement(element));
|
|
306 if (v->IsNull())
|
|
307 {
|
|
308 target[tagName] = Json::nullValue;
|
|
309 }
|
|
310 else
|
|
311 {
|
|
312 std::string s = v->AsString();
|
|
313 if (maxStringLength == 0 ||
|
|
314 s.size() <= maxStringLength)
|
|
315 {
|
|
316 target[tagName] = s;
|
|
317 }
|
|
318 else
|
|
319 {
|
|
320 // An integer value of 0 in JSON indicates too long field
|
|
321 target[tagName] = 0;
|
|
322 }
|
|
323 }
|
|
324 }
|
|
325 else
|
|
326 {
|
|
327 target[tagName] = Json::Value(Json::arrayValue);
|
|
328
|
|
329 // "All subclasses of DcmElement except for DcmSequenceOfItems
|
|
330 // are leaf nodes, while DcmSequenceOfItems, DcmItem, DcmDataset
|
|
331 // etc. are not." The following cast is thus OK.
|
|
332 DcmSequenceOfItems& sequence = dynamic_cast<DcmSequenceOfItems&>(element);
|
|
333
|
|
334 for (unsigned long i = 0; i < sequence.card(); i++)
|
|
335 {
|
|
336 DcmItem* child = sequence.getItem(i);
|
|
337 Json::Value& v = target[tagName].append(Json::objectValue);
|
|
338 StoreItem(v, *child, maxStringLength);
|
|
339 }
|
|
340 }
|
|
341 }
|
|
342
|
|
343
|
|
344 void FromDcmtkBridge::ToJson(Json::Value& root,
|
|
345 DcmDataset& dataset,
|
|
346 unsigned int maxStringLength)
|
|
347 {
|
|
348 StoreItem(root, dataset, maxStringLength);
|
|
349 }
|
|
350
|
|
351
|
|
352
|
|
353 void FromDcmtkBridge::ToJson(Json::Value& target,
|
|
354 const std::string& path,
|
|
355 unsigned int maxStringLength)
|
|
356 {
|
|
357 DcmFileFormat dicom;
|
|
358 if (!dicom.loadFile(path.c_str()).good())
|
|
359 {
|
|
360 throw PalantirException(ErrorCode_BadFileFormat);
|
|
361 }
|
|
362 else
|
|
363 {
|
|
364 FromDcmtkBridge::ToJson(target, *dicom.getDataset(), maxStringLength);
|
|
365 }
|
|
366 }
|
|
367
|
|
368
|
|
369 void FromDcmtkBridge::ExtractNormalizedImage(std::string& result,
|
|
370 DcmDataset& dataset)
|
|
371 {
|
|
372 // See also: http://support.dcmtk.org/wiki/dcmtk/howto/accessing-compressed-data
|
|
373
|
|
374 PngWriter w;
|
|
375 std::auto_ptr<DicomIntegerPixelAccessor> accessor;
|
|
376
|
|
377 DicomMap m;
|
|
378 FromDcmtkBridge::Convert(m, dataset);
|
|
379
|
|
380 DcmElement* e;
|
|
381 if (dataset.findAndGetElement(ToDcmtkBridge::Convert(DicomTag::PIXEL_DATA), e).good() &&
|
|
382 e != NULL)
|
|
383 {
|
|
384 Uint8* pixData = NULL;
|
|
385 if (e->getUint8Array(pixData) == EC_Normal)
|
|
386 {
|
|
387 accessor.reset(new DicomIntegerPixelAccessor(m, pixData, e->getLength()));
|
|
388 }
|
|
389 }
|
|
390
|
|
391 if (accessor.get() == NULL ||
|
|
392 accessor->GetWidth() == 0 ||
|
|
393 accessor->GetHeight() == 0)
|
|
394 {
|
|
395 w.WriteToMemory(result, 0, 0, 0, PixelFormat_Grayscale8, NULL);
|
|
396 }
|
|
397 else
|
|
398 {
|
|
399 int32_t min, max;
|
|
400 accessor->GetExtremeValues(min, max);
|
|
401
|
|
402 std::vector<uint8_t> image(accessor->GetWidth() * accessor->GetHeight(), 0);
|
|
403 if (min != max)
|
|
404 {
|
|
405 uint8_t* result = &image[0];
|
|
406 for (unsigned int y = 0; y < accessor->GetHeight(); y++)
|
|
407 {
|
|
408 for (unsigned int x = 0; x < accessor->GetWidth(); x++, result++)
|
|
409 {
|
|
410 int32_t v = accessor->GetValue(x, y);
|
2
|
411 *result = static_cast<uint8_t>(boost::math::lround(static_cast<float>(v - min) / static_cast<float>(max - min) * 255.0f));
|
0
|
412 }
|
|
413 }
|
|
414 }
|
|
415
|
|
416 w.WriteToMemory(result, accessor->GetWidth(), accessor->GetHeight(),
|
|
417 accessor->GetWidth(), PixelFormat_Grayscale8, &image[0]);
|
|
418 }
|
|
419 }
|
|
420
|
|
421
|
|
422 void FromDcmtkBridge::ExtractNormalizedImage(std::string& result,
|
|
423 const std::string& dicomContent)
|
|
424 {
|
|
425 DcmInputBufferStream is;
|
|
426 if (dicomContent.size() > 0)
|
|
427 {
|
|
428 is.setBuffer(&dicomContent[0], dicomContent.size());
|
|
429 }
|
|
430 is.setEos();
|
|
431
|
|
432 DcmFileFormat dicom;
|
|
433 if (dicom.read(is).good())
|
|
434 {
|
|
435 ExtractNormalizedImage(result, *dicom.getDataset());
|
|
436 }
|
|
437 else
|
|
438 {
|
|
439 throw PalantirException(ErrorCode_BadFileFormat);
|
|
440 }
|
|
441 }
|
|
442
|
|
443
|
|
444
|
|
445 std::string FromDcmtkBridge::GetName(const DicomTag& t)
|
|
446 {
|
|
447 DcmTagKey tag(t.GetGroup(), t.GetElement());
|
|
448 const DcmDataDictionary& dict = dcmDataDict.rdlock();
|
|
449 const DcmDictEntry* entry = dict.findEntry(tag, NULL);
|
|
450
|
|
451 std::string s("Unknown");
|
|
452 if (entry != NULL)
|
|
453 {
|
|
454 s = std::string(entry->getTagName());
|
|
455 }
|
|
456
|
|
457 dcmDataDict.unlock();
|
|
458 return s;
|
|
459 }
|
|
460
|
|
461
|
|
462 DicomTag FromDcmtkBridge::FindTag(const char* name)
|
|
463 {
|
|
464 const DcmDataDictionary& dict = dcmDataDict.rdlock();
|
|
465 const DcmDictEntry* entry = dict.findEntry(name);
|
|
466
|
|
467 if (entry == NULL)
|
|
468 {
|
|
469 dcmDataDict.unlock();
|
|
470 throw PalantirException("Unknown DICOM tag");
|
|
471 }
|
|
472 else
|
|
473 {
|
|
474 DcmTagKey key = entry->getKey();
|
|
475 DicomTag tag(key.getGroup(), key.getElement());
|
|
476 dcmDataDict.unlock();
|
|
477 return tag;
|
|
478 }
|
|
479 }
|
|
480
|
|
481
|
|
482 void FromDcmtkBridge::Print(FILE* fp, const DicomMap& m)
|
|
483 {
|
|
484 for (DicomMap::Map::const_iterator
|
|
485 it = m.map_.begin(); it != m.map_.end(); it++)
|
|
486 {
|
|
487 DicomTag t = it->first;
|
|
488 std::string s = it->second->AsString();
|
|
489 printf("0x%04x 0x%04x (%s) [%s]\n", t.GetGroup(), t.GetElement(), GetName(t).c_str(), s.c_str());
|
|
490 }
|
|
491 }
|
|
492
|
|
493
|
|
494 void FromDcmtkBridge::ToJson(Json::Value& result,
|
|
495 const DicomMap& values)
|
|
496 {
|
|
497 if (result.type() != Json::objectValue)
|
|
498 {
|
|
499 throw PalantirException(ErrorCode_BadParameterType);
|
|
500 }
|
|
501
|
|
502 result.clear();
|
|
503
|
|
504 for (DicomMap::Map::const_iterator
|
|
505 it = values.map_.begin(); it != values.map_.end(); it++)
|
|
506 {
|
|
507 result[GetName(it->first)] = it->second->AsString();
|
|
508 }
|
|
509 }
|
|
510 }
|