Skip to main content

A dose-volume histogram based decision-support system for dosimetric comparison of radiotherapy treatment plans



The choice of any radiotherapy treatment plan is usually made after the evaluation of a few preliminary isodose distributions obtained from different beam configurations. Despite considerable advances in planning techniques, such final decision remains a challenging task that would greatly benefit from efficient and reliable assessment tools.


For any dosimetric plan considered, data on dose-volume histograms supplied by treatment planning systems are used to provide estimates on planning target coverage as well as on sparing of organs at risk and the remaining healthy tissue. These partial metrics are then combined into a dose distribution index (DDI), which provides a unified, easy-to-read score for each competing radiotherapy plan. To assess the performance of the proposed scoring system, DDI figures for fifty brain cancer patients were retrospectively evaluated. Patients were divided in three groups depending on tumor location and malignancy. For each patient, three tentative plans were designed and recorded during planning, one of which was eventually selected for treatment. We thus were able to compare the plans with better DDI scores and those actually delivered.


When planning target coverage and organs at risk sparing are considered as equally important, the tentative plan with the highest DDI score is shown to coincide with that actually delivered in 32 of the 50 patients considered. In 15 (respectively 3) of the remaining 18 cases, the plan with highest DDI value still coincides with that actually selected, provided that organs at risk sparing is given higher priority (respectively, lower priority) than target coverage.


DDI provides a straightforward and non-subjective tool for dosimetric comparison of tentative radiotherapy plans. In particular, DDI readily quantifies differences among competing plans with similar-looking dose-volume histograms and can be easily implemented for any tumor type and localization, irrespective of the planning system and irradiation technique considered. Moreover, DDI permits to estimate the dosimetry impact of different priorities being assigned to sparing of organs at risk or to better target coverage.


A radiotherapy treatment plan is expected to uniformly irradiate a selected planning target volume (PTV), while at the same time minimizing radiation-induced damage to organs at risk (OARs) and the remaining volume at risk (RVR) [13]. To that end, a prescription dose (D p ) to be delivered at the PTV should be defined, whereas radiation doses at OARs and RVRs should be kept as low as possible. To address that issue, computerized treatment planning systems (TPS) are used to simulate a number of tentative plans. For any such plan, the TPS provides dose-volume histograms (DVHs) for the PTV and each OAR delineated, as well as the isodose curves over the whole treatment domain. The choice of the treatment plan is then usually made, regardless of further physical and radiobiological considerations, mainly by inspection of the DVHs and isodose curves obtained for each tentative plan considered.

It is widely accepted that the decision-making process just recalled is highly subjective. Actually, considerable skills are required to pick up the most suitable tentative plan from a mere inspection of the DVHs and isodose curves of competing plans, since these may look quite similar to the naked eye. During the last decades, a number of approaches have been proposed to improve the choice of treatment plans. However, and in spite of the large efforts made so far, no widely accepted quantitative decision-support tool seems to be currently in clinical use [3, 4]. A standard approach towards deriving such decision-support tools has consisted in proposing physical-based quantitative metrics. Such metrics assign a score to a dosimetric plan according to its compliance with specified physical criteria [3]. Paramount among them, the conventionally used homogeneity index (HI) is defined as the ratio between the maximum radiation dose in the PTV and the prescription dose [5, 6]. In addition, alternative indices of dose homogeneity on the target have been proposed [710], which however fail often to accurately account for dose homogeneity within the PTV [7, 11]. On the other hand, a conformity index (CI) was first proposed in the Radiation Therapy Oncology Group (RTOG) radiosurgery guidelines [5], where it was defined as the ratio of the prescription isodose volume to the target volume. However, CI can yield false positives in extreme cases of nonconcordance between the target volume and isodose curves [6, 12]. Other conformity indices have been proposed [1318], but practical limitations to their use have been pointed out in ICRU reports [13]. A major drawback of most of such indices is that they do not consider the overall dosimetric information provided by DVHs. As a matter of fact, a consensus has been built to consider DVHs as key indicators of the compliance with clinical requirements, and for that reason DVHs are routinely provided by the TPS.

The previous figures of merit provide a measure of target coverage and conformity only, without assessing organs at risk and healthy tissue sparing. This is a serious drawback, since normal tissue morbidity often represents the main limiting factor to radiation delivery and treatment outcomes [19, 20]. For that reason, several volume-based indices have been proposed to estimate such effects. These include, the healthy tissue overdosage factor (HTOF) [6], healthy tissue conformity index (HTCI) [6, 15] and critical organ-scoring index (COSI) [21]. Moreover, dose conformality at the PTV and OARs has also been simultaneously quantified [6, 2224]. However, questions have been raised about the suitability of such indices to account for cold and hot spots (regions receiving particularly low or high radiation doses) in the PTV and OARs respectively [6, 7, 24].

It follows from our previous remarks that a single assessment index, accounting for a particular clinical aspect during planning, may not be suitable to compare tentative plans since it might ignore other key dosimetric issues. For this reason, indices have been proposed that provide combined scores for several treatment requirements. For instance, an overall quality factor (QF) and a unified dosimetry index (UDI) have been defined as weighted combinations of some of the indices mentioned before [25, 26]. However, the difficulty in interpreting the overall score thus obtained sets serious limitations to their practical use. Moreover, the figures of merit proposed often depend on many parameters [21, 24, 27], provide ambiguous scoring due to averaging effects [23, 25] and are not easy to implement in clinical routine [26, 28].

The indices recalled above are purely physical quantities based on dose-volume constraints. Radiobiological indices, such as tumor control probability (TCP) and normal tissue complication probability (NTCP), aim instead at estimating the effects of radiation on pathological and healthy tissues resulting from a given dose distribution. A large range of TCP [2931] and NTCP [3236] models have been reported, but recommendations in the ICRU reports suggest that they should be used with caution in clinical practice [13]. In particular, a major limitation to their use stems from the difficulty to estimate tissue-dependent parameter values to characterize the radiobiological response of tumors and organs involved. In many situations of clinical interest, such parameters have only been estimated in vitro or remain unknown.


Dose distribution index (DDI)

A dosimetric comparison system, the dose distribution index (DDI), is proposed to assist in the choice of radiotherapy treatment plans for any given cancer patient. DDI simultaneously takes into account key dosimetric issues as dose coverage, conformity and homogeneity over the PTV, as well as sparing of OARs and the RVR. Only standard information provided by the TPS, namely the DVHs and the prescription dose on the PTV, are required to grade and evaluate tentative plans. Moreover, one can readily estimate the separate impact on the PTV coverage and sparing of OARs and RVRs of any such modification, and thus use this information as a guidance to improve a tentative treatment plan.

The proposed decision-support system is defined as a weighted sum of three components, which quantify the overall quality of a dosimetric plan. The first of those components assesses the dose coverage, conformity and homogeneity on the planning target, and it is defined as follows:

$$ I_{T} = \left(1 - \left| 1- \left(\frac{\int\limits_{0}^{D_{M}} V_{T}(D) dD}{D_{p} \cdot PTV} \right) \cdot \left(\frac{D_{m}}{D_{M}} \right) \right| \right), $$

where PTV is the planning target volume, D p is the prescription dose on the PTV, D m is the maximum dose received at least by 100 % of the PTV, D M is the maximum dose in the PTV and V T (D) represents the DVH curve corresponding to the target. More precisely, this index computes the ratio between the areas under the DVH curve for any tentative plan considered and that corresponding to an ideal plan of constant D p uniformly delivered on the PTV, see Fig. 1(a). The ratio (D m /D M ) is included in (1) to account for dose homogeneity in the PTV. Alternatively, and depending on clinical requirements, this ratio might be taken as (D 95/D 5), where D 5 and D 95 are the dose received by the 5 and 95 % of the PTV respectively [7, 11]. Notice that as far as the PTV is considered, an ideal plan ensuring perfect dose coverage would yield a score I T =1.

Fig. 1
figure 1

Schematic representation of dose-volume histograms. a A DVH corresponding to a PTV and b an OAR or the RVR. The factors involved in the formulation of I T , I O and I R are also indicated, see (1), (2) and (3) respectively

We argue in a similar way to estimate the merit of a tentative treatment plan with respect to sparing of OARs and the RVR. In this case, a good plan would make close to zero the area under the DVH corresponding to such critical structures, see Fig. 1(b). The contributions to the DDI of OARs and RVR components are respectively formulated as follows:

$$ I_{O} = \frac{1}{N} \sum\limits_{i=1}^{N} w_{Oi} \left(1- \left(\frac{\int\limits_{0}^{D_{Mi}} V_{Oi}(D) dD}{D_{p} \cdot OARV_{i}} \right) \cdot \left(\frac{D_{M}}{D_{p}} \right) \right) $$


$$ I_{R} = \left(1- \left(\frac{\int\limits_{0}^{D_{M}} V_{R}(D) dD}{D_{p} \cdot RVR} \right) \cdot \left(\frac{D_{M}}{D_{p}} \right) \right), $$

where D M is the maximum dose in each structure involved, N is the number of OARs considered, OARV is the volume of a given OAR, V O (D) and V R (D) represent the DVH curves of a given OAR and the RVR respectively. The ratio (D M /D p ) is included to weight the deviation of the maximum dose received by OARs and the RVR from the prescription dose D p on the PTV, the latter remaining constant for each competing plan considered for a same patient. The nonnegative weighting factors w Oi ≤1.0 in (2) account for the relative clinical importance assigned to each OAR considered. In what concerns OARs and the RVR, an ideal plan would yield I O =I R =1, which would result in no irradiation at all on such critical volumes.

The DDI decision-support system is given by:

$$ DDI = \frac{1}{3} \left(w_{T} I_{T} + w_{O} I_{O} + w_{R} I_{R} \right), $$

where w T , w O and w R are nonnegative weighting parameters, each ranging between zero and one, that make precise the relative clinical priority assigned to the indices (1), (2) and (3) during the treatment planning process. These parameters allow to account for the relative importance assigned to target coverage and organs at risk sparing, from highly important (weight factor equal to 1) to less important (weight factor close to 0). In particular, if weights are selected equal 1 then all the structures involved are considered as equally important. In that manner, the DDI encodes the contributions from all three components into a unified score as defined by equation (4). Notice that an ideal plan would give D D I=1, which corresponds to uniform delivery of the D p to the PTV and no irradiation at all in OARs and the RVR. The previous formulation can be easily extended to include several PTVs, similarly to what has been done for OARs in (2). In fact, this may be the case when a heterogeneous planning target is considered, and non-homogeneous radiation doses are prescribed [10, 28, 37]. For convenience of the reader, we summarize in Table 1 the notation used in the formulation of this decision-support system.

Table 1 Description of the notation used in the formulation of the DDI decision-support system

To illustrate how the DDI would work in clinical practice, we have performed a retrospective study over fifty brain cancer patients already treated at Hospital Universitario Puerta de Hierro in Madrid (HUPH). For each patient, three tentative plans were designed with the same TPS during the planning process. Out of the three competing plans considered, one was the selected by clinicians and delivered to the patient. We have then compared in each case the applied treatment plan with the tentative plan achieving the highest DDI score.

Patients’ characteristics

Fifty brain cancer patients randomly selected and previously treated by external beam radiotherapy at HUPH between the years 2012 and 2014 were considered for the present study. In particular, patients with brain tumors were selected due the large amount of OARs involved, and the consequent difficulty of global evaluation of PTV coverage and OARs sparing. The patients were divided in three different groups, based on tumor location and malignancy, as follows: (i) 25 patients with centrally located benign tumors (13 Meningiomas, 6 Pituitary Adenomas, 5 Acoustic Neuromas and 1 Facial Nerve Neuroma) representing 50 % of the sample; (ii) 10 patients with centrally located malignant tumors (4 Malignant Craniopharyngiomas, 2 Cavum Carcinomas, 3 Chondrosarcomas and 1 Clear Cell Carcinoma (CCC) of the Base of the Tongue) making up for 20 % of the total, and (iii) 15 patients with peripheral benign and malignant tumors (12 Brain Metastasis, 2 Astrocytomas and 1 Meningioma) which correspond to 30 % of the cases considered.

Volume definitions

PTVs, OARs and RVRs were defined as recommended by ICRU reports [13]. The gross tumour volume (GTV) was delineated on all relevant slices of a planning CT scan. The margins were then expanded isotropically to allow for microscopic invasion, thus giving raise to the clinical target volume (CTV). This last is further enlarged to allow for organ motion and set-up error. The resulting structure of interest is then defined as the planning target volume (PTV). The organs selected as OARs including the Brainstem, Eyes, Crystalline Lens, Lacrimal Glands, Retinas, Optic Chiasm, Hypothalamus, Cochleas, Spinal Cord, Parotid Glands, Pituitary Gland, Optic Nerves, Optic Tracts and Cranial Nerves (Trigeminal, Facial and Vestibulocochlear Nerves) were contoured on the relevant transverse CT slices of the planning system.

Treatment planning

For each patient, three different tentative plans using three-dimensional conformal radiation therapy (3D-CRT) or sliding windows intensity-modulated radiation therapy (IMRT) were simulated. The dose-volume planning prescription on the target required that at least 95 % of the PTV had to receive no less than 95 % of the prescription dose. This last was established at 50.4 Gy in 28 daily fractions (5 days a week) of 1.8 Gy for conventional fractionation protocols (in twenty-nine out of fifty patients, 58 % of total), and between 12 Gy–18 Gy in a single session for stereotactic radiosurgery (in the remaining twenty-one patients, 42 % of total). The dose planification system was an iPLAN RT Dose 4.1.1 (BrainLAB AG, Germany). Sets of three alternative treatment plans, with the same radiotherapy technique, structures of interest and clinical requirements, were prepared for each patient. Out of the three tentative plans considered, one was selected and then delivered to the patient. For all tentative plans the DVHs corresponding to each volume delineated were extracted from the treatment planning system for computation of the DDI.


In Table 2 the resulting mean I T , I O , I R and DDI scores, as well as standard deviations for the overall set of tentative plans are provided. Weighting factors in components (2) and (4) were set equal to one, so that the PTV coverage and sparing of OARs and RVRs are considered as equally important. The lowest and highest DDI scores, corresponding to the minimum and maximum deviations from an ideal dosimetric plan, were 0.70 and 0.98 respectively. A comparison of mean I T , I O , I R and DDI scores, as well as standard deviations for all tentative plans corresponding to patients with centrally and peripherally located tumors is also shown in Table 2. From the results displayed there, we observe that tentative plans for peripheral targets are in average closer to an ideal dosimetric plan than those for centrally located targets with respect to PTV coverage and sparing of OARs and RVRs. Figure 2 shows the DDI scores for all tentative plans for each group of patients considered, as well as the corresponding means and standard deviations. The small differences observed on the DDI values between competing plans in some patients are nearly imperceptible to the naked eye and thus difficult to detect and quantify. However, the DDI was able to significantly discriminate between similar-looking DVHs when scoring the corresponding tentative plans.

Fig. 2
figure 2

Plots of the DDI scores for each set of three tentative plans corresponding to the fifty brain cancer patients considered. Applied treatment plan in red, and discarded tentative plans in green and purple. a 25 centrally located benign tumors. b 15 peripherally located benign and malignant tumors. c 10 centrally located malignant tumors. Mean DDI scores (black solid lines) and standard deviations (blue solid lines) are also represented

Table 2 Scoring results corresponding to all tentative plans for the fifty brain cancer patients considered and differences between centrally and peripherally located tumors

Concerning the degree of coincidence between the applied treatment plan (ATP) and the tentative plan (DDIP) achieving the highest DDI score, we separately analyzed each group of patients classified by their tumor location and malignancy. In the cases of (i) centrally located benign tumors, (ii) centrally located malignant tumors and (iii) peripherally located benign and malignant tumors, we obtained that ATP and DDIP agree in 17 out of the 25, 6 out of the 10, and 9 out of the 15 patients respectively, see Fig. 2. Altogether, this results in an agreement between ATP and DDIP for 32 out of the 50 patients in this study under the assumption that PTV coverage and sparing of OARs and the RVR are considered as equally important.

To shed some light on the clinical choices made for the patients where ATP and DDIP did not coincide, we separately analyzed the DDI components for the planning target and critical structures in such cases. Figure 3 shows a comparison of the mean I T and \(\left (I_{O} + I_{R}\right)/2\) scores, which respectively estimate the quality of PTV coverage and critical structures sparing, for ATPs and the remaining competing plans for patients where ATP and DDIP did not coincide. We observe that, although the APTs were not in average the choices providing better PTV coverage, they yielded lower irradiation of critical structures. More precisely, in 15 out of the 18 cases of no coincidence \(\left (I_{O} + I_{R}\right)/2\) was higher in the APTs compared to other competing plans, while I T was higher in the ATPs of the remaining 3 patients.

Fig. 3
figure 3

Comparison of scoring results for the cases where ATP and DDIP did not coincide. The mean I T and joint I O and I R defined as (I O +I R )/2 scores and standard deviations are provided for the ATPs (label AP) and remaining tentative plans (label RP)

To gain further insight in one of the discrepancy cases, we selected a patient diagnosed with a centrally located Meningioma, where the ATP and DDIP did not coincide. Figure 4 shows the DVHs for the PTV and nearby OARs corresponding to the three different tentative plans prepared for that patient. OARs receiving low radiation doses as Eyes, Lacrimal Glands, Crystalline Lens and Retinas have not been represented for visualization simplicity, but all of them were considered in the calculation of the DDI scores. A quick glance at Fig. 4 reveals that in all tentative plans the target coverage is significantly high. However, what is unclear is which plan best provides target coverage and OARs sparing. In fact, when DVHs for all competing plans are superposed as shown in Fig. 4, the high number of OARs involved makes difficult to select the most suitable plan by mere inspection.

Fig. 4
figure 4

DVHs corresponding to the PTV (a), OARs (b) and RVR (c) for the three tentative plans of a patient diagnosed with a centrally located Meningioma. The D p on the PTV was set equal to 50.4 Gy in 28 daily fractions (5 days per week) of 1.8 Gy

In Table 3 the I T , I O , I R and DDI scores for such tentative plans are provided. We see that if PTV coverage is given maximum priority, the best choice corresponds to Plan 3, followed by Plans 2 and respectively 1. On the contrary, when lowering radiation doses on OARs and the RVR are selected as the primary clinical goal, the most suitable choice would be Plan 1, followed by Plans 2 and 3 respectively. Thus, consideration of partial DDI scores reveals that increasing (respectively decreasing) the target coverage was detrimental (respectively favorable) to sparing of OARs and the RVR. Although based on the DDI scores the Plan 2 provided the best balance between target coverage and sparing of critical structures, the Plan 1 was actually selected for treatment. This is consistent with highest priority having been assigned to sparing of OARs and the RVR at that time. For instance, a choice of relative weighting factors in (4) as w T =0.5, w O =1.0 and w R =1.0 would have suggested that choice.

Table 3 Scoring results for a patient diagnosed with a centrally located Meningioma


The availability of commercial TPS allows radiophysicists and radiation oncologists to virtually explore different tentative plans, and to perform successive modifications of a starting plan, for each radiotherapy patient. However, to this day the choice of treatment plans is mostly made by side-by-side comparison of DVHs and isodose curves for each tentative plan considered. Since differences between them are usually small, any selection made by inspection, as it is generally the case, may not be easy to justify on clinical terms. Therefore, the development of clinically based quantitative scoring indices could be quite useful as decision-support tools to assist in making clinical choices. For that purpose, several indices for scoring dosimetric plans have been proposed and shortly recalled at the ‘Background’ Section.

Bearing these facts in mind, we have introduced here a decision-support system, termed dose distribution index (DDI), to compare different tentative plans and to estimate the consequences of modifications on any treatment plan under consideration. The DDI makes use of the overall dosimetric information provided by the DVHs to simultaneously evaluate the target coverage and sparing of critical structures. Moreover, this evaluation tool can be easily extended to include several PTVs, as well as multiple dose coverage objectives on heterogeneous tumors. In our view, such features represent a comparative advantage with respect to previous combined indices. For instance, the unified dosimetry index (UDI) [25] makes use of four dosimetric objectives of dose coverage, conformity, homogeneity, and dose gradient on the PTV. Although the UDI provides a useful tool for ranking tentative plans, this figure of merit relies in specific dose-volume constraints instead to consider the overall dosimetric information provided by the DVHs. Moreover, UDI does not directly account for sparing of critical structures and does not seem easy to use when multiple target dose coverage objectives are specified. To consider a further example, a histogram analysis in radiation therapy (HART) computational tool has been proposed in [26]. This software incorporates various physically and biologically based indices for plan evaluation. Whereas HART is a complete and suitable program to compare tentative plans, the results are prone to the intrinsic limitations of the indices therein considered, as reported in the ICRU reports [13].

We point out that the DDI is very sensitive to small dosimetric variations between tentative plans and distinctly discriminates situations where DVHs look quite similar to the naked eye. Furthermore, it may be also used to compare tentative plans obtained from any TPS irrespective of tumor type and location, and can be formulated for any irradiation technique that provides DVHs. At the practical level, computing the DDI is straightforward, and that decision-support system can be easily displayed on TPS simultaneously to the visualization of the corresponding DVHs. In addition, it may readily provide the separate values I T , I O and I R , as well as those for each single OAR, which could be a useful complement to the DVH analysis.

Although the DDI can be used to globally assess the quality of any treatment plan, its main practical interest could be to compare tentative plans for a same patient. In fact, depending on diverse factors as tumor malignancy and location, curative or palliative purposes or overall patient’s medical history, a dosimetric plan can have a lower DDI score compared to other treatment plans for different patients, and yet be the best possible choice that satisfies concrete clinical requirements and treatment goals in that particular case. On the other hand, for a given patient different planners may well pursue different, and occasionally conflicting, goals. For instance, dose coverage and conformity on the PTV may be given highest priority in some cases, while palliative considerations may require of minimal damage to OARs and the RVR in others. DDI is not designed to change such priorities, but rather to help clinicians to precisely assess the consequences of their choices during treatment planning.

To give an idea of what DDI can accomplish, we have presented a retrospective study, corresponding to fifty patients treated at HUPH of centrally and peripherally located benign and malignant brain tumors. In each case, three tentative plans were simulated and recorded, out of which one was eventually delivered to the patient. We then compared the applied treatment plan (ATP) with that which provides the highest DDI score (DDIP). We observed that, under the assumption that PTV covering and OARs and RVRs sparing are equally important, ATP and DDIP coincided in 32 of the 50 cases. Discrepancies between ATP and DDIP in the remaining 18 cases may well result from different clinical priorities being attributed to critical structures during planning. In fact, in those cases ATP and DDIP also coincide, provided that different orders of priorities (for instance, deeming OARs sparing as the most important goal) are established instead.

We finally remark on the modular character of the quantitative index provided. In particular, DDI can be combined with multi-parameter optimization algorithms and included in objective functions during radiotherapy planning processes. In fact, considerable attention is being paid to the design of optimization methods to obtain dose distributions satisfying specific criteria as high tumor coverage and low radiation dose deposition on OARs [3, 3840]. However, no method has been found as yet to computationally find a plan that satisfies current clinical and logistical requirements in an optimal way. We expect that implementing the DDI decision-support system in current TPS would improve reliability standards, by precisely quantifying the impact of planner subjectivity in assessing and selecting radiotherapy treatment plans.


In this work a decision-support system has been developed to compare and evaluate different tentative plans during radiotherapy treatment planning. The proposed DDI is a straightforward and non-subjective system to assess and compare competing tentative plans. The DDI takes into account the overall dosimetric information provided by DVHs and not only that given by selected dose-volume constraints. Moreover, DDI can be easily computed irrespective of the TPS in clinical use, tumor type and localization, and irradiation technique. This decision-support tool is not meant to be a replacement for planners’ clinical decisions. Instead, it is intended to provide clinical personnel with a precise quantitative assessment of the impact of dosimetric changes over any starting tentative plan that they might consider.



three-dimensional conformal radiation therapy


applied treatment plan


clear cell carcinoma


conformity index


critical organ scoring index


clinical target volume


dose distribution index


tentative plan achieving the highest DDI score


dose-volume histogram


gross tumour volume


computed tomography


graphical user interface


homogeneity index


healthy tissue conformity index


healthy tissue overdosage factor


Hospital Universitario Puerta de Hierro, Madrid


International Commission on Radiation Units and Measurements


intensity-modulated radiation therapy


normal tissue complication probability


planning target volume


organ at risk


quality factor


Radiation Therapy Oncology Group


remaining volume at risk


treatment planning system


tumor control probability


unified dosimetry index


  1. Prescribing, recording, and reporting photon beam therapy. Bethesda, MD, USA: ICRU Publications; International Commission on Radiation Units and Measurements; 1993. ICRU Report 50.

  2. Prescribing recording and reporting photon-beam therapy. (Supplement to ICRU Report 50) Bethesda, MD, USA: ICRU Publications; International Commission on Radiation Units and Measurements; 1999. ICRU Report 62.

  3. International Commission of Radiation Units and Measurements. ICRU report 83: Prescribing, Recording, and Reporting Photon-Beam Intensity-Modulated Radiation Therapy (IMRT). J ICRU. 2010; 10(1):1–106. ISSN 1473-6691

    Article  Google Scholar 

  4. Moore KL, Brame RS, Low DA, Mutic S. Quantitative metrics for assessing plan quality. In: Seminars in radiation oncology, vol. 22. Elsevier: 2012. p. 62–9, doi:10.1016/j.semradonc.2011.09.005.

  5. Shaw E, Kline R, Gillin M, Souhami L, Hirschfeld A, Dinapoli R, et al. Radiation therapy oncology group: radiosurgery quality assurance guidelines. Int J Radiat Oncol Biol Phys. 1993; 27(5):1231–9.

    Article  PubMed  CAS  Google Scholar 

  6. Feuvret L, Noël G, Mazeron JJ, Bey P. Conformity index: a review. Int J Radiat Oncol Biol Phys. 2006; 64(2):333–42.

    Article  PubMed  Google Scholar 

  7. Yoon M. A new homogeneityx index based on the statistical analysis of dose volume histogram. J Appl Clin Med Phys. 2007; 8(2):9–17.

    PubMed  Google Scholar 

  8. Thilmann C, Zabel A, Milker-Zabel S, Schlegel W, Wannenmacher M, Debus J. Number and orientation of beams in inversely planned intensity-modulated radiotherapy of the female breast and the parasternal lymph nodes. Am J Clin Oncol. 2003; 26(5):136–43.

    Article  Google Scholar 

  9. Wu V, Sham J, Kwong D. Inverse planning in three-dimensional conformal and intensity-modulated radiotherapy of mid-thoracic oesophageal cancer. Br J Radiol. 2014; 77(919):568–72.

    Article  Google Scholar 

  10. Park Y-K, Park S, Wu H-G, Kim S. A new plan quality index for dose painting radiotherapy. J Appl Clin Med Phys. 2014; 15(4):316–5.

    Google Scholar 

  11. Kataria T, Sharma K, Subramani V, Karrthick K, Bisht SS. Homogeneity index: An objective tool for assessment of conformal radiation treatments. J medical physics/Association of Med Physicists of India. 2012; 37(4):207.

    Google Scholar 

  12. Leung LHT, Kan MWK, Cheng ACK, Wong WKH, Yau CC. A new dose–volume-based plan quality index for imrt plan comparison. Radiother Oncol. 2007; 85(3):407–17.

    Article  PubMed  Google Scholar 

  13. Knöös T, Kristensen I, Nilsson P. Volumetric and dosimetric evaluation of radiation treatment plans: radiation conformity index. Int J Radiat Oncol Biol Phys. 1998; 42(5):1169–76.

    Article  PubMed  Google Scholar 

  14. Wagner TH, Bova FJ, Friedman WA, Buatti JM, Bouchet LG, Meeks SL. A simple and reliable index for scoring rival stereotactic radiosurgery plans. Int J Radiat Oncol Biol Phys. 2003; 57(4):1141–9.

    Article  PubMed  Google Scholar 

  15. Lomax NJ, Scheib SG. Quantifying the degree of conformity in radiosurgery treatment planning. Int J Radiat Oncol Biol Phys. 2003; 55(5):1409–19.

    Article  PubMed  Google Scholar 

  16. Paddick I. A simple scoring ratio to index the conformity of radiosurgical treatment plans: technical note. J Neurosurg. 2000; 93(Supplement 3):219–22.

    PubMed  Google Scholar 

  17. Wu Q-RJ, Wessels B, Einstein D, Maciunas R, Kim E, Kinsella T. Quality of coverage: conformity measures for stereotactic radiosurgery. J Appl Clin Med Phys. 2003; 4(4):374–81.

    Article  PubMed  Google Scholar 

  18. Cheung FW, Law MY. A novel conformity index for intensity modulated radiation therapy plan evaluation. Med Phys. 2012; 39(9):5740–56.

    Article  PubMed  Google Scholar 

  19. Grimm J, LaCouture T, Croce R, Yeo I, Zhu Y, Xue J. Dose tolerance limits and dose volume histogram evaluation for stereotactic body radiotherapy. J Appl Clin Med Phys. 2011; 12(2):267–92.

    Google Scholar 

  20. Bentzen SM, Constine LS, Deasy JO, Eisbruch A, Jackson A, Marks LB, et al. Quantitative analyses of normal tissue effects in the clinic (quantec): an introduction to the scientific issues. Int J Radiat Oncol Biol Phys. 2010; 76(3):3–9.

    Article  Google Scholar 

  21. Menhel J, Levin D, Alezra D, Symon Z, Pfeffer R. Assessing the quality of conformal treatment planning: a new tool for quantitative comparison. Phys Med Biol. 2006; 51(20):5363.

    Article  PubMed  CAS  Google Scholar 

  22. Van’t Riet A, Mak AC, Moerland MA, Elders LH, van der Zee W. A conformation number to quantify the degree of conformality in brachytherapy and external beam irradiation: application to the prostate. Int J Radiat Oncol Biol Phys. 1997; 37(3):731–6.

    Article  PubMed  Google Scholar 

  23. Baltas D, Kolotas C, Geramani K, Mould RF, Ioannidis G, Kekchidi M, et al. A conformal index (coin) to evaluate implant quality and dose specification in brachytherapy. Int J Radiat Oncol Biol Phys. 1998; 40(2):515–24.

    Article  PubMed  CAS  Google Scholar 

  24. Miften MM, Das SK, Su M, Marks LB. A dose-volume based tool for evaluating and ranking imrt treatment plans. J Appl Clin Med Phys. 2004; 5(4):1–14.

    Article  PubMed  Google Scholar 

  25. Akpati HC, Kim C, Kim B, Park T, Meek A. Unified dosimetry index (udi): a figure of merit for ranking treatment plans. J Appl Clin Med Phys. 2008; 9(3):99–108.

    Article  Google Scholar 

  26. Pyakuryal A, Myint WK, Gopalakrishnan M, Jang S, Logemann JA, Mittal BB. A computational tool for the efficient analysis of dose-volume histograms for radiation therapy treatment plans. J Appl Clin Med Phys/American College of Medical Physics. 2010; 11(1):3013.

    Google Scholar 

  27. Meyer RR, Zhang HH, Goadrich L, Nazareth DP, Shi L, D’Souza WD. A multiplan treatment-planning framework: a paradigm shift for intensity-modulated radiotherapy. Int J Radiat Oncol Biol Phys. 2007; 68(4):1178–89.

    Article  PubMed  Google Scholar 

  28. Tomé WA, Fowler JF. Selective boosting of tumor subvolumes. Int J Radiat Oncol Biol Phys. 2000; 48(2):593–9.

    Article  PubMed  Google Scholar 

  29. Bentzen M, Sl Tucker S. Quantifying the position and steepness of radiation dose-response curves. Int J Radiat Biol. 1997; 71(5):531–42.

    Article  PubMed  CAS  Google Scholar 

  30. Zaider M, Minerbo G. Tumour control probability: a formulation applicable to any temporal protocol of dose delivery. Phys Med Biol. 2000; 45(2):279.

    Article  PubMed  CAS  Google Scholar 

  31. Dawson A, Hillen T. Derivation of the tumour control probability (tcp) from a cell cycle model. Comput Math Methods Med. 2006; 7(2–3):121–41.

    Article  Google Scholar 

  32. Niemierko A, Goitein M. Modeling of normal tissue response to radiation: the critical volume model. Int J Radiat Oncol Biol Phys. 1993; 25(1):135–45.

    Article  PubMed  CAS  Google Scholar 

  33. Kehwar T. Analytical approach to estimate normal tissue complication probability using best fit of normal tissue tolerance doses into the ntcp equation of the linear quadratic model. J Cancer Res Ther. 2005; 1(3):168–79.

    Article  PubMed  CAS  Google Scholar 

  34. Lyman JT. Complication probability as assessed from dose-volume histograms. Radiat Res. 1985; 104(2s):13–19.

    Google Scholar 

  35. Burman C, Kutcher G, Emami B, Goitein M. Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys. 1991; 21(1):123–35.

    Article  PubMed  CAS  Google Scholar 

  36. Luxton G, Keall PJ, King CR. A new formula for normal tissue complication probability (ntcp) as a function of equivalent uniform dose (eud). Phys Med Biol. 2008; 53(1):23.

    Article  PubMed  Google Scholar 

  37. López Alfonso JC, Jagiella N, Núñez L, Herrero MA, Drasdo D. Estimating dose painting effects in radiotherapy: a mathematical model. PloS ONE. 2014; 9(2):89380.

    Article  Google Scholar 

  38. Webb S. Intensity-Modulated Radiation Therapy, Series in Medical Physics and Biomedical Engineering. Bristol (UK): CRC Press; 2001.

    Book  Google Scholar 

  39. Alfonso JCL, Buttazzo G, García-Archilla B, Herrero MA, Núñez L. A class of optimization problems in radiotherapy dosimetry planning. Discr Cont Dyn Systems B. 2012; 17:1651–72.

    Article  Google Scholar 

  40. Alfonso J, Buttazzo G, García-Archilla B, Herrero M, Núñez L. Selecting radiotherapy dose distributions by means of constrained optimization problems. Bull Math Biol. 2014; 76(5):1017–44.

    Article  PubMed  CAS  Google Scholar 

Download references


JCLA, MAH and LN have been partially supported by the Spanish MINECO Grant MTM 2014-53156-P. JCLA gratefully acknowledges the support of the German Excellence Initiative via the Cluster of Excellence EXC 1056 Center for Advancing Electronics Dresden (cfAED) at the Technische Universität Dresden and the German Federal Ministry of Education and Research (BMBF) for the eMED project SYSIMIT (01ZX1308D). We acknowledge support by the German Research Foundation and the Open Access Publication Funds of the TU Dresden.

Author information

Authors and Affiliations


Corresponding author

Correspondence to J. C. L. Alfonso.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JCLA and LN conceived the study and were responsible for data acquisition. JCLA, MAH and LN carried out the study, analyzed the data and drafted the manuscript. All authors read and approved the final manuscript.


This material has never been published and is not currently under evaluation in any other peer-reviewed publication or journal.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

L. Alfonso, J.C., Herrero, M.A. & Núñez, L. A dose-volume histogram based decision-support system for dosimetric comparison of radiotherapy treatment plans. Radiat Oncol 10, 263 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: