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:
- 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.
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
# Multiply the entire dose distribution before or after conversion:
factor = 8
a.multiply(factor)
Python 3 Standard Library with the additional packages Matplotlib, NumPy, pydicom.