X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=Example%2FWriteDicomAsJPEG.cxx;h=b8d52fbc2371ec7142346c44624329428eb244cd;hb=5f2a9e86acc7c14559f0c692af86f4e92ccf8a87;hp=c2d74a02ba3367482a6650159b0f6d05b57ac785;hpb=a632cce0349ca495574f98c5b338a34f581943e4;p=gdcm.git diff --git a/Example/WriteDicomAsJPEG.cxx b/Example/WriteDicomAsJPEG.cxx index c2d74a02..b8d52fbc 100644 --- a/Example/WriteDicomAsJPEG.cxx +++ b/Example/WriteDicomAsJPEG.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: WriteDicomAsJPEG.cxx,v $ Language: C++ - Date: $Date: 2007/08/21 15:10:04 $ - Version: $Revision: 1.15 $ + Date: $Date: 2007/08/29 08:13:40 $ + Version: $Revision: 1.20 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -19,13 +19,15 @@ #include "gdcmFile.h" #include "gdcmFileHelper.h" #include "gdcmUtil.h" +#include "gdcmDebug.h" // Open a dicom file and compress it as JPEG stream int main(int argc, char *argv[]) { if( argc < 2) { - std::cerr << argv[0] << " inputfilename.dcm\n"; + std::cerr << argv[0] << " inputfilename.dcm [ outputfilename.dcm" + << " quality debug]\n"; return 1; } @@ -37,8 +39,11 @@ int main(int argc, char *argv[]) if( argc >= 4 ) quality = atoi(argv[3]); std::cerr << "Using quality: " << quality << std::endl; - -// Step 1 : Create the header of the image + + if (argc > 4) + GDCM_NAME_SPACE::Debug::DebugOn(); + +// Step 1 : Read the image GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New(); f->SetLoadMode ( GDCM_NAME_SPACE::LD_ALL ); // Load everything f->SetFileName( filename ); @@ -52,12 +57,14 @@ int main(int argc, char *argv[]) //tested->Print( std::cout ); int samplesPerPixel = f->GetSamplesPerPixel(); - size_t testedDataSize = tested->GetImageDataSize(); - std::cerr << "testedDataSize:" << testedDataSize << std::endl; - uint8_t *testedImageData = tested->GetImageData(); - -// Step 1 : Create the header of the image - + size_t testedDataSize = tested->GetImageDataRawSize(); // Raw : Don't convert gray pixels+LUT to RBG pixels + uint8_t *testedImageData = tested->GetImageDataRaw(); + + if( GDCM_NAME_SPACE::Debug::GetDebugFlag() ) { + tested->Print( std::cout ); + std::cout << "-------------------------------------------------------------------------------" << std::endl; + } +// Step 1 : Create the header of the new file GDCM_NAME_SPACE::File *fileToBuild = GDCM_NAME_SPACE::File::New(); std::ostringstream str; @@ -68,7 +75,6 @@ int main(int argc, char *argv[]) str.str(""); str << ysize; fileToBuild->InsertEntryString(str.str(),0x0028,0x0010,"US"); // Rows - if(zsize>1) { @@ -76,9 +82,10 @@ int main(int argc, char *argv[]) str << zsize; fileToBuild->InsertEntryString(str.str(),0x0028,0x0008,"IS"); // Number of Frames } + int bitsallocated = f->GetBitsAllocated(); - int bitsstored = f->GetBitsStored(); - int highbit = f->GetHighBitPosition(); + int bitsstored = f->GetBitsStored(); + int highbit = f->GetHighBitPosition(); //std::string pixtype = f->GetPixelType(); int sign = f->IsSignedPixelData(); @@ -104,19 +111,56 @@ int main(int argc, char *argv[]) str << samplesPerPixel; //img.components; fileToBuild->InsertEntryString(str.str(),0x0028,0x0002,"US"); // Samples per Pixel +// The image may be displayed uncorectly if these fields are missing + + // Set the Pixel Aspect Ratio, if any + std::string par = f->GetEntryString(0x0028,0x0034); + std::cerr <<"Pixel Aspect Ratio [" << par << "]" << std::endl; + if ( par != GDCM_NAME_SPACE::GDCM_UNFOUND ) + fileToBuild->InsertEntryString(par,0x0028,0x0034,"IS"); // Pixel Aspect Ratio + + // Set the Modality, if any + std::string modality = f->GetEntryString(0x0008,0x0060); + std::cerr <<"Modality [" << modality << "]" << std::endl; + if ( modality != GDCM_NAME_SPACE::GDCM_UNFOUND ) + fileToBuild->InsertEntryString(modality,0x0008,0x0060,"CS"); // Modality + + // Set the Media Storage SOP Class UID, if any + std::string mssop = f->GetEntryString(0x0002,0x0002); + std::cerr <<"Media Storage SOP Class UID [" << mssop << "]" << std::endl; + if ( mssop != GDCM_NAME_SPACE::GDCM_UNFOUND ) + fileToBuild->InsertEntryString(mssop,0x0002,0x0002,"UI"); // Media Storage SOP Class UID + + // This one is mandatory to deal with Pixel Aspect Ratio, in ultrasound images ! + // Set the SOP Class UID, if any + std::string sop = f->GetEntryString(0x0008,0x0016); + std::cerr <<"SOP Class UID [" << sop << "]" << std::endl; + if ( sop != GDCM_NAME_SPACE::GDCM_UNFOUND ) + fileToBuild->InsertEntryString(sop,0x0008,0x0016,"UI"); // SOP Class UID + // Step 2 : Create the output image size_t size = xsize * ysize * zsize * samplesPerPixel * bitsallocated / 8; GDCM_NAME_SPACE::FileHelper *fileH = GDCM_NAME_SPACE::FileHelper::New(fileToBuild); - assert( size == testedDataSize ); + + // Consider that pixels are unmodified + fileH->SetContentType(GDCM_NAME_SPACE::UNMODIFIED_PIXELS_IMAGE); + std::cerr << "xsize " << xsize << " ysize " << ysize << " zsize " << zsize << " samplesPerPixel " << samplesPerPixel + << " bitsallocated " << bitsallocated << std::endl; + std::cerr << "size " << size << " testedDataSize " << testedDataSize << + std::endl; + assert(abs (size-testedDataSize) <= 1 ); fileH->SetWriteTypeToJPEG( ); - fileH->SetImageData(testedImageData, testedDataSize); + //fileH->SetImageData(testedImageData, testedDataSize); + + // SetUserData will ensure the compression + fileH->SetUserData(testedImageData, testedDataSize); if( !fileH->Write(outfilename) ) { std::cerr << "write fails" << std::endl; } - + f->Delete(); tested->Delete(); fileToBuild->Delete();