Developing a dose-volume histogram computation program for brachytherapy

Methods, software and datasets to verify DVH calculations against analytical values: Twenty years late(r)

PURPOSE

The authors designed data, methods, and metrics that can serve as a standard, independent of any software package, to evaluate dose-volume histogram (DVH) calculation accuracy and detect limitations. The authors use simple geometrical objects at different orientations combined with dose grids of varying spatial resolution with linear 1D dose gradients; when combined, ground truth DVH curves can be calculated analytically in closed form to serve as the absolute standards.

METHODS

dicom RT structure sets containing a small sphere, cylinder, and cone were created programmatically with axial plane spacing varying from 0.2 to 3 mm. Cylinders and cones were modeled in two different orientations with respect to the IEC 1217 Y axis. The contours were designed to stringently but methodically test voxelation methods required for DVH. Synthetic RT dose files were generated with 1D linear dose gradient and with grid resolution varying from 0.4 to 3 mm. Two commercial DVH algorithms-pinnacle (Philips Radiation Oncology Systems) and PlanIQ (Sun Nuclear Corp.)-were tested against analytical values using custom, noncommercial analysis software. In Test 1, axial contour spacing was constant at 0.2 mm while dose grid resolution varied. In Tests 2 and 3, the dose grid resolution was matched to varying subsampled axial contours with spacing of 1, 2, and 3 mm, and difference analysis and metrics were employed: (1) histograms of the accuracy of various DVH parameters (total volume, Dmax, Dmin, and doses to % volume: D99, D95, D5, D1, D0.03 cm(3)) and (2) volume errors extracted along the DVH curves were generated and summarized in tabular and graphical forms.

RESULTS

In Test 1, pinnacle produced 52 deviations (15%) while PlanIQ produced 5 (1.5%). In Test 2, pinnacle and PlanIQ differed from analytical by >3% in 93 (36%) and 18 (7%) times, respectively. Excluding Dmin and Dmax as least clinically relevant would result in 32 (15%) vs 5 (2%) scored deviations for pinnacle vs PlanIQ in Test 1, while Test 2 would yield 53 (25%) vs 17 (8%). In Test 3, statistical analyses of volume errors extracted continuously along the curves show pinnacle to have more errors and higher variability (relative to PlanIQ), primarily due to pinnacle's lack of sufficient 3D grid supersampling. Another major driver for pinnacle errors is an inconsistency in implementation of the "end-capping"; the additional volume resulting from expanding superior and inferior contours halfway to the next slice is included in the total volume calculation, but dose voxels in this expanded volume are excluded from the DVH. PlanIQ had fewer deviations, and most were associated with a rotated cylinder modeled by rectangular axial contours; for coarser axial spacing, the limited number of cross-sectional rectangles hinders the ability to render the true structure volume.

CONCLUSIONS

The method is applicable to any DVH-calculating software capable of importing dicom RT structure set and dose objects (the authors' examples are available for download). It includes a collection of tests that probe the design of the DVH algorithm, measure its accuracy, and identify failure modes. Merits and applicability of each test are discussed.

Quality of interstitial PDR-brachytherapy-implants of head-and-neck-cancers: predictive factors for local control and late toxicity?

BACKGROUND AND PURPOSE

Parameters and indices related to the implant geometry in use for describing the quality of volume implants in interstitial brachytherapy were developed on the basis of LDR-brachytherapy. The aim of our study was to evaluate their usefulness for predicting late toxicity and local control in the PDR-brachytherapy of head-and-neck-tumors.

PATIENTS AND METHODS

Between January 2000 and October 2004, 210 patients were treated with PDR-brachytherapy which was administered either postoperatively or as definitive treatment. Brachytherapy was used as sole treatment in some cases while in others a combination with EBRT was used. For assessment of quality of implants we analyzed the following indices and parameters using the univariate chi2 test and multivariate logistic regression analysis: V85, V120 and V150 (volume enclosed by the surface of the 85%-, 120%- and 150%-isodose), UI (uniformity index), QI (quality index), HI (homogeneity index), VGR (volume gradient ratio), DNR (dose non-uniformity ratio), LD (low dose), HD (high dose), PD (peak dose) and the intersource spacing.

RESULTS

After a median follow-up of 24 months (4-50) the rate of - usually transient - soft tissue necrosis (STN) was 11%, osteoradionecrosis (ORN) was seen in 7.6% of cases and local relapse occurred in 7% of cases. Univariate analysis shows a significant influence on the development of soft tissue necrosis for V85, and on osteoradionecrosis for HD and PD. In the multivariate analysis a correlation between soft tissue necrosis and QI was found. For local control a correlation with QI, VGR and minimal tube distance was found using univariate analysis.

CONCLUSIONS

Using interstitial PDR-brachytherapy in head-and-neck-tumors the probability of local control and of the development of soft tissue necrosis or osteoradionecrosis is dependent on dose and volume parameter like the volume of the reference isodose, the high and peak dose values, on the homogeneity of the dose distribution, quantified by the quality index or the volume gradient ratio as well on the minimal tube distance.

The application of a radiotherapy-planning system: consideration from a case of intraoral squamous cell carcinoma that lacked a control

OBJECTIVE

The purpose of this case study was to present an applied example of a radiotherapy-planning system used during the treatment of a malignant tumor in the oral and maxillofacial region.

STUDY DESIGN

The various radiotherapy modalities were performed on a patient with oral squamous cell carcinoma of multiple metachronous recurrences. Estimated radiation dose-distribution curves for each radiotherapy modality were computed by using a commercially available radiotherapy-planning system. A personal computer was used to make the superimposed radiation dose-distribution curve.

RESULTS

The 3-dimensional dose-distribution curves were determined with the radiotherapy-planning system. In addition, with the use of superimposition of the dose-distribution curves from all sources of radiotherapy, it was possible to estimate regions receiving extremely high radiation dosages. This information serves as a road map to potential postradiation complication sites at follow-up examinations.

CONCLUSION

The radiotherapy-planning system is very useful for the evaluation of a radiotherapy dosage to treat malignant tumors in the oral and maxillofacial region.

Comparison of brachytherapy strategies based on dose-volume histograms derived from quantitative intravascular ultrasound

PURPOSE

We present in this paper the comparison, by simulation, of different treatment strategies based either on beta- or gamma-sources, both with and without a centering device. Ionizing radiation to prevent restenosis is an emerging modality in interventional cardiology. Numerous clinical studies are presently being performed or planned, but there is variability in dose prescription, and both gamma- and beta-emitters are used, leading to a wide range of possible dose distributions over the arterial vessel wall. This paper discusses the potential merits of dose-volume histograms (DVH) based on three-dimensional (3-D) reconstruction of electrocardiogram (ECG)-gated intravascular ultrasound (IVUS) to compare brachytherapy treatment strategies.

MATERIALS AND METHODS

DVH describe the cumulative distribution of dose over three specific volumes: (1) at the level of the luminal surface, a volume was defined with a thickness of 0.1 mm from the automatically detected contour of the highly echogenic blood-vessel interface; (2) at the level of the IVUS echogenic media-adventitia interface (external elastic lamina [EEL]), an adventitial volume was computed considering a 0.5-mm thickness from EEL; and (3) the volume encompassed between the luminal surface and the EEL (plaque + media). The IVUS data used were recorded in 23 of 31 patients during the Beta Energy Restenosis Trial (BERT) conducted in our institution.

RESULTS

On average, the minimal dose in 90% of the adventitial volume was 37 +/- 16% of the prescribed dose; the minimal dose in 90% of the plaque + media volume was 58 +/- 24% and of the luminal surface volume was 67 +/- 31%. The minimal dose in the 10% most exposed luminal surface volume was 296 +/- 42%. Simulations of the use of a gamma-emitter and/or a radioactive source train centered in the lumen are reported, with a comparison of the homogeneity of the dose distribution.

CONCLUSIONS

It is possible to derive DVH from IVUS, to evaluate the dose delivered to different parts of the coronary wall. This process should improve our understanding of the mechanisms of action of brachytherapy.