/*========================================================================= Program: gdcm Module: $RCSfile: exColorToRGB.cxx,v $ Language: C++ Date: $Date: 2007/05/23 14:18:05 $ Version: $Revision: 1.10 $ 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 "gdcmFileHelper.h" #include "gdcmDocument.h" #include "gdcmDataEntry.h" #include "gdcmSeqEntry.h" #include // for exit int main(int argc, char *argv[]) { std::cout << "------------------------------------------------" << std::endl; std::cout << "Tranforms a full gdcm-readable 'color' Dicom image " << " (e.g Palette Color, YBR, 3-Planes RGB) " << " into an 'RGB_Pixels' Dicom Image " << " " << std::endl; if (argc < 3) { std::cerr << "Usage :" << std::endl << argv[0] << " input_dicom output_dicom" << std::endl; return 1; } std::string fileName = argv[1]; std::string output = argv[2]; // ============================================================ // Read the input image. // ============================================================ // a GDCM_NAME_SPACE::File contains all the Dicom Field but the Pixels Element std::cout << argv[1] << std::endl; GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New(); f->SetLoadMode( GDCM_NAME_SPACE::LD_ALL); f->SetFileName( fileName ); bool res = f->Load(); if (!res) { std::cerr << "Sorry, " << fileName <<" not a gdcm-readable " << "DICOM / ACR File" <Delete(); return 0; } std::cout << " ... is readable " << std::endl; /* if (!f->IsMonochrome()) { std::cerr << "Sorry, " << fileName <<" not a 'color' File " << " " <Delete(); return 0; } */ // ============================================================ // Load the pixels in memory. // ============================================================ // We need a GDCM_NAME_SPACE::FileHelper, since we want to load the pixels GDCM_NAME_SPACE::FileHelper *fh = GDCM_NAME_SPACE::FileHelper::New(f); // uint8_t DOESN'T mean it's mandatory for the image to be a 8 bits one ! // It's just for prototyping. // Feel free to cast it. uint8_t *imageData = fh->GetImageData(); if ( imageData == 0 ) { std::cerr << "Sorry, Pixels of" << fileName <<" are not " << " gdcm-readable." << std::endl; f->Delete(); return 0; } // ------ User wants write a new image without shadow groups ------------- // ------ without Sequences ------------- GDCM_NAME_SPACE::FileHelper *copy = GDCM_NAME_SPACE::FileHelper::New( ); copy->SetFileName( output ); copy->Load(); GDCM_NAME_SPACE::DocEntry *d = f->GetFirstEntry(); while(d) { // We skip SeqEntries, since user cannot do much with them if ( !(dynamic_cast(d)) // We skip Shadow Groups, since nobody knows what they mean && !( d->GetGroup()%2 ) ) { if ( GDCM_NAME_SPACE::DataEntry *de = dynamic_cast(d) ) { copy->GetFile()->InsertEntryBinArea( de->GetBinArea(),de->GetLength(), de->GetGroup(),de->GetElement(), de->GetVR() ); } else { // We skip GDCM_NAME_SPACE::SeqEntries } } d = f->GetNextEntry(); } // User knows the image is a 'color' one -RGB, YBR, Palette Color- // and wants to write it as RGB copy->SetImageData(imageData, fh->GetImageDataSize()); copy->SetWriteModeToRGB(); copy->WriteDcmExplVR( output ); f->Delete(); fh->Delete(); copy->Delete(); exit (0); }