src/gdcmSerieHelper.h

   1 /*=========================================================================
   2
   3   Program:   gdcm
   4   Module:    $RCSfile: gdcmSerieHelper.h,v $
   5   Language:  C++
   6   Date:      $Date: 2010/04/09 15:38:18 $
   7   Version:   $Revision: 1.46 $
   8
   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.
  12
  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.
  16
  17 =========================================================================*/
  18
  19 #ifndef _GDCMSERIEHELPER_H_
  20 #define _GDCMSERIEHELPER_H_
  21
  22 #include "gdcmRefCounter.h"
  23 #include "gdcmTagKey.h"
  24 #include "gdcmDebug.h"  // for LEGACY
  25
  26 #include "gdcmDicomDirSerie.h"
  27
  28 #include <vector>
  29 #include <iostream>
  30 #include <map>
  31
  32 namespace GDCM_NAME_SPACE
  33 {
  34 class File;
  35
  36    typedef std::vector<File* > FileList;
  37 #ifndef GDCM_LEGACY_REMOVE
  38   typedef std::vector<File* > GdcmFileList;
  39 #endif
  40
  41    /// \brief XCoherent stands for 'Extra Coherent',
  42    ///        (The name 'Coherent' would be enough but it was used before;
  43    ///        I don't want to put a bomb in the code)
  44    ///        Any 'better name' is welcome !
  45    typedef std::map<std::string, FileList *> XCoherentFileSetmap;
  46
  47    typedef bool (*BOOL_FUNCTION_PFILE_PFILE_POINTER)(File *, File *);
  48
  49 //-----------------------------------------------------------------------------
  50 /**
  51  * \brief
  52  *  This class should be used for a stack of 2D dicom images.
  53  *
  54  *   - It allows to explore (recursively or not) a directory and
  55  *   makes a set of 'Single SerieUID Filesets'
  56  *   - It allows :
  57  *   - - to sort any of the 'Single SerieUID Fileset' on the image position.
  58  *   - - to split any of the Single SerieUID Filesets (better use this name than
  59  *   'CoherentFileList' : it's NOT a std::list, files are NOT coherent ...)
  60  *    into several XCoherent Filesets
  61  *   XCoherent stands for 'Extra Coherent' (same orientation, or same position)
  62  */
  63
  64 class GDCM_EXPORT SerieHelper  : public RefCounter
  65 {
  66    gdcmTypeMacro(SerieHelper);
  67
  68 public:
  69    /// SingleSerieUIDFileSetmap replaces the former CoherentFileListmap
  70    /// (List were actually std::vectors, and where no coherent at all :
  71    ///   They were only Single SeriesInstanceUID File sets)
  72    typedef std::map<std::string, FileList *> SingleSerieUIDFileSetmap;
  73
  74    typedef std::vector<File* > FileVector;
  75
  76 /// \brief Constructs a SerieHelper with a RefCounter
  77     static SerieHelper *New() {return new SerieHelper();}
  78
  79    virtual ~SerieHelper();
  80    void Print(std::ostream &os = std::cout, std::string const &indent = "" );
  81
  82    /// \todo should return bool or throw error ?
  83    void AddFileName(std::string const &filename);
  84    bool AddFile(File *header);
  85  //  GDCM_LEGACY(bool AddGdcmFile(File* header))
  86
  87    void SetDirectory(std::string const &dir, bool recursive=false);
  88    void SetDicomDirSerie(DicomDirSerie *se);
  89
  90    bool IsCoherent(FileList *fileSet);
  91    void OrderFileList(FileList *fileSet);
  92    void Clear() { ClearAll(); }
  93
  94    /// \brief Gets the FIRST Single SerieUID Fileset.
  95    ///        Deprecated; kept not to break the API
  96    /// \note Caller must call OrderFileList first
  97    /// @return the (first) Single SerieUID Fileset
  98    const FileList &GetFileList()
  99                            { return *SingleSerieUIDFileSetHT.begin()->second; }
 100
 101  //  GDCM_LEGACY(   FileList *GetFirstCoherentFileList()  )
 102  //  GDCM_LEGACY(   FileList *GetNextCoherentFileList()   )
 103  //  GDCM_LEGACY(   FileList *GetCoherentFileList(std::string serieUID)  )
 104
 105    FileList *GetFirstSingleSerieUIDFileSet();
 106    FileList *GetNextSingleSerieUIDFileSet();
 107    FileList *GetSingleSerieUIDFileSet(std::string serieUID);
 108    /// brief returns the 'Series Instance UID' Single SerieUID FileSet
 109    std::string GetCurrentSerieUIDFileSetUID()
 110                              { return  (*ItFileSetHt).first; }
 111
 112    /// \brief returns the distance between the 2 first -adjacent- slices,
 113    ///        along the slice normal; -1.0 if process failed
 114    ///        Computed during  ImagePositionPatientOrdering() process.
 115    ///        To be use *immediately after* ImagePositionPatientOrdering()
 116    ///        (it's a member of the FileHelper, *not* of the FileList !)
 117    double GetZSpacing() { return ZSpacing; }
 118
 119    /// All the following allow user to restrict DICOM file to be part
 120    /// of a particular serie
 121    /// \todo : find a trick to allow user to say the restrictions are ored
 122    ///         (not only anded)
 123    ///         ex : keep the images whose SerieNumber is 101 or 102 or 103.
 124    void AddRestriction(TagKey const &key);
 125    void AddRestriction(TagKey const &key, std::string const &value, int op);
 126    void AddRestriction(uint16_t group, uint16_t elem, std::string const &value,
 127                                                                     int op);
 128
 129    /// \brief Use additional series information such as ProtocolName
 130    ///        and SeriesName to identify when a single SeriesUID contains
 131    ///        multiple 3D volumes - as can occur with perfusion and DTI imaging
 132    void SetUseSeriesDetails( bool useSeriesDetails )
 133      {
 134      m_UseSeriesDetails = useSeriesDetails;
 135      }
 136    bool GetUseSeriesDetails()
 137      {
 138      return m_UseSeriesDetails;
 139      }
 140    /// \brief This function will add the following DICOM tag as being part of a
 141    /// 'fake' uid. This is usefull when the Serie UID is not enough to disseminate
 142    /// into multiple sub serie when needed:
 143    /// 0020 0011 Series Number
 144    /// 0018 0024 Sequence Name
 145    /// 0018 0050 Slice Thickness
 146    /// 0028 0010 Rows
 147    /// 0028 0011 Columns
 148    void CreateDefaultUniqueSeriesIdentifier();
 149
 150    void AddSeriesDetail(uint16_t group, uint16_t elem, bool convert);
 151
 152    std::string CreateUserDefinedFileIdentifier( File * inFile );
 153
 154    /// \brief Create a string that uniquely identifies a series.   By default
 155    //         uses the SeriesUID.   If UseSeriesDetails(true) has been called,
 156    //         then additional identifying information is used.
 157    std::string CreateUniqueSeriesIdentifier( File * inFile );
 158
 159 /**
 160  * \brief Sets the LoadMode as a boolean string.
 161  *        LD_NOSEQ, LD_NOSHADOW, LD_NOSHADOWSEQ
 162  *        ... (nothing more, right now)
 163  *        WARNING : before using LD_NOSHADOW, be sure *all* your files
 164  *        contain accurate values in the 0x0000 element (if any)
 165  *        of *each* Shadow Group. The parser will fail if the size is wrong !
 166  * @param   mode Load mode to be used
 167  */
 168    void SetLoadMode (int mode) { LoadMode = mode; }
 169
 170 /**
 171  * \brief Sets the DropDuplicatePositions as a boolean.
 172  * @param   drop  DropDuplicatePositions mode to be used
 173  */
 174    void SetDropDuplicatePositions (bool drop) { DropDuplicatePositions = drop; }
 175
 176 /// Brief User wants the files to be sorted Direct Order (default value)
 177    void SetSortOrderToDirect()  { DirectOrder = true;  }
 178
 179 /// Brief User wants the files to be sorted Reverse Order
 180    void SetSortOrderToReverse() { DirectOrder = false; }
 181
 182    /// to allow user to give is own comparison function
 183    void SetUserLessThanFunction( BOOL_FUNCTION_PFILE_PFILE_POINTER userFunc )
 184                         { UserLessThanFunction = userFunc; }
 185
 186    XCoherentFileSetmap SplitOnOrientation(FileList *fileSet);
 187    XCoherentFileSetmap SplitOnPosition(FileList *fileSet);
 188    XCoherentFileSetmap SplitOnTagValue(FileList *fileSet,
 189                                                uint16_t group, uint16_t element);
 190 protected :
 191    SerieHelper();
 192
 193 private:
 194    void ClearAll();
 195    bool UserOrdering(FileList *fileSet);
 196    bool ImagePositionPatientOrdering(FileList *fileSet);
 197    bool ImageNumberOrdering(FileList *fileSet);
 198    bool FileNameOrdering(FileList *fileSet);
 199
 200    static bool ImageNumberLessThan(File *file1, File *file2);
 201    static bool ImageNumberGreaterThan(File *file1, File *file2);
 202    static bool FileNameLessThan(File *file1, File *file2);
 203    static bool FileNameGreaterThan(File *file1, File *file2);
 204
 205 //Attributes:
 206
 207    SingleSerieUIDFileSetmap SingleSerieUIDFileSetHT;
 208    SingleSerieUIDFileSetmap::iterator ItFileSetHt;
 209
 210 #ifndef GDCM_LEGACY_REMOVE
 211    typedef std::pair<TagKey, std::string> Rule;
 212    typedef std::vector<Rule> SerieRestrictions;
 213    SerieRestrictions Restrictions;
 214 #endif
 215
 216    // New style for (extented) Rules
 217    typedef struct {
 218       uint16_t group;
 219       uint16_t elem;
 220       std::string value;
 221       int op;
 222    } ExRule;
 223    typedef std::vector<ExRule> SerieExRestrictions;
 224    SerieExRestrictions ExRestrictions;
 225    SerieExRestrictions ExRefine;
 226
 227    typedef struct {
 228       uint16_t group;
 229       uint16_t elem;
 230       bool convert;
 231    } ExDetail;
 232    typedef std::vector<ExDetail> SeriesExDetails;
 233    SeriesExDetails ExDetails;
 234
 235
 236    /// \brief Bit string integer (each one considered as a boolean)
 237    ///        Bit 0 : Skip Sequences,    if possible
 238    ///        Bit 1 : Skip Shadow Groups if possible
 239    ///        Probabely, some more to add
 240    int LoadMode;
 241
 242    /// \brief whether we want to sort in direct order or not (reverse order).
 243    ///        To be used by aware user only
 244    bool DirectOrder;
 245
 246    /// \brief If user knows more about his images than gdcm does,
 247    ///        he may supply his own comparison function.
 248    BOOL_FUNCTION_PFILE_PFILE_POINTER UserLessThanFunction;
 249
 250    bool DropDuplicatePositions;
 251
 252    void Sort(FileList *fileList, bool (*pt2Func)( File *file1, File *file2) );
 253
 254    bool m_UseSeriesDetails;
 255    /// \brief distance between the 2 first -adjancent- slices, along the slice normal
 256    ///        Computed during  ImagePositionPatientOrdering() process.
 257    ///        (we need to sort the whole image set to do that)
 258    ///        To be used *immediately after* ImagePositionPatientOrdering()
 259    ///        (it's a member of the FileHelper, *not* of the FileList !)
 260    double ZSpacing;
 261 };
 262
 263 } // end namespace gdcm
 264
 265 //-----------------------------------------------------------------------------
 266 #endif