Skip to content

Convert DICOM-RT dose matrices to EQD2 via Python

Notifications You must be signed in to change notification settings

chris201706/PyEQD2

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

PyEQD2

Use Python to convert a DICOM-RT dose distribution to its Linear Quadratic Equivalent Dose in 2-Gy Fractions (LQED2, EQD2), voxel by voxel, according to the equation:

Important Facts

  • The input dose distribution may be the total dose or a fraction dose. Use the 'multiply' method if necessary (see examples).
  • Tested only for use with RayStation 6.

Example (Main Functions)

import pyeqd2 as eqd

dicompath = 'my_directory/testfile.dcm'

a = eqd.RTDose(dicompath)   # load dicom-rt file
a.plot()                    # get a glimpse
a.plot_legend()             # interpret the glimpse

fractions = 8                       # number of fractions
abratio = 3                         # alpha/beta ratio
a.make_eqd2(fractions, abratio)     # conversion
a.export()                          # save as new dicom-rt file

Example (Other Commands)

# Multiply the entire dose distribution before or after conversion:
factor = 8
a.multiply(factor)

Requirements

Python 3 Standard Library with the additional packages Matplotlib, NumPy, pydicom.