1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
############################################################################
#
# Program: GDCM (Grassroots DICOM). A DICOM library
#
# Copyright (c) 2006-2011 Mathieu Malaterre
# All rights reserved.
# See Copyright.txt or http://gdcm.sourceforge.net/Copyright.html for details.
#
# This software is distributed WITHOUT ANY WARRANTY; without even
# the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
# PURPOSE. See the above copyright notice for more information.
#
############################################################################
"""
Usage:
python ManipulateFile.py input.dcm output.dcm
Footnote:
GDCM 1.2.x would create incorrect Multiframes MR Image Storage file. Try to recover from
the issues to recreate a MultiframeGrayscaleByteSecondaryCaptureImageStorage file.
e.g:
python ManipulateFile.py Insight/Testing/Temporary/itkGDCMImageIOTest5-j2k.dcm manipulated.dcm
"""
import sys
import gdcm
if __name__ == "__main__":
file1 = sys.argv[1]
file2 = sys.argv[2]
r = gdcm.Reader()
r.SetFileName( file1 )
if not r.Read():
sys.exit(1)
ano = gdcm.Anonymizer()
ano.SetFile( r.GetFile() )
ano.RemovePrivateTags()
ano.Remove( gdcm.Tag(0x0032,0x1030) )
ano.Remove( gdcm.Tag(0x008,0x14) )
ano.Remove( gdcm.Tag(0x008,0x1111) )
ano.Remove( gdcm.Tag(0x008,0x1120) )
ano.Remove( gdcm.Tag(0x008,0x1140) )
ano.Remove( gdcm.Tag(0x10,0x21b0) )
ano.Empty( gdcm.Tag(0x10,0x10) )
ano.Empty( gdcm.Tag(0x10,0x20) )
ano.Empty( gdcm.Tag(0x10,0x30) )
ano.Empty( gdcm.Tag(0x20,0x10) )
ano.Empty( gdcm.Tag(0x32,0x1032) )
ano.Empty( gdcm.Tag(0x32,0x1033) )
ano.Empty( gdcm.Tag(0x40,0x241) )
ano.Empty( gdcm.Tag(0x40,0x254) )
ano.Empty( gdcm.Tag(0x40,0x253) )
ano.Empty( gdcm.Tag(0x40,0x1001) )
ano.Empty( gdcm.Tag(0x8,0x80) )
ano.Empty( gdcm.Tag(0x8,0x50) )
ano.Empty( gdcm.Tag(0x8,0x1030) )
ano.Empty( gdcm.Tag(0x8,0x103e) )
ano.Empty( gdcm.Tag(0x18,0x1030) )
ano.Empty( gdcm.Tag(0x38,0x300) )
g = gdcm.UIDGenerator()
ano.Replace( gdcm.Tag(0x0008,0x0018), g.Generate() )
ano.Replace( gdcm.Tag(0x0020,0x00d), g.Generate() )
ano.Replace( gdcm.Tag(0x0020,0x00e), g.Generate() )
ano.Replace( gdcm.Tag(0x0020,0x052), g.Generate() )
#ano.Replace( gdcm.Tag(0x0008,0x0016), "1.2.840.10008.5.1.4.1.1.7.2" )
"""
ano.Remove( gdcm.Tag(0x0018,0x0020) ) # ScanningSequence
ano.Remove( gdcm.Tag(0x0018,0x0021) ) # SequenceVariant
ano.Remove( gdcm.Tag(0x0018,0x0022) ) # ScanOptions
ano.Remove( gdcm.Tag(0x0018,0x0023) ) # MRAcquisitionType
ano.Remove( gdcm.Tag(0x0018,0x0050) ) # SliceThickness
ano.Remove( gdcm.Tag(0x0018,0x0080) ) # RepetitionTime
ano.Remove( gdcm.Tag(0x0018,0x0081) ) # EchoTime
ano.Remove( gdcm.Tag(0x0018,0x0088) ) # SpacingBetweenSlices
ano.Remove( gdcm.Tag(0x0018,0x0091) ) # EchoTrainLength
ano.Remove( gdcm.Tag(0x0018,0x1164) ) # ImagerPixelSpacing
ano.Remove( gdcm.Tag(0x0020,0x0032) ) # Image Position (Patient)
ano.Remove( gdcm.Tag(0x0020,0x0037) ) # Image Orientation (Patient)
ano.Remove( gdcm.Tag(0x0020,0x0052) ) # Frame of Reference UID
ano.Remove( gdcm.Tag(0x0020,0x1040) ) # Position Reference Indicator
ano.Replace( gdcm.Tag(0x0028,0x0301), "NO" ) # Burned In Annotation
ano.Empty( gdcm.Tag(0x0020,0x0020) )
ano.Remove( gdcm.Tag(0x7fe0,0x0000) )
#ano.Empty( gdcm.Tag(0x0028,0x0009) ) # Frame Increment Pointer
#ano.Empty( gdcm.Tag(0x0028,0x1052) ) #<entry group="0028" element="1052" vr="DS" vm="1" name="Rescale Intercept"/>
#ano.Empty( gdcm.Tag(0x0028,0x1053) ) #<entry group="0028" element="1053" vr="DS" vm="1" name="Rescale Slope"/>
#ano.Replace( gdcm.Tag(0x0028,0x1054), "US" ) #<entry group="0028" element="1054" vr="LO" vm="1" name="Rescale Type"/>
ano.Replace( gdcm.Tag(0x2050, 0x0020), "IDENTITY")
"""
w = gdcm.Writer()
w.SetFile( ano.GetFile() )
w.SetFileName( file2 )
if not w.Write():
sys.exit(1)
|
