790
|
1 /**
|
|
2 * Orthanc - A Lightweight, RESTful DICOM Store
|
|
3 * Copyright (C) 2012-2014 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
|
|
33
|
|
34 /*=========================================================================
|
|
35
|
|
36 This file is based on portions of the following project:
|
|
37
|
|
38 Program: GDCM (Grassroots DICOM). A DICOM library
|
|
39 Module: http://gdcm.sourceforge.net/Copyright.html
|
|
40
|
|
41 Copyright (c) 2006-2011 Mathieu Malaterre
|
|
42 Copyright (c) 1993-2005 CREATIS
|
|
43 (CREATIS = Centre de Recherche et d'Applications en Traitement de l'Image)
|
|
44 All rights reserved.
|
|
45
|
|
46 Redistribution and use in source and binary forms, with or without
|
|
47 modification, are permitted provided that the following conditions are met:
|
|
48
|
|
49 * Redistributions of source code must retain the above copyright notice,
|
|
50 this list of conditions and the following disclaimer.
|
|
51
|
|
52 * Redistributions in binary form must reproduce the above copyright notice,
|
|
53 this list of conditions and the following disclaimer in the documentation
|
|
54 and/or other materials provided with the distribution.
|
|
55
|
|
56 * Neither name of Mathieu Malaterre, or CREATIS, nor the names of any
|
|
57 contributors (CNRS, INSERM, UCB, Universite Lyon I), may be used to
|
|
58 endorse or promote products derived from this software without specific
|
|
59 prior written permission.
|
|
60
|
|
61 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS ``AS IS''
|
|
62 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
63 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
64 ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHORS OR CONTRIBUTORS BE LIABLE FOR
|
|
65 ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
|
|
66 DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
|
|
67 SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
|
|
68 CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
|
|
69 OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
|
|
70 OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
71
|
|
72 =========================================================================*/
|
|
73
|
|
74
|
|
75 #ifndef NOMINMAX
|
|
76 #define NOMINMAX
|
|
77 #endif
|
|
78
|
|
79 #include "ParsedDicomFile.h"
|
|
80
|
|
81 #include "FromDcmtkBridge.h"
|
|
82 #include "ToDcmtkBridge.h"
|
|
83 #include "../Core/Toolbox.h"
|
|
84 #include "../Core/OrthancException.h"
|
|
85 #include "../Core/FileFormats/PngWriter.h"
|
|
86 #include "../Core/Uuid.h"
|
|
87 #include "../Core/DicomFormat/DicomString.h"
|
|
88 #include "../Core/DicomFormat/DicomNullValue.h"
|
|
89 #include "../Core/DicomFormat/DicomIntegerPixelAccessor.h"
|
|
90
|
|
91 #include <list>
|
|
92 #include <limits>
|
|
93
|
|
94 #include <boost/lexical_cast.hpp>
|
|
95
|
|
96 #include <dcmtk/dcmdata/dcchrstr.h>
|
|
97 #include <dcmtk/dcmdata/dcdicent.h>
|
|
98 #include <dcmtk/dcmdata/dcdict.h>
|
|
99 #include <dcmtk/dcmdata/dcfilefo.h>
|
|
100 #include <dcmtk/dcmdata/dcistrmb.h>
|
|
101 #include <dcmtk/dcmdata/dcuid.h>
|
|
102 #include <dcmtk/dcmdata/dcmetinf.h>
|
|
103
|
|
104 #include <dcmtk/dcmdata/dcvrae.h>
|
|
105 #include <dcmtk/dcmdata/dcvras.h>
|
|
106 #include <dcmtk/dcmdata/dcvrcs.h>
|
|
107 #include <dcmtk/dcmdata/dcvrda.h>
|
|
108 #include <dcmtk/dcmdata/dcvrds.h>
|
|
109 #include <dcmtk/dcmdata/dcvrdt.h>
|
|
110 #include <dcmtk/dcmdata/dcvrfd.h>
|
|
111 #include <dcmtk/dcmdata/dcvrfl.h>
|
|
112 #include <dcmtk/dcmdata/dcvris.h>
|
|
113 #include <dcmtk/dcmdata/dcvrlo.h>
|
|
114 #include <dcmtk/dcmdata/dcvrlt.h>
|
|
115 #include <dcmtk/dcmdata/dcvrpn.h>
|
|
116 #include <dcmtk/dcmdata/dcvrsh.h>
|
|
117 #include <dcmtk/dcmdata/dcvrsl.h>
|
|
118 #include <dcmtk/dcmdata/dcvrss.h>
|
|
119 #include <dcmtk/dcmdata/dcvrst.h>
|
|
120 #include <dcmtk/dcmdata/dcvrtm.h>
|
|
121 #include <dcmtk/dcmdata/dcvrui.h>
|
|
122 #include <dcmtk/dcmdata/dcvrul.h>
|
|
123 #include <dcmtk/dcmdata/dcvrus.h>
|
|
124 #include <dcmtk/dcmdata/dcvrut.h>
|
|
125 #include <dcmtk/dcmdata/dcpixel.h>
|
|
126 #include <dcmtk/dcmdata/dcpixseq.h>
|
|
127 #include <dcmtk/dcmdata/dcpxitem.h>
|
|
128
|
|
129
|
|
130 #include <boost/math/special_functions/round.hpp>
|
|
131 #include <glog/logging.h>
|
|
132 #include <dcmtk/dcmdata/dcostrmb.h>
|
|
133
|
|
134
|
|
135 static const char* CONTENT_TYPE_OCTET_STREAM = "application/octet-stream";
|
|
136
|
|
137
|
|
138
|
|
139 namespace Orthanc
|
|
140 {
|
794
|
141 struct ParsedDicomFile::PImpl
|
|
142 {
|
|
143 std::auto_ptr<DcmFileFormat> file_;
|
|
144 };
|
|
145
|
|
146
|
|
147 // This method can only be called from the constructors!
|
790
|
148 void ParsedDicomFile::Setup(const char* buffer, size_t size)
|
|
149 {
|
|
150 DcmInputBufferStream is;
|
|
151 if (size > 0)
|
|
152 {
|
|
153 is.setBuffer(buffer, size);
|
|
154 }
|
|
155 is.setEos();
|
|
156
|
794
|
157 pimpl_->file_.reset(new DcmFileFormat);
|
|
158 pimpl_->file_->transferInit();
|
|
159 if (!pimpl_->file_->read(is).good())
|
790
|
160 {
|
794
|
161 delete pimpl_; // Avoid a memory leak due to exception
|
|
162 // throwing, as we are in the constructor
|
|
163
|
790
|
164 throw OrthancException(ErrorCode_BadFileFormat);
|
|
165 }
|
794
|
166 pimpl_->file_->loadAllDataIntoMemory();
|
|
167 pimpl_->file_->transferEnd();
|
790
|
168 }
|
|
169
|
|
170
|
|
171 static void SendPathValueForDictionary(RestApiOutput& output,
|
|
172 DcmItem& dicom)
|
|
173 {
|
|
174 Json::Value v = Json::arrayValue;
|
|
175
|
|
176 for (unsigned long i = 0; i < dicom.card(); i++)
|
|
177 {
|
|
178 DcmElement* element = dicom.getElement(i);
|
|
179 if (element)
|
|
180 {
|
|
181 char buf[16];
|
|
182 sprintf(buf, "%04x-%04x", element->getTag().getGTag(), element->getTag().getETag());
|
|
183 v.append(buf);
|
|
184 }
|
|
185 }
|
|
186
|
|
187 output.AnswerJson(v);
|
|
188 }
|
|
189
|
|
190 static inline uint16_t GetCharValue(char c)
|
|
191 {
|
|
192 if (c >= '0' && c <= '9')
|
|
193 return c - '0';
|
|
194 else if (c >= 'a' && c <= 'f')
|
|
195 return c - 'a' + 10;
|
|
196 else if (c >= 'A' && c <= 'F')
|
|
197 return c - 'A' + 10;
|
|
198 else
|
|
199 return 0;
|
|
200 }
|
|
201
|
|
202 static inline uint16_t GetTagValue(const char* c)
|
|
203 {
|
|
204 return ((GetCharValue(c[0]) << 12) +
|
|
205 (GetCharValue(c[1]) << 8) +
|
|
206 (GetCharValue(c[2]) << 4) +
|
|
207 GetCharValue(c[3]));
|
|
208 }
|
|
209
|
|
210 static void ParseTagAndGroup(DcmTagKey& key,
|
|
211 const std::string& tag)
|
|
212 {
|
|
213 DicomTag t = FromDcmtkBridge::ParseTag(tag);
|
|
214 key = DcmTagKey(t.GetGroup(), t.GetElement());
|
|
215 }
|
|
216
|
|
217
|
|
218 static void SendSequence(RestApiOutput& output,
|
|
219 DcmSequenceOfItems& sequence)
|
|
220 {
|
|
221 // This element is a sequence
|
|
222 Json::Value v = Json::arrayValue;
|
|
223
|
|
224 for (unsigned long i = 0; i < sequence.card(); i++)
|
|
225 {
|
|
226 v.append(boost::lexical_cast<std::string>(i));
|
|
227 }
|
|
228
|
|
229 output.AnswerJson(v);
|
|
230 }
|
|
231
|
|
232
|
|
233 static unsigned int GetPixelDataBlockCount(DcmPixelData& pixelData,
|
|
234 E_TransferSyntax transferSyntax)
|
|
235 {
|
|
236 DcmPixelSequence* pixelSequence = NULL;
|
|
237 if (pixelData.getEncapsulatedRepresentation
|
|
238 (transferSyntax, NULL, pixelSequence).good() && pixelSequence)
|
|
239 {
|
|
240 return pixelSequence->card();
|
|
241 }
|
|
242 else
|
|
243 {
|
|
244 return 1;
|
|
245 }
|
|
246 }
|
|
247
|
|
248
|
|
249 static void AnswerDicomField(RestApiOutput& output,
|
|
250 DcmElement& element,
|
|
251 E_TransferSyntax transferSyntax)
|
|
252 {
|
|
253 // This element is nor a sequence, neither a pixel-data
|
|
254 std::string buffer;
|
|
255 buffer.resize(65536);
|
|
256 Uint32 length = element.getLength(transferSyntax);
|
|
257 Uint32 offset = 0;
|
|
258
|
|
259 output.GetLowLevelOutput().SendOkHeader(CONTENT_TYPE_OCTET_STREAM, true, length, NULL);
|
|
260
|
|
261 while (offset < length)
|
|
262 {
|
|
263 Uint32 nbytes;
|
|
264 if (length - offset < buffer.size())
|
|
265 {
|
|
266 nbytes = length - offset;
|
|
267 }
|
|
268 else
|
|
269 {
|
|
270 nbytes = buffer.size();
|
|
271 }
|
|
272
|
|
273 OFCondition cond = element.getPartialValue(&buffer[0], offset, nbytes);
|
|
274
|
|
275 if (cond.good())
|
|
276 {
|
|
277 output.GetLowLevelOutput().Send(&buffer[0], nbytes);
|
|
278 offset += nbytes;
|
|
279 }
|
|
280 else
|
|
281 {
|
|
282 LOG(ERROR) << "Error while sending a DICOM field: " << cond.text();
|
|
283 return;
|
|
284 }
|
|
285 }
|
|
286
|
|
287 output.MarkLowLevelOutputDone();
|
|
288 }
|
|
289
|
|
290
|
|
291 static bool AnswerPixelData(RestApiOutput& output,
|
|
292 DcmItem& dicom,
|
|
293 E_TransferSyntax transferSyntax,
|
|
294 const std::string* blockUri)
|
|
295 {
|
|
296 DcmTag k(DICOM_TAG_PIXEL_DATA.GetGroup(),
|
|
297 DICOM_TAG_PIXEL_DATA.GetElement());
|
|
298
|
|
299 DcmElement *element = NULL;
|
|
300 if (!dicom.findAndGetElement(k, element).good() ||
|
|
301 element == NULL)
|
|
302 {
|
|
303 return false;
|
|
304 }
|
|
305
|
|
306 try
|
|
307 {
|
|
308 DcmPixelData& pixelData = dynamic_cast<DcmPixelData&>(*element);
|
|
309 if (blockUri == NULL)
|
|
310 {
|
|
311 // The user asks how many blocks are presents in this pixel data
|
|
312 unsigned int blocks = GetPixelDataBlockCount(pixelData, transferSyntax);
|
|
313
|
|
314 Json::Value result(Json::arrayValue);
|
|
315 for (unsigned int i = 0; i < blocks; i++)
|
|
316 {
|
|
317 result.append(boost::lexical_cast<std::string>(i));
|
|
318 }
|
|
319
|
|
320 output.AnswerJson(result);
|
|
321 return true;
|
|
322 }
|
|
323
|
|
324
|
|
325 unsigned int block = boost::lexical_cast<unsigned int>(*blockUri);
|
|
326
|
|
327 if (block < GetPixelDataBlockCount(pixelData, transferSyntax))
|
|
328 {
|
|
329 DcmPixelSequence* pixelSequence = NULL;
|
|
330 if (pixelData.getEncapsulatedRepresentation
|
|
331 (transferSyntax, NULL, pixelSequence).good() && pixelSequence)
|
|
332 {
|
|
333 // This is the case for JPEG transfer syntaxes
|
|
334 if (block < pixelSequence->card())
|
|
335 {
|
|
336 DcmPixelItem* pixelItem = NULL;
|
|
337 if (pixelSequence->getItem(pixelItem, block).good() && pixelItem)
|
|
338 {
|
|
339 if (pixelItem->getLength() == 0)
|
|
340 {
|
|
341 output.AnswerBuffer(NULL, 0, CONTENT_TYPE_OCTET_STREAM);
|
|
342 return true;
|
|
343 }
|
|
344
|
|
345 Uint8* buffer = NULL;
|
|
346 if (pixelItem->getUint8Array(buffer).good() && buffer)
|
|
347 {
|
|
348 output.AnswerBuffer(buffer, pixelItem->getLength(), CONTENT_TYPE_OCTET_STREAM);
|
|
349 return true;
|
|
350 }
|
|
351 }
|
|
352 }
|
|
353 }
|
|
354 else
|
|
355 {
|
|
356 // This is the case for raw, uncompressed image buffers
|
|
357 assert(*blockUri == "0");
|
|
358 AnswerDicomField(output, *element, transferSyntax);
|
|
359 }
|
|
360 }
|
|
361 }
|
|
362 catch (boost::bad_lexical_cast&)
|
|
363 {
|
|
364 // The URI entered by the user is not a number
|
|
365 }
|
|
366 catch (std::bad_cast&)
|
|
367 {
|
|
368 // This should never happen
|
|
369 }
|
|
370
|
|
371 return false;
|
|
372 }
|
|
373
|
|
374
|
|
375
|
|
376 static void SendPathValueForLeaf(RestApiOutput& output,
|
|
377 const std::string& tag,
|
|
378 DcmItem& dicom,
|
|
379 E_TransferSyntax transferSyntax)
|
|
380 {
|
|
381 DcmTagKey k;
|
|
382 ParseTagAndGroup(k, tag);
|
|
383
|
|
384 DcmSequenceOfItems* sequence = NULL;
|
|
385 if (dicom.findAndGetSequence(k, sequence).good() &&
|
|
386 sequence != NULL &&
|
|
387 sequence->getVR() == EVR_SQ)
|
|
388 {
|
|
389 SendSequence(output, *sequence);
|
|
390 return;
|
|
391 }
|
|
392
|
|
393 DcmElement* element = NULL;
|
|
394 if (dicom.findAndGetElement(k, element).good() &&
|
|
395 element != NULL &&
|
|
396 //element->getVR() != EVR_UNKNOWN && // This would forbid private tags
|
|
397 element->getVR() != EVR_SQ)
|
|
398 {
|
|
399 AnswerDicomField(output, *element, transferSyntax);
|
|
400 }
|
|
401 }
|
|
402
|
|
403 void ParsedDicomFile::SendPathValue(RestApiOutput& output,
|
|
404 const UriComponents& uri)
|
|
405 {
|
794
|
406 DcmItem* dicom = pimpl_->file_->getDataset();
|
|
407 E_TransferSyntax transferSyntax = pimpl_->file_->getDataset()->getOriginalXfer();
|
790
|
408
|
|
409 // Special case: Accessing the pixel data
|
|
410 if (uri.size() == 1 ||
|
|
411 uri.size() == 2)
|
|
412 {
|
|
413 DcmTagKey tag;
|
|
414 ParseTagAndGroup(tag, uri[0]);
|
|
415
|
|
416 if (tag.getGroup() == DICOM_TAG_PIXEL_DATA.GetGroup() &&
|
|
417 tag.getElement() == DICOM_TAG_PIXEL_DATA.GetElement())
|
|
418 {
|
|
419 AnswerPixelData(output, *dicom, transferSyntax, uri.size() == 1 ? NULL : &uri[1]);
|
|
420 return;
|
|
421 }
|
|
422 }
|
|
423
|
|
424 // Go down in the tag hierarchy according to the URI
|
|
425 for (size_t pos = 0; pos < uri.size() / 2; pos++)
|
|
426 {
|
|
427 size_t index;
|
|
428 try
|
|
429 {
|
|
430 index = boost::lexical_cast<size_t>(uri[2 * pos + 1]);
|
|
431 }
|
|
432 catch (boost::bad_lexical_cast&)
|
|
433 {
|
|
434 return;
|
|
435 }
|
|
436
|
|
437 DcmTagKey k;
|
|
438 DcmItem *child = NULL;
|
|
439 ParseTagAndGroup(k, uri[2 * pos]);
|
|
440 if (!dicom->findAndGetSequenceItem(k, child, index).good() ||
|
|
441 child == NULL)
|
|
442 {
|
|
443 return;
|
|
444 }
|
|
445
|
|
446 dicom = child;
|
|
447 }
|
|
448
|
|
449 // We have reached the end of the URI
|
|
450 if (uri.size() % 2 == 0)
|
|
451 {
|
|
452 SendPathValueForDictionary(output, *dicom);
|
|
453 }
|
|
454 else
|
|
455 {
|
|
456 SendPathValueForLeaf(output, uri.back(), *dicom, transferSyntax);
|
|
457 }
|
|
458 }
|
|
459
|
|
460
|
|
461
|
|
462
|
|
463
|
|
464 static DcmElement* CreateElementForTag(const DicomTag& tag)
|
|
465 {
|
|
466 DcmTag key(tag.GetGroup(), tag.GetElement());
|
|
467
|
|
468 switch (key.getEVR())
|
|
469 {
|
|
470 // http://support.dcmtk.org/docs/dcvr_8h-source.html
|
|
471
|
|
472 /**
|
|
473 * TODO.
|
|
474 **/
|
|
475
|
|
476 case EVR_OB: // other byte
|
|
477 case EVR_OF: // other float
|
|
478 case EVR_OW: // other word
|
|
479 case EVR_AT: // attribute tag
|
|
480 throw OrthancException(ErrorCode_NotImplemented);
|
|
481
|
|
482 case EVR_UN: // unknown value representation
|
|
483 throw OrthancException(ErrorCode_ParameterOutOfRange);
|
|
484
|
|
485
|
|
486 /**
|
|
487 * String types.
|
|
488 * http://support.dcmtk.org/docs/classDcmByteString.html
|
|
489 **/
|
|
490
|
|
491 case EVR_AS: // age string
|
|
492 return new DcmAgeString(key);
|
|
493
|
|
494 case EVR_AE: // application entity title
|
|
495 return new DcmApplicationEntity(key);
|
|
496
|
|
497 case EVR_CS: // code string
|
|
498 return new DcmCodeString(key);
|
|
499
|
|
500 case EVR_DA: // date string
|
|
501 return new DcmDate(key);
|
|
502
|
|
503 case EVR_DT: // date time string
|
|
504 return new DcmDateTime(key);
|
|
505
|
|
506 case EVR_DS: // decimal string
|
|
507 return new DcmDecimalString(key);
|
|
508
|
|
509 case EVR_IS: // integer string
|
|
510 return new DcmIntegerString(key);
|
|
511
|
|
512 case EVR_TM: // time string
|
|
513 return new DcmTime(key);
|
|
514
|
|
515 case EVR_UI: // unique identifier
|
|
516 return new DcmUniqueIdentifier(key);
|
|
517
|
|
518 case EVR_ST: // short text
|
|
519 return new DcmShortText(key);
|
|
520
|
|
521 case EVR_LO: // long string
|
|
522 return new DcmLongString(key);
|
|
523
|
|
524 case EVR_LT: // long text
|
|
525 return new DcmLongText(key);
|
|
526
|
|
527 case EVR_UT: // unlimited text
|
|
528 return new DcmUnlimitedText(key);
|
|
529
|
|
530 case EVR_SH: // short string
|
|
531 return new DcmShortString(key);
|
|
532
|
|
533 case EVR_PN: // person name
|
|
534 return new DcmPersonName(key);
|
|
535
|
|
536
|
|
537 /**
|
|
538 * Numerical types
|
|
539 **/
|
|
540
|
|
541 case EVR_SL: // signed long
|
|
542 return new DcmSignedLong(key);
|
|
543
|
|
544 case EVR_SS: // signed short
|
|
545 return new DcmSignedShort(key);
|
|
546
|
|
547 case EVR_UL: // unsigned long
|
|
548 return new DcmUnsignedLong(key);
|
|
549
|
|
550 case EVR_US: // unsigned short
|
|
551 return new DcmUnsignedShort(key);
|
|
552
|
|
553 case EVR_FL: // float single-precision
|
|
554 return new DcmFloatingPointSingle(key);
|
|
555
|
|
556 case EVR_FD: // float double-precision
|
|
557 return new DcmFloatingPointDouble(key);
|
|
558
|
|
559
|
|
560 /**
|
|
561 * Sequence types, should never occur at this point.
|
|
562 **/
|
|
563
|
|
564 case EVR_SQ: // sequence of items
|
|
565 throw OrthancException(ErrorCode_ParameterOutOfRange);
|
|
566
|
|
567
|
|
568 /**
|
|
569 * Internal to DCMTK.
|
|
570 **/
|
|
571
|
|
572 case EVR_ox: // OB or OW depending on context
|
|
573 case EVR_xs: // SS or US depending on context
|
|
574 case EVR_lt: // US, SS or OW depending on context, used for LUT Data (thus the name)
|
|
575 case EVR_na: // na="not applicable", for data which has no VR
|
|
576 case EVR_up: // up="unsigned pointer", used internally for DICOMDIR suppor
|
|
577 case EVR_item: // used internally for items
|
|
578 case EVR_metainfo: // used internally for meta info datasets
|
|
579 case EVR_dataset: // used internally for datasets
|
|
580 case EVR_fileFormat: // used internally for DICOM files
|
|
581 case EVR_dicomDir: // used internally for DICOMDIR objects
|
|
582 case EVR_dirRecord: // used internally for DICOMDIR records
|
|
583 case EVR_pixelSQ: // used internally for pixel sequences in a compressed image
|
|
584 case EVR_pixelItem: // used internally for pixel items in a compressed image
|
|
585 case EVR_UNKNOWN: // used internally for elements with unknown VR (encoded with 4-byte length field in explicit VR)
|
|
586 case EVR_PixelData: // used internally for uncompressed pixeld data
|
|
587 case EVR_OverlayData: // used internally for overlay data
|
|
588 case EVR_UNKNOWN2B: // used internally for elements with unknown VR with 2-byte length field in explicit VR
|
|
589 default:
|
|
590 break;
|
|
591 }
|
|
592
|
|
593 throw OrthancException(ErrorCode_InternalError);
|
|
594 }
|
|
595
|
|
596
|
|
597
|
|
598 static void FillElementWithString(DcmElement& element,
|
|
599 const DicomTag& tag,
|
|
600 const std::string& value)
|
|
601 {
|
|
602 DcmTag key(tag.GetGroup(), tag.GetElement());
|
|
603 bool ok = false;
|
|
604
|
|
605 try
|
|
606 {
|
|
607 switch (key.getEVR())
|
|
608 {
|
|
609 // http://support.dcmtk.org/docs/dcvr_8h-source.html
|
|
610
|
|
611 /**
|
|
612 * TODO.
|
|
613 **/
|
|
614
|
|
615 case EVR_OB: // other byte
|
|
616 case EVR_OF: // other float
|
|
617 case EVR_OW: // other word
|
|
618 case EVR_AT: // attribute tag
|
|
619 throw OrthancException(ErrorCode_NotImplemented);
|
|
620
|
|
621 case EVR_UN: // unknown value representation
|
|
622 throw OrthancException(ErrorCode_ParameterOutOfRange);
|
|
623
|
|
624
|
|
625 /**
|
|
626 * String types.
|
|
627 **/
|
|
628
|
|
629 case EVR_DS: // decimal string
|
|
630 case EVR_IS: // integer string
|
|
631 case EVR_AS: // age string
|
|
632 case EVR_DA: // date string
|
|
633 case EVR_DT: // date time string
|
|
634 case EVR_TM: // time string
|
|
635 case EVR_AE: // application entity title
|
|
636 case EVR_CS: // code string
|
|
637 case EVR_SH: // short string
|
|
638 case EVR_LO: // long string
|
|
639 case EVR_ST: // short text
|
|
640 case EVR_LT: // long text
|
|
641 case EVR_UT: // unlimited text
|
|
642 case EVR_PN: // person name
|
|
643 case EVR_UI: // unique identifier
|
|
644 {
|
|
645 ok = element.putString(value.c_str()).good();
|
|
646 break;
|
|
647 }
|
|
648
|
|
649
|
|
650 /**
|
|
651 * Numerical types
|
|
652 **/
|
|
653
|
|
654 case EVR_SL: // signed long
|
|
655 {
|
|
656 ok = element.putSint32(boost::lexical_cast<Sint32>(value)).good();
|
|
657 break;
|
|
658 }
|
|
659
|
|
660 case EVR_SS: // signed short
|
|
661 {
|
|
662 ok = element.putSint16(boost::lexical_cast<Sint16>(value)).good();
|
|
663 break;
|
|
664 }
|
|
665
|
|
666 case EVR_UL: // unsigned long
|
|
667 {
|
|
668 ok = element.putUint32(boost::lexical_cast<Uint32>(value)).good();
|
|
669 break;
|
|
670 }
|
|
671
|
|
672 case EVR_US: // unsigned short
|
|
673 {
|
|
674 ok = element.putUint16(boost::lexical_cast<Uint16>(value)).good();
|
|
675 break;
|
|
676 }
|
|
677
|
|
678 case EVR_FL: // float single-precision
|
|
679 {
|
|
680 ok = element.putFloat32(boost::lexical_cast<float>(value)).good();
|
|
681 break;
|
|
682 }
|
|
683
|
|
684 case EVR_FD: // float double-precision
|
|
685 {
|
|
686 ok = element.putFloat64(boost::lexical_cast<double>(value)).good();
|
|
687 break;
|
|
688 }
|
|
689
|
|
690
|
|
691 /**
|
|
692 * Sequence types, should never occur at this point.
|
|
693 **/
|
|
694
|
|
695 case EVR_SQ: // sequence of items
|
|
696 {
|
|
697 ok = false;
|
|
698 break;
|
|
699 }
|
|
700
|
|
701
|
|
702 /**
|
|
703 * Internal to DCMTK.
|
|
704 **/
|
|
705
|
|
706 case EVR_ox: // OB or OW depending on context
|
|
707 case EVR_xs: // SS or US depending on context
|
|
708 case EVR_lt: // US, SS or OW depending on context, used for LUT Data (thus the name)
|
|
709 case EVR_na: // na="not applicable", for data which has no VR
|
|
710 case EVR_up: // up="unsigned pointer", used internally for DICOMDIR suppor
|
|
711 case EVR_item: // used internally for items
|
|
712 case EVR_metainfo: // used internally for meta info datasets
|
|
713 case EVR_dataset: // used internally for datasets
|
|
714 case EVR_fileFormat: // used internally for DICOM files
|
|
715 case EVR_dicomDir: // used internally for DICOMDIR objects
|
|
716 case EVR_dirRecord: // used internally for DICOMDIR records
|
|
717 case EVR_pixelSQ: // used internally for pixel sequences in a compressed image
|
|
718 case EVR_pixelItem: // used internally for pixel items in a compressed image
|
|
719 case EVR_UNKNOWN: // used internally for elements with unknown VR (encoded with 4-byte length field in explicit VR)
|
|
720 case EVR_PixelData: // used internally for uncompressed pixeld data
|
|
721 case EVR_OverlayData: // used internally for overlay data
|
|
722 case EVR_UNKNOWN2B: // used internally for elements with unknown VR with 2-byte length field in explicit VR
|
|
723 default:
|
|
724 break;
|
|
725 }
|
|
726 }
|
|
727 catch (boost::bad_lexical_cast&)
|
|
728 {
|
|
729 ok = false;
|
|
730 }
|
|
731
|
|
732 if (!ok)
|
|
733 {
|
|
734 throw OrthancException(ErrorCode_InternalError);
|
|
735 }
|
|
736 }
|
|
737
|
|
738
|
|
739 void ParsedDicomFile::Remove(const DicomTag& tag)
|
|
740 {
|
|
741 DcmTagKey key(tag.GetGroup(), tag.GetElement());
|
794
|
742 DcmElement* element = pimpl_->file_->getDataset()->remove(key);
|
790
|
743 if (element != NULL)
|
|
744 {
|
|
745 delete element;
|
|
746 }
|
|
747 }
|
|
748
|
|
749
|
|
750
|
|
751 void ParsedDicomFile::RemovePrivateTags()
|
|
752 {
|
|
753 typedef std::list<DcmElement*> Tags;
|
|
754
|
|
755 Tags privateTags;
|
|
756
|
794
|
757 DcmDataset& dataset = *pimpl_->file_->getDataset();
|
790
|
758 for (unsigned long i = 0; i < dataset.card(); i++)
|
|
759 {
|
|
760 DcmElement* element = dataset.getElement(i);
|
|
761 DcmTag tag(element->getTag());
|
|
762 if (!strcmp("PrivateCreator", tag.getTagName()) || // TODO - This may change with future versions of DCMTK
|
|
763 tag.getPrivateCreator() != NULL)
|
|
764 {
|
|
765 privateTags.push_back(element);
|
|
766 }
|
|
767 }
|
|
768
|
|
769 for (Tags::iterator it = privateTags.begin();
|
|
770 it != privateTags.end(); ++it)
|
|
771 {
|
|
772 DcmElement* tmp = dataset.remove(*it);
|
|
773 if (tmp != NULL)
|
|
774 {
|
|
775 delete tmp;
|
|
776 }
|
|
777 }
|
|
778 }
|
|
779
|
|
780
|
|
781
|
|
782 void ParsedDicomFile::Insert(const DicomTag& tag,
|
|
783 const std::string& value)
|
|
784 {
|
|
785 std::auto_ptr<DcmElement> element(CreateElementForTag(tag));
|
|
786 FillElementWithString(*element, tag, value);
|
|
787
|
794
|
788 if (!pimpl_->file_->getDataset()->insert(element.release(), false, false).good())
|
790
|
789 {
|
|
790 // This field already exists
|
|
791 throw OrthancException(ErrorCode_InternalError);
|
|
792 }
|
|
793 }
|
|
794
|
|
795
|
|
796 void ParsedDicomFile::Replace(const DicomTag& tag,
|
|
797 const std::string& value,
|
|
798 DicomReplaceMode mode)
|
|
799 {
|
|
800 DcmTagKey key(tag.GetGroup(), tag.GetElement());
|
|
801 DcmElement* element = NULL;
|
|
802
|
794
|
803 if (!pimpl_->file_->getDataset()->findAndGetElement(key, element).good() ||
|
790
|
804 element == NULL)
|
|
805 {
|
|
806 // This field does not exist, act wrt. the specified "mode"
|
|
807 switch (mode)
|
|
808 {
|
|
809 case DicomReplaceMode_InsertIfAbsent:
|
|
810 Insert(tag, value);
|
|
811 break;
|
|
812
|
|
813 case DicomReplaceMode_ThrowIfAbsent:
|
|
814 throw OrthancException(ErrorCode_InexistentItem);
|
|
815
|
|
816 case DicomReplaceMode_IgnoreIfAbsent:
|
|
817 return;
|
|
818 }
|
|
819 }
|
|
820 else
|
|
821 {
|
|
822 FillElementWithString(*element, tag, value);
|
|
823 }
|
|
824
|
|
825
|
|
826 /**
|
|
827 * dcmodify will automatically correct 'Media Storage SOP Class
|
|
828 * UID' and 'Media Storage SOP Instance UID' in the metaheader, if
|
|
829 * you make changes to the related tags in the dataset ('SOP Class
|
|
830 * UID' and 'SOP Instance UID') via insert or modify mode
|
|
831 * options. You can disable this behaviour by using the -nmu
|
|
832 * option.
|
|
833 **/
|
|
834
|
|
835 if (tag == DICOM_TAG_SOP_CLASS_UID)
|
|
836 {
|
|
837 Replace(DICOM_TAG_MEDIA_STORAGE_SOP_CLASS_UID, value, DicomReplaceMode_InsertIfAbsent);
|
|
838 }
|
|
839
|
|
840 if (tag == DICOM_TAG_SOP_INSTANCE_UID)
|
|
841 {
|
|
842 Replace(DICOM_TAG_MEDIA_STORAGE_SOP_INSTANCE_UID, value, DicomReplaceMode_InsertIfAbsent);
|
|
843 }
|
|
844 }
|
|
845
|
|
846
|
|
847 void ParsedDicomFile::Answer(RestApiOutput& output)
|
|
848 {
|
|
849 std::string serialized;
|
794
|
850 if (FromDcmtkBridge::SaveToMemoryBuffer(serialized, pimpl_->file_->getDataset()))
|
790
|
851 {
|
|
852 output.AnswerBuffer(serialized, CONTENT_TYPE_OCTET_STREAM);
|
|
853 }
|
|
854 }
|
|
855
|
|
856
|
|
857
|
|
858 bool ParsedDicomFile::GetTagValue(std::string& value,
|
|
859 const DicomTag& tag)
|
|
860 {
|
|
861 DcmTagKey k(tag.GetGroup(), tag.GetElement());
|
794
|
862 DcmDataset& dataset = *pimpl_->file_->getDataset();
|
790
|
863 DcmElement* element = NULL;
|
|
864 if (!dataset.findAndGetElement(k, element).good() ||
|
|
865 element == NULL)
|
|
866 {
|
|
867 return false;
|
|
868 }
|
|
869
|
|
870 std::auto_ptr<DicomValue> v(FromDcmtkBridge::ConvertLeafElement(*element));
|
|
871
|
|
872 if (v.get() == NULL)
|
|
873 {
|
|
874 value = "";
|
|
875 }
|
|
876 else
|
|
877 {
|
|
878 value = v->AsString();
|
|
879 }
|
|
880
|
|
881 return true;
|
|
882 }
|
|
883
|
|
884
|
|
885
|
|
886 DicomInstanceHasher ParsedDicomFile::GetHasher()
|
|
887 {
|
|
888 std::string patientId, studyUid, seriesUid, instanceUid;
|
|
889
|
|
890 if (!GetTagValue(patientId, DICOM_TAG_PATIENT_ID) ||
|
|
891 !GetTagValue(studyUid, DICOM_TAG_STUDY_INSTANCE_UID) ||
|
|
892 !GetTagValue(seriesUid, DICOM_TAG_SERIES_INSTANCE_UID) ||
|
|
893 !GetTagValue(instanceUid, DICOM_TAG_SOP_INSTANCE_UID))
|
|
894 {
|
|
895 throw OrthancException(ErrorCode_BadFileFormat);
|
|
896 }
|
|
897
|
|
898 return DicomInstanceHasher(patientId, studyUid, seriesUid, instanceUid);
|
|
899 }
|
|
900
|
|
901
|
|
902 static void StoreElement(Json::Value& target,
|
|
903 DcmElement& element,
|
|
904 unsigned int maxStringLength);
|
|
905
|
|
906 static void StoreItem(Json::Value& target,
|
|
907 DcmItem& item,
|
|
908 unsigned int maxStringLength)
|
|
909 {
|
|
910 target = Json::Value(Json::objectValue);
|
|
911
|
|
912 for (unsigned long i = 0; i < item.card(); i++)
|
|
913 {
|
|
914 DcmElement* element = item.getElement(i);
|
|
915 StoreElement(target, *element, maxStringLength);
|
|
916 }
|
|
917 }
|
|
918
|
|
919
|
|
920 static void StoreElement(Json::Value& target,
|
|
921 DcmElement& element,
|
|
922 unsigned int maxStringLength)
|
|
923 {
|
|
924 assert(target.type() == Json::objectValue);
|
|
925
|
|
926 DicomTag tag(FromDcmtkBridge::GetTag(element));
|
|
927 const std::string formattedTag = tag.Format();
|
|
928
|
|
929 #if 0
|
|
930 const std::string tagName = FromDcmtkBridge::GetName(tag);
|
|
931 #else
|
|
932 // This version of the code gives access to the name of the private tags
|
|
933 DcmTag tagbis(element.getTag());
|
|
934 const std::string tagName(tagbis.getTagName());
|
|
935 #endif
|
|
936
|
|
937 if (element.isLeaf())
|
|
938 {
|
|
939 Json::Value value(Json::objectValue);
|
|
940 value["Name"] = tagName;
|
|
941
|
|
942 if (tagbis.getPrivateCreator() != NULL)
|
|
943 {
|
|
944 value["PrivateCreator"] = tagbis.getPrivateCreator();
|
|
945 }
|
|
946
|
|
947 std::auto_ptr<DicomValue> v(FromDcmtkBridge::ConvertLeafElement(element));
|
|
948 if (v->IsNull())
|
|
949 {
|
|
950 value["Type"] = "Null";
|
|
951 value["Value"] = Json::nullValue;
|
|
952 }
|
|
953 else
|
|
954 {
|
|
955 std::string s = v->AsString();
|
|
956 if (maxStringLength == 0 ||
|
|
957 s.size() <= maxStringLength)
|
|
958 {
|
|
959 value["Type"] = "String";
|
|
960 value["Value"] = s;
|
|
961 }
|
|
962 else
|
|
963 {
|
|
964 value["Type"] = "TooLong";
|
|
965 value["Value"] = Json::nullValue;
|
|
966 }
|
|
967 }
|
|
968
|
|
969 target[formattedTag] = value;
|
|
970 }
|
|
971 else
|
|
972 {
|
|
973 Json::Value children(Json::arrayValue);
|
|
974
|
|
975 // "All subclasses of DcmElement except for DcmSequenceOfItems
|
|
976 // are leaf nodes, while DcmSequenceOfItems, DcmItem, DcmDataset
|
|
977 // etc. are not." The following cast is thus OK.
|
|
978 DcmSequenceOfItems& sequence = dynamic_cast<DcmSequenceOfItems&>(element);
|
|
979
|
|
980 for (unsigned long i = 0; i < sequence.card(); i++)
|
|
981 {
|
|
982 DcmItem* child = sequence.getItem(i);
|
|
983 Json::Value& v = children.append(Json::objectValue);
|
|
984 StoreItem(v, *child, maxStringLength);
|
|
985 }
|
|
986
|
|
987 target[formattedTag]["Name"] = tagName;
|
|
988 target[formattedTag]["Type"] = "Sequence";
|
|
989 target[formattedTag]["Value"] = children;
|
|
990 }
|
|
991 }
|
|
992
|
|
993
|
|
994 template <typename T>
|
|
995 static void ExtractPngImageTruncate(std::string& result,
|
|
996 DicomIntegerPixelAccessor& accessor,
|
|
997 PixelFormat format)
|
|
998 {
|
|
999 assert(accessor.GetChannelCount() == 1);
|
|
1000
|
|
1001 PngWriter w;
|
|
1002
|
|
1003 std::vector<T> image(accessor.GetWidth() * accessor.GetHeight(), 0);
|
|
1004 T* pixel = &image[0];
|
|
1005 for (unsigned int y = 0; y < accessor.GetHeight(); y++)
|
|
1006 {
|
|
1007 for (unsigned int x = 0; x < accessor.GetWidth(); x++, pixel++)
|
|
1008 {
|
|
1009 int32_t v = accessor.GetValue(x, y);
|
|
1010 if (v < static_cast<int32_t>(std::numeric_limits<T>::min()))
|
|
1011 *pixel = std::numeric_limits<T>::min();
|
|
1012 else if (v > static_cast<int32_t>(std::numeric_limits<T>::max()))
|
|
1013 *pixel = std::numeric_limits<T>::max();
|
|
1014 else
|
|
1015 *pixel = static_cast<T>(v);
|
|
1016 }
|
|
1017 }
|
|
1018
|
|
1019 w.WriteToMemory(result, accessor.GetWidth(), accessor.GetHeight(),
|
|
1020 accessor.GetWidth() * sizeof(T), format, &image[0]);
|
|
1021 }
|
|
1022
|
|
1023
|
|
1024 void ParsedDicomFile::SaveToMemoryBuffer(std::string& buffer)
|
|
1025 {
|
794
|
1026 FromDcmtkBridge::SaveToMemoryBuffer(buffer, pimpl_->file_->getDataset());
|
790
|
1027 }
|
|
1028
|
|
1029
|
|
1030 void ParsedDicomFile::SaveToFile(const std::string& path)
|
|
1031 {
|
|
1032 // TODO Avoid using a temporary memory buffer, write directly on disk
|
|
1033 std::string content;
|
|
1034 SaveToMemoryBuffer(content);
|
|
1035 Toolbox::WriteFile(content, path);
|
|
1036 }
|
|
1037
|
|
1038
|
794
|
1039 ParsedDicomFile::ParsedDicomFile() : pimpl_(new PImpl)
|
790
|
1040 {
|
794
|
1041 pimpl_->file_.reset(new DcmFileFormat);
|
790
|
1042 Replace(DICOM_TAG_PATIENT_ID, FromDcmtkBridge::GenerateUniqueIdentifier(ResourceType_Patient));
|
|
1043 Replace(DICOM_TAG_STUDY_INSTANCE_UID, FromDcmtkBridge::GenerateUniqueIdentifier(ResourceType_Study));
|
|
1044 Replace(DICOM_TAG_SERIES_INSTANCE_UID, FromDcmtkBridge::GenerateUniqueIdentifier(ResourceType_Series));
|
|
1045 Replace(DICOM_TAG_SOP_INSTANCE_UID, FromDcmtkBridge::GenerateUniqueIdentifier(ResourceType_Instance));
|
|
1046 }
|
|
1047
|
791
|
1048
|
794
|
1049 ParsedDicomFile::ParsedDicomFile(const char* content, size_t size) : pimpl_(new PImpl)
|
791
|
1050 {
|
|
1051 Setup(content, size);
|
|
1052 }
|
|
1053
|
794
|
1054 ParsedDicomFile::ParsedDicomFile(const std::string& content) : pimpl_(new PImpl)
|
791
|
1055 {
|
|
1056 if (content.size() == 0)
|
|
1057 {
|
|
1058 Setup(NULL, 0);
|
|
1059 }
|
|
1060 else
|
|
1061 {
|
|
1062 Setup(&content[0], content.size());
|
|
1063 }
|
|
1064 }
|
792
|
1065
|
|
1066
|
794
|
1067 ParsedDicomFile::ParsedDicomFile(ParsedDicomFile& other) :
|
|
1068 pimpl_(new PImpl)
|
792
|
1069 {
|
794
|
1070 pimpl_->file_.reset(dynamic_cast<DcmFileFormat*>(other.pimpl_->file_->clone()));
|
792
|
1071 }
|
|
1072
|
|
1073
|
|
1074 ParsedDicomFile::~ParsedDicomFile()
|
|
1075 {
|
794
|
1076 delete pimpl_;
|
792
|
1077 }
|
|
1078
|
793
|
1079
|
|
1080 void* ParsedDicomFile::GetDcmtkObject()
|
|
1081 {
|
794
|
1082 return pimpl_->file_.get();
|
793
|
1083 }
|
|
1084
|
|
1085
|
|
1086 ParsedDicomFile* ParsedDicomFile::Clone()
|
|
1087 {
|
794
|
1088 return new ParsedDicomFile(*this);
|
793
|
1089 }
|
790
|
1090 }
|