X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmTS.cxx;h=93e9ef3b1b3bf317f9c81b9441db48799e1e9c93;hb=76ac46b16fd92d0fb444f786d4946424d9029315;hp=f0f2049c98b396da2975088157dc1df339c6db42;hpb=31c63a6cb25b49a119e0132201c051ade2e17704;p=gdcm.git diff --git a/src/gdcmTS.cxx b/src/gdcmTS.cxx index f0f2049c..93e9ef3b 100644 --- a/src/gdcmTS.cxx +++ b/src/gdcmTS.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: gdcmTS.cxx,v $ Language: C++ - Date: $Date: 2005/01/11 16:44:43 $ - Version: $Revision: 1.35 $ + Date: $Date: 2009/03/04 08:57:42 $ + Version: $Revision: 1.56 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -24,20 +24,26 @@ #include #include #include +#include // for isdigit // TODO -// a lot of troubles expected with TS : 1.2.840.113619.5.2 +// troubles expected with TS : 1.2.840.113619.5.2 // Implicit VR - Big Endian -// see : http://www.gemedicalsystemseurope.com/euen/it_solutions/pdf/lsqxi_rev2.pdf -// +// http://www.gemedicalsystemseurope.com/euen/it_solutions/pdf/lsqxi_rev2.pdf +// G.E. deliberately violated a lot of Dicom rules are +// (probabely to avoid other people to read their images) +// Just try and error on new images : +// PrintFile debug filein=... +// and fix the bugs -namespace gdcm +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 DLX (G.E Private) + // Implicit VR Big Endian (G.E Private) "1.2.840.113619.5.2", // Explicit VR Little Endian "1.2.840.10008.1.2.1", @@ -57,28 +63,46 @@ static const char *SpecialStrings[] = { "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]) + // 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" + "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 TS::TS() { + std::string filename = DictSet::BuildDictPath() + DICT_TS; std::ifstream from(filename.c_str()); - if( !from ) + if ( !from ) { - gdcmVerboseMacro("Can't open dictionary" << filename.c_str()); + gdcmWarningMacro("Can't open dictionary" << filename.c_str()); FillDefaultTSDict( TsMap ); } else @@ -92,67 +116,71 @@ TS::TS() from >> std::ws; std::getline(from, name); - if(key != "") + if (key != "") { TsMap[key] = name; } } + from.close(); } } -//----------------------------------------------------------------------------- TS::~TS() { TsMap.clear(); } -//----------------------------------------------------------------------------- -// Print -/** - * \ingroup VR - * \brief Print all - * @param os The output stream to be written to. - */ -void TS::Print(std::ostream &os) -{ - std::ostringstream s; - - for (TSHT::const_iterator it = TsMap.begin(); it != TsMap.end(); ++it) - { - s << "TS : " << it->first << " = " << it->second << std::endl; - } - os << s.str(); -} - //----------------------------------------------------------------------------- // Public + +/// \brief returns occurence number of the given key int TS::Count(TSKey const &key) { return TsMap.count(key); } -TSAtr const & TS::GetValue(TSKey const &key) +/// \brief returns the human readable value of a Transfer Syntax string +TSAtr const &TS::GetValue(TSKey const &key) { - TSHT::const_iterator it = TsMap.find(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]) ) + { + 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(); } +/** + * \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 ( IsTransferSyntax(key) ) { if ( key == SpecialStrings[RLELossless] ) { @@ -162,11 +190,17 @@ bool TS::IsRLELossless(TSKey const &key) return r; } +/** + * \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 ( IsTransferSyntax(key) ) { if ( key == SpecialStrings[JPEGFullProgressionProcess10_12] || key == SpecialStrings[JPEGLosslessProcess14] @@ -178,17 +212,40 @@ bool TS::IsJPEGLossless(TSKey const &key) return r; } +/** + * \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) +{ + 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; +} + /** * \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. + * @return True when JPEG2000 (Lossly or LossLess) found. + * False in all other cases. */ bool TS::IsJPEG2000(TSKey const &key) { bool r = false; // First check this is an actual transfer syntax - if( IsTransferSyntax(key) ) + if ( IsTransferSyntax(key) ) { if ( key == SpecialStrings[JPEG2000Lossless] || key == SpecialStrings[JPEG2000] ) @@ -200,22 +257,41 @@ bool TS::IsJPEG2000(TSKey const &key) } /** - * \brief Determines if the Transfer Syntax corresponds to any form - * of Jpeg encoded Pixel data. + * \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) ) + if ( IsTransferSyntax(key) ) { - if ( key == SpecialStrings[JPEGBaselineProcess1] - || key == SpecialStrings[JPEGExtendedProcess2_4] - || key == SpecialStrings[JPEGExtendedProcess3_5] - || key == SpecialStrings[JPEGSpectralSelectionProcess6_8] - || IsJPEGLossless( key ) - || IsJPEG2000( key ) ) + if ( IsJPEGLossy( key ) + || IsJPEGLossless( key ) + || IsJPEG2000( key ) + || IsJPEGLS( key ) + ) + { + r = true; + } + } + return r; +} + +/** + * \brief Determines if the Transfer Syntax corresponds to any form + * of Jpeg-LS encoded Pixel data. + * @return True when any form of JPEG-LS found. False otherwise. + */ +bool TS::IsJPEGLS(TSKey const &key) +{ + bool r = false; + // First check this is an actual transfer syntax + if ( IsTransferSyntax(key) ) + { + if ( key == SpecialStrings[JPEGLSLossless] + || key == SpecialStrings[JPEGLSNearLossless] ) { r = true; } @@ -224,18 +300,17 @@ bool TS::IsJPEG(TSKey const &key) } /** - * \brief Determines if the Transfer Syntax corresponds to encapsulated - * of encoded Pixel Data (as opposed to native). - * @return True when encapsulated. False when native. + * \brief Determines if the Transfer Syntax corresponds to any form + * of MPEG encoded Pixel data. + * @return True when any form of MPEG found. False otherwise. */ -bool TS::IsEncapsulate(TSKey const &key) +bool TS::IsMPEG(TSKey const &key) { bool r = false; // First check this is an actual transfer syntax - if( IsTransferSyntax(key) ) + if ( IsTransferSyntax(key) ) { - if ( key == SpecialStrings[RLELossless] - || IsJPEG(key) ) + if ( key == SpecialStrings[MPEG2MainProfile] ) { r = true; } @@ -243,6 +318,29 @@ bool TS::IsEncapsulate(TSKey const &key) return r; } +/** + * \brief Determines if the SOP id corresponds to any form + * of UltrasoundImageStorage_Retired. + * @return True when Ultrasound Image Storage Retired. False otherwise. + */ +bool TS::IsUltrasoundImageStorage_Retired(TSKey const &key) +{ + bool r = false; + // First check this is an actual SOP id + if ( IsTransferSyntax(key) ) + { + if ( key == SpecialStrings[UltrasoundImageStorage_Retired] ) + { + r = true; + } + } + return r; +} +/** + * \brief GetSpecialTransferSyntax ?? + * @param key TSKey const &key ?? + * @return TS::SpecialType ??. + */ TS::SpecialType TS::GetSpecialTransferSyntax(TSKey const &key) { for (int i = 0; SpecialStrings[i] != NULL; i++) @@ -252,11 +350,15 @@ TS::SpecialType TS::GetSpecialTransferSyntax(TSKey const &key) return SpecialType(i); } } - return UnknownTS; } -const char* TS::GetSpecialTransferSyntax(SpecialType t) +/** + * \brief GetSpecialTransferSyntax ?? + * @param t SpecialType t ?? + * @return char* TS : SpecialStrings[t] ??. + */ +const char *TS::GetSpecialTransferSyntax(SpecialType t) { return SpecialStrings[t]; } @@ -268,5 +370,21 @@ const char* TS::GetSpecialTransferSyntax(SpecialType t) // Private //----------------------------------------------------------------------------- +// Print +/** + * \brief Print all + * @param os The output stream to be written to. + */ +void TS::Print(std::ostream &os,std::string const &) +{ + std::ostringstream s; + + for (TSHT::const_iterator it = TsMap.begin(); it != TsMap.end(); ++it) + { + s << "TS : " << it->first << " = " << it->second << std::endl; + } + os << s.str(); +} +//----------------------------------------------------------------------------- } // end namespace gdcm