[Dcmlib] [Fwd: [Insight-developers] anonymize DICOM files]

[Jean-Pierre Roux]

Mathieu Malaterre wrote:

> Je fais suivre un example de la ML ITK.
>
> Il lui manquait la ligne super importante dans gdcm 0.6:
>
> gdcm::File f(&h);
> f.GetImageData(); //super important
> f.WriteDcmExplVR(OutputFilename);
>
> Sinon c'est realisable un vrai 'anonymiser' ? C'est a dire ne modifier 
> que l'entete sans essayer de lire/analyser l'image data ? (=a la 
> hexedit).


Il y a la methode gdcm::File::AnonymizeFile()

/**
* \brief anonymize a File (removes Patient's personal info)
* (read the code to see which ones ...)
*/
bool File::AnonymizeFile()
{
// If exist, replace by spaces
SetValEntry (" ",0x0010, 0x2154); // Telephone
SetValEntry (" ",0x0010, 0x1040); // Adress
SetValEntry (" ",0x0010, 0x0020); // Patient ID

DocEntry* patientNameHE = GetDocEntry (0x0010, 0x0010);

if ( patientNameHE ) // we replace it by Study Instance UID (why not)
{
std::string studyInstanceUID = GetEntryValue (0x0020, 0x000d);
if ( studyInstanceUID != GDCM_UNFOUND )
{
InsertValEntry(studyInstanceUID, 0x0010, 0x0010);
}
else
{
InsertValEntry("anonymised", 0x0010, 0x0010);
}
}

// Just for fun :-(
// (if any) remove or replace all the stuff that contains a Date

//0008 0012 DA ID Instance Creation Date
//0008 0020 DA ID Study Date
//0008 0021 DA ID Series Date
//0008 0022 DA ID Acquisition Date
//0008 0023 DA ID Content Date
//0008 0024 DA ID Overlay Date
//0008 0025 DA ID Curve Date
//0008 002a DT ID Acquisition Datetime
//0018 9074 DT ACQ Frame Acquisition Datetime
//0018 9151 DT ACQ Frame Reference Datetime
//0018 a002 DT ACQ Contribution Date Time
//0020 3403 SH REL Modified Image Date (RET)
//0032 0032 DA SDY Study Verified Date
//0032 0034 DA SDY Study Read Date
//0032 1000 DA SDY Scheduled Study Start Date
//0032 1010 DA SDY Scheduled Study Stop Date
//0032 1040 DA SDY Study Arrival Date
//0032 1050 DA SDY Study Completion Date
//0038 001a DA VIS Scheduled Admission Date
//0038 001c DA VIS Scheduled Discharge Date
//0038 0020 DA VIS Admitting Date
//0038 0030 DA VIS Discharge Date
//0040 0002 DA PRC Scheduled Procedure Step Start Date
//0040 0004 DA PRC Scheduled Procedure Step End Date
//0040 0244 DA PRC Performed Procedure Step Start Date
//0040 0250 DA PRC Performed Procedure Step End Date
//0040 2004 DA PRC Issue Date of Imaging Service Request
//0040 4005 DT PRC Scheduled Procedure Step Start Date and Time
//0040 4011 DT PRC Expected Completion Date and Time
//0040 a030 DT PRC Verification Date Time
//0040 a032 DT PRC Observation Date Time
//0040 a120 DT PRC DateTime
//0040 a121 DA PRC Date
//0040 a13a DT PRC Referenced Datetime
//0070 0082 DA ??? Presentation Creation Date
//0100 0420 DT ??? SOP Autorization Date and Time
//0400 0105 DT ??? Digital Signature DateTime
//2100 0040 DA PJ Creation Date
//3006 0008 DA SSET Structure Set Date
//3008 0024 DA ??? Treatment Control Point Date
//3008 0054 DA ??? First Treatment Date
//3008 0056 DA ??? Most Recent Treatment Date
//3008 0162 DA ??? Safe Position Exit Date
//3008 0166 DA ??? Safe Position Return Date
//3008 0250 DA ??? Treatment Date
//300a 0006 DA RT RT Plan Date
//300a 022c DA RT Air Kerma Rate Reference Date
//300e 0004 DA RT Review Date

return true;
}

>
> Mathieu
>
>
> -------- Original Message --------
> Subject: [Insight-developers] anonymize DICOM files
> Date: Thu, 03 Feb 2005 17:26:17 -0600
> From: Kent Williams <norman-k-williams at uiowa.edu>
> To: insight-developers at itk.org
>
> I have several DICOM data sets that for one reason or another are not
> working properly with the gdcm Dicom reader. The problems we're having
> aren't a big deal -- and I can deal with it off the list with the GDCM
> guy -- but I need to anonymize the datasets before sending them offsite.
>
> My idea was to write a simple app that anonymizes the header and dumps a
> new file out. gdcm::Header has a function called 'AnonymizeHeader' that
> I copied to blank out the problem fields, but it's not clear to me how
> to read in the whole dicom slice and write it back out along with the
> pixel data -- using GDCM.
>
> What's happening is that it appears the header is getting written out OK
> but the pixel data is not.
>
> Below is the program I wrote -- clearly I'm missing the 'read in the
> data from the input data, attach it to the output file stage.'
>
> I would have used an anonymizer someone else wrote -- if there was one I
> could trust, that was free, and didn't have a stupid user interface!
> Google didn't find anything like that!
>
>
> #include "gdcm.h"
> #include "unistd.h"
>
> int
> main(int argc, char **argv)
> {
> char *inputfilename(0),
> *outputfilename(0);
> int c;
> bool quitopt = false;
> while(!quitopt && (c = getopt(argc,argv,"i:o:")) != -1)
> {
> switch(c)
> {
> case 'i':
> inputfilename = optarg;
> break;
> case 'o':
> outputfilename = optarg;
> break;
> default:
> quitopt = true;
> break;
> }
> }
> if(inputfilename == 0)
> {
> std::cerr << "No input file given" << std::endl;
> exit(-1);
> }
> if(outputfilename == 0)
> {
> std::cerr << "No output file given" << std::endl;
> exit(-1);
> }
> std::string InputFilename(inputfilename),
> OutputFilename(outputfilename);
>
> gdcm::Header h(InputFilename);
> if(!h.IsReadable())
> {
> std::cerr << "Can't open input file " << inputfilename << std::endl;
> exit(-1);
> }
> // If exist, replace by spaces
> h.SetEntryByNumber (" ",0x0010, 0x2154); // Telephone
> h.SetEntryByNumber (" ",0x0010, 0x1040); // Adress
> h.SetEntryByNumber (" ",0x0010, 0x0020); // Patient ID
>
> gdcm::DocEntry* patientNameHE = h.GetDocEntryByNumber (0x0010, 0x0010);
>
> if ( patientNameHE ) // we replace it with bogus name
> {
> h.ReplaceOrCreateByNumber("Ernest_T_Bass", 0x0010, 0x0010);
> }
> gdcm::File f(&h);
> f.WriteDcmExplVR(OutputFilename);
> }
>
> _______________________________________________
> Insight-developers mailing list
> Insight-developers at itk.org
> http://www.itk.org/mailman/listinfo/insight-developers
>
> _______________________________________________
> Dcmlib mailing list
> Dcmlib at creatis.insa-lyon.fr
> http://www.creatis.insa-lyon.fr/mailman/listinfo/dcmlib
>