X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FPrintFile.cxx;h=771f07c7058bde5c5909f535899515a3f08b7c57;hb=f45f440f9ba41ad0a37c23a93c1c8f5705c511a6;hp=ad890128ab1b86909b3487b072a98ad865fbab8b;hpb=ab62b1a706df8e2ab1fbf26fa5037d6c22293655;p=gdcm.git diff --git a/Example/PrintFile.cxx b/Example/PrintFile.cxx index ad890128..771f07c7 100644 --- a/Example/PrintFile.cxx +++ b/Example/PrintFile.cxx @@ -1,80 +1,328 @@ +/*========================================================================= + + Program: gdcm + Module: $RCSfile: PrintFile.cxx,v $ + Language: C++ + 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 + http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details. + + This software is distributed WITHOUT ANY WARRANTY; without even + the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR + PURPOSE. See the above copyright notices for more information. + +=========================================================================*/ +#include "gdcmFile.h" +#include "gdcmSeqEntry.h" +#include "gdcmSQItem.h" +#include "gdcmBinEntry.h" + +#include "gdcmFileHelper.h" +#include "gdcmDebug.h" + +#include "gdcmArgMgr.h" + #include -#include "gdcm.h" -int main(int argc, char* argv[]) + +void ShowLutData(gdcm::File *f); + +void ShowLutData(gdcm::File *f) { - - gdcmHeader *e1; - gdcmFile *f1; - std::string fileName; - if (argc != 2) { - std::cout << " usage : PrintDocument fileName" << std::endl; + // 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; + } - if (argc > 1) { - fileName=argv[1]; - } else { - fileName += GDCM_DATA_ROOT; - fileName += "/test.acr"; +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] [noshadowseq][noshadow][noseq][debug] ", + " level = 0,1,2 : depending on the amount of details user wants to see", + " 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 || am->ArgMgrDefined("usage") ) + { + am->ArgMgrUsage(usage); // Display 'usage' + delete am; + return 0; } - - e1= new gdcmHeader( fileName.c_str() ); - f1 = new gdcmFile(e1); + char *fileName = am->ArgMgrWantString("filein",usage); + + 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(); + + /* if unused Param we give up */ + if ( am->ArgMgrPrintUnusedLabels() ) + { + am->ArgMgrUsage(usage); + delete am; + return 0; + } + + delete am; // we don't need Argument Manager any longer + + // ----------- End Arguments Manager --------- + + // gdcm::File::IsReadable() is no usable here, because we deal with + // any kind of gdcm-Parsable *document* + // not only gdcm::File (as opposed to gdcm::DicomDir) + + gdcm::File *f = new gdcm::File(); + f->SetLoadMode(loadMode); + f->SetFileName( fileName ); + bool res = f->Load(); // just to see + + if ( !res ) + { + delete f; + return 0; + } + + gdcm::FileHelper *fh = new gdcm::FileHelper(f); + fh->SetPrintLevel( level ); + + fh->Print(); - e1->SetPrintLevel(2); - - e1->Print(); - std::cout << "\n\n" << std::endl; - int dataSize = f1->GetImageDataSize(); - std::cout <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 - << " SampleserPixel=" << sPP - << " PlanarConfiguration=" << planarConfig - << std::endl - << " PhotometricInterpretation=" - << e1->GetEntryByNumber(0x0028,0x0004) + std::cout << " pixelType= [" << pixelType + << "] SamplesPerPixel= [" << sPP + << "] PlanarConfiguration= [" << planarConfig + << "] "<< std::endl + << " PhotometricInterpretation= [" + << f->GetEntryValue(0x0028,0x0004) + << "] "<< std::endl; + + int numberOfScalarComponents=f->GetNumberOfScalarComponents(); + std::cout << " NumberOfScalarComponents = " << numberOfScalarComponents + <GetEntryValue(0x0002,0x0010) == gdcm::GDCM_NOTLOADED ) + { + std::cout << "Transfer Syntax not loaded. " << std::endl + << "Better you increase MAX_SIZE_LOAD_ELEMENT_VALUE" + << std::endl; + return 0; + } + + std::string transferSyntaxName = f->GetTransferSyntaxName(); + std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" << std::endl; + std::cout << " SwapCode= " << f->GetSwapCode() << std::endl; - int numberOfScalarComponents=e1->GetNumberOfScalarComponents(); - std::cout << " NumberOfScalarComponents " << numberOfScalarComponents <GetXSpacing() << std::endl; + //std::cout << "Y spacing " << f->GetYSpacing() << std::endl; + //std::cout << "Z spacing " << f->GetZSpacing() << std::endl; - std::string transferSyntaxName = e1->GetTransfertSyntaxName(); - std::cout << " TransferSyntaxName= [" << transferSyntaxName << "]" << std::endl; - - if ( transferSyntaxName != "Implicit VR - Little Endian" - && transferSyntaxName != "Explicit VR - Little Endian" - && transferSyntaxName != "Deflated Explicit VR - Little Endian" - && transferSyntaxName != "Explicit VR - Big Endian" - && transferSyntaxName != "Uncompressed ACR-NEMA" ) - { - std::cout << std::endl << "===========================================" - << std::endl; - f1->ParsePixelData(); - std::cout << std::endl << "===========================================" - << std::endl; - } - - if(e1->IsReadable()) + // Display the LUT as an int array (for debugging purpose) + if ( f->HasLUT() && showlut ) + { + uint8_t* lutrgba = fh->GetLutRGBA(); + if ( lutrgba == 0 ) + { + 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 ( fh->GetLutItemSize() == 8 ) + { + 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 << 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 (f->IsReadable()) std::cout <