]> Creatis software - gdcm.git/blob - Example/TestCopyDicom.cxx
BUG: Comitting patch from JPR, to support IRAD Data Set
[gdcm.git] / Example / TestCopyDicom.cxx
1 /*=========================================================================
2                                                                                 
3   Program:   gdcm
4   Module:    $RCSfile: TestCopyDicom.cxx,v $
5   Language:  C++
6   Date:      $Date: 2004/11/16 04:26:18 $
7   Version:   $Revision: 1.12 $
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 "gdcmHeader.h"
19 #include "gdcmFile.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 *original = new gdcm::File( filename );
84    
85       std::cout << "--- Original ----------------------" << std::endl;
86       //original->GetHeader()->Print();
87    
88       gdcm::File *copy = new gdcm::File( output );
89
90       const gdcm::TagDocEntryHT & Ht = original->GetHeader()->GetTagHT();
91
92       size_t dataSize = original->GetImageDataSize();
93       uint8_t* imageData = original->GetImageData();
94   
95       //First of all copy the header field by field
96   
97       // Warning :Accessor gdcmElementSet::GetEntry() should not exist 
98       // It was commented out by Mathieu, that was a *good* idea
99       // (the user does NOT have to know the way we implemented the Header !)
100       // Waiting for a 'clean' solution, I keep the method ...JPRx
101
102       gdcm::DocEntry* d;
103
104       for (gdcm::TagDocEntryHT::const_iterator tag = Ht.begin(); tag != Ht.end(); ++tag)
105       {
106          d = tag->second;
107          d->Print(); std::cout << std::endl;
108          if ( gdcm::BinEntry* b = dynamic_cast<gdcm::BinEntry*>(d) )
109          {              
110             copy->GetHeader()->ReplaceOrCreateByNumber( 
111                                  b->GetBinArea(),
112                                  b->GetLength(),
113                                  b->GetGroup(), 
114                                  b->GetElement(),
115                                  b->GetVR() );
116          }
117          else if ( gdcm::ValEntry* v = dynamic_cast<gdcm::ValEntry*>(d) )
118          {   
119             copy->GetHeader()->ReplaceOrCreateByNumber( 
120                                  v->GetValue(),
121                                  v->GetGroup(), 
122                                  v->GetElement(),
123                                  v->GetVR() ); 
124          }
125          else
126          {
127           // We skip pb of SQ recursive exploration
128           //std::cout << "Skipped Sequence " 
129           //          << "------------- " << d->GetVR() << " "<< std::hex
130           //          << d->GetGroup() << " " << d->GetElement()
131           //  << std::endl;    
132          }
133       }
134
135
136
137
138
139       //copy->GetImageData();
140       copy->SetImageData(imageData, dataSize);
141
142       std::cout << "--- Copy ----------------------" << std::endl;
143       std::cout <<std::endl << "DO NOT care about Offset"  <<std::endl<<std::endl;; 
144       copy->GetHeader()->Print();
145       std::cout << "--- ---- ----------------------" << std::endl;
146    
147       copy->WriteDcmExplVR( output );
148
149       return 0;
150 }
151
152
153