1 /*=========================================================================
4 Module: $RCSfile: gdcmDataEntry.cxx,v $
6 Date: $Date: 2005/10/25 14:27:37 $
7 Version: $Revision: 1.11 $
9 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
10 l'Image). All rights reserved. See Doc/License.txt or
11 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
13 This software is distributed WITHOUT ANY WARRANTY; without even
14 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 PURPOSE. See the above copyright notices for more information.
17 =========================================================================*/
19 #include "gdcmDataEntry.h"
22 #include "gdcmGlobal.h"
24 #include "gdcmDebug.h"
30 //-----------------------------------------------------------------------------
31 #define MAX_SIZE_PRINT_ELEMENT_VALUE 0x7fffffff
32 uint32_t DataEntry::MaxSizePrintEntry = MAX_SIZE_PRINT_ELEMENT_VALUE;
34 //-----------------------------------------------------------------------------
35 // Constructor / Destructor
37 * \brief Constructor for a given DictEntry
38 * @param e Pointer to existing dictionary entry
40 DataEntry::DataEntry(DictEntry *e)
51 * \brief Constructor for a given DocEntry
52 * @param e Pointer to existing Doc entry
54 DataEntry::DataEntry(DocEntry *e)
55 : DocEntry(e->GetDictEntry())
65 * \brief Canonical destructor.
67 DataEntry::~DataEntry ()
72 //-----------------------------------------------------------------------------
75 //-----------------------------------------------------------------------------
78 * \brief Sets the value (non string) of the current Dicom Header Entry
80 void DataEntry::SetBinArea( uint8_t *area, bool self )
90 * \brief Inserts the value (non string) into the current Dicom Header Entry
92 void DataEntry::CopyBinArea( uint8_t *area, uint32_t length )
96 uint32_t lgh = length + length%2;
99 if( area && length > 0 )
102 memcpy(BinArea,area,length);
106 State = STATE_LOADED;
110 void DataEntry::SetValue(const uint32_t &id, const double &val)
114 State = STATE_LOADED;
116 if( id > GetValueCount() )
118 gdcmErrorMacro("Index (" << id << ")is greater than the data size");
122 const VRKey &vr = GetVR();
123 if( vr == "US" || vr == "SS" )
125 uint16_t *data = (uint16_t *)BinArea;
126 data[id] = (uint16_t)val;
128 else if( vr == "UL" || vr == "SL" )
130 uint32_t *data = (uint32_t *)BinArea;
131 data[id] = (uint32_t)val;
133 else if( vr == "FL" )
135 float *data = (float *)BinArea;
136 data[id] = (float)val;
138 else if( vr == "FD" )
140 double *data = (double *)BinArea;
141 data[id] = (double)val;
143 else if( Global::GetVR()->IsVROfStringRepresentable(vr) )
145 gdcmErrorMacro("SetValue on String representable not implemented yet");
149 BinArea[id] = (uint8_t)val;
153 * \brief returns, as a double (?!?) one of the values
154 // (when entry is multivaluated), identified by its index.
155 // Returns 0.0 if index is wrong
156 // FIXME : warn the user there was a problem !
158 double DataEntry::GetValue(const uint32_t &id) const
162 gdcmErrorMacro("BinArea not set. Can't get the value");
166 uint32_t count = GetValueCount();
169 gdcmErrorMacro("Index (" << id << ")is greater than the data size");
173 // FIX the API : user *knows* that entry contains a US
174 // and he receives a double ?!?
176 const VRKey &vr = GetVR();
177 if( vr == "US" || vr == "SS" )
178 return ((uint16_t *)BinArea)[id];
179 else if( vr == "UL" || vr == "SL" )
180 return ((uint32_t *)BinArea)[id];
181 else if( vr == "FL" )
182 return ((float *)BinArea)[id];
183 else if( vr == "FD" )
184 return ((double *)BinArea)[id];
185 else if( Global::GetVR()->IsVROfStringRepresentable(vr) )
189 // Don't use std::string to accelerate processing
191 char *tmp = new char[GetLength()+1];
192 memcpy(tmp,BinArea,GetLength());
203 for(uint32_t i=0;i<GetLength();i++)
233 * \brief Checks if the multiplicity of the value follows Dictionary VM
235 bool DataEntry::IsValueCountValid() const
239 const std::string &strVM = GetVM();
240 uint32_t vc = GetValueCount();
242 // FIXME : what shall we do with VM = "2-n", "3-n", etc
246 // make sure it is at least one ??? FIXME
247 valid = vc >= 1 || vc == 0;
251 std::istringstream os;
255 // vm respect the one from the dict
256 // vm is 0 (we need to check if this element is allowed to be empty) FIXME
260 // Entries whose type is 1 are mandatory, with a mandatory value.
261 // Entries whose type is 1c are mandatory-inside-a-Sequence,
262 // with a mandatory value.
263 // Entries whose type is 2 are mandatory, with an optional value.
264 // Entries whose type is 2c are mandatory-inside-a-Sequence,
265 // with an optional value.
266 // Entries whose type is 3 are optional.
268 // case vc == 0 is only applicable for 'type 2' entries.
269 // Problem : entry type may depend on the modality and/or the Sequence
270 // it's embedded in !
271 // (Get the information in the 'Conformance Statements' ...)
272 valid = vc == vm || vc == 0;
278 * \brief returns the number of elementary values
280 uint32_t DataEntry::GetValueCount( ) const
282 const VRKey &vr = GetVR();
283 if( vr == "US" || vr == "SS" )
284 return GetLength()/sizeof(uint16_t);
285 else if( vr == "UL" || vr == "SL" )
286 return GetLength()/sizeof(uint32_t);
287 else if( vr == "FL" )
288 return GetLength()/sizeof(float);
289 else if( vr == "FD" )
290 return GetLength()/sizeof(double);
291 else if( Global::GetVR()->IsVROfStringRepresentable(vr) )
293 // Some element in DICOM are allowed to be empty
294 if( !GetLength() ) return 0;
295 // Don't use std::string to accelerate processing
297 for(uint32_t i=0;i<GetLength();i++)
299 if( BinArea[i] == '\\')
308 * \brief Sets the 'value' of an Entry, passed as a std::string
309 * @param value string representation of the value to be set
311 void DataEntry::SetString(std::string const &value)
315 const VRKey &vr = GetVR();
316 if ( vr == "US" || vr == "SS" )
318 std::vector<std::string> tokens;
319 Util::Tokenize (value, tokens, "\\");
320 SetLength(tokens.size()*sizeof(uint16_t));
323 uint16_t *data = (uint16_t *)BinArea;
324 for (unsigned int i=0; i<tokens.size();i++)
325 data[i] = atoi(tokens[i].c_str());
328 else if ( vr == "UL" || vr == "SL" )
330 std::vector<std::string> tokens;
331 Util::Tokenize (value, tokens, "\\");
332 SetLength(tokens.size()*sizeof(uint32_t));
335 uint32_t *data = (uint32_t *)BinArea;
336 for (unsigned int i=0; i<tokens.size();i++)
337 data[i] = atoi(tokens[i].c_str());
340 else if ( vr == "FL" )
342 std::vector<std::string> tokens;
343 Util::Tokenize (value, tokens, "\\");
344 SetLength(tokens.size()*sizeof(float));
347 float *data = (float *)BinArea;
348 for (unsigned int i=0; i<tokens.size();i++)
349 data[i] = (float)atof(tokens[i].c_str());
352 else if ( vr == "FD" )
354 std::vector<std::string> tokens;
355 Util::Tokenize (value, tokens, "\\");
356 SetLength(tokens.size()*sizeof(double));
359 double *data = (double *)BinArea;
360 for (unsigned int i=0; i<tokens.size();i++)
361 data[i] = atof(tokens[i].c_str());
366 if( value.size() > 0 )
368 // FIXME : should be quicker if we don't create one more std::string
369 // just to make even the length of a char array ...
372 std::string finalVal = Util::DicomString( value.c_str() );
373 SetLength(finalVal.size());
376 memcpy(BinArea, &(finalVal[0]), finalVal.size());
379 size_t l = value.size();
382 memcpy(BinArea, &(value[0]), l);
387 State = STATE_LOADED;
390 * \brief returns as a string (when possible) the value of the DataEntry
392 std::string const &DataEntry::GetString() const
394 static std::ostringstream s;
395 const VRKey &vr = GetVR();
403 // When short integer(s) are stored, convert the following (n * 2) characters
404 // as a displayable string, the values being separated by a back-slash
406 if( vr == "US" || vr == "SS" )
408 uint16_t *data=(uint16_t *)BinArea;
410 for (unsigned int i=0; i < GetValueCount(); i++)
418 // See above comment on multiple short integers (mutatis mutandis).
419 else if( vr == "UL" || vr == "SL" )
421 uint32_t *data=(uint32_t *)BinArea;
423 for (unsigned int i=0; i < GetValueCount(); i++)
431 else if( vr == "FL" )
433 float *data=(float *)BinArea;
435 for (unsigned int i=0; i < GetValueCount(); i++)
443 else if( vr == "FD" )
445 double *data=(double *)BinArea;
447 for (unsigned int i=0; i < GetValueCount(); i++)
456 StrArea.append((const char *)BinArea,GetLength());
461 * \brief Copies all the attributes from an other DocEntry
462 * @param doc entry to copy from
464 void DataEntry::Copy(DocEntry *doc)
468 DataEntry *entry = dynamic_cast<DataEntry *>(doc);
471 State = entry->State;
473 CopyBinArea(entry->BinArea,entry->GetLength());
477 * \brief Writes the value of a DataEntry
478 * @param fp already open ofstream pointer
479 * @param filetype type of the file (ACR, ImplicitVR, ExplicitVR, ...)
481 void DataEntry::WriteContent(std::ofstream *fp, FileType filetype)
483 DocEntry::WriteContent(fp, filetype);
485 if ( GetGroup() == 0xfffe )
487 return; //delimitors have NO value
490 uint8_t *binArea8 = BinArea; //safe notation
491 size_t lgr = GetLength();
492 if (BinArea) // the binArea was *actually* loaded
495 // The same operation should be done if we wanted
496 // to write image with Big Endian Transfer Syntax,
497 // while working on Little Endian Processor
498 // --> forget Big Endian Transfer Syntax writting!
499 // Next DICOM version will give it up ...
502 // The stuff looks nice, but it's probably bugged,
503 // since troubles occur on big endian processors (SunSparc, Motorola)
504 // while reading the pixels of a
505 // gdcm-written Little-Endian 16 bits per pixel image
507 #if defined(GDCM_WORDS_BIGENDIAN) || defined(GDCM_FORCE_BIGENDIAN_EMULATION)
509 /// \todo FIXME : Right now, we only care of Pixels element
510 /// we should deal with *all* the BinEntries
511 /// Well, not really since we are not interpreting values read...
513 // 8 Bits Pixels *are* OB, 16 Bits Pixels *are* OW
514 // -value forced while Reading process-
517 // --> the following lines *looked* very clever,
518 // --> but they don't work on big endian processors.
519 // --> since I've no access for the moment to a big endian proc :-(
520 // --> I comment them out, to see the result on the dash board
523 // --> Revert to initial code : TestWriteSimple hangs on Darwin :-(
524 if (GetGroup() == 0x7fe0 && GetVR() == "OW")
526 uint16_t *binArea16 = (uint16_t*)binArea8;
527 binary_write (*fp, binArea16, lgr );
531 // For any other VR, DataEntry is re-written as-is
532 binary_write (*fp, binArea8, lgr );
536 // --> remove the following line, an uncomment the previous ones,
537 // --> if it doesn't work better
539 /*binary_write ( *fp, binArea8, lgr ); // Elem value*/
542 binary_write ( *fp, binArea8, lgr ); // Elem value
543 #endif //GDCM_WORDS_BIGENDIAN
547 // nothing was loaded, but we need to skip space on disc
549 // --> WARNING : nothing is written;
550 // --> the initial data (on the the source image) is lost
551 // --> user is *not* informed !
553 fp->seekp(lgr, std::ios::cur);
557 //-----------------------------------------------------------------------------
559 void DataEntry::NewBinArea(void)
562 if( GetLength() > 0 )
563 BinArea = new uint8_t[GetLength()];
567 void DataEntry::DeleteBinArea(void)
569 if (BinArea && SelfArea)
576 //-----------------------------------------------------------------------------
579 //-----------------------------------------------------------------------------
582 * \brief Prints a DataEntry (Dicom entry)
583 * @param os ostream we want to print in
584 * @param indent Indentation string to be prepended during printing
586 void DataEntry::Print(std::ostream &os, std::string const & )
591 uint16_t g = GetGroup();
592 if (g == 0xfffe) // delimiters have NO value
594 return; // just to avoid identing all the remaining code
597 std::ostringstream s;
603 const VRKey &vr = GetVR();
605 if( vr == "US" || vr == "SS" )
606 s << " [" << GetString() << "]";
607 else if( vr == "UL" || vr == "SL" )
608 s << " [" << GetString() << "]";
609 else if ( vr == "FL" )
610 s << " [" << GetString() << "]";
611 else if ( vr == "FD" )
612 s << " [" << GetString() << "]";
615 if(Global::GetVR()->IsVROfStringRepresentable(vr))
617 // replace non printable characters by '.'
618 std::string cleanString = Util::CreateCleanString(v);
619 if ( cleanString.length() <= GetMaxSizePrintEntry()
622 // FIXME : when IsNotLoaded(), you create a Clean String ?!?
623 // FIXME : PrintLevel<2 *does* print the values
624 // (3 is only for extra offsets printing)
625 // What do you wanted to do ? JPR
627 s << " [" << cleanString << "]";
631 s << " [gdcm::too long for print (" << cleanString.length() << ") ]";
636 // A lot of Private elements (with no VR) contain actually
637 // only printable characters;
638 // Let's deal with them as is they were VR std::string representable
640 if ( Util::IsCleanArea( GetBinArea(), GetLength() ) )
642 // FIXME : since the 'Area' *is* clean, just use
643 // a 'CreateString' method, to save CPU time.
644 std::string cleanString =
645 Util::CreateCleanString( BinArea,GetLength() );
646 s << " [" << cleanString << "]";
650 s << " [" << GDCM_BINLOADED << ";"
651 << "length = " << GetLength() << "]";
659 s << " [" << GDCM_NOTLOADED << "]";
660 else if( IsUnfound() )
661 s << " [" << GDCM_UNFOUND << "]";
662 else if( IsUnread() )
663 s << " [" << GDCM_UNREAD << "]";
664 else if ( GetLength() == 0 )
669 s << " (" << GDCM_PIXELDATA << ")";
671 // Display the UID value (instead of displaying only the rough code)
672 // First 'clean' trailing character (space or zero)
675 const uint16_t &gr = GetGroup();
676 const uint16_t &elt = GetElement();
677 TS *ts = Global::GetTS();
681 // Any more to be displayed ?
682 if ( elt == 0x0010 || elt == 0x0002 )
684 if ( v.length() != 0 ) // for brain damaged headers
686 if ( ! isdigit((unsigned char)v[v.length()-1]) )
688 v.erase(v.length()-1, 1);
691 s << " ==>\t[" << ts->GetValue(v) << "]";
694 else if (gr == 0x0008)
696 if ( elt == 0x0016 || elt == 0x1150 )
698 if ( v.length() != 0 ) // for brain damaged headers
700 if ( ! isdigit((unsigned char)v[v.length()-1]) )
702 v.erase(v.length()-1, 1);
705 s << " ==>\t[" << ts->GetValue(v) << "]";
708 else if (gr == 0x0004)
710 if ( elt == 0x1510 || elt == 0x1512 )
712 if ( v.length() != 0 ) // for brain damaged headers
714 if ( ! isdigit((unsigned char)v[v.length()-1]) )
716 v.erase(v.length()-1, 1);
719 s << " ==>\t[" << ts->GetValue(v) << "]";
727 //-----------------------------------------------------------------------------
728 } // end namespace gdcm