1 /*=========================================================================
4 Module: $RCSfile: TestCopyDicom.cxx,v $
6 Date: $Date: 2005/10/17 10:41:59 $
7 Version: $Revision: 1.30 $
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 "gdcmValEntry.h"
22 #include "gdcmBinEntry.h"
25 #include <unistd.h> //for access, unlink
27 #include <io.h> //for _access
30 // return true if the file exists
31 bool FileExists(const char *filename)
34 # define access _access
39 if ( access(filename, R_OK) != 0 )
49 bool RemoveFile(const char *source)
52 #define _unlink unlink
54 return unlink(source) != 0 ? false : true;
57 // Here we load a gdcmFile and then try to create from scratch a copy of it,
58 // copying field by field the dicom image
60 int main(int argc, char *argv[])
64 std::cerr << "Usage :" << std::endl <<
65 argv[0] << " input_dicom output_dicom" << std::endl;
69 // don't modify identation in order to let this source xdiffable with ../Test
71 std::string filename = argv[1];
72 std::string output = argv[2];
74 if( FileExists( output.c_str() ) )
76 std::cerr << "Don't try to cheat, I am removing the file anyway"
78 if( !RemoveFile( output.c_str() ) )
80 std::cerr << "Ouch, the file exist, but I cannot remove it"
85 gdcm::File *fileOr = new gdcm::File();
86 fileOr->SetFileName( filename );
88 gdcm::FileHelper *original = new gdcm::FileHelper( fileOr );
90 std::cout << "--- Original ----------------------" << std::endl;
91 //original->GetFile()->Print();
93 gdcm::FileHelper *copy = new gdcm::FileHelper( );
94 copy->SetFileName( output );
99 //dataSize = original->GetImageDataSize();// just an accesor :useless here
101 imageData = original->GetImageData(); // VERY important :
102 // brings pixels into memory !
103 //(void)imageData; // not enough to avoid warning with icc compiler
104 //(void)dataSize; // not enough to avoid warning on 'Golgot'
106 std::cout << imageData << std::endl; // to avoid warning ?
108 //First of all copy the header field by field
110 gdcm::DocEntry *d = original->GetFile()->GetFirstEntry();
113 if ( gdcm::BinEntry *b = dynamic_cast<gdcm::BinEntry*>(d) )
115 copy->GetFile()->InsertBinEntry( b->GetBinArea(),b->GetLength(),
116 b->GetGroup(),b->GetElement(),
119 else if ( gdcm::ValEntry *v = dynamic_cast<gdcm::ValEntry*>(d) )
121 copy->GetFile()->InsertValEntry( v->GetValue(),
122 v->GetGroup(),v->GetElement(),
127 // We skip pb of SQ recursive exploration
128 std::cout << "Skipped Sequence "
129 << "------------- " << d->GetVR() << " "<< std::hex
130 << d->GetGroup() << "," << d->GetElement()
134 d=original->GetFile()->GetNextEntry();
137 //copy->GetImageData();
138 //copy->SetImageData(imageData, dataSize);
140 std::cout << "--- Copy ----------------------" << std::endl;
141 std::cout <<std::endl << "DO NOT care about Offset"
142 <<std::endl << std::endl;;
143 copy->GetFile()->Print();
144 std::cout << "--- ---- ----------------------" << std::endl;
146 copy->WriteDcmExplVR( output );
149 delete fileOr; // File
150 delete original; // FileHelper
151 delete copy; // FileHelper