X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2FgdcmHeaderHelper.cxx;h=54f9b4ec897fcf5ec4890e5e0fdab234e62d9526;hb=c5618e6b8276e1b21bbea354fe8bcba3c78eb4bd;hp=16d1157690fbea3cfd32b810b94a39fc5735500b;hpb=0574f31372f78bf8399e42ee04bd0010e684c440;p=gdcm.git diff --git a/src/gdcmHeaderHelper.cxx b/src/gdcmHeaderHelper.cxx index 16d11576..54f9b4ec 100644 --- a/src/gdcmHeaderHelper.cxx +++ b/src/gdcmHeaderHelper.cxx @@ -3,64 +3,13 @@ #include "gdcmHeaderHelper.h" #include "gdcmDirList.h" -#include "gdcmUtil.h" //for debug +#include "gdcmDebug.h" #include #include +#include -/*#ifdef _MSC_VER - #include - - int GetDir(std::string dPath, std::list &filenames) - { - //For now dPath should have an ending "\" - WIN32_FIND_DATA FileData; - HANDLE hFile; - hFile = FindFirstFile((dPath+"*").c_str(), &FileData); - if ( hFile == INVALID_HANDLE_VALUE ) - { - //No files ! - return false; - } - - if( strncmp(FileData.cFileName, ".", 1) != 0 ) - filenames.push_back( dPath+FileData.cFileName ); - while( FindNextFile(hFile, &FileData ) != 0) - { - if( strncmp(FileData.cFileName, ".", 1) != 0 ) - filenames.push_back( dPath+FileData.cFileName ); - } - return true; - } - -#else - #include - - int GetDir(std::string dPath, std::list &filenames) - { - DIR *dir = opendir( dPath.c_str() ); - if (dir == NULL) - return false; - struct dirent *entry; - while((entry = readdir(dir)) != NULL) - { - // if( strncmp(entry->d_name, ".", 1) != 0 && strncmp(entry->d_name, "..", 2) != 0) - if( strncmp(entry->d_name, ".", 1) != 0 ) - { - filenames.push_back( dPath + "/" + entry->d_name ); - } - } - closedir(dir); - return true; - } - -#endif*/ - -//----------------------------------------------------------------------------- -// gdcmHeaderHelper //----------------------------------------------------------------------------- -// Constructor / Destructor /** - * \ingroup gdcmHeaderHelper * \brief constructor */ gdcmHeaderHelper::gdcmHeaderHelper() : gdcmHeader( ) { @@ -68,7 +17,6 @@ gdcmHeaderHelper::gdcmHeaderHelper() : gdcmHeader( ) { } /** - * \ingroup gdcmHeaderHelper * \brief constructor * @param InFilename Name of the file to deal with * @param exception_on_error @@ -95,7 +43,6 @@ gdcmHeaderHelper::gdcmHeaderHelper(const char *InFilename, //----------------------------------------------------------------------------- // Public /** - * \ingroup gdcmHeaderHelper * \brief Returns the size (in bytes) of a single pixel of data. * @return The size in bytes of a single pixel of data. * @@ -170,7 +117,7 @@ std::string gdcmHeaderHelper::GetPixelType() { /** * \ingroup gdcmHeaderHelper * \brief gets the info from 0028,0030 : Pixel Spacing - * \ else 1. + * else 1.0 * @return X dimension of a pixel */ float gdcmHeaderHelper::GetXSpacing() { @@ -181,8 +128,11 @@ float gdcmHeaderHelper::GetXSpacing() { dbg.Verbose(0, "gdcmHeader::GetXSpacing: unfound Pixel Spacing (0028,0030)"); return 1.; } - if( sscanf( StrSpacing.c_str(), "%f\\%f", &yspacing, &xspacing) != 2) - return 0.; + int nbValues; + if( (nbValues = sscanf( StrSpacing.c_str(), "%f\\%f", &yspacing, &xspacing)) != 2) { + if (nbValues==1) // if single value is found, xspacing is defaulted to yspacing + return yspacing; + } if (xspacing == 0.) { dbg.Verbose(0, "gdcmHeader::GetYSpacing: gdcmData/CT-MONO2-8-abdo.dcm problem"); // seems to be a bug in the header ... @@ -194,32 +144,26 @@ float gdcmHeaderHelper::GetXSpacing() { /** * \ingroup gdcmHeaderHelper * \brief gets the info from 0028,0030 : Pixel Spacing - * \ else 1. + * else 1.0 * @return Y dimension of a pixel */ float gdcmHeaderHelper::GetYSpacing() { - float xspacing, yspacing; + float yspacing; std::string StrSpacing = GetEntryByNumber(0x0028,0x0030); if (StrSpacing == GDCM_UNFOUND) { dbg.Verbose(0, "gdcmHeader::GetYSpacing: unfound Pixel Spacing (0028,0030)"); return 1.; } - if( sscanf( StrSpacing.c_str(), "%f\\%f", &yspacing, &xspacing) != 2) - return 0.; - if (xspacing == 0.) { - dbg.Verbose(0, "gdcmHeader::GetYSpacing: gdcmData/CT-MONO2-8-abdo.dcm problem"); - // seems to be a bug in the header ... - sscanf( StrSpacing.c_str(), "%f\\0\\%f", &yspacing, &xspacing); - } + sscanf( StrSpacing.c_str(), "%f", &yspacing); return yspacing; } /** *\ingroup gdcmHeaderHelper *\brief gets the info from 0018,0088 : Space Between Slices - *\ else from 0018,0050 : Slice Thickness - *\ else 1. + * else from 0018,0050 : Slice Thickness + * else 1.0 * @return Z dimension of a voxel-to be */ float gdcmHeaderHelper::GetZSpacing() { @@ -288,8 +232,8 @@ float gdcmHeaderHelper::GetRescaleIntercept() { /** * \ingroup gdcmHeaderHelper * \brief This function is intended to user who doesn't want - * \ to have to manage a LUT and expects to get an RBG Pixel image - * \ (or a monochrome one ...) + * to have to manage a LUT and expects to get an RBG Pixel image + * (or a monochrome one ...) * \warning to be used with GetImagePixels() * @return 1 if Gray level, 3 if Color (RGB, YBR or PALETTE COLOR) */ @@ -328,8 +272,8 @@ int gdcmHeaderHelper::GetNumberOfScalarComponents() { /** * \ingroup gdcmHeaderHelper * \brief This function is intended to user that DOESN'T want - * \to get RGB pixels image when it's stored as a PALETTE COLOR image - * \ - the (vtk) user is supposed to know how deal with LUTs - + * to get RGB pixels image when it's stored as a PALETTE COLOR image + * - the (vtk) user is supposed to know how deal with LUTs - * \warning to be used with GetImagePixelsRaw() * @return 1 if Gray level, 3 if Color (RGB or YBR - NOT 'PALETTE COLOR' -) */ @@ -375,7 +319,6 @@ int gdcmHeaderHelper::GetNumberOfScalarComponentsRaw() { } /** - *\ingroup gdcmHeaderHelper *\brief gets the info from 0008,0018 : SOP Instance UID *\todo ? : return the ACR-NEMA element value if DICOM one is not found * @return SOP Instance UID @@ -393,17 +336,16 @@ int gdcmHeaderHelper::GetNumberOfScalarComponentsRaw() { // as the Z coordinate, // 0. for all the coordinates if nothing is found -// TODO : find a way to inform the caller nothing was found -// TODO : How to tell the caller a wrong number of values was found? +// \todo find a way to inform the caller nothing was found +// \todo How to tell the caller a wrong number of values was found? // // --------------------------------------------------------------- // /** - * \ingroup gdcmHeaderHelper * \brief gets the info from 0020,0032 : Image Position Patient - *\ else from 0020,0030 : Image Position (RET) - *\ else 0. + * else from 0020,0030 : Image Position (RET) + * else 0. * @return up-left image corner X position */ @@ -416,7 +358,7 @@ float gdcmHeaderHelper::GetXOrigin() { StrImPos = GetEntryByNumber(0x0020,0x0030); // For ACR-NEMA images if (StrImPos == GDCM_UNFOUND) { dbg.Verbose(0, "gdcmHeader::GetXImagePosition: unfound Image Position (RET) (0020,0030)"); - // How to tell the caller nothing was found ? + /// \todo How to tell the caller nothing was found ? return 0.; } } @@ -426,10 +368,9 @@ float gdcmHeaderHelper::GetXOrigin() { } /** - * \ingroup gdcmHeaderHelper * \brief gets the info from 0020,0032 : Image Position Patient - * \ else from 0020,0030 : Image Position (RET) - * \ else 0. + * else from 0020,0030 : Image Position (RET) + * else 0. * @return up-left image corner Y position */ float gdcmHeaderHelper::GetYOrigin() { @@ -441,7 +382,7 @@ float gdcmHeaderHelper::GetYOrigin() { StrImPos = GetEntryByNumber(0x0020,0x0030); // For ACR-NEMA images if (StrImPos == GDCM_UNFOUND) { dbg.Verbose(0, "gdcmHeader::GetYImagePosition: unfound Image Position (RET) (0020,0030)"); - // How to tell the caller nothing was found ? + /// \todo How to tell the caller nothing was found ? return 0.; } } @@ -451,7 +392,6 @@ float gdcmHeaderHelper::GetYOrigin() { } /** - * \ingroup gdcmHeaderHelper * \brief gets the info from 0020,0032 : Image Position Patient * \ else from 0020,0030 : Image Position (RET) * \ else from 0020,1041 : Slice Location @@ -503,17 +443,17 @@ float gdcmHeaderHelper::GetZOrigin() { } /** - * \ingroup gdcmHeaderHelper * \brief gets the info from 0020,0013 : Image Number * \ else 0. * @return image number */ int gdcmHeaderHelper::GetImageNumber() { - //The function i atoi() takes the address of an area of memory as parameter and converts - //the string stored at that location to an integer using the external decimal to internal - //binary conversion rules. This may be preferable to sscanf() since atoi() is a much smaller, - // simpler and faster function. sscanf() can do all possible conversions whereas atoi() can - //only do single decimal integer conversions. + // The function i atoi() takes the address of an area of memory as + // parameter and converts the string stored at that location to an integer + // using the external decimal to internal binary conversion rules. This may + // be preferable to sscanf() since atoi() is a much smaller, simpler and + // faster function. sscanf() can do all possible conversions whereas + // atoi() can only do single decimal integer conversions. std::string StrImNumber = GetEntryByNumber(0x0020,0x0013); //0020 0013 IS REL Image Number if (StrImNumber != GDCM_UNFOUND) { return atoi( StrImNumber.c_str() ); @@ -522,12 +462,12 @@ int gdcmHeaderHelper::GetImageNumber() { } /** - * \ingroup gdcmHeaderHelper * \brief gets the info from 0008,0060 : Modality * @return Modality Type */ ModalityType gdcmHeaderHelper::GetModality(void) { - std::string StrModality = GetEntryByNumber(0x0008,0x0060); //0008 0060 CS ID Modality + // 0008 0060 CS ID Modality + std::string StrModality = GetEntryByNumber(0x0008,0x0060); if (StrModality != GDCM_UNFOUND) { if ( StrModality.find("AU") < StrModality.length()) return AU; else if ( StrModality.find("AS") < StrModality.length()) return AS; @@ -572,8 +512,8 @@ ModalityType gdcmHeaderHelper::GetModality(void) { else { - //throw error return value ??? - // specified <> unknow in our database + /// \todo throw error return value ??? + /// specified <> unknow in our database return Unknow; } } @@ -581,7 +521,6 @@ ModalityType gdcmHeaderHelper::GetModality(void) { } /** - * \ingroup gdcmHeaderHelper * \brief gets the info from 0020,0037 : Image Orientation Patient * @param iop adress of the (6)float aray to receive values * @return cosines of image orientation patient @@ -591,7 +530,8 @@ void gdcmHeaderHelper::GetImageOrientationPatient( float* iop ) { //iop is supposed to be float[6] iop[0] = iop[1] = iop[2] = iop[3] = iop[4] = iop[5] = 0; - std::string StrImOriPat = GetEntryByNumber(0x0020,0x0037); // 0020 0037 DS REL Image Orientation (Patient) + // 0020 0037 DS REL Image Orientation (Patient) + std::string StrImOriPat = GetEntryByNumber(0x0020,0x0037); if (StrImOriPat != GDCM_UNFOUND) { if( sscanf( StrImOriPat.c_str(), "%f\\%f\\%f\\%f\\%f\\%f", &iop[0], &iop[1], &iop[2], &iop[3], &iop[4], &iop[5]) != 6) { @@ -603,7 +543,8 @@ void gdcmHeaderHelper::GetImageOrientationPatient( float* iop ) { } //For ACR-NEMA - StrImOriPat = GetEntryByNumber(0x0020,0x0035); //0020 0035 DS REL Image Orientation (RET) + // 0020 0035 DS REL Image Orientation (RET) + StrImOriPat = GetEntryByNumber(0x0020,0x0035); if (StrImOriPat != GDCM_UNFOUND) { if( sscanf( StrImOriPat.c_str(), "%f\\%f\\%f\\%f\\%f\\%f", &iop[0], &iop[1], &iop[2], &iop[3], &iop[4], &iop[5]) != 6) { @@ -625,12 +566,10 @@ void gdcmHeaderHelper::GetImageOrientationPatient( float* iop ) { -//----------------------------------------------------------------------------- -// gdcmSerieHeaderHelper //----------------------------------------------------------------------------- // Constructor / Destructor gdcmSerieHeaderHelper::~gdcmSerieHeaderHelper(){ - //! \todo + /// \todo for (std::list::iterator it = CoherentGdcmFileList.begin(); it != CoherentGdcmFileList.end(); it++) { @@ -645,7 +584,6 @@ gdcmSerieHeaderHelper::~gdcmSerieHeaderHelper(){ //----------------------------------------------------------------------------- // Public /** - * \ingroup gdcmHeaderHelper * \brief add a gdcmFile to the list based on file name * @param filename Name of the file to deal with */ @@ -655,7 +593,6 @@ void gdcmSerieHeaderHelper::AddFileName(std::string filename) { } /** - * \ingroup gdcmHeaderHelper * \brief add a gdcmFile to the list * @param file gdcmHeaderHelper to add */ @@ -664,7 +601,6 @@ void gdcmSerieHeaderHelper::AddGdcmFile(gdcmHeaderHelper *file){ } /** - * \ingroup gdcmHeaderHelper * \brief Sets the Directory * @param dir Name of the directory to deal with */ @@ -680,7 +616,6 @@ void gdcmSerieHeaderHelper::SetDirectory(std::string dir){ } /** - * \ingroup gdcmHeaderHelper * \brief Sorts the File List * \warning This could be implemented in a 'Strategy Pattern' approach * But as I don't know how to do it, I leave it this way @@ -698,7 +633,6 @@ void gdcmSerieHeaderHelper::OrderGdcmFileList(){ } /** - * \ingroup gdcmHeaderHelper * \brief Gets the *coherent* File List * @return the *coherent* File List */ @@ -836,7 +770,7 @@ bool gdcmSerieHeaderHelper::ImagePositionPatientOrdering() bool gdcmSerieHeaderHelper::ImageNumberOrdering() { int min, max, pos; - int n = 0;//CoherentGdcmFileList.size(); //O(N) operation !! + int n = 0;//CoherentGdcmFileList.size() is a O(N) operation !! unsigned char *partition; std::list::iterator it = CoherentGdcmFileList.begin();