1 /*=========================================================================
4 Module: $RCSfile: TestCopyDicom.cxx,v $
6 Date: $Date: 2005/10/18 08:35:43 $
7 Version: $Revision: 1.31 $
9 Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
10 l'Image). All rights reserved. See Doc/License.txt or
11 http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
13 This software is distributed WITHOUT ANY WARRANTY; without even
14 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 PURPOSE. See the above copyright notices for more information.
17 =========================================================================*/
19 #include "gdcmFileHelper.h"
20 #include "gdcmDocument.h"
21 #include "gdcmDataEntry.h"
24 #include <unistd.h> //for access, unlink
26 #include <io.h> //for _access
29 // return true if the file exists
30 bool FileExists(const char *filename)
33 # define access _access
38 if ( access(filename, R_OK) != 0 )
48 bool RemoveFile(const char *source)
51 #define _unlink unlink
53 return unlink(source) != 0 ? false : true;
56 // Here we load a gdcmFile and then try to create from scratch a copy of it,
57 // copying field by field the dicom image
59 int main(int argc, char *argv[])
63 std::cerr << "Usage :" << std::endl <<
64 argv[0] << " input_dicom output_dicom" << std::endl;
68 // don't modify identation in order to let this source xdiffable with ../Test
70 std::string filename = argv[1];
71 std::string output = argv[2];
73 if( FileExists( output.c_str() ) )
75 std::cerr << "Don't try to cheat, I am removing the file anyway"
77 if( !RemoveFile( output.c_str() ) )
79 std::cerr << "Ouch, the file exist, but I cannot remove it"
84 gdcm::File *fileOr = new gdcm::File();
85 fileOr->SetFileName( filename );
87 gdcm::FileHelper *original = new gdcm::FileHelper( fileOr );
89 std::cout << "--- Original ----------------------" << std::endl;
90 //original->GetFile()->Print();
92 gdcm::FileHelper *copy = new gdcm::FileHelper( );
93 copy->SetFileName( output );
98 //dataSize = original->GetImageDataSize();// just an accesor :useless here
100 imageData = original->GetImageData(); // VERY important :
101 // brings pixels into memory !
102 //(void)imageData; // not enough to avoid warning with icc compiler
103 //(void)dataSize; // not enough to avoid warning on 'Golgot'
105 std::cout << imageData << std::endl; // to avoid warning ?
107 //First of all copy the header field by field
109 gdcm::DocEntry *d = original->GetFile()->GetFirstEntry();
112 if ( gdcm::DataEntry *de = dynamic_cast<gdcm::DataEntry *>(d) )
114 copy->GetFile()->InsertEntryBinArea( de->GetBinArea(),de->GetLength(),
115 de->GetGroup(),de->GetElement(),
120 // We skip pb of SQ recursive exploration
121 std::cout << "Skipped Sequence "
122 << "------------- " << d->GetVR() << " "<< std::hex
123 << d->GetGroup() << "," << d->GetElement()
127 d=original->GetFile()->GetNextEntry();
130 //copy->GetImageData();
131 //copy->SetImageData(imageData, dataSize);
133 std::cout << "--- Copy ----------------------" << std::endl;
134 std::cout <<std::endl << "DO NOT care about Offset"
135 <<std::endl << std::endl;;
136 copy->GetFile()->Print();
137 std::cout << "--- ---- ----------------------" << std::endl;
139 copy->WriteDcmExplVR( output );
142 delete fileOr; // File
143 delete original; // FileHelper
144 delete copy; // FileHelper