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
|
############################################################################
#
# 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 MergeFile.py input1.dcm input2.dcm
It will produce a 'merge.dcm' output file, which contains all meta information from input1.dcm
and copy the Stored Pixel values from input2.dcm
This script even works when input2.dcm is a Secondary Capture and does not contains information
such as IOP and IPP...
"""
import sys
import gdcm
if __name__ == "__main__":
file1 = sys.argv[1]
file2 = sys.argv[2]
r1 = gdcm.ImageReader()
r1.SetFileName( file1 )
if not r1.Read():
sys.exit(1)
r2 = gdcm.ImageReader()
r2.SetFileName( file2 )
if not r2.Read():
sys.exit(1)
# Image from r2 could be Secondary Capture and thus would not contains neither IPP nor IOP
# Instead always prefer to only copy the Raw Data Element.
# Warning ! Image need to be identical ! Only the value of Stored Pixel can be different.
r1.GetImage().SetDataElement( r2.GetImage().GetDataElement() )
w = gdcm.ImageWriter()
w.SetFile( r1.GetFile() )
#w.SetImage( r2.GetImage() ) # See comment above
w.SetImage( r1.GetImage() )
w.SetFileName( "merge.dcm" )
if not w.Write():
sys.exit(1)
sys.exit(0)
|
