]> Creatis software - gdcm.git/blob - Testing/TestAnonymize.cxx
ENH: Removing the forbidden tags (now in NIH.dic)
[gdcm.git] / Testing / TestAnonymize.cxx
1 /*=========================================================================
2                                                                                 
3   Program:   gdcm
4   Module:    $RCSfile: TestAnonymize.cxx,v $
5   Language:  C++
6   Date:      $Date: 2005/02/09 16:41:51 $
7   Version:   $Revision: 1.1 $
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 "gdcmValEntry.h"
21 #include "gdcmBinEntry.h"
22 #include "gdcmGlobal.h"
23 #include "gdcmDebug.h"
24
25 //Generated file:
26 #include "gdcmDataImages.h"
27
28 #ifndef _WIN32
29 #include <unistd.h> //for access, unlink
30 #else
31 #include <io.h> //for _access on Win32
32 #endif
33
34
35 bool FileExists(const char *filename);
36 bool RemoveFile(const char *source);
37
38 // ---------------------------------------------------------
39
40 int Anonymize(std::string const &filename, 
41               std::string const &output )
42 {
43       std::cout << "   Testing: " << filename << std::endl;
44
45       if( FileExists( output.c_str() ) )
46       {
47          if( !RemoveFile( output.c_str() ) )
48          {
49             std::cout << "Ouch, the file exist, but I cannot remove it" 
50                       << std::endl;
51             return 1;
52          }
53       }
54
55       //////////////// Step 1:
56       std::cout << "      1...";
57       std::cout << std::endl;
58       gdcm::File *f1 = new gdcm::File( filename );
59       std::cout << " ... Read !" << std::endl;   
60       // ============================================================
61       //   Load the pixels in memory.
62       //   Write a new file
63       // ============================================================
64
65       // We need a gdcm::FileHelper, since we want to load the pixels        
66       gdcm::FileHelper *fh1 = new gdcm::FileHelper(f1);
67
68       // --- Don't forget to load the Pixels ...
69       // We shall not use them, but we have to load them
70   
71       fh1->GetImageData();
72
73       std::cout << " Image Data... Got ! " << std::endl;
74
75       // Institution name 
76       f1->AddAnonymizeElement( 0x0008, 0x0080, "Xanadoo" ); 
77       // Patient's name 
78       f1->AddAnonymizeElement( 0x0010, 0x0010, "Fantomas" );   
79       // Patient's ID
80       f1->AddAnonymizeElement( 0x0010, 0x0020,"1515" );   
81       // Study Instance UID
82       f1->AddAnonymizeElement( 0x0020, 0x000d, "9.99.999.9999" );
83       // Telephone
84       f1->AddAnonymizeElement( 0x0010, 0x2154, "3615" );
85
86       std::cout << " Anonymize list... Done ! " << std::endl;
87
88       f1->AnonymizeFile();
89
90       std::cout << " Anonymize File... Done ! " << std::endl;
91  
92       fh1->WriteDcmExplVR(output);  
93       std::cout << " Anonymized File... Written ! " << std::endl;
94
95       f1->ClearAnonymizeList();
96       std::cout << " Anonymize list... Cleared ! " << std::endl;
97     
98       delete f1;
99       delete fh1;
100
101       // Read the file we just wrote
102       f1 = new gdcm::File( output );
103
104       std::cout << " Anonymized File... Re-Read ! " << std::endl;
105
106       std::string v;
107       bool plouf = false;
108
109       // Compare and abort if different.
110       v = f1->GetEntryValue(0x0008, 0x0080);
111       if ( v != gdcm::GDCM_UNFOUND ) 
112          if (v.find("Xanadoo") >= v.length() )
113             plouf = true;
114
115       v = f1->GetEntryValue(0x0010, 0x0010);
116       if ( v != gdcm::GDCM_UNFOUND ) 
117          if (v.find("Fantomas") >= v.length() )
118             plouf = true;
119
120       if ( v != gdcm::GDCM_UNFOUND )
121          v = f1->GetEntryValue(0x0010, 0x0020);
122          if (v.find("1515") >= v.length() )
123             plouf = true;
124
125        if ( v != gdcm::GDCM_UNFOUND )
126          v = f1->GetEntryValue(0x0010, 0x000d);
127          if (v.find("9.99.999.9999") >= v.length() )
128             plouf = true;
129       
130       delete f1;
131
132       if ( !plouf)
133       { 
134          return 1;
135       } 
136
137  
138    // ============================================================
139    //   Don't load the pixels in memory.
140    //   Overwrite the file
141    // ============================================================
142
143       // Read the file we just anonymize and check
144
145       f1 = new gdcm::File( output );
146
147       // First, we set values to replace the ones we want to hide
148    
149       // Patient's name 
150       f1->AddAnonymizeElement(0x0010, 0x0010, "XXL");  
151       std::cout << " replace Patient's Name " << std::endl;
152       // Patient's ID
153       f1->AddAnonymizeElement( 0x0010, 0x0020,"007" );
154       std::cout << " replace Patient's ID " << std::endl;
155       // Study Instance UID
156       f1->AddAnonymizeElement(0x0020, 0x000d, "6.66.666.6666" );
157       std::cout << " replace Study ID " << std::endl;
158
159       // --------------------- we overwrite the file
160
161       // No need to load the pixels.
162       // The gdcm::File remains untouched in memory
163
164       std::cout <<"Let's AnonymizeNoLoad " << std::endl;;
165       f1->AnonymizeNoLoad();
166
167       std::cout <<"End AnonymizeNoLoad" << std::endl;
168       // No need to write the File : modif were done on disc !
169  
170       delete f1;
171       f1 = new gdcm::File( output );
172
173       std::string val;
174       plouf = false;
175
176       val = f1->GetEntryValue(0x0010, 0x0010);
177       if ( val != gdcm::GDCM_UNFOUND ) 
178          if (val.find("XXL") >= v.length() )
179             plouf = true;
180      
181       if ( val != gdcm::GDCM_UNFOUND )
182          val = f1->GetEntryValue(0x0010, 0x0020);
183          if (val.find("007") >= v.length() )
184       plouf = true;
185
186       if ( val != gdcm::GDCM_UNFOUND )
187          val = f1->GetEntryValue(0x0010, 0x000d);
188          if (val.find("6.66.666.6666") >= v.length() )
189       plouf = true;
190       
191       delete f1;
192
193       if ( !plouf)
194       { 
195          return 1;
196       }                       
197
198       return 0; // no fail
199 }
200
201 // Here we load a gdcmFile we anonymize it with and without loading the pixels,
202
203 int TestAnonymize(int argc, char *argv[])
204 {
205    //gdcm::Debug::DebugOn();
206
207    if ( argc == 3 )
208    {
209       // The test is specified a specific filename, use it instead of looping
210       // over all images
211       const std::string input = argv[1];
212       const std::string reference = argv[2];
213       return Anonymize( input, reference );
214    }
215    else if ( argc > 3 || argc == 2 )
216    {
217       std::cout << "   Usage: " << argv[0]
218                 << " (no arguments needed)." << std::endl;
219       std::cout << "or   Usage: " << argv[0]
220                 << " filename.dcm reference.dcm" << std::endl;
221       return 1;
222    }
223    // else other cases:
224
225    std::cout << "   Description (Test::TestAnonymize): "
226              << std::endl;
227    std::cout << "   For all images in gdcmData (and not blacklisted in "
228                 "Test/CMakeLists.txt)"
229              << std::endl;
230    std::cout << "   apply the following to each filename.xxx: "
231              << std::endl;
232    std::cout << "   - Create a new file " 
233              << std::endl;
234    std::cout << "   - Anonymize the file (creates a new file) " 
235              << std::endl;
236    std::cout << "   - Read this new file " 
237              << std::endl;
238    std::cout << "   - Check if the written values are ok" 
239              << std::endl;
240    std::cout << "   - AnonymizeNoLoad the new file (it s overwritten)" 
241              << std::endl;
242    std::cout << "   - Check if the written values are ok" 
243              << std::endl;
244    std::cout << std::endl;
245
246    int i =0;
247    int retVal = 0;  //by default this is an error
248    while( gdcmDataImages[i] != 0 )
249    {
250       std::string filename = GDCM_DATA_ROOT;
251       filename += "/";  //doh!
252       filename += gdcmDataImages[i];
253
254       std::string output = "output.dcm";
255
256       if( Anonymize( filename, output ) != 0 )
257       {
258          retVal++;
259       }
260       i++;
261    }
262    return retVal;
263 }