
ORIGINAL ARTICLE 



Year : 2010  Volume
: 35
 Issue : 4  Page : 197201 

Dosevolume histogram quality assurance for linacbased treatment planning systems
Michael S Gossman^{1}, Morris I Bank^{2}
^{1} TriState Regional Cancer Center, Medical Physics Section, 706, 23rd Street, Ashland, Kentucky, USA ^{2} St. Vincent Oncology Center, 8301 Harcourt Road, Indianapolis, Indiana, USA
Date of Submission  26Oct2009 
Date of Decision  02Dec2009 
Date of Acceptance  03Apr2010 
Date of Web Publication  19Oct2010 
Correspondence Address: Michael S Gossman TriState Regional Cancer Center, Medical Physics Section, 706 23rd Street Ashland, Kentucky 41101 USA
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.71759
Abstract   
Dosevolume histograms provide key information to radiation oncologists when they assess the adequacy of a patient treatment plan in radiation therapy. It is important therefore that all clinically relevant data be accurate. In this article we present the first quality assurance routine involving a direct comparison of planning system results with the results obtained from independent hand calculations. Given a known threedimensional (3D) structure such as a parallelepiped, a simple beam arrangement, and known physics beam data, a timeefficient and reproducible method for verifying the accuracy of volumetric statistics (DVH) from a radiation therapy treatment planning system (TPS) can be employed rapidly, satisfying the QA requirements for (TPS) commissioning, upgrades, and annual checks. Using this method, the maximum disagreement was only 1.7% for 6 MV and 1.3% for 18 MV photon energies. The average accuracy was within 0.6% for 6 MV and 0.4% for 18 MV for all depthdose results. A 2% disagreement was observed with the treatment planning system DVH from defined volume comparison to the known structure dimensions.
Keywords: Dose, dosevolume, dosevolume histogram, histogram, quality assurance, treatment planning system, Volume
How to cite this article: Gossman MS, Bank MI. Dosevolume histogram quality assurance for linacbased treatment planning systems. J Med Phys 2010;35:197201 
How to cite this URL: Gossman MS, Bank MI. Dosevolume histogram quality assurance for linacbased treatment planning systems. J Med Phys [serial online] 2010 [cited 2019 Feb 19];35:197201. Available from: http://www.jmp.org.in/text.asp?2010/35/4/197/71759 
Introduction   
A qualitative isodose distribution superpositioned over computed tomographic (CT) data is often insufficient to provide qualitative data for the Radiation Oncologist to determine the adequacy of a patient treatment plan. A more quantitative result is necessary which modern treatment planning systems provide in statistical analysis plots known as dosevolume histograms (DVHs). These plots describe the dose throughout the volume of each structure contoured within the plan. Since it is important to have accurate statistics for clinical use, it is important to test the system against an independent source. Concise summary data verifications for the 3D dose distributions are needed. ^{[1]} A literature search shows a lack of any quality control procedures performed and described to ensure that computerized DVH analyses are correct. ^{[2],[3]} Our research provides this methodology and discussion on the accuracy of its use.
Often the statistics of the dosetovolume relationship are the deciding factor when one plan is chosen over another. ^{[4]} However, this is not always true since the DVH has no spatial resolution. ^{[5]} The criterion of dosevolume constraints may also be associated with the initial prescription. In order to verify dosevolume (DVH) statistics from a treatment planning system, a known volume must be used. A complete understanding of the mathematical components involved, the 'untimely' cost of analysis, and the methodic conceptual construction of a valid quality assurance test has kept such research from being published. ^{[1]} Quality assurance of the DVH has been strongly recommended though the authors do not explain exactly how this is to be done. ^{[3],[6],[7],[8]} This article presents a reproducible and efficient DVH quality assurance process, with a clear stepbystep description of the methodology.
Using a commercial radiation therapy treatment planning system and given an artificial CT data set representing a homogeneous water phantom, one may construct a single treatment plan that can be utilized annually or when the TPS is updated. In the method discussed in this article the majority of effort is involved with the initial contouring of a structure with a known volume and the associated singlefield setup. Here, we chose a rectangular parallelepiped or 'block.' The results of volumetric dose are mathematically calculated from depthdose equations and then directly compared to the results from the treatment planning system DVH plots. Being spreadsheet based, this methodology is repeatable when needed, merely requiring a recalculation of the computer treatment plan. From each new calculation, DVH data can be extracted and input for direct comparison. The method described complements the focus of attention on DVH errors, by avoiding the effects due to grid size and sampling frequency problems. Further, it provides volumetric and geometric resolution, compares absolute and relative dose results, and tests statistics generated for dose maxima, minima, and partial volume doses within the contoured structure.
Materials and Methods   
Object and beam definitions
The Eclipse Treatment Planning System (Varian Medical Systems, Inc.) was used to define the volumes of interest and for dose delivery simulation. A homogeneous phantom set to density of 1.0 resembling tissue equivalent media was chosen for use with all treatment planning experiments. The phantom consisted of an artificial data set created within the treatment planning software. The dimensions of the phantom data set were suitable for radiation fields of all sizes at 40 Χ 40 cm ^{2} . A slice spacing of 1.25 mm was chosen for data usage.
The Eclipse software was commissioned to compute modeled dose from a 21EX Clinac (Varian Medical Systems, Inc.) using the Anisotropic Analytical Algorithm (AAA) version 8.6 in External Beam Planning Software build 8.6.17. Since the 21EX particle accelerator was engineered to generate photon beams at 6 MV and 18 MV, both energy modalities were used. Each beam was calibrated in accordance with the clinical reference dosimetry standard: the American Association of Physicists in Medicine (AAPM) TG51 formalism. ^{[9]} For all photon beams at this facility, a calibration was performed such that 1.00 cGy/MU was achieved at the sourceaxis distance (SAD) of the LINAC and at the depth where the maximum dose was achieved.
Using the contouring workspace in the software, an object contour was created on one superior slice of the CT data of the water phantom. The contour was given dimensions of 2 cm Χ 10 cm. The square object drawing was copied and reapplied to slices consecutively inferior in the phantom material through an additional 2 cm. A total of 17 slices were contoured using the 'copy and paste' tool, since (1 slice/0.125 cm) Χ (2 cm − 0 cm) + 1 = 17 slices. The resulting structure is a rectangular solid known as a parallelepiped, with dimensions of 2 Χ 2 Χ 10 cm ^{3} . This 40 cm ^{3} known structural volume was used for dose analysis. The phantom was assigned a Hounsfield unit (HU) value of 0, making it a unit density homogeneous phantom material.
After completing the contouring required, a treatment field was assigned to the CT data set. The isocenter of the field was centered at a depth of 10 cm in the water phantom (5 cm deep in the parallelepiped). The calculation reference point was at the same point. The prescribed dose to that point was 100 cGy. The field size of the beam was 30 Χ 30 cm ^{2} . For Varian standard scaling, the gantry angle was 180, with the collimator angle and the couch angle each set at 180. The field was copied so that both 6 MV and 18 MV results could be tested independently. The volume calculation grid at 0.15 cm was chosen to be small, since large voxel sizes result in less accurate results. ^{[10]} Once all parameters in the plan were correctly verified, the dose throughout the 3D volume was calculated..
Hand calculations
The resulting computation of dose is qualitatively exhibited in [Figure 1]. It is again shown in the following illustration with a focus on the region comprising the contoured structural volume [Figure 2]. Isodose lines ranging from 75% to 125% (in 5% increments) are displayed in the [Figure 1] and [Figure 2] as seen from the planning software for 18 MV photons. The depth of dose occurrence is analyzed for accuracy against the known dose data commissioned for the planning system. The accuracy of the analysis is strictly a function of knowing the depthdose for that point. Further knowledge of the precise volume of the contour receiving such dose at each depth permits taking into account field edge horn effects created by the flattening filter of the accelerator. The calculation of the number of monitor units necessary to arrive at this prescription dose is shown in Equation 1.
 Figure 1: Axial view isodose distribution indicating the position of the rectangular volume within the water phantom. The isodose levels represent the percentage of dose, relative to 100 cGy at 10 cm phantom depth (5 cm structure depth), for an 18 MV beam having normal geometry and a field size of 30 × 30 cm2 at the SAD.
Click here to view 
 Figure 2: As in the previous illustration [Figure 1], this is a zoomin on the distribution of dose that is seen inside the structure in the axial view
Click here to view 
The number of monitor units (MU) is determined from the quotient of the prescribed tumor dose (TD) with the product of four variables: the machines calibration doserate (Output), the tumor maximum dose ratio (TMR_{d} , _{FSd}), the collimator head scatter factor (Sc_{FScol} ), and the phantom scatter factor (Sp_{FSd} ). For this exercise, the tumor is the contoured block. The calculation depth is denoted d, with FS_{col} representing the collimator field size and FS_{d} representing the field size at that depth. For the depth of 10 cm in the water phantom (5 cm within the block contour) the dose to the point of interest is 100 cGy. For all other depths (d'), the percentage depthdose (%DD) along the central ray may be calculated from Equation 2:
The sourceaxis distance (SAD) defined as the distance from the source of radiation to the center of gantry rotation on the accelerator is precisely 100 cm. The sourcesurface distance (SSD) is the distance from the source to the phantom surface. The distance was fixed in this study at 90 cm.
It is expected that the 100% isodose line will pass through the centroid of the contoured parallelepiped. As shown in [Figure 2], the horn effects of the planning field clearly present a distribution of dose that is curved, such that the amount of volume receiving 100% of the dose is more than half of the object volume. Although the field generally exhibits flatness and symmetry within industry standards, this is not to be confused with the intent of analytically accounting for the change in volume receiving the dose inside of the isodose line associated at that depth. Offaxis ratios (OAR) from scanned data were used to more accurately quantify the contour volume receiving the dose at each depth reviewed. These were used to correct the geometric volume of the object for accuracy. From the dimensions of the contour described, it was possible to study the accuracy of dose mathematically from scanned data vs the treatment planning system's DVH results. Dose depth studies here range from a phantom depth of 6 cm (1 cm within the block contour) to a phantom depth of 13 cm (8 cm within the contour).
Computer calculations
The statistical DVH computed by the treatment planning software was created through a lengthy interpolation process. First, the shape of the parallelepiped was registered in 3D coordinates space along with the dose values for each voxel. By binning each voxel, a statistical account of the dose received to each segment of the structure can be determined by weighting the dose of each bin with neighboring bins throughout the entire object space. This iterative weighting is referred to as the dose matrix coordinates. It is within the dose matrix that differences in the binned object volume calculation are observed as a source of error. Imported data from CT acquisition is a common cause of resampling errors, where the dose matrix and the dose grid chosen for computation are considerably different. For quality assurance studies such as this, we recommend the creation of an artificial data set to avoid the significant levels of error that could be otherwise introduced. Finally, the cumulative DVH is generated as an integral of the sampled dose over the interpolated structures. ^{[11]}
Results   
Using a spreadsheet for calculations, all of the parameters discussed here were imbedded. The dose was determined at each depth using Equations 1 and 2. Since 100% of the dose is normalized to a single point in the center of the block, half of it can be expected to receive the full prescribed dose along with a subtle increase due to horn effects that is accounted for using OARs. Thus, the dose was renormalized at each depth in accordance to the prescription. At 10 cm in the phantom (5 cm in the block contour) the dose should be 100 cGy. All other depthdoses are relative to the dose at this point. Therefore, each percentage depthdose was renormalized with respect to the weight point. Points anterior to this resulted in doses greater than 100 cGy, whereas points downstream resulted in doses less than 100 cGy.
The fractional volume known was mathematically determined at each depth. Dose to each tenth fraction of volume was calculated from knowledge of the isodose line traversed at each depth. Then, while taking this result in product with a percentage increase in the dose attributed to the horn effect within the beam profile, the resulting percentage of volume enveloping the isodose line at each depth is evident.
The percentage of volume encompassing the dose at each depth as well as the stated renormalized percentage depthdose at each volumetric depth was calculated and directly compared to the results from the planning system DVH plots. Both Xray energy modalities were reviewed independently and identically in the planning system. The treatment planning system DVH statistics for 6 MV and 18 MV are presented in [Figure 3]. A sample of the resulting calculation routine screen is presented in [Table 1].  Figure 3: Dosevolume histogram statistical analysis from the treatment planning system for 6 MV (dashed) and 18 MV (solid)
Click here to view 
 Table 1: Spreadsheet of independent dosevolume calculation correlation to treatment planning system DVH results. The independent calculations appear on the left of the spreadsheet, with DVH results on the right along with their overall agreement. Both 6 MV and 18 MV results are presented.
Click here to view 
Data presented in the table on the left are utilized in the calculation of dose at the prescription depth point of interest, which is shown in the middle of the spreadsheet. From Equation 2, the ratio of the calculated point dose to the 100 cGy prescribed dose is indicated as %DDd ^{calc} . The percentage of geometric volume receiving that dose is given as %Vd ^{Calc} . The object volume and beam assignment were chosen such that the increase in volume is in the anterior to posterior direction. This is the same directional path of the radiation beam. Therefore, %Vd ^{Calc} is determined knowing the geometric parallelepiped volume encompassed at each depth along with a residual amount to account for the curvature of the isodose line. The additional volume correction can be determined by averaging the OARs from the central axis through the halfwidth of the object diameter. The object is only 5 cm wide, so offaxis factors are only considered up to 2.5 cm. At 6 MV and 18 MV, the average OAR was determined to be 1.001 and 1.000, respectively.
Data analyzed from the treatment planning software are revealed in [Table 1] on the right. The DVH computed percentage depthdose at the calculation point is denoted %DDd ^{DVH} . The percentage of reconstructed volume receiving that dose is given as %Vd ^{DVH} . The calculated percentage depthdoses and the software computed percentage depthdoses were compared relative to the volume encompassed at each depth in the object. The results yield a maximum deviation of −1.7% for 6 MV and −1.3% for 18 MV. Equation 3 below indicates the calculation performed.
As compared to the manually created known volume of 40.0 cm ^{3} , the DVH indicated a total volume of 40.8 cm ^{3} . This is within tolerance at 2.0%. All data statistics were seen to satisfy requirements of analysis for the treatment planning software utilized here. In general, if independent calculations differ by more than 5.0%, the disparity should be investigated further and resolved. ^{[12]}
Discussion and Conclusions   
For quality assurance in medical physics practice, it is often advantageous to employ reproducible routines for calculations. Such is the case for software treatment planning system annual quality assurance testing. Here, it is shown how a routine can be created to test the treatment planning system on its ability to accurately calculate dose to a manually created parallelepiped volume and properly quantify the volumetric statistics from its computations. Dosevolume histograms are readily available in most clinically available TPS software. As important as statistical data are to radiation oncologists in this era, it is beneficial to understand the accuracy of such quantitative results upon which clinical decisions are based. The methods described here are reproducible, accurate, and serves as an effective tool to perform periodic quality assurance of the treatment planning system. Depending on the accuracy of the scanning data employed into the treatment planning system and the algorithms to which the DVHs are computed, the accuracy of the system could ideally be verified to within a few percent.
The planning DVH ties together the percentage of total structure volume to the percentage of relative dose it receives. The ratio of the two can be directly compared to the results from independent calculations. Here, a maximum deviation of 1.7% was observed and accepted for both 6 MV and 18 MV Xrays in this quantitative study. The accuracy average difference was 0.6% for 6 MV and 0.4% for 18 MV. The total volume difference for the 40 cm ^{3} block structure was 2%.
This research describes a method which is ideal for validating the accuracy of the DVH statistical analysis in a commercially available radiation therapy treatment planning system. It is not necessary to recreate the entirety of the plan for the next test scheduled. It is only necessary to recalculate the dose using the system algorithm of choice. Further, the implementation of a spreadsheet to expedite the presented calculations can permit complete dosevolume histogram quality assurance in a timely manner. This research provides a sound basis for performing quality assurance tests for treatment planning system dosevolume histograms with successful results. Such quality assurance should be performed at the time of commissioning, annually, and after any major algorithm upgrades.
References   
1.  Kessler ML, Ten Haken RK, Fraass BA, McShan DL. Expanding the use and effectiveness of dosevolume histograms for 3D treatment planning. I: Integration of 3D dosedisplay. Int J Radiat Oncol Biol Phys 1994;29:112531. 
2.  Panitsa E, Rosenwald JC, Kappas C. Quality control of dose volume histogram computation characteristics of 3D treatment planning systems. Phys Med Biol 1998;43:280716. [PUBMED] [FULLTEXT] 
3.  Drzymala RE, Mohan R, Brewster L, Chu J, Goitein M, Harms W, et al. Dosevolume histograms. Int J Radiat Oncol Biol Phys 1991;21:718. [PUBMED] [FULLTEXT] 
4.  Van Dyk J, Barnett RB, Cygler JE, Shragge PC. Commissioning and quality assurance of treatment planning computers. Int J Radiat Oncol Biol Phys 1993;26:26173. [PUBMED] 
5.  Ezzell GA, Galvin JM, Low D, Palta JR, Rosen I, Sharpe MB, et al. Guidance document on delivery, treatment planning, and clinical implementation of IMRT: Report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee. Med Phys 2003;30:2089115. [PUBMED] 
6.  Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall G, Stern R, et al. American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: Quality assurance for clinical radiotherapy treatment planning. Med Phys 1998;25:1773829. [PUBMED] 
7.  Purdy JA. Acceptance testing and quality assurance of radiation treatment planning systems Seminar: Evolving Strategies in Radiation Oncology. Rome, Italy; 2004 . 
8.  International Atomic Energy Agency (IAEA). Commissioning and Quality Assurance of Computerized Treatment Planning Systems for Radiation Treatment of Cancer. Technical Reports Series (TRS) No. 430. Vienna, Austria; 2004. 
9.  Almond PR, Biggs PJ, Coursey BM, Hanson WF, Huq MS, Nath R, et al. AAPM's TG51 protocol for clinical reference dosimetry of highenergy photon and electron beams. Med Phys 1999;26:184770. [PUBMED] 
10.  Corbett JF, Jezioranski J, Crook J, Yeung I. The effect of voxel size on the accuracy of dosevolume histograms of prostate ^{125} I seed implants. Med Phys 2002;29:10036. [PUBMED] 
11.  Varian Medical Systems, Inc. Eclipse Algorithms Reference Guide. Palo Alto, CA; 2007 . 
12.  Kutcher GJ, Coia L, Gillin M, Hanson WF, Leibel S, Morton RJ, et al. Comprehensive QA for radiation oncology: Report of AAPM Radiation Therapy Committee Task Group 40. Med Phys 1994;21:581618. [PUBMED] 
[Figure 1], [Figure 2], [Figure 3]
[Table 1]
