]> Creatis software - gdcm.git/blob - Example/TestCopyDicom.cxx
2dd66ff34f0fa87e632f87a1b6e0e62660ec169a
[gdcm.git] / Example / TestCopyDicom.cxx
1 /*=========================================================================
2                                                                                 
3   Program:   gdcm
4   Module:    $RCSfile: TestCopyDicom.cxx,v $
5   Language:  C++
6   Date:      $Date: 2005/07/19 15:19:25 $
7   Version:   $Revision: 1.26 $
8                                                                                 
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.
12                                                                                 
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.
16                                                                                 
17 =========================================================================*/
18 #include "gdcmFile.h"
19 #include "gdcmFileHelper.h"
20 #include "gdcmDocument.h"
21 #include "gdcmValEntry.h"
22 #include "gdcmBinEntry.h"
23
24 #ifndef _WIN32
25 #include <unistd.h> //for access, unlink
26 #else
27 #include <io.h> //for _access
28 #endif
29
30 // return true if the file exists
31 bool FileExists(const char *filename)
32 {
33 #ifdef _MSC_VER
34 # define access _access
35 #endif
36 #ifndef R_OK
37 # define R_OK 04
38 #endif
39   if ( access(filename, R_OK) != 0 )
40     {
41     return false;
42     }
43   else
44     {
45     return true;
46     }
47 }
48
49 bool RemoveFile(const char *source)
50 {
51 #ifdef _MSC_VER
52 #define _unlink unlink
53 #endif
54   return unlink(source) != 0 ? false : true;
55 }
56
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
59
60 int main(int argc, char *argv[])
61 {
62    if (argc < 3)
63    {
64       std::cerr << "Usage :" << std::endl << 
65       argv[0] << " input_dicom output_dicom" << std::endl;
66       return 1;
67    }
68
69 // don't modify identation in order to let this source xdiffable with ../Test
70
71       std::string filename = argv[1];
72       std::string output = argv[2];
73
74       if( FileExists( output.c_str() ) )
75       {
76          std::cerr << "Don't try to cheat, I am removing the file anyway" << std::endl;
77          if( !RemoveFile( output.c_str() ) )
78          {
79             std::cerr << "Ouch, the file exist, but I cannot remove it" << std::endl;
80             return 1;
81          }
82       }
83       gdcm::File *fileOr = new gdcm::File();
84       fileOr->SetFileName( filename );
85       fileOr->Load();
86       gdcm::FileHelper *original = new gdcm::FileHelper( fileOr );
87    
88       std::cout << "--- Original ----------------------" << std::endl;
89       //original->GetFile()->Print();
90    
91       gdcm::FileHelper *copy = new gdcm::FileHelper( );
92       copy->SetFileName( output );
93       copy->Load();
94
95       size_t dataSize = original->GetImageDataSize();
96       uint8_t *imageData = original->GetImageData();
97       (void)imageData;
98       (void)dataSize;
99   
100       //First of all copy the header field by field
101
102       gdcm::DocEntry *d = original->GetFile()->GetFirstEntry();
103       while(d)
104       {
105          if ( gdcm::BinEntry *b = dynamic_cast<gdcm::BinEntry*>(d) )
106          {              
107             copy->GetFile()->InsertBinEntry( b->GetBinArea(),b->GetLength(),
108                                              b->GetGroup(),b->GetElement(),
109                                              b->GetVR() );
110          }
111          else if ( gdcm::ValEntry *v = dynamic_cast<gdcm::ValEntry*>(d) )
112          {   
113             copy->GetFile()->InsertValEntry( v->GetValue(),
114                                              v->GetGroup(),v->GetElement(),
115                                              v->GetVR() ); 
116          }
117          else
118          {
119           // We skip pb of SQ recursive exploration
120           std::cout << "Skipped Sequence " 
121                     << "------------- " << d->GetVR() << " "<< std::hex
122                     << d->GetGroup() << "," << d->GetElement()
123                     << std::endl;    
124          }
125
126          d=original->GetFile()->GetNextEntry();
127       }
128
129       //copy->GetImageData();
130       //copy->SetImageData(imageData, dataSize);
131
132       std::cout << "--- Copy ----------------------" << std::endl;
133       std::cout <<std::endl << "DO NOT care about Offset"  <<std::endl<<std::endl;; 
134       copy->GetFile()->Print();
135       std::cout << "--- ---- ----------------------" << std::endl;
136    
137       copy->WriteDcmExplVR( output );
138
139       return 0;
140 }
141
142
143