1 /*=========================================================================
4 Module: $RCSfile: gdcmFileHelper.cxx,v $
7 Date: $Date: 2005/04/27 09:52:28 $
8 Version: $Revision: 1.35 $
10 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
11 l'Image). All rights reserved. See Doc/License.txt or
12 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
14 This software is distributed WITHOUT ANY WARRANTY; without even
15 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
16 PURPOSE. See the above copyright notices for more information.
18 =========================================================================*/
20 #include "gdcmFileHelper.h"
21 #include "gdcmGlobal.h"
23 #include "gdcmDocument.h"
24 #include "gdcmDebug.h"
26 #include "gdcmBinEntry.h"
27 #include "gdcmValEntry.h"
28 #include "gdcmSeqEntry.h"
29 #include "gdcmSQItem.h"
30 #include "gdcmContentEntry.h"
32 #include "gdcmPixelReadConvert.h"
33 #include "gdcmPixelWriteConvert.h"
34 #include "gdcmDocEntryArchive.h"
35 #include "gdcmDictSet.h"
40 // ----------------------------- WARNING -------------------------
42 These lines will be moved to the document-to-be 'User's Guide'
44 // To read an image, user needs a gdcm::File
45 gdcm::File *f1 = new gdcm::File(fileName);
46 // user can now check some values
47 std::string v = f1->GetEntryValue(groupNb,ElementNb);
48 // to get the pixels, user needs a gdcm::FileHelper
49 gdcm::FileHelper *fh1 = new gdcm::FileHelper(f1);
50 // user may ask not to convert Palette to RGB
51 uint8_t *pixels = fh1->GetImageDataRaw();
52 int imageLength = fh1->GetImageDataRawSize();
53 // He can now use the pixels, create a new image, ...
54 uint8_t *userPixels = ...
56 To re-write the image, user re-uses the gdcm::FileHelper
58 fh1->SetImageData( userPixels, userPixelsLength);
59 fh1->SetTypeToRaw(); // Even if it was possible to convert Palette to RGB
62 fh1->SetWriteTypeToDcmExpl(); // he wants Explicit Value Representation
63 // Little Endian is the default
64 // no other value is allowed
65 (-->SetWriteType(ExplicitVR);)
66 -->WriteType = ExplicitVR;
67 fh1->Write(newFileName); // overwrites the file, if any
70 fh1->WriteDcmExplVR(newFileName);
73 // ----------------------------- WARNING -------------------------
75 These lines will be moved to the document-to-be 'Developer's Guide'
77 WriteMode : WMODE_RAW / WMODE_RGB
78 WriteType : ImplicitVR, ExplicitVR, ACR, ACR_LIBIDO
80 fh1->Write(newFileName);
81 SetWriteFileTypeToImplicitVR() / SetWriteFileTypeToExplicitVR();
82 (modifies TransferSyntax)
83 SetWriteToRaw(); / SetWriteToRGB();
84 (modifies, when necessary : photochromatic interpretation,
85 samples per pixel, Planar configuration,
86 bits allocated, bits stored, high bit -ACR 24 bits-
87 Pixels element VR, pushes out the LUT )
88 CheckWriteIntegrity();
89 (checks user given pixels length)
90 FileInternal->Write(fileName,WriteType)
91 fp = opens file(fileName);
92 ComputeGroup0002Length(writetype);
94 RemoveEntryNoDestroy(palettes, etc)
95 Document::WriteContent(fp, writetype);
97 (moves back to the File all the archived elements)
98 RestoreWriteFileType();
99 (pushes back group 0002, with TransferSyntax)
107 //-------------------------------------------------------------------------
108 // Constructor / Destructor
110 * \brief Constructor dedicated to deal with the *pixels* area of a ACR/DICOMV3
111 * file (gdcm::File only deals with the ... header)
112 * Opens (in read only and when possible) an existing file and checks
113 * for DICOM compliance. Returns NULL on failure.
114 * It will be up to the user to load the pixels into memory
115 * \note the in-memory representation of all available tags found in
116 * the DICOM header is post-poned to first header information access.
117 * This avoid a double parsing of public part of the header when
118 * one sets an a posteriori shadow dictionary (efficiency can be
119 * seen as a side effect).
121 FileHelper::FileHelper( )
123 FileInternal = new File( );
129 * \brief Constructor dedicated to deal with the *pixels* area of a ACR/DICOMV3
130 * file (File only deals with the ... header)
131 * Opens (in read only and when possible) an existing file and checks
132 * for DICOM compliance. Returns NULL on failure.
133 * It will be up to the user to load the pixels into memory
134 * \note the in-memory representation of all available tags found in
135 * the DICOM header is post-poned to first header information access.
136 * This avoid a double parsing of public part of the header when
137 * user sets an a posteriori shadow dictionary (efficiency can be
138 * seen as a side effect).
139 * @param header already built Header
141 FileHelper::FileHelper(File *header)
143 FileInternal = header;
149 * \brief Constructor dedicated to deal with the *pixels* area of a ACR/DICOMV3
150 * file (gdcm::File only deals with the ... header)
151 * Opens (in read only and when possible) an existing file and checks
152 * for DICOM compliance. Returns NULL on failure.
153 * It will be up to the user to load the pixels into memory
154 * \note the in-memory representation of all available tags found in
155 * the DICOM header is post-poned to first header information access.
156 * This avoid a double parsing of public part of the header when
157 * one sets an a posteriori shadow dictionary (efficiency can be
158 * seen as a side effect).
159 * @param filename file to be opened for parsing
161 FileHelper::FileHelper(std::string const &filename )
163 FileInternal = new File( filename );
169 * \brief canonical destructor
170 * \note If the header (gdcm::File) was created by the FileHelper constructor,
171 * it is destroyed by the FileHelper
173 FileHelper::~FileHelper()
175 if( PixelReadConverter )
177 delete PixelReadConverter;
179 if( PixelWriteConverter )
181 delete PixelWriteConverter;
195 //-----------------------------------------------------------------------------
198 * \brief Accesses an existing DocEntry (i.e. a Dicom Element)
199 * through it's (group, element) and modifies it's content with
201 * @param content new value (string) to substitute with
202 * @param group group number of the Dicom Element to modify
203 * @param elem element number of the Dicom Element to modify
204 * \return false if DocEntry not found
206 bool FileHelper::SetValEntry(std::string const &content,
207 uint16_t group, uint16_t elem)
209 return FileInternal->SetValEntry(content, group, elem);
214 * \brief Accesses an existing DocEntry (i.e. a Dicom Element)
215 * through it's (group, element) and modifies it's content with
217 * @param content new value (void* -> uint8_t*) to substitute with
218 * @param lgth new value length
219 * @param group group number of the Dicom Element to modify
220 * @param elem element number of the Dicom Element to modify
221 * \return false if DocEntry not found
223 bool FileHelper::SetBinEntry(uint8_t *content, int lgth,
224 uint16_t group, uint16_t elem)
226 return FileInternal->SetBinEntry(content, lgth, group, elem);
230 * \brief Modifies the value of a given DocEntry (Dicom entry)
231 * when it exists. Creates it with the given value when unexistant.
232 * @param content (string)value to be set
233 * @param group Group number of the Entry
234 * @param elem Element number of the Entry
235 * \return pointer to the modified/created Dicom entry (NULL when creation
238 ValEntry *FileHelper::InsertValEntry(std::string const &content,
239 uint16_t group, uint16_t elem)
241 return FileInternal->InsertValEntry(content,group,elem);
245 * \brief Modifies the value of a given DocEntry (Dicom entry)
246 * when it exists. Creates it with the given value when unexistant.
247 * A copy of the binArea is made to be kept in the Document.
248 * @param binArea (binary)value to be set
249 * @param lgth new value length
250 * @param group Group number of the Entry
251 * @param elem Element number of the Entry
252 * \return pointer to the modified/created Dicom entry (NULL when creation
255 BinEntry *FileHelper::InsertBinEntry(uint8_t *binArea, int lgth,
256 uint16_t group, uint16_t elem)
258 return FileInternal->InsertBinEntry(binArea, lgth, group, elem);
262 * \brief Modifies the value of a given DocEntry (Dicom entry)
263 * when it exists. Creates it, empty (?!) when unexistant.
264 * @param group Group number of the Entry
265 * @param elem Element number of the Entry
266 * \return pointer to the modified/created Dicom entry (NULL when creation
269 SeqEntry *FileHelper::InsertSeqEntry(uint16_t group, uint16_t elem)
271 return FileInternal->InsertSeqEntry(group, elem);
275 * \brief Get the size of the image data
276 * If the image can be RGB (with a lut or by default), the size
277 * corresponds to the RGB image
278 * (use GetImageDataRawSize if you want to be sure to get *only*
279 * the size of the pixels)
280 * @return The image size
282 size_t FileHelper::GetImageDataSize()
284 if ( PixelWriteConverter->GetUserData() )
286 return PixelWriteConverter->GetUserDataSize();
289 return PixelReadConverter->GetRGBSize();
293 * \brief Get the size of the image data
294 * If the image could be converted to RGB using a LUT,
295 * this transformation is not taken into account by GetImageDataRawSize
296 * (use GetImageDataSize if you wish)
297 * @return The raw image size
299 size_t FileHelper::GetImageDataRawSize()
301 if ( PixelWriteConverter->GetUserData() )
303 return PixelWriteConverter->GetUserDataSize();
306 return PixelReadConverter->GetRawSize();
310 * \brief - Allocates necessary memory,
311 * - Reads the pixels from disk (uncompress if necessary),
312 * - Transforms YBR pixels, if any, into RGB pixels,
313 * - Transforms 3 planes R, G, B, if any, into a single RGB Plane
314 * - Transforms single Grey plane + 3 Palettes into a RGB Plane
315 * - Copies the pixel data (image[s]/volume[s]) to newly allocated zone.
316 * @return Pointer to newly allocated pixel data.
317 * NULL if alloc fails
319 uint8_t *FileHelper::GetImageData()
321 if ( PixelWriteConverter->GetUserData() )
323 return PixelWriteConverter->GetUserData();
328 // If the decompression failed nothing can be done.
332 if ( FileInternal->HasLUT() && PixelReadConverter->BuildRGBImage() )
334 return PixelReadConverter->GetRGB();
338 // When no LUT or LUT conversion fails, return the Raw
339 return PixelReadConverter->GetRaw();
344 * \brief Allocates necessary memory,
345 * Transforms YBR pixels (if any) into RGB pixels
346 * Transforms 3 planes R, G, B (if any) into a single RGB Plane
347 * Copies the pixel data (image[s]/volume[s]) to newly allocated zone.
348 * DOES NOT transform Grey plane + 3 Palettes into a RGB Plane
349 * @return Pointer to newly allocated pixel data.
350 * NULL if alloc fails
352 uint8_t *FileHelper::GetImageDataRaw ()
359 * Read the pixels from disk (uncompress if necessary),
360 * Transforms YBR pixels, if any, into RGB pixels
361 * Transforms 3 planes R, G, B, if any, into a single RGB Plane
362 * Transforms single Grey plane + 3 Palettes into a RGB Plane
363 * Copies at most MaxSize bytes of pixel data to caller allocated
365 * \warning This function allows people that want to build a volume
366 * from an image stack *not to* have, first to get the image pixels,
367 * and then move them to the volume area.
368 * It's absolutely useless for any VTK user since vtk chooses
369 * to invert the lines of an image, that is the last line comes first
370 * (for some axis related reasons?). Hence he will have
371 * to load the image line by line, starting from the end.
372 * VTK users have to call GetImageData
374 * @param destination Address (in caller's memory space) at which the
375 * pixel data should be copied
376 * @param maxSize Maximum number of bytes to be copied. When MaxSize
377 * is not sufficient to hold the pixel data the copy is not
378 * executed (i.e. no partial copy).
379 * @return On success, the number of bytes actually copied. Zero on
380 * failure e.g. MaxSize is lower than necessary.
382 size_t FileHelper::GetImageDataIntoVector (void *destination, size_t maxSize)
386 // If the decompression failed nothing can be done.
390 if ( FileInternal->HasLUT() && PixelReadConverter->BuildRGBImage() )
392 if ( PixelReadConverter->GetRGBSize() > maxSize )
394 gdcmWarningMacro( "Pixel data bigger than caller's expected MaxSize");
398 (void*)PixelReadConverter->GetRGB(),
399 PixelReadConverter->GetRGBSize() );
400 return PixelReadConverter->GetRGBSize();
403 // Either no LUT conversion necessary or LUT conversion failed
404 if ( PixelReadConverter->GetRawSize() > maxSize )
406 gdcmWarningMacro( "Pixel data bigger than caller's expected MaxSize");
410 (void*)PixelReadConverter->GetRaw(),
411 PixelReadConverter->GetRawSize() );
412 return PixelReadConverter->GetRawSize();
416 * \brief Points the internal pointer to the callers inData
417 * image representation, BUT WITHOUT COPYING THE DATA.
418 * 'image' Pixels are presented as C-like 2D arrays : line per line.
419 * 'volume'Pixels are presented as C-like 3D arrays : plane per plane
420 * \warning Since the pixels are not copied, it is the caller's responsability
421 * not to deallocate its data before gdcm uses them (e.g. with
422 * the Write() method )
423 * @param inData user supplied pixel area (uint8_t* is just for the compiler.
424 * user is allowed to pass any kind of pixelsn since the size is
426 * @param expectedSize total image size, *in Bytes*
428 void FileHelper::SetImageData(uint8_t *inData, size_t expectedSize)
430 SetUserData(inData, expectedSize);
434 * \brief Set the image data defined by the user
435 * \warning When writting the file, this data are get as default data to write
436 * @param inData user supplied pixel area (uint8_t* is just for the compiler.
437 * user is allowed to pass any kind of pixels since the size is
439 * @param expectedSize total image size, *in Bytes*
441 void FileHelper::SetUserData(uint8_t *inData, size_t expectedSize)
443 PixelWriteConverter->SetUserData(inData, expectedSize);
447 * \brief Get the image data defined by the user
448 * \warning When writting the file, this data are get as default data to write
450 uint8_t *FileHelper::GetUserData()
452 return PixelWriteConverter->GetUserData();
456 * \brief Get the image data size defined by the user
457 * \warning When writting the file, this data are get as default data to write
459 size_t FileHelper::GetUserDataSize()
461 return PixelWriteConverter->GetUserDataSize();
465 * \brief Get the image data from the file.
466 * If a LUT is found, the data are expanded to be RGB
468 uint8_t *FileHelper::GetRGBData()
470 return PixelReadConverter->GetRGB();
474 * \brief Get the image data size from the file.
475 * If a LUT is found, the data are expanded to be RGB
477 size_t FileHelper::GetRGBDataSize()
479 return PixelReadConverter->GetRGBSize();
483 * \brief Get the image data from the file.
484 * Even when a LUT is found, the data are not expanded to RGB!
486 uint8_t *FileHelper::GetRawData()
488 return PixelReadConverter->GetRaw();
492 * \brief Get the image data size from the file.
493 * Even when a LUT is found, the data are not expanded to RGB!
495 size_t FileHelper::GetRawDataSize()
497 return PixelReadConverter->GetRawSize();
501 * \brief Access to the underlying \ref PixelReadConverter RGBA LUT
503 uint8_t* FileHelper::GetLutRGBA()
505 return PixelReadConverter->GetLutRGBA();
509 * \brief Writes on disk A SINGLE Dicom file
510 * NO test is performed on processor "Endiannity".
511 * It's up to the user to call his Reader properly
512 * @param fileName name of the file to be created
513 * (any already existing file is over written)
514 * @return false if write fails
516 bool FileHelper::WriteRawData(std::string const &fileName)
518 std::ofstream fp1(fileName.c_str(), std::ios::out | std::ios::binary );
521 gdcmWarningMacro( "Fail to open (write) file:" << fileName.c_str());
525 if( PixelWriteConverter->GetUserData() )
527 fp1.write( (char*)PixelWriteConverter->GetUserData(),
528 PixelWriteConverter->GetUserDataSize() );
530 else if ( PixelReadConverter->GetRGB() )
532 fp1.write( (char*)PixelReadConverter->GetRGB(),
533 PixelReadConverter->GetRGBSize());
535 else if ( PixelReadConverter->GetRaw() )
537 fp1.write( (char*)PixelReadConverter->GetRaw(),
538 PixelReadConverter->GetRawSize());
542 gdcmErrorMacro( "Nothing written." );
551 * \brief Writes on disk A SINGLE Dicom file,
552 * using the Implicit Value Representation convention
553 * NO test is performed on processor "Endiannity".
554 * @param fileName name of the file to be created
555 * (any already existing file is overwritten)
556 * @return false if write fails
559 bool FileHelper::WriteDcmImplVR (std::string const &fileName)
561 SetWriteTypeToDcmImplVR();
562 return Write(fileName);
566 * \brief Writes on disk A SINGLE Dicom file,
567 * using the Explicit Value Representation convention
568 * NO test is performed on processor "Endiannity".
569 * @param fileName name of the file to be created
570 * (any already existing file is overwritten)
571 * @return false if write fails
574 bool FileHelper::WriteDcmExplVR (std::string const &fileName)
576 SetWriteTypeToDcmExplVR();
577 return Write(fileName);
581 * \brief Writes on disk A SINGLE Dicom file,
582 * using the ACR-NEMA convention
583 * NO test is performed on processor "Endiannity".
584 * (a l'attention des logiciels cliniques
585 * qui ne prennent en entrée QUE des images ACR ...
586 * \warning if a DICOM_V3 header is supplied,
587 * groups < 0x0008 and shadow groups are ignored
588 * \warning NO TEST is performed on processor "Endiannity".
589 * @param fileName name of the file to be created
590 * (any already existing file is overwritten)
591 * @return false if write fails
594 bool FileHelper::WriteAcr (std::string const &fileName)
597 return Write(fileName);
601 * \brief Writes on disk A SINGLE Dicom file,
602 * @param fileName name of the file to be created
603 * (any already existing file is overwritten)
604 * @return false if write fails
606 bool FileHelper::Write(std::string const &fileName)
611 SetWriteFileTypeToImplicitVR();
612 CheckMandatoryElements();
614 case Unknown: // should never happen; ExplicitVR is the default value
616 SetWriteFileTypeToExplicitVR();
617 CheckMandatoryElements();
621 // Just to avoid further trouble if user create a file ex-nihilo,
622 // wants to write it as an ACR-NEMA file,
623 // and forget to create any Entry belonging to group 0008
625 // We add Recognition Code (RET)
626 if ( ! FileInternal->GetValEntry(0x0008, 0x0010) )
627 FileInternal->InsertValEntry("", 0x0008, 0x0010);
628 SetWriteFileTypeToACR();
629 SetWriteFileTypeToImplicitVR();
630 CheckMandatoryElements();
635 // --------------------------------------------------------------
636 // Special Patch to allow gdcm to re-write ACR-LibIDO formated images
638 // if recognition code tells us we dealt with a LibIDO image
639 // we reproduce on disk the switch between lineNumber and columnNumber
640 // just before writting ...
641 /// \todo the best trick would be *change* the recognition code
642 /// but pb expected if user deals with, e.g. COMPLEX images
644 if( WriteType == ACR_LIBIDO )
650 SetWriteToNoLibido();
652 // ----------------- End of Special Patch ----------------
657 SetWriteToRaw(); // modifies and pushes to the archive, when necessary
660 SetWriteToRGB(); // modifies and pushes to the archive, when necessary
664 bool check = CheckWriteIntegrity(); // verifies length
667 check = FileInternal->Write(fileName,WriteType);
671 RestoreWriteFileType();
672 RestoreWriteMandatory();
674 // --------------------------------------------------------------
675 // Special Patch to allow gdcm to re-write ACR-LibIDO formated images
677 // ...and we restore the header to be Dicom Compliant again
678 // just after writting
679 RestoreWriteOfLibido();
680 // ----------------- End of Special Patch ----------------
685 //-----------------------------------------------------------------------------
688 * \brief Checks the write integrity
690 * The tests made are :
691 * - verify the size of the image to write with the possible write
692 * when the user set an image data
693 * @return true if check is successfull
695 bool FileHelper::CheckWriteIntegrity()
697 if(PixelWriteConverter->GetUserData())
699 int numberBitsAllocated = FileInternal->GetBitsAllocated();
700 if ( numberBitsAllocated == 0 || numberBitsAllocated == 12 )
702 numberBitsAllocated = 16;
705 size_t decSize = FileInternal->GetXSize()
706 * FileInternal->GetYSize()
707 * FileInternal->GetZSize()
708 * ( numberBitsAllocated / 8 )
709 * FileInternal->GetSamplesPerPixel();
710 size_t rgbSize = decSize;
711 if( FileInternal->HasLUT() )
712 rgbSize = decSize * 3;
717 if( decSize!=PixelWriteConverter->GetUserDataSize() )
719 gdcmWarningMacro( "Data size (Raw) is incorrect. Should be "
720 << decSize << " / Found :"
721 << PixelWriteConverter->GetUserDataSize() );
726 if( rgbSize!=PixelWriteConverter->GetUserDataSize() )
728 gdcmWarningMacro( "Data size (RGB) is incorrect. Should be "
729 << decSize << " / Found "
730 << PixelWriteConverter->GetUserDataSize() );
741 * \brief Updates the File to write RAW data (as opposed to RGB data)
742 * (modifies, when necessary, photochromatic interpretation,
743 * bits allocated, Pixels element VR)
745 void FileHelper::SetWriteToRaw()
747 if( FileInternal->GetNumberOfScalarComponents() == 3
748 && !FileInternal->HasLUT())
754 ValEntry *photInt = CopyValEntry(0x0028,0x0004);
755 if(FileInternal->HasLUT())
757 photInt->SetValue("PALETTE COLOR ");
761 photInt->SetValue("MONOCHROME2 ");
764 PixelWriteConverter->SetReadData(PixelReadConverter->GetRaw(),
765 PixelReadConverter->GetRawSize());
767 std::string vr = "OB";
768 if( FileInternal->GetBitsAllocated()>8 )
770 if( FileInternal->GetBitsAllocated()==24 ) // For RGB ACR files
773 CopyBinEntry(GetFile()->GetGrPixel(),GetFile()->GetNumPixel(),vr);
774 pixel->SetValue(GDCM_BINLOADED);
775 pixel->SetBinArea(PixelWriteConverter->GetData(),false);
776 pixel->SetLength(PixelWriteConverter->GetDataSize());
778 Archive->Push(photInt);
779 Archive->Push(pixel);
784 * \brief Updates the File to write RGB data (as opposed to RAW data)
785 * (modifies, when necessary, photochromatic interpretation,
786 * samples per pixel, Planar configuration,
787 * bits allocated, bits stored, high bit -ACR 24 bits-
788 * Pixels element VR, pushes out the LUT, )
790 void FileHelper::SetWriteToRGB()
792 if(FileInternal->GetNumberOfScalarComponents()==3)
794 PixelReadConverter->BuildRGBImage();
796 ValEntry *spp = CopyValEntry(0x0028,0x0002);
799 ValEntry *planConfig = CopyValEntry(0x0028,0x0006);
800 planConfig->SetValue("0 ");
802 ValEntry *photInt = CopyValEntry(0x0028,0x0004);
803 photInt->SetValue("RGB ");
805 if(PixelReadConverter->GetRGB())
807 PixelWriteConverter->SetReadData(PixelReadConverter->GetRGB(),
808 PixelReadConverter->GetRGBSize());
812 PixelWriteConverter->SetReadData(PixelReadConverter->GetRaw(),
813 PixelReadConverter->GetRawSize());
816 std::string vr = "OB";
817 if( FileInternal->GetBitsAllocated()>8 )
819 if( FileInternal->GetBitsAllocated()==24 ) // For RGB ACR files
822 CopyBinEntry(GetFile()->GetGrPixel(),GetFile()->GetNumPixel(),vr);
823 pixel->SetValue(GDCM_BINLOADED);
824 pixel->SetBinArea(PixelWriteConverter->GetData(),false);
825 pixel->SetLength(PixelWriteConverter->GetDataSize());
828 Archive->Push(planConfig);
829 Archive->Push(photInt);
830 Archive->Push(pixel);
833 Archive->Push(0x0028,0x1101);
834 Archive->Push(0x0028,0x1102);
835 Archive->Push(0x0028,0x1103);
836 Archive->Push(0x0028,0x1201);
837 Archive->Push(0x0028,0x1202);
838 Archive->Push(0x0028,0x1203);
840 // For old '24 Bits' ACR-NEMA
841 // Thus, we have a RGB image and the bits allocated = 24 and
842 // samples per pixels = 1 (in the read file)
843 if(FileInternal->GetBitsAllocated()==24)
845 ValEntry *bitsAlloc = CopyValEntry(0x0028,0x0100);
846 bitsAlloc->SetValue("8 ");
848 ValEntry *bitsStored = CopyValEntry(0x0028,0x0101);
849 bitsStored->SetValue("8 ");
851 ValEntry *highBit = CopyValEntry(0x0028,0x0102);
852 highBit->SetValue("7 ");
854 Archive->Push(bitsAlloc);
855 Archive->Push(bitsStored);
856 Archive->Push(highBit);
866 * \brief Restore the File write mode
868 void FileHelper::RestoreWrite()
870 Archive->Restore(0x0028,0x0002);
871 Archive->Restore(0x0028,0x0004);
872 Archive->Restore(0x0028,0x0006);
873 Archive->Restore(GetFile()->GetGrPixel(),GetFile()->GetNumPixel());
875 // For old ACR-NEMA (24 bits problem)
876 Archive->Restore(0x0028,0x0100);
877 Archive->Restore(0x0028,0x0101);
878 Archive->Restore(0x0028,0x0102);
881 Archive->Restore(0x0028,0x1101);
882 Archive->Restore(0x0028,0x1102);
883 Archive->Restore(0x0028,0x1103);
884 Archive->Restore(0x0028,0x1201);
885 Archive->Restore(0x0028,0x1202);
886 Archive->Restore(0x0028,0x1203);
888 // group 0002 may be pushed out for ACR-NEMA writting purposes
889 Archive->Restore(0x0002,0x0000);
890 Archive->Restore(0x0002,0x0001);
891 Archive->Restore(0x0002,0x0002);
892 Archive->Restore(0x0002,0x0003);
893 Archive->Restore(0x0002,0x0010);
894 Archive->Restore(0x0002,0x0012);
895 Archive->Restore(0x0002,0x0013);
896 Archive->Restore(0x0002,0x0016);
897 Archive->Restore(0x0002,0x0100);
898 Archive->Restore(0x0002,0x0102);
902 * \brief Pushes out the whole group 0002
903 * FIXME : better, set a flag to tell the writer not to write it ...
904 * FIXME : method should probably have an other name !
905 * SetWriteFileTypeToACR is NOT opposed to
906 * SetWriteFileTypeToExplicitVR and SetWriteFileTypeToImplicitVR
908 void FileHelper::SetWriteFileTypeToACR()
910 Archive->Push(0x0002,0x0000);
911 Archive->Push(0x0002,0x0001);
912 Archive->Push(0x0002,0x0002);
913 Archive->Push(0x0002,0x0003);
914 Archive->Push(0x0002,0x0010);
915 Archive->Push(0x0002,0x0012);
916 Archive->Push(0x0002,0x0013);
917 Archive->Push(0x0002,0x0016);
918 Archive->Push(0x0002,0x0100);
919 Archive->Push(0x0002,0x0102);
923 * \brief Sets in the File the TransferSyntax to 'Explicit VR Little Endian"
925 void FileHelper::SetWriteFileTypeToExplicitVR()
927 std::string ts = Util::DicomString(
928 Global::GetTS()->GetSpecialTransferSyntax(TS::ExplicitVRLittleEndian) );
930 ValEntry *tss = CopyValEntry(0x0002,0x0010);
937 * \brief Sets in the File the TransferSyntax to 'Implicit VR Little Endian"
939 void FileHelper::SetWriteFileTypeToImplicitVR()
941 std::string ts = Util::DicomString(
942 Global::GetTS()->GetSpecialTransferSyntax(TS::ImplicitVRLittleEndian) );
944 ValEntry *tss = CopyValEntry(0x0002,0x0010);
952 * \brief Restore in the File the initial group 0002
954 void FileHelper::RestoreWriteFileType()
959 * \brief Set the Write not to Libido format
961 void FileHelper::SetWriteToLibido()
963 ValEntry *oldRow = dynamic_cast<ValEntry *>
964 (FileInternal->GetDocEntry(0x0028, 0x0010));
965 ValEntry *oldCol = dynamic_cast<ValEntry *>
966 (FileInternal->GetDocEntry(0x0028, 0x0011));
968 if( oldRow && oldCol )
970 std::string rows, columns;
972 ValEntry *newRow=new ValEntry(oldRow->GetDictEntry());
973 ValEntry *newCol=new ValEntry(oldCol->GetDictEntry());
975 newRow->Copy(oldCol);
976 newCol->Copy(oldRow);
978 newRow->SetValue(oldCol->GetValue());
979 newCol->SetValue(oldRow->GetValue());
981 Archive->Push(newRow);
982 Archive->Push(newCol);
985 ValEntry *libidoCode = CopyValEntry(0x0008,0x0010);
986 libidoCode->SetValue("ACRNEMA_LIBIDO_1.1");
987 Archive->Push(libidoCode);
991 * \brief Set the Write not to No Libido format
993 void FileHelper::SetWriteToNoLibido()
995 ValEntry *recCode = dynamic_cast<ValEntry *>
996 (FileInternal->GetDocEntry(0x0008,0x0010));
999 if( recCode->GetValue() == "ACRNEMA_LIBIDO_1.1" )
1001 ValEntry *libidoCode = CopyValEntry(0x0008,0x0010);
1002 libidoCode->SetValue("");
1003 Archive->Push(libidoCode);
1009 * \brief Restore the Write format
1011 void FileHelper::RestoreWriteOfLibido()
1013 Archive->Restore(0x0028,0x0010);
1014 Archive->Restore(0x0028,0x0011);
1015 Archive->Restore(0x0008,0x0010);
1017 // Restore 'LibIDO-special' entries, if any
1018 Archive->Restore(0x0028,0x0015);
1019 Archive->Restore(0x0028,0x0016);
1020 Archive->Restore(0x0028,0x0017);
1021 Archive->Restore(0x0028,0x00199);
1025 * \brief Duplicates a ValEntry or creates it.
1026 * @param group Group number of the Entry
1027 * @param elem Element number of the Entry
1028 * \return pointer to the new Val Entry (NULL when creation failed).
1030 ValEntry *FileHelper::CopyValEntry(uint16_t group, uint16_t elem)
1032 DocEntry *oldE = FileInternal->GetDocEntry(group, elem);
1037 newE = new ValEntry(oldE->GetDictEntry());
1042 newE = GetFile()->NewValEntry(group, elem);
1049 * \brief Duplicates a BinEntry or creates it.
1050 * @param group Group number of the Entry
1051 * @param elem Element number of the Entry
1052 * @param vr Value Representation of the Entry
1053 * FIXME : what is it used for?
1054 * \return pointer to the new Bin Entry (NULL when creation failed).
1056 BinEntry *FileHelper::CopyBinEntry(uint16_t group, uint16_t elem,
1057 const std::string &vr)
1059 DocEntry *oldE = FileInternal->GetDocEntry(group, elem);
1063 if( oldE->GetVR()!=vr )
1068 newE = new BinEntry(oldE->GetDictEntry());
1073 newE = GetFile()->NewBinEntry(group, elem, vr);
1080 * \brief This method is called automatically, just before writting
1081 * in order to produce a 'True Dicom V3' image
1082 * We cannot know *how* the user made the File (reading an old ACR-NEMA
1083 * file or a not very clean DICOM file ...)
1085 * Just before writting :
1086 * - we check the Entries
1087 * - we create the mandatory entries if they are missing
1088 * - we modify the values if necessary
1089 * - we push the sensitive entries to the Archive
1090 * The writing process will restore the entries as they where before
1091 * entering FileHelper::CheckMandatoryElements, so the user will always
1092 * see the entries just as he left them.
1094 * \todo : - warn the user if we had to add some entries :
1095 * even if a mandatory entry is missing, we add it, with a default value
1096 * (we don't want to give up the writting process if user forgot to
1097 * specify Lena's Patient ID, for instance ...)
1098 * - read the whole PS 3.3 Part of DICOM (890 pages)
1099 * and write a *full* checker (probably one method per Modality ...)
1100 * Any contribution is welcome.
1101 * - write a user callable full checker, to allow post reading
1102 * and/or pre writting image consistency check.
1105 void FileHelper::CheckMandatoryElements()
1107 // just to remember : 'official' 0002 group
1109 //0002 0000 UL 1 Meta Group Length
1110 //0002 0001 OB 1 File Meta Information Version
1111 //0002 0002 UI 1 Media Stored SOP Class UID
1112 //0002 0003 UI 1 Media Stored SOP Instance UID
1113 //0002 0010 UI 1 Transfer Syntax UID
1114 //0002 0012 UI 1 Implementation Class UID
1115 //0002 0013 SH 1 Implementation Version Name
1116 //0002 0016 AE 1 Source Application Entity Title
1117 //0002 0100 UI 1 Private Information Creator
1118 //0002 0102 OB 1 Private Information
1120 // Create them if not found
1121 // Always modify the value
1122 // Push the entries to the archive.
1123 ValEntry *e_0002_0000 = CopyValEntry(0x0002,0x0000);
1124 e_0002_0000->SetValue("0"); // for the moment
1125 Archive->Push(e_0002_0000);
1127 BinEntry *e_0002_0001 = CopyBinEntry(0x0002,0x0001, "OB");
1128 e_0002_0001->SetBinArea((uint8_t*)Util::GetFileMetaInformationVersion(),
1130 e_0002_0001->SetLength(2);
1131 Archive->Push(e_0002_0001);
1133 // 'Media Stored SOP Class UID'
1134 ValEntry *e_0002_0002 = CopyValEntry(0x0002,0x0002);
1135 // [Secondary Capture Image Storage]
1136 e_0002_0002->SetValue("1.2.840.10008.5.1.4.1.1.7");
1137 Archive->Push(e_0002_0002);
1139 // 'Media Stored SOP Instance UID'
1140 ValEntry *e_0002_0003 = CopyValEntry(0x0002,0x0003);
1141 e_0002_0003->SetValue(Util::CreateUniqueUID());
1142 Archive->Push(e_0002_0003);
1144 // 'Implementation Class UID'
1145 ValEntry *e_0002_0012 = CopyValEntry(0x0002,0x0012);
1146 e_0002_0012->SetValue(Util::CreateUniqueUID());
1147 Archive->Push(e_0002_0012);
1149 // 'Implementation Version Name'
1150 ValEntry *e_0002_0013 = CopyValEntry(0x0002,0x0013);
1151 e_0002_0013->SetValue("GDCM 1.0");
1152 Archive->Push(e_0002_0013);
1154 //'Source Application Entity Title' Not Mandatory
1155 //ValEntry *e_0002_0016 = CopyValEntry(0x0002,0x0016);
1156 // e_0002_0016->SetValue("1.2.840.10008.5.1.4.1.1.7");
1157 // Archive->Push(e_0002_0016);
1159 // Push out 'LibIDO-special' entries, if any
1160 Archive->Push(0x0028,0x0015);
1161 Archive->Push(0x0028,0x0016);
1162 Archive->Push(0x0028,0x0017);
1163 Archive->Push(0x0028,0x00199);
1165 // Deal with the pb of (Bits Stored = 12)
1166 // - we're gonna write the image as Bits Stored = 16
1167 if ( FileInternal->GetEntryValue(0x0028,0x0100) == "12")
1169 ValEntry *e_0028_0100 = CopyValEntry(0x0028,0x0100);
1170 e_0028_0100->SetValue("16");
1171 Archive->Push(e_0028_0100);
1174 // --- Check UID-related Entries ---
1176 // If 'SOP Class UID' exists ('true DICOM' image)
1177 // we create the 'Source Image Sequence' SeqEntry
1178 // to hold informations about the Source Image
1180 ValEntry *e_0008_0016 = FileInternal->GetValEntry(0x0008, 0x0016);
1181 if ( e_0008_0016 != 0 )
1183 // Create 'Source Image Sequence' SeqEntry
1184 SeqEntry *s = new SeqEntry (
1185 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x2112) );
1186 SQItem *sqi = new SQItem(1);
1187 // (we assume 'SOP Instance UID' exists too)
1188 // create 'Referenced SOP Class UID'
1189 ValEntry *e_0008_1150 = new ValEntry(
1190 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x1150) );
1191 e_0008_1150->SetValue( e_0008_0016->GetValue());
1192 sqi->AddEntry(e_0008_1150);
1194 // create 'Referenced SOP Instance UID'
1195 ValEntry *e_0008_0018 = FileInternal->GetValEntry(0x0008, 0x0018);
1196 ValEntry *e_0008_1155 = new ValEntry(
1197 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x1155) );
1198 e_0008_1155->SetValue( e_0008_0018->GetValue());
1199 sqi->AddEntry(e_0008_1155);
1201 s->AddSQItem(sqi,1);
1202 // temporarily replaces any previous 'Source Image Sequence'
1205 // 'Image Type' (The written image is no longer an 'ORIGINAL' one)
1206 ValEntry *e_0008_0008 = CopyValEntry(0x0008,0x0008);
1207 e_0008_0008->SetValue("DERIVED\\PRIMARY");
1208 Archive->Push(e_0008_0008);
1212 // There was no 'SOP Class UID'.
1213 // the source image was NOT a true Dicom one.
1214 // We consider the image is a 'Secondary Capture' one
1216 e_0008_0016 = new ValEntry(
1217 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x0016) );
1218 // [Secondary Capture Image Storage]
1219 e_0008_0016 ->SetValue("1.2.840.10008.5.1.4.1.1.7");
1220 Archive->Push(e_0008_0016);
1223 // ---- The user will never have to take any action on the following ----.
1225 // new value for 'SOP Instance UID'
1226 ValEntry *e_0008_0018 = new ValEntry(
1227 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x0018) );
1228 e_0008_0018->SetValue( Util::CreateUniqueUID() );
1229 Archive->Push(e_0008_0018);
1231 // Instance Creation Date
1232 ValEntry *e_0008_0012 = CopyValEntry(0x0008,0x0012);
1233 std::string date = Util::GetCurrentDate();
1234 e_0008_0012->SetValue(date.c_str());
1235 Archive->Push(e_0008_0012);
1237 // Instance Creation Time
1238 ValEntry *e_0008_0013 = CopyValEntry(0x0008,0x0013);
1239 std::string time = Util::GetCurrentTime();
1240 e_0008_0013->SetValue(time.c_str());
1241 Archive->Push(e_0008_0013);
1243 // ----- Add Mandatory Entries if missing ---
1245 // Entries whose type is 1 are mandatory, with a mandatory value
1246 // Entries whose type is 1c are mandatory-inside-a-Sequence
1247 // Entries whose type is 2 are mandatory, with a optional value
1248 // Entries whose type is 2c are mandatory-inside-a-Sequence
1249 // Entries whose type is 3 are optional
1251 // 'Serie Instance UID'
1252 // Keep the value if exists
1253 // The user is allowed to create his own Series,
1254 // keeping the same 'Serie Instance UID' for various images
1255 // The user shouldn't add any image to a 'Manufacturer Serie'
1256 // but there is no way no to allowed him to do that
1257 ValEntry *e_0020_000e = FileInternal->GetValEntry(0x0020, 0x000e);
1260 e_0020_000e = new ValEntry(
1261 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0020, 0x000e) );
1262 e_0020_000e->SetValue(Util::CreateUniqueUID() );
1263 Archive->Push(e_0020_000e);
1266 // 'Study Instance UID'
1267 // Keep the value if exists
1268 // The user is allowed to create his own Study,
1269 // keeping the same 'Study Instance UID' for various images
1270 // The user may add images to a 'Manufacturer Study',
1271 // adding new series to an already existing Study
1272 ValEntry *e_0020_000d = FileInternal->GetValEntry(0x0020, 0x000d);
1275 e_0020_000d = new ValEntry(
1276 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0020, 0x000d) );
1277 e_0020_000d->SetValue(Util::CreateUniqueUID() );
1278 Archive->Push(e_0020_000d);
1281 // Modality : if missing we set it to 'OTher'
1282 ValEntry *e_0008_0060 = FileInternal->GetValEntry(0x0008, 0x0060);
1285 e_0008_0060 = new ValEntry(
1286 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x0060) );
1287 e_0008_0060->SetValue("OT");
1288 Archive->Push(e_0008_0060);
1291 // Manufacturer : if missing we set it to 'GDCM Factory'
1292 ValEntry *e_0008_0070 = FileInternal->GetValEntry(0x0008, 0x0070);
1295 e_0008_0070 = new ValEntry(
1296 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x0070) );
1297 e_0008_0070->SetValue("GDCM Factory");
1298 Archive->Push(e_0008_0070);
1301 // Institution Name : if missing we set it to 'GDCM Hospital'
1302 ValEntry *e_0008_0080 = FileInternal->GetValEntry(0x0008, 0x0080);
1305 e_0008_0080 = new ValEntry(
1306 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x0080) );
1307 e_0008_0080->SetValue("GDCM Hospital");
1308 Archive->Push(e_0008_0080);
1311 // Patient's Name : if missing, we set it to 'GDCM^Patient'
1312 ValEntry *e_0010_0010 = FileInternal->GetValEntry(0x0010, 0x0010);
1315 e_0010_0010 = new ValEntry(
1316 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0010, 0x0010) );
1317 e_0010_0010->SetValue("GDCM^Patient");
1318 Archive->Push(e_0010_0010);
1321 // Patient's Birth Date : 'type 2' entry -> must exist, value not mandatory
1322 ValEntry *e_0010_0030 = FileInternal->GetValEntry(0x0010, 0x0030);
1325 e_0010_0030 = new ValEntry(
1326 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0010, 0x0030) );
1327 e_0010_0030->SetValue("");
1328 Archive->Push(e_0010_0030);
1331 // Patient's Sex :'type 2' entry -> must exist, value not mandatory
1332 ValEntry *e_0010_0040 = FileInternal->GetValEntry(0x0010, 0x0040);
1335 e_0010_0040 = new ValEntry(
1336 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0010, 0x0040) );
1337 e_0010_0040->SetValue("");
1338 Archive->Push(e_0010_0040);
1341 // Referring Physician's Name :'type 2' entry -> must exist, value not mandatory
1342 ValEntry *e_0008_0090 = FileInternal->GetValEntry(0x0008, 0x0090);
1345 e_0008_0090 = new ValEntry(
1346 Global::GetDicts()->GetDefaultPubDict()->GetEntry(0x0008, 0x0090) );
1347 e_0008_0090->SetValue("");
1348 Archive->Push(e_0008_0090);
1351 // Remove some inconstencies (probably some more will be added)
1353 // if (0028 0008)Number of Frames exists
1354 // Push out (0020 0052),Frame of Reference UID
1355 // (only meaningfull within a Serie)
1356 ValEntry *e_0028_0008 = FileInternal->GetValEntry(0x0028, 0x0008);
1359 Archive->Push(0x0020, 0X0052);
1364 * \brief Restore in the File the initial group 0002
1366 void FileHelper::RestoreWriteMandatory()
1368 // group 0002 may be pushed out for ACR-NEMA writting purposes
1369 Archive->Restore(0x0002,0x0000);
1370 Archive->Restore(0x0002,0x0001);
1371 Archive->Restore(0x0002,0x0002);
1372 Archive->Restore(0x0002,0x0003);
1373 Archive->Restore(0x0002,0x0012);
1374 Archive->Restore(0x0002,0x0013);
1375 Archive->Restore(0x0002,0x0016);
1376 Archive->Restore(0x0002,0x0100);
1377 Archive->Restore(0x0002,0x0102);
1380 //-----------------------------------------------------------------------------
1383 * \brief Factorization for various forms of constructors.
1385 void FileHelper::Initialize()
1387 WriteMode = WMODE_RAW;
1388 WriteType = ExplicitVR;
1390 PixelReadConverter = new PixelReadConvert;
1391 PixelWriteConverter = new PixelWriteConvert;
1392 Archive = new DocEntryArchive( FileInternal );
1394 if ( FileInternal->IsReadable() )
1396 PixelReadConverter->GrabInformationsFromFile( FileInternal );
1401 * \brief Reads/[decompresses] the pixels,
1402 * *without* making RGB from Palette Colors
1403 * @return the pixels area, whatever its type
1404 * (uint8_t is just for prototyping : feel free to Cast it)
1406 uint8_t *FileHelper::GetRaw()
1408 uint8_t *raw = PixelReadConverter->GetRaw();
1411 // The Raw image migth not be loaded yet:
1412 std::ifstream *fp = FileInternal->OpenFile();
1413 PixelReadConverter->ReadAndDecompressPixelData( fp );
1415 FileInternal->CloseFile();
1417 raw = PixelReadConverter->GetRaw();
1420 gdcmWarningMacro( "Read/decompress of pixel data apparently went wrong.");
1428 //-----------------------------------------------------------------------------
1430 * \brief Prints the common part of ValEntry, BinEntry, SeqEntry
1431 * @param os ostream we want to print in
1432 * @param indent (unused)
1434 void FileHelper::Print(std::ostream &os, std::string const &)
1436 FileInternal->SetPrintLevel(PrintLevel);
1437 FileInternal->Print(os);
1439 PixelReadConverter->SetPrintLevel(PrintLevel);
1440 PixelReadConverter->Print(os);
1443 //-----------------------------------------------------------------------------
1444 } // end namespace gdcm