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
|
from __future__ import print_function
import argparse
import odil
def add_subparser(subparsers):
parser = subparsers.add_parser(
"store", help="DICOM store (C-STORE)",
formatter_class=argparse.ArgumentDefaultsHelpFormatter)
parser.add_argument("host", help="Remote host address")
parser.add_argument("port", type=int, help="Remote host port")
parser.add_argument(
"calling_ae_title", help="AE title of the calling application")
parser.add_argument(
"called_ae_title", help="AE title of the called application")
parser.add_argument("filenames", nargs="+", help="File names")
parser.set_defaults(function=store)
return parser
def store(host, port, calling_ae_title, called_ae_title, filenames):
transfer_syntaxes = [
odil.registry.ExplicitVRLittleEndian,
odil.registry.ImplicitVRLittleEndian,
]
# Find all SOP classes to negotiate at association time. We don't need to
# read the whole data set for this
sop_classes = set()
for filename in filenames:
with odil.open(filename) as stream:
_, data_set = odil.Reader.read_file(
stream,
halt_condition=lambda tag: tag>odil.registry.SOPClassUID)
sop_classes.update(data_set.as_string("SOPClassUID"))
presentation_contexts = [
odil.AssociationParameters.PresentationContext(
2*i+1, sop_class, transfer_syntaxes,
odil.AssociationParameters.PresentationContext.Role.SCU)
for i, sop_class in enumerate(sop_classes)
]
# Create the association and the Store SCU
association = odil.Association()
association.set_peer_host(host)
association.set_peer_port(port)
association.update_parameters()\
.set_calling_ae_title(calling_ae_title)\
.set_called_ae_title(called_ae_title)\
.set_presentation_contexts(presentation_contexts)
association.associate()
store = odil.StoreSCU(association)
for filename in filenames:
with odil.open(filename) as stream:
_, data_set = odil.Reader.read_file(stream)
try:
store.set_affected_sop_class(data_set)
store.store(data_set)
except Exception as e:
print("Could not store {}: {}".format(filename, e))
association.release()
|
