X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FPrintFile.cxx;h=771f07c7058bde5c5909f535899515a3f08b7c57;hb=f45f440f9ba41ad0a37c23a93c1c8f5705c511a6;hp=9ff3a681d6ebd695db02b19f5dc24b2a9665d768;hpb=4d688a4874380bbde4dff31775ab4d56d49f7f9e;p=gdcm.git diff --git a/Example/PrintFile.cxx b/Example/PrintFile.cxx index 9ff3a681..771f07c7 100644 --- a/Example/PrintFile.cxx +++ b/Example/PrintFile.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: PrintFile.cxx,v $ Language: C++ - Date: $Date: 2005/06/17 14:36:19 $ - Version: $Revision: 1.40 $ + Date: $Date: 2005/07/07 17:31:53 $ + Version: $Revision: 1.45 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -16,6 +16,10 @@ =========================================================================*/ #include "gdcmFile.h" +#include "gdcmSeqEntry.h" +#include "gdcmSQItem.h" +#include "gdcmBinEntry.h" + #include "gdcmFileHelper.h" #include "gdcmDebug.h" @@ -23,24 +27,115 @@ #include + +void ShowLutData(gdcm::File *f); + +void ShowLutData(gdcm::File *f) +{ + // Nothing is written yet to get LUT Data user friendly + // The following is to be moved into a PixelReadConvert method + // Let here, waiting for a clever idea on the way to do it. + + gdcm::SeqEntry *modLutSeq = f->GetSeqEntry(0x0028,0x3000); + if ( modLutSeq !=0 ) + { + gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem(); + if ( sqi != 0 ) + { + std::string lutDescriptor = sqi->GetEntryValue(0x0028,0x3002); + if ( /*lutDescriptor == GDCM_UNFOUND*/ 0 ) + { + //gdcmWarningMacro( "LUT Descriptor is missing" ); + std::cout << "LUT Descriptor is missing" << std::endl; + return; + } + int length; // LUT length in Bytes + int deb; // Subscript of the first Lut Value + int nbits; // Lut item size (in Bits) + + int nbRead; // nb of items in LUT descriptor (must be = 3) + + nbRead = sscanf( lutDescriptor.c_str(), + "%d\\%d\\%d", + &length, &deb, &nbits ); + std::cout << "length " << length + << " deb " << deb + << " nbits " << nbits + << std::endl; + if ( nbRead != 3 ) + { + //gdcmWarningMacro( "Wrong LUT descriptor" ); + std::cout << "Wrong LUT descriptor" << std::endl; + } + //LUT Data (CTX dependent) + gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006); + if ( b != 0 ) + { + int BitsAllocated = f->GetBitsAllocated(); + if ( BitsAllocated <= 8 ) + { + int mult; + if ( ( nbits == 16 ) && ( BitsAllocated == 8 ) ) + { + // when LUT item size is different than pixel size + mult = 2; // high byte must be = low byte + } + else + { + // See PS 3.3-2003 C.11.1.1.2 p 619 + mult = 1; + } + uint8_t *lut = b->GetBinArea(); + for( int i=0; i < length; ++i ) + { + std::cout << i+deb << " : \t" + << (int) (lut[i*mult + 1]) << std::endl; + } + } + else + { + uint16_t *lut = (uint16_t *)(b->GetBinArea()); + for( int i=0; i < length; ++i ) + { + std::cout << i+deb << " : \t" + << (int) (((uint16_t *)lut)[i]) + << std::endl; + } + } + } + else + std::cout << "No LUT Data BinEntry (0x0028,0x3006) found?!? " + << std::endl; + } + else + std::cout << "No First SQ Item within (0x0028,0x3000) ?!? " + << std::endl; + } + else + std::cout << "No LUT Data SeqEntry (0x0028,0x3000) found " + << std::endl; + + } + int main(int argc, char *argv[]) { START_USAGE(usage) " \n PrintFile : \n", " Display the header of a ACR-NEMA/PAPYRUS/DICOM File", - " usage: PrintFile filein=fileName [level=n] [noshadow] [noseq] [debug] ", + " usage: PrintFile filein=fileName [level=n] [noshadowseq][noshadow][noseq][debug] ", " level = 0,1,2 : depending on the amount of details user wants to see", - " noshadow : user doesn't want to load Private groups (odd number)", - " noseq : user doesn't want to load Sequences ", - " debug : user wants to run the program in 'debug mode' ", - " showlut : user wants to display the Palette Color (as an int array)", + " noshadowseq: user doesn't want to load Private Sequences", + " noshadow : user doesn't want to load Private groups (odd number)", + " noseq : user doesn't want to load Sequences ", + " debug : user wants to run the program in 'debug mode' ", + " showlut :user wants to display the Palette Color (as an int array)", FINISH_USAGE // Initialize Arguments Manager gdcm::ArgMgr *am= new gdcm::ArgMgr(argc, argv); - if (argc == 1) + if (argc == 1 || am->ArgMgrDefined("usage") ) { am->ArgMgrUsage(usage); // Display 'usage' delete am; @@ -49,26 +144,27 @@ int main(int argc, char *argv[]) char *fileName = am->ArgMgrWantString("filein",usage); - int loadMode; - if ( am->ArgMgrDefined("noshadow") && am->ArgMgrDefined("noseq") ) - loadMode = NO_SEQ | NO_SHADOW; - else if ( am->ArgMgrDefined("noshadow") ) - loadMode = NO_SHADOW; - else if ( am->ArgMgrDefined("noseq") ) - loadMode = NO_SEQ; - else - loadMode = 0; + int loadMode = 0x00000000; + if ( am->ArgMgrDefined("noshadowseq") ) + loadMode |= NO_SHADOWSEQ; + else + { + if ( am->ArgMgrDefined("noshadow") ) + loadMode |= NO_SHADOW; + if ( am->ArgMgrDefined("noseq") ) + loadMode |= NO_SEQ; + } int level = am->ArgMgrGetInt("level", 2); + bool showlut = ( 0 != am->ArgMgrDefined("SHOWLUT") ); + if (am->ArgMgrDefined("debug")) gdcm::Debug::DebugOn(); - - bool showlut = ( 0 != am->ArgMgrDefined("SHOWLUT") ); /* if unused Param we give up */ if ( am->ArgMgrPrintUnusedLabels() ) - { + { am->ArgMgrUsage(usage); delete am; return 0; @@ -82,53 +178,54 @@ int main(int argc, char *argv[]) // any kind of gdcm-Parsable *document* // not only gdcm::File (as opposed to gdcm::DicomDir) - gdcm::File *e1 = new gdcm::File(); - e1->SetLoadMode(loadMode); + gdcm::File *f = new gdcm::File(); + f->SetLoadMode(loadMode); + f->SetFileName( fileName ); + bool res = f->Load(); // just to see - bool res = e1->Load( fileName ); if ( !res ) { - delete e1; + delete f; return 0; } - gdcm::FileHelper *f1 = new gdcm::FileHelper(e1); - f1->SetPrintLevel( level ); + gdcm::FileHelper *fh = new gdcm::FileHelper(f); + fh->SetPrintLevel( level ); - f1->Print(); + fh->Print(); std::cout << "\n\n" << std::endl; std::cout <GetImageDataSize() << std::endl; - std::cout <<" dataSizeRaw " << f1->GetImageDataRawSize() << std::endl; + std::cout <<" dataSize " << fh->GetImageDataSize() << std::endl; + std::cout <<" dataSizeRaw " << fh->GetImageDataRawSize() << std::endl; int nX,nY,nZ,sPP,planarConfig; std::string pixelType; - nX=e1->GetXSize(); - nY=e1->GetYSize(); - nZ=e1->GetZSize(); + nX=f->GetXSize(); + nY=f->GetYSize(); + nZ=f->GetZSize(); std::cout << " DIMX=" << nX << " DIMY=" << nY << " DIMZ=" << nZ << std::endl; - pixelType = e1->GetPixelType(); - sPP = e1->GetSamplesPerPixel(); - planarConfig = e1->GetPlanarConfiguration(); + pixelType = f->GetPixelType(); + sPP = f->GetSamplesPerPixel(); + planarConfig = f->GetPlanarConfiguration(); std::cout << " pixelType= [" << pixelType << "] SamplesPerPixel= [" << sPP << "] PlanarConfiguration= [" << planarConfig << "] "<< std::endl << " PhotometricInterpretation= [" - << e1->GetEntryValue(0x0028,0x0004) + << f->GetEntryValue(0x0028,0x0004) << "] "<< std::endl; - int numberOfScalarComponents=e1->GetNumberOfScalarComponents(); + int numberOfScalarComponents=f->GetNumberOfScalarComponents(); std::cout << " NumberOfScalarComponents = " << numberOfScalarComponents <GetEntryValue(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) + if ( f->GetEntryValue(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) { std::cout << "Transfer Syntax not loaded. " << std::endl << "Better you increase MAX_SIZE_LOAD_ELEMENT_VALUE" @@ -136,45 +233,96 @@ int main(int argc, char *argv[]) return 0; } - std::string transferSyntaxName = e1->GetTransferSyntaxName(); + std::string transferSyntaxName = f->GetTransferSyntaxName(); std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" << std::endl; - std::cout << " SwapCode= " << e1->GetSwapCode() << std::endl; - - // Display the LUT as an int array (for degugging purpose) - if ( e1->HasLUT() && showlut ) + std::cout << " SwapCode= " << f->GetSwapCode() << std::endl; + + //std::cout << "\n\n" << std::endl; + //std::cout << "X spacing " << f->GetXSpacing() << std::endl; + //std::cout << "Y spacing " << f->GetYSpacing() << std::endl; + //std::cout << "Z spacing " << f->GetZSpacing() << std::endl; + + // Display the LUT as an int array (for debugging purpose) + if ( f->HasLUT() && showlut ) { - uint8_t* lutrgba = f1->GetLutRGBA(); + uint8_t* lutrgba = fh->GetLutRGBA(); if ( lutrgba == 0 ) { - std::cout << "Lut RGBA not built ?!?" << std::endl; - } + std::cout << "Lut RGBA (Palette Color) not built " << std::endl; + + // Nothing is written yet to get LUT Data user friendly + // The following is to be moved into a PixelRedaConvert method + + gdcm::SeqEntry *modLutSeq = f->GetSeqEntry(0x0028,0x3000); + if ( modLutSeq !=0 ) + { + gdcm::SQItem *sqi= modLutSeq->GetFirstSQItem(); + if ( !sqi ) + { + std::string lutDescriptor = sqi->GetEntryValue(0x0028,0x3002); + int length; // LUT length in Bytes + int deb; // Subscript of the first Lut Value + int nbits; // Lut item size (in Bits) + int nbRead; // nb of items in LUT descriptor (must be = 3) + + nbRead = sscanf( lutDescriptor.c_str(), + "%d\\%d\\%d", + &length, &deb, &nbits ); + if ( nbRead != 3 ) + { + //gdcmWarningMacro( "Wrong LUT descriptor" ); + std::cout << "Wrong LUT descriptor" << std::endl; + } + gdcm::BinEntry *b = sqi->GetBinEntry(0x0028,0x3006); + if ( b != 0 ) + { + if ( b->GetLength() != 0 ) + { + std::cout << "---------------------------------------" + << " We should never reach this point " + << std::endl; + //LoadEntryBinArea(b); //LUT Data (CTX dependent) + } + } + } + } + else + std::cout << "No LUT Data (0x0028,0x3000) found " << std::endl; + } else { - if ( f1->GetLutItemSize() == 8 ) + if ( fh->GetLutItemSize() == 8 ) { - for (int i=0;iGetLutItemNumber();i++) - std::cout << (int)(lutrgba[i*4]) << " " + for (int i=0;iGetLutItemNumber();i++) + std::cout << i << " : \t" + << (int)(lutrgba[i*4]) << " " << (int)(lutrgba[i*4+1]) << " " << (int)(lutrgba[i*4+2]) << std::endl; } else // LutItemSize assumed to be = 16 { uint16_t* lutrgba16 = (uint16_t*)lutrgba; - for (int i=0;iGetLutItemNumber();i++) - std::cout << lutrgba16[i*4] << " " - << lutrgba16[i*4+1] << " " - << lutrgba16[i*4+2] << std::endl; + for (int i=0;iGetLutItemNumber();i++) + std::cout << i << " : \t" + << (int)(lutrgba16[i*4]) << " " + << (int)(lutrgba16[i*4+1]) << " " + << (int)(lutrgba16[i*4+2]) << std::endl; } } } + else if (showlut) + { + std::cout << "Try LUT Data "<< std::endl; + ShowLutData(f); + } - if(e1->IsReadable()) + if (f->IsReadable()) std::cout <