1 /*=========================================================================
4 Module: $RCSfile: AnonymizeDicomDir.cxx,v $
6 Date: $Date: 2005/07/07 17:31:53 $
7 Version: $Revision: 1.4 $
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 =========================================================================*/
18 #include "gdcmGlobal.h"
19 #include "gdcmCommon.h"
20 #include "gdcmDebug.h"
23 #include "gdcmSQItem.h"
24 #include "gdcmSeqEntry.h"
25 #include "gdcmValEntry.h"
27 #include "gdcmDocument.h"
30 #include "gdcmArgMgr.h"
35 * \brief AnonymizeDicomDir
38 void AnoNoLoad(gdcm::SQItem *s, std::fstream *fp,
39 uint16_t group, uint16_t elem,
42 void AnoNoLoad(gdcm::SQItem *s, std::fstream *fp,
43 uint16_t group, uint16_t elem,
53 d = s->GetDocEntry( group, elem);
58 if ( ! dynamic_cast<gdcm::ValEntry *>(d) )
61 offset = d->GetOffset();
62 lgth = d->GetLength();
65 spaces = new std::string( lgth-valLgth, ' ');
69 fp->seekp( offset, std::ios::beg );
70 fp->write( v.c_str(), lgth );
73 int main(int argc, char *argv[])
77 " \n AnonymizeDicomDir :\n",
78 " Anonymize a gdcm-readable DICOMDIR ",
79 " even when some 'Objects' are not yet taken into account",
80 " Warning : the DICOMDIR is overwritten",
81 " usage: AnonymizeDicomDir filein=dicomDirName [debug] ",
82 " debug : user wants to run the program in 'debug mode' ",
85 // ----- Initialize Arguments Manager ------
86 gdcm::ArgMgr *am = new gdcm::ArgMgr(argc, argv);
88 if (am->ArgMgrDefined("usage"))
90 am->ArgMgrUsage(usage); // Display 'usage'
95 char *fileName = am->ArgMgrWantString("filein",usage);
97 delete am; // we don't need Argument Manager any longer
99 // ============================================================
100 // Read the input DICOMDIR
101 // ============================================================
104 f = new gdcm::File( );
106 f->SetFileName( fileName );
107 bool res = f->Load();
110 std::cerr << "Sorry, " << fileName <<" not a gdcm-readable "
111 << "file" <<std::endl;
113 std::cout << " ... is readable " << std::endl;
115 // Directory record sequence
116 gdcm::DocEntry *e = f->GetDocEntry(0x0004, 0x1220);
119 std::cout << "No Directory Record Sequence (0004,1220) found" <<std::endl;;
125 gdcm::SeqEntry *s = dynamic_cast<gdcm::SeqEntry *>(e);
128 std::cout << "Element (0004,1220) is not a Sequence ?!?" <<std::endl;
134 // Open the file LTTG (aka ALAP)
135 std::fstream *fp = new std::fstream(fileName,
136 std::ios::in | std::ios::out | std::ios::binary);
140 int patientNumber = 0;
141 std::ostringstream oss;
143 gdcm::SQItem *tmpSI=s->GetFirstSQItem(); // For all the SQItems
146 d = tmpSI->GetDocEntry(0x0004, 0x1430); // Directory Record Type
147 if ( gdcm::ValEntry* valEntry = dynamic_cast<gdcm::ValEntry *>(d) )
149 v = valEntry->GetValue();
153 std::cout << "(0004,1430) not a ValEntry ?!?" << std::endl;
157 if( v != "PATIENT " )
159 continue; // Work only on PATIENT
162 oss << patientNumber;
164 // Overwrite the sensitive Entries
167 AnoNoLoad(tmpSI, fp, 0x0010, 0x0010, oss.str());
169 AnoNoLoad(tmpSI, fp, 0x0010, 0x0020, oss.str());
170 // Patient's Birth Date
171 AnoNoLoad(tmpSI, fp, 0x0010, 0x0030, oss.str());
173 AnoNoLoad(tmpSI, fp, 0x0010, 0x2154, oss.str());
175 // Aware use will add more Entries he wants to rubb out here
179 tmpSI=s->GetNextSQItem();
182 // Close the file ASAP