X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmTS.cxx;h=d3751af010960f019bbcb709a2a12d97ef3a2692;hb=fc4f7706c78baf3f60474c06d871b82db221f37b;hp=82a25780a35fc20d3cf3cc8d0cc40e83b3331f31;hpb=f5575c1671810f3eeecae9967861dff3f1f75f05;p=gdcm.git diff --git a/src/gdcmTS.cxx b/src/gdcmTS.cxx index 82a25780..d3751af0 100644 --- a/src/gdcmTS.cxx +++ b/src/gdcmTS.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: gdcmTS.cxx,v $ Language: C++ - Date: $Date: 2004/10/08 04:43:38 $ - Version: $Revision: 1.25 $ + Date: $Date: 2007/05/23 14:18:11 $ + Version: $Revision: 1.55 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -24,71 +24,343 @@ #include #include #include +#include // for isdigit + +// TODO +// troubles expected with TS : 1.2.840.113619.5.2 +// Implicit VR - Big Endian +// http://www.gemedicalsystemseurope.com/euen/it_solutions/pdf/lsqxi_rev2.pdf +// G.E. deliberately violated a lot of Dicom rules are +// (probabely to to avoid other people to read their images) +// Just try and error on new images : +// PrintFile debug filein=... +// and fix the bugs + +namespace GDCM_NAME_SPACE +{ +//----------------------------------------------------------------------------- +/// \brief Transfer Syntaxes gdcm deals with (internal use only) +static const char *SpecialStrings[] = { + // Implicit VR Little Endian + "1.2.840.10008.1.2", + // Implicit VR Big Endian (G.E Private) + "1.2.840.113619.5.2", + // Explicit VR Little Endian + "1.2.840.10008.1.2.1", + // Deflated Explicit VR Little Endian + "1.2.840.10008.1.2.1.99", + // Explicit VR Big Endian + "1.2.840.10008.1.2.2", + // JPEG Baseline (Process 1) + "1.2.840.10008.1.2.4.50", + // JPEG Extended (Process 2 & 4) + "1.2.840.10008.1.2.4.51", + // JPEG Extended (Process 3 & 5) + "1.2.840.10008.1.2.4.52", + // JPEG Spectral Selection, Non-Hierarchical (Process 6 & 8) + "1.2.840.10008.1.2.4.53", + // JPEG Full Progression, Non-Hierarchical (Process 10 & 12) + "1.2.840.10008.1.2.4.55", + // JPEG Lossless, Non-Hierarchical (Process 14) + "1.2.840.10008.1.2.4.57", + // JPEG Lossless, Hierarchical, First-Order Prediction (Process 14, + // [Selection Value 1]) + "1.2.840.10008.1.2.4.70", + // JPEG-LS Lossless Image Compression + "1.2.840.10008.1.2.4.80", + // JPEG-LS Lossy (Near-Lossless) Image Compression + "1.2.840.10008.1.2.4.81", + // JPEG 2000 Lossless + "1.2.840.10008.1.2.4.90", + // JPEG 2000 + "1.2.840.10008.1.2.4.91", + // RLE Lossless + "1.2.840.10008.1.2.5", + // MPEG2 Main Profile @ Main Level + "1.2.840.10008.1.2.4.100", + + // The following are *not* t.s. but SOP uid + // Ultrasound Image Storage (Retired) + "1.2.840.10008.5.1.4.1.1.6", + + // Unknown + "Unknown Transfer Syntax", // Pretty sure we never use this case... + NULL // Compilers have no obligation to finish by NULL, do it ourself +}; + +//----------------------------------------------------------------------------- +/// \brief auto generated function, to fill up the Dicom Dictionnary, +/// if relevant file is not found on user's disk +void FillDefaultTSDict(TSHT &ts); //----------------------------------------------------------------------------- // Constructor / Destructor -gdcmTS::gdcmTS() +TS::TS() { - std::string filename=gdcmDictSet::BuildDictPath() + std::string(DICT_TS); + + std::string filename = DictSet::BuildDictPath() + DICT_TS; std::ifstream from(filename.c_str()); - dbg.Error(!from, "gdcmTS::gdcmTS: can't open dictionary",filename.c_str()); + if ( !from ) + { + gdcmWarningMacro("Can't open dictionary" << filename.c_str()); + FillDefaultTSDict( TsMap ); + } + else + { + TSKey key; + TSAtr name; + + while (!from.eof()) + { + from >> key; + from >> std::ws; + std::getline(from, name); + + if (key != "") + { + TsMap[key] = name; + } + } + + from.close(); + } +} + +TS::~TS() +{ + TsMap.clear(); +} - std::string key; - std::string name; +//----------------------------------------------------------------------------- +// Public - while (!from.eof()) +/// \brief returns occurence number of the given key +int TS::Count(TSKey const &key) +{ + return TsMap.count(key); +} + +/// \brief returns the human readable value of a Transfer Syntax string +TSAtr const &TS::GetValue(TSKey const &key) +{ + // First thing clean up the string + // (sometimes the transfer syntax is padded with spaces) + std::string copy = key; + while ( copy.size() && !isdigit((unsigned char)copy[copy.size()-1]) ) { - from >> key; + copy.erase(copy.size()-1, 1); + } + + TSHT::const_iterator it = TsMap.find(copy); + if (it == TsMap.end()) + { + return GDCM_UNFOUND; + } + return it->second; +} +/** + * \brief Determines if the key passed corresponds to a 'Transfer Syntax' + * as defined in DICOM (and stored in gdcm::TS class) + * @return True when key is an actual 'Transfer Syntax'. + * False in all other cases. + */ +bool TS::IsTransferSyntax(TSKey const &key) +{ + TSHT::const_iterator it = TsMap.find(key); + return it != TsMap.end(); +} - from >> std::ws; // used to be eatwhite(from); - std::getline(from, name); /// MEMORY LEAK +/** + * \brief Determines if the Transfer Syntax was already encountered + * and if it corresponds to a Run Length Encoding Lossless one + * @return True when Run Length Encoding Lossless found. + * False in all other cases. + */ +bool TS::IsRLELossless(TSKey const &key) +{ + bool r = false; + // First check this is an actual transfer syntax + if ( IsTransferSyntax(key) ) + { + if ( key == SpecialStrings[RLELossless] ) + { + r = true; + } + } + return r; +} - if(key!="") +/** + * \brief Determines if the Transfer Syntax was already encountered + * and if it corresponds to a 'classical' JPEG Lossless one + * @return True when 'classical' Lossless found. + * False in all other cases. + */ +bool TS::IsJPEGLossless(TSKey const &key) +{ + bool r = false; + // First check this is an actual transfer syntax + if ( IsTransferSyntax(key) ) + { + if ( key == SpecialStrings[JPEGFullProgressionProcess10_12] + || key == SpecialStrings[JPEGLosslessProcess14] + || key == SpecialStrings[JPEGLosslessProcess14_1] ) { - ts[key]=name; + r = true; } } - from.close(); + return r; } -//----------------------------------------------------------------------------- -gdcmTS::~gdcmTS() +/** + * \brief Determines if the Transfer Syntax was already encountered + * and if it corresponds to a 'classical' JPEG Lossy one + * @return True when 'classical' Lossy found. + * False in all other cases. + */ +bool TS::IsJPEGLossy(TSKey const &key) { - ts.clear(); + bool r = false; + // First check this is an actual transfer syntax + if ( IsTransferSyntax(key) ) + { + if ( key == SpecialStrings[JPEGBaselineProcess1] + || key == SpecialStrings[JPEGExtendedProcess2_4] + || key == SpecialStrings[JPEGExtendedProcess3_5] + || key == SpecialStrings[JPEGSpectralSelectionProcess6_8] ) + { + r = true; + } + } + return r; } -//----------------------------------------------------------------------------- -// Print /** - * \ingroup gdcmVR - * \brief Print all - * @param os The output stream to be written to. + * \brief Determines if the Transfer Syntax was already encountered + * and if it corresponds to a JPEG2000 one + * @return True when JPEG2000 (Lossly or LossLess) found. + * False in all other cases. */ -void gdcmTS::Print(std::ostream &os) +bool TS::IsJPEG2000(TSKey const &key) { - std::ostringstream s; + bool r = false; + // First check this is an actual transfer syntax + if ( IsTransferSyntax(key) ) + { + if ( key == SpecialStrings[JPEG2000Lossless] + || key == SpecialStrings[JPEG2000] ) + { + r = true; + } + } + return r; +} - for (TSHT::iterator it = ts.begin(); it != ts.end(); ++it) +/** + * \brief Determines if the Transfer Syntax corresponds to + * 'classical' Jpeg Lossless or Jpeg lossy. + * @return True when any form of JPEG found. False otherwise. + */ +bool TS::IsJPEG(TSKey const &key) +{ + bool r = false; + // First check this is an actual transfer syntax + if ( IsTransferSyntax(key) ) { - s << "TS : "<first<<" = "<second<first << " = " << it->second << std::endl; + } + os << s.str(); +} + +//----------------------------------------------------------------------------- +} // end namespace gdcm