X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmSerieHelper.cxx;h=16616963b698c1c38ca51dad6d38e43ef2c62afa;hb=f2dd0107ba49eda3386b34715cdc0c046d12397f;hp=c8c1cbba86453cc3d502085159704e8562811264;hpb=fc71b6d6df72ee061f0a88e3de407a83511f8ed2;p=gdcm.git diff --git a/src/gdcmSerieHelper.cxx b/src/gdcmSerieHelper.cxx index c8c1cbba..16616963 100644 --- a/src/gdcmSerieHelper.cxx +++ b/src/gdcmSerieHelper.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: gdcmSerieHelper.cxx,v $ Language: C++ - Date: $Date: 2006/03/30 16:51:29 $ - Version: $Revision: 1.49 $ + Date: $Date: 2007/07/05 09:34:12 $ + Version: $Revision: 1.59 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -19,7 +19,7 @@ #include "gdcmSerieHelper.h" #include "gdcmDirList.h" #include "gdcmFile.h" -#include "gdcmDictEntry.h" // for TranslateToKey +//#include "gdcmDictEntry.h" // for TranslateToKey : no more ! #include "gdcmDebug.h" #include "gdcmUtil.h" @@ -29,7 +29,7 @@ #include #include //for sscanf -namespace gdcm +namespace GDCM_NAME_SPACE { //----------------------------------------------------------------------------- @@ -43,6 +43,7 @@ SerieHelper::SerieHelper() ClearAll(); UserLessThanFunction = 0; DirectOrder = true; + } /** @@ -63,8 +64,8 @@ void SerieHelper::ClearAll() FileList *l = GetFirstSingleSerieUIDFileSet(); while (l) { - // For all the gdcm::File of a File set - for (gdcm::FileList::iterator it = l->begin(); + // For all the GDCM_NAME_SPACE::File of a File set + for (GDCM_NAME_SPACE::FileList::iterator it = l->begin(); it != l->end(); ++it) { @@ -84,7 +85,7 @@ void SerieHelper::ClearAll() // Public /** - * \brief add a gdcm::File to the Fileset corresponding to its Serie UID + * \brief add a GDCM_NAME_SPACE::File to the Fileset corresponding to its Serie UID * @param filename Name of the file to deal with */ void SerieHelper::AddFileName(std::string const &filename) @@ -111,10 +112,10 @@ void SerieHelper::AddFileName(std::string const &filename) } /** - * \brief add a gdcm::File to the first (and supposed to be unique) file set - * of the gdcm::SerieHelper. + * \brief add a GDCM_NAME_SPACE::File to the first (and supposed to be unique) file set + * of the GDCM_NAME_SPACE::SerieHelper. * \warning : this method should be used by aware users only! - * Passing a gdcm::File* has the same effect than passing a file name! + * Passing a GDCM_NAME_SPACE::File* has the same effect than passing a file name! * \todo : decide which one is wrong (the method, or the commentary)! * the following comment doesn't match the method :-( * User is supposed to know the files he want to deal with @@ -126,7 +127,7 @@ void SerieHelper::AddFileName(std::string const &filename) * vtkGdcmReader parsing twice the same files. * *no* coherence check is performed, but those specified * by SerieHelper::AddRestriction() - * @param header gdcm::File* of the file to deal with + * @param header GDCM_NAME_SPACE::File* of the file to deal with * @return true if file was added, false if file was rejected */ bool SerieHelper::AddFile(File *header) @@ -163,7 +164,7 @@ bool SerieHelper::AddFile(File *header) if ( SingleSerieUIDFileSetHT.count(id) == 0 ) { - gdcmDebugMacro(" New Serie UID :[" << id << "]"); + gdcmDebugMacro(" New/gdcmSerieHelper.cxx Serie UID :[" << id << "]"); // create a std::list in 'id' position SingleSerieUIDFileSetHT[id] = new FileList; } @@ -209,7 +210,7 @@ void SerieHelper::AddRestriction(TagKey const &key) ExRefine.push_back( r ); } -#ifndef GDCM_LEGACY_REMOVE +//#ifndef GDCM_LEGACY_REMOVE /** * \brief add a rule for restricting a DICOM file to be in the serie we are * trying to find. For example you can select only the DICOM files from a @@ -223,13 +224,15 @@ void SerieHelper::AddRestriction(TagKey const &key) * @deprecated use : AddRestriction(TagKey const &key, * std::string const &value, int op); */ + void SerieHelper::AddRestriction(uint16_t group, uint16_t elem, std::string const &value, int op) { TagKey t(group, elem); AddRestriction(t, value, op); } -#endif + +//#endif /** * \brief add an extra 'SerieDetail' for building a 'Serie Identifier' @@ -240,8 +243,7 @@ void SerieHelper::AddRestriction(uint16_t group, uint16_t elem, * e.g : 100 would be *before* 20; 000020.00 vs 00100.00 : OK */ void SerieHelper::AddSeriesDetail(uint16_t group, uint16_t elem, bool convert) -{ - +{ ExDetail d; d.group = group; d.elem = elem; @@ -265,6 +267,22 @@ void SerieHelper::SetDirectory(std::string const &dir, bool recursive) } } +/** + * \brief Sets the DicomDirSerie + * @param se DicomDirSerie to deal with + */ +void SerieHelper::SetDicomDirSerie(DicomDirSerie *se) +{ + DirList dirList(se); + + DirListType filenames_list = dirList.GetFilenames(); + for( DirListType::const_iterator it = filenames_list.begin(); + it != filenames_list.end(); ++it) + { + AddFileName( *it ); + } +} + /** * \brief Sorts the given Fileset * \warning This could be implemented in a 'Strategy Pattern' approach @@ -274,7 +292,10 @@ void SerieHelper::SetDirectory(std::string const &dir, bool recursive) */ void SerieHelper::OrderFileList(FileList *fileSet) { - + // Only computed during ImagePositionPatientOrdering + // (need to sort the FileList using IPP and IOP !) + ZSpacing = -1.0; + if ( SerieHelper::UserLessThanFunction ) { UserOrdering( fileSet ); @@ -327,12 +348,13 @@ bool SerieHelper::IsCoherent(FileList *fileSet) return true; } -#ifndef GDCM_LEGACY_REMOVE +//#ifndef GDCM_LEGACY_REMOVE /** * \brief accessor (DEPRECATED : use GetFirstSingleSerieUIDFileSet ) * Warning : 'coherent' means here they have the same Serie UID * @return The first FileList if found, otherwhise NULL */ + /* FileList *SerieHelper::GetFirstCoherentFileList() { ItFileSetHt = SingleSerieUIDFileSetHT.begin(); @@ -340,13 +362,14 @@ FileList *SerieHelper::GetFirstCoherentFileList() return ItFileSetHt->second; return NULL; } - +*/ /** * \brief accessor (DEPRECATED : use GetNextSingleSerieUIDFileSet ) * Warning : 'coherent' means here they have the same Serie UID * \note : meaningfull only if GetFirstCoherentFileList() already called * @return The next FileList if found, otherwhise NULL */ + /* FileList *SerieHelper::GetNextCoherentFileList() { gdcmAssertMacro (ItFileSetHt != SingleSerieUIDFileSetHT.end()); @@ -356,6 +379,7 @@ FileList *SerieHelper::GetNextCoherentFileList() return ItFileSetHt->second; return NULL; } +*/ /** * \brief accessor (DEPRECATED : use GetSingleSerieUIDFileSet ) @@ -363,13 +387,15 @@ FileList *SerieHelper::GetNextCoherentFileList() * @param SerieUID SerieUID * \return pointer to the FileList if found, otherwhise NULL */ + /* FileList *SerieHelper::GetCoherentFileList(std::string SerieUID) { if ( SingleSerieUIDFileSetHT.count(SerieUID) == 0 ) return 0; return SingleSerieUIDFileSetHT[SerieUID]; } -#endif +*/ +//#endif /** @@ -429,7 +455,7 @@ XCoherentFileSetmap SerieHelper::SplitOnOrientation(FileList *fileSet) std::ostringstream ossOrient; FileList::const_iterator it = fileSet->begin(); - it ++; + //it ++; for ( ; it != fileSet->end(); ++it) @@ -479,7 +505,7 @@ XCoherentFileSetmap SerieHelper::SplitOnPosition(FileList *fileSet) std::ostringstream ossPosition; std::string strPosition; // re computed FileList::const_iterator it = fileSet->begin(); - it ++; + //it ++; for ( ; it != fileSet->end(); ++it) @@ -539,7 +565,7 @@ XCoherentFileSetmap SerieHelper::SplitOnPosition(FileList *fileSet) * value of a given Tag * @param fileSet File Set to be splitted * @param group group number of the target Element - * @param elem element number of the target Element + * @param element element number of the target Element * \return std::map of 'Xcoherent' File sets */ @@ -555,7 +581,7 @@ XCoherentFileSetmap SerieHelper::SplitOnTagValue(FileList *fileSet, std::string strTagValue; // read on disc FileList::const_iterator it = fileSet->begin(); - it ++; + //it ++; for ( ; it != fileSet->end(); ++it) @@ -684,8 +710,9 @@ bool SerieHelper::ImagePositionPatientOrdering( FileList *fileList ) // Find out if min/max are coherent if ( min == max ) { - gdcmWarningMacro("Looks like all images have the exact same image position" - << ". No PositionPatientOrdering sort performed" ); + gdcmWarningMacro("Looks like all images have the exact same image position. " + << "No PositionPatientOrdering sort performed. " + << "No 'ZSpacing' calculated! "); return false; } @@ -697,18 +724,21 @@ bool SerieHelper::ImagePositionPatientOrdering( FileList *fileList ) { if (distmultimap.count((*it2).first) != 1) { - gdcmErrorMacro("File: " + gdcmErrorMacro("File: [" << ((*it2).second->GetFileName()) - << " Distance: " + << "] : more than ONE file at distance: '" << (*it2).first - << " position is not unique"); + << " (position is not unique!) " + << "No PositionPatientOrdering sort performed. " + << "No 'ZSpacing' calculated! "); ok = false; } } if (!ok) { - return false; + if (! DropDuplicatePositions) + return false; } // Now, we could calculate Z Spacing as the difference @@ -735,6 +765,12 @@ bool SerieHelper::ImagePositionPatientOrdering( FileList *fileList ) ++it3) { fileList->push_back( (*it3).second ); + if (DropDuplicatePositions) + { + it3 = distmultimap.upper_bound((*it3).first); // skip all duplicates + if (it3 == distmultimap.end() ) // if last image, stop iterate + break; + } } } else // user asked for reverse order @@ -745,6 +781,12 @@ bool SerieHelper::ImagePositionPatientOrdering( FileList *fileList ) { it4--; fileList->push_back( (*it4).second ); + if (DropDuplicatePositions) // skip all duplicates + { + it4 = distmultimap.upper_bound((*it4).first); + if (it4 == distmultimap.begin() ) // if first image, stop iterate + break; + } } while (it4 != distmultimap.begin() ); } @@ -847,7 +889,7 @@ bool SerieHelper::UserOrdering(FileList *fileList) */ void SerieHelper::Print(std::ostream &os, std::string const &indent) { - // For all the Coherent File lists of the gdcm::Serie + // For all the Coherent File lists of the GDCM_NAME_SPACE::Serie SingleSerieUIDFileSetmap::iterator itl = SingleSerieUIDFileSetHT.begin(); if ( itl == SingleSerieUIDFileSetHT.end() ) { @@ -879,19 +921,23 @@ void SerieHelper::CreateDefaultUniqueSeriesIdentifier() // A scout scan prior to a CT volume scan can share the same // SeriesUID, but they will sometimes have a different Series Number AddRestriction( TagKey(0x0020, 0x0011) ); + // 0018 0024 Sequence Name // For T1-map and phase-contrast MRA, the different flip angles and // directions are only distinguished by the Sequence Name AddRestriction( TagKey(0x0018, 0x0024) ); + // 0018 0050 Slice Thickness // On some CT systems, scout scans and subsequence volume scans will // have the same SeriesUID and Series Number - YET the slice // thickness will differ from the scout slice and the volume slices. AddRestriction( TagKey(0x0018, 0x0050)); + // 0028 0010 Rows // If the 2D images in a sequence don't have the same number of rows, // then it is difficult to reconstruct them into a 3D volume. AddRestriction( TagKey(0x0028, 0x0010)); + // 0028 0011 Columns // If the 2D images in a sequence don't have the same number of columns, // then it is difficult to reconstruct them into a 3D volume. @@ -908,7 +954,7 @@ void SerieHelper::CreateDefaultUniqueSeriesIdentifier() * (he knows more than we do about his images!) * ex : in tagging series, the only pertnent tag is * 0018|1312 [In-plane Phase Encoding Direction] value : ROW/COLUMN - * @param inFile gdcm::File we want to build a Serie Identifier for. + * @param inFile GDCM_NAME_SPACE::File we want to build a Serie Identifier for. * @return the SeriesIdentifier */ std::string SerieHelper::CreateUniqueSeriesIdentifier( File *inFile ) @@ -926,7 +972,7 @@ std::string SerieHelper::CreateUniqueSeriesIdentifier( File *inFile ) { const ExRule &r = *it2; std::string s = inFile->GetEntryString( r.group, r.elem ); - if( s == gdcm::GDCM_UNFOUND ) + if( s == GDCM_UNFOUND ) { s = ""; } @@ -955,7 +1001,7 @@ std::string SerieHelper::CreateUniqueSeriesIdentifier( File *inFile ) else // Could not open inFile { gdcmWarningMacro("Could not parse series info."); - std::string id = gdcm::GDCM_UNFOUND; + std::string id = GDCM_UNFOUND; return id; } } @@ -985,7 +1031,7 @@ std::string SerieHelper::CreateUserDefinedFileIdentifier( File * inFile ) const ExDetail &r = *it2; s = inFile->GetEntryString( r.group, r.elem ); - // User is allowed to ask 'convertion', to allow further ordering + // User is allowed to ask for 'convertion', to allow further ordering // e.g : 100 would be *before* 20; 000020.00 vs 00100.00 : OK if (it2->convert) { @@ -1003,10 +1049,10 @@ std::string SerieHelper::CreateUserDefinedFileIdentifier( File * inFile ) for(unsigned int i=0; i= 'a' && s[i] <= 'z') - || (s[i] >= '0' && s[i] <= '9') - || (s[i] >= 'A' && s[i] <= 'Z'))) + && !( s[i] == '.' || s[i] == '%' || s[i] == '_' + || (s[i] >= 'a' && s[i] <= 'z') + || (s[i] >= '0' && s[i] <= '9') + || (s[i] >= 'A' && s[i] <= 'Z'))) { s.erase(i, 1); } @@ -1015,7 +1061,8 @@ std::string SerieHelper::CreateUserDefinedFileIdentifier( File * inFile ) id += s.c_str(); id += "%%%"; // make the FileIdentifier Tokenizable } - + id += inFile->GetFileName(); + id += "%%%"; return id; } @@ -1029,5 +1076,14 @@ void SerieHelper::Sort(FileList *fileList, bool (*pt2Func)( File *file1, File *f std::sort(fileList->begin(), fileList->end(), pt2Func ); } +/* +#ifndef GDCM_LEGACY_REMOVE +bool SerieHelper::AddGdcmFile(File* header) +{ + return AddFile(header); +} +#endif +*/ + //----------------------------------------------------------------------------- } // end namespace gdcm