# gdcmSequenceOfItems.h

Go to the documentation of this file.

    00001 /*=========================================================================
    00002
    00003   Program: GDCM (Grassroots DICOM). A DICOM library
    00004   Module:  $URL$
    00005
    00006   Copyright (c) 2006-2010 Mathieu Malaterre
    00007   All rights reserved.
    00008   See Copyright.txt or http://gdcm.sourceforge.net/Copyright.html for details.
    00009
    00010      This software is distributed WITHOUT ANY WARRANTY; without even
    00011      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
    00012      PURPOSE.  See the above copyright notice for more information.
    00013
    00014 =========================================================================*/
    00015
    00016 #ifndef GDCMSEQUENCEOFITEMS_H
    00017 #define GDCMSEQUENCEOFITEMS_H
    00018
    00019 #include "gdcmValue.h"
    00020 #include "gdcmItem.h"
    00021
    00022 #include <vector>
    00023
    00024 namespace gdcm
    00025 {
    00026
    00039 class GDCM_EXPORT SequenceOfItems : public Value
    00040 {
    00041 public:
    00042   // Typdefs:
    00043   typedef std::vector< Item > ItemVector;
    00044   typedef ItemVector::iterator Iterator;
    00045   typedef ItemVector::const_iterator ConstIterator;
    00046   Iterator Begin() { return Items.begin(); }
    00047   Iterator End() { return Items.end(); }
    00048   ConstIterator Begin() const { return Items.begin(); }
    00049   ConstIterator End() const { return Items.end(); }
    00050
    00052   SequenceOfItems():SequenceLengthField(0xFFFFFFFF) { }
    00053   //SequenceOfItems(VL const &vl = 0xFFFFFFFF):SequenceLengthField(vl),NType(type) { }
    00054
    00056   VL GetLength() const { return SequenceLengthField; }
    00058   void SetLength(VL length) {
    00059     SequenceLengthField = length;
    00060   }
    00061   void SetLengthToUndefined() {
    00062     SequenceLengthField = 0xFFFFFFFF;
    00063   }
    00065   bool IsUndefinedLength() const {
    00066     return SequenceLengthField.IsUndefined();
    00067   }
    00068
    00069   template <typename TDE>
    00070   VL ComputeLength() const;
    00071   void Clear() {}
    00072
    00074   void AddItem(Item const &item);
    00075
    00076   unsigned int GetNumberOfItems() const {  return Items.size(); }
    00077   void SetNumberOfItems(unsigned int n) {  Items.resize(n); }
    00078
    00079   /* WARNING: first item is #1 (see DICOM standard)
    00080    *  Each Item shall be implicitly assigned an ordinal position starting with the value 1 for the
    00081    * first Item in the Sequence, and incremented by 1 with each subsequent Item. The last Item in the
    00082    * Sequence shall have an ordinal position equal to the number of Items in the Sequence.
    00083    */
    00084   const Item &GetItem(unsigned int position) const;
    00085   Item &GetItem(unsigned int position);
    00086
    00087   SequenceOfItems &operator=(const SequenceOfItems &val) {
    00088     SequenceLengthField = val.SequenceLengthField;
    00089     Items = val.Items;
    00090     return *this;
    00091     }
    00092
    00093   template <typename TDE, typename TSwap>
    00094   std::istream &Read(std::istream &is)
    00095     {
    00096     const Tag seqDelItem(0xfffe,0xe0dd);
    00097     if( SequenceLengthField.IsUndefined() )
    00098       {
    00099       Item item;
    00100       while( item.Read<TDE,TSwap>(is) && item.GetTag() != seqDelItem )
    00101         {
    00102         //gdcmDebugMacro( "Item: " << item );
    00103         assert( item.GetTag() != seqDelItem );
    00104         Items.push_back( item );
    00105         item.Clear();
    00106         }
    00107       //assert( item.GetTag() == seqDelItem && item.GetVL() == 0 );
    00108       }
    00109     else
    00110       {
    00111       Item item;
    00112       VL l = 0;
    00113       //is.seekg( SequenceLengthField, std::ios::cur ); return is;
    00114       while( l != SequenceLengthField )
    00115         {
    00116         try
    00117           {
    00118           item.Read<TDE,TSwap>(is);
    00119           }
    00120         catch( Exception &ex )
    00121           {
    00122           if( strcmp( ex.GetDescription(), "Changed Length" ) == 0 )
    00123             {
    00124             VL newlength = l + item.template GetLength<TDE>();
    00125             if( newlength > SequenceLengthField )
    00126               {
    00127               // BogugsItemAndSequenceLength.dcm
    00128               gdcmWarningMacro( "SQ length is wrong" );
    00129               SequenceLengthField = newlength;
    00130               }
    00131             }
    00132           else
    00133             {
    00134             throw ex;
    00135             }
    00136           }
    00137         if( item.GetTag() == seqDelItem )
    00138           {
    00139           gdcmWarningMacro( "SegDelItem found in defined length Sequence" );
    00140           }
    00141         //assert( item.GetTag() == Tag(0xfffe,0xe000) );
    00142         Items.push_back( item );
    00143         l += item.template GetLength<TDE>();
    00144         if( l > SequenceLengthField ) throw "Length of Item larger than expected";
    00145         assert( l <= SequenceLengthField );
    00146 #ifdef GDCM_SUPPORT_BROKEN_IMPLEMENTATION
    00147         // MR_Philips_Intera_No_PrivateSequenceImplicitVR.dcm
    00148         // (0x2005, 0x1080): for some reason computation of length fails...
    00149         if( SequenceLengthField == 778 && l == 774 )
    00150           {
    00151           gdcmWarningMacro( "PMS: Super bad hack" );
    00152           SequenceLengthField = l;
    00153           throw Exception( "Wrong Length" );
    00154           //l = SequenceLengthField;
    00155           }
    00156         // Bug_Philips_ItemTag_3F3F
    00157         // (0x2005, 0x1080): Because we do not handle fully the bug at the item
    00158         // level we need to check here too
    00159         else if ( SequenceLengthField == 444 && l == 3*71 )
    00160           {
    00161           // This one is a double bug. Item length is wrong and impact SQ length
    00162           gdcmWarningMacro( "PMS: Super bad hack" );
    00163           l = SequenceLengthField;
    00164           }
    00165 #endif
    00166         }
    00167       assert( l == SequenceLengthField );
    00168       }
    00169     return is;
    00170     }
    00171
    00172   template <typename TDE,typename TSwap>
    00173   std::ostream const &Write(std::ostream &os) const
    00174     {
    00175     typename ItemVector::const_iterator it = Items.begin();
    00176     for(;it != Items.end(); ++it)
    00177       {
    00178       it->Write<TDE,TSwap>(os);
    00179       }
    00180     if( SequenceLengthField.IsUndefined() )
    00181       {
    00182       // seq del item is not stored, write it !
    00183       const Tag seqDelItem(0xfffe,0xe0dd);
    00184       seqDelItem.Write<TSwap>(os);
    00185       VL zero = 0;
    00186       zero.Write<TSwap>(os);
    00187       }
    00188
    00189     return os;
    00190     }
    00191
    00192 //protected:
    00193   void Print(std::ostream &os) const {
    00194     os << "\t(" << SequenceLengthField << ")\n";
    00195     ItemVector::const_iterator it =
    00196       Items.begin();
    00197     for(;it != Items.end(); ++it)
    00198       {
    00199       os << "  " << *it;
    00200       }
    00201     if( SequenceLengthField.IsUndefined() )
    00202       {
    00203       const Tag seqDelItem(0xfffe,0xe0dd);
    00204       VL zero = 0;
    00205       os << seqDelItem;
    00206       os << "\t" << zero;
    00207       }
    00208   }
    00209
    00210   static SmartPointer<SequenceOfItems> New()
    00211   {
    00212      return new SequenceOfItems;
    00213   }
    00214   bool FindDataElement(const Tag &t) const;
    00215
    00216   bool operator==(const Value &val) const
    00217     {
    00218     const SequenceOfItems &sqi = dynamic_cast<const SequenceOfItems&>(val);
    00219     return SequenceLengthField == sqi.SequenceLengthField &&
    00220       Items == sqi.Items;
    00221     }
    00222
    00223 private:
    00224 public:
    00226   VL SequenceLengthField;
    00228   ItemVector Items;
    00229 };
    00230
    00231 } // end namespace gdcm
    00232
    00233 #include "gdcmSequenceOfItems.txx"
    00234
    00235 #endif //GDCMSEQUENCEOFITEMS_H
    00236