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
|
"""
======================================
Read RTPLAN DICOM and list information
======================================
Given an RTPLAN DICOM file, list basic info for the beams in it.
"""
# authors : Guillaume Lemaitre <g.lemaitre58@gmail.com>
# license : MIT
import pydicom
from pydicom.data import get_testdata_file
print(__doc__)
def list_beams(plan_dataset):
"""Summarizes the RTPLAN beam information in the dataset."""
lines = ["{name:^13s} {num:^8s} {gantry:^8s} {ssd:^11s}".format(
name="Beam name", num="Number", gantry="Gantry", ssd="SSD (cm)")]
for beam in plan_dataset.BeamSequence:
cp0 = beam.ControlPointSequence[0]
SSD = float(cp0.SourceToSurfaceDistance / 10)
lines.append("{b.BeamName:^13s} {b.BeamNumber:8d} "
"{gantry:8.1f} {ssd:8.1f}".format(b=beam,
gantry=cp0.GantryAngle,
ssd=SSD))
return "\n".join(lines)
filename = get_testdata_file('rtplan.dcm')
dataset = pydicom.dcmread(filename)
print(list_beams(dataset))
|
