A comparison of conventional 'forward planning' with inverse planning for 3D conformal radiotherapy of the prostate

Novel hybrid treatment planning approach for irradiation a pediatric craniospinal axis

This study presents a new treatment planning approach merging 3D-CRT and VMAT fields into a hybrid treatment plan (HybTP), in order to achieve an optimum dose coverage of the planning target volume (PTV) and protection of OAR. Craniospinal axis irradiation (CSI) treated with 3D conformal radiotherapy (3D-CRT) is associated with high doses to the heart and eye lenses but provides better sparing of lungs and kidneys compared to volumetric modulated arc therapy (VMAT). VMAT treatment spares eye lenses and the heart, but lungs and kidneys are not as effective as 3D-CRT. Thus, a combination of both techniques (HybTP) may be optimal in sparing all these organs at risk (OAR). The results of HybTP are compared with helical tomotherapy (HT), intensity modulated radio therapy (IMRT), VMAT, and 3D-CRT plans. Hybrid, HT, VMAT, IMRT, and 3D-CRT treatment plans for a male child (age 6 years) with medulloblastoma were created and compared. A total dose of 35.2 Gy (PTV) with a dose per fraction of 1.6 Gy was prescribed. The following dose acceptance criteria were defined: The plans were compared regarding dose homogeneity index (HI) and conformity index (CI), PTV coverage, (particularly at cribriform plate) and doses at OARs. Best conformity was achieved with HT (CI = 0.98) followed by VMAT (CI = 0.96), IMRT (CI = 0.91), HybTP (CI = 0.86), and 3D-CRT (CI = 0.83). The homogeneity index varied marginally. For both HT and IMRT the HI was 0.07, and for 3D-CRT, VMAT and HybTP the HI was between 0.13 and 0.15. The cribriform plate was sufficiently covered by HybTP, VMAT, and 3D-CRT. The dose acceptance criteria for OARs were met by HT and HybTP. VMAT did not meet the criteria for lung (Dmean = right 10.4 Gy/left 10.2 Gy), 3D-CRT did not meet the criteria for eye lenses (Dmax = right 32.3 Gy/left 33.1), and heart (V25≈44%) and IMRT did not meet the criteria for lung (Dmean = right 11.1 Gy/left 11.2 Gy) and eye lenses (Dmax = right 12.2 Gy/left 13.1). HybTP meets all defined acceptance criteria and has proved to be a reasonable alternative for CSI. With HybTP that combines VMAT at the brain and heart with 3D-CRT posterior spinal fields (to spare lungs and kidneys), both appropriate coverage of the PTV and sparing of OAR can be achieved.

Past, present, and future of radiotherapy for the benefit of patients

Radiotherapy has been driven by constant technological advances since the discovery of X-rays in 1895. Radiotherapy aims to sculpt the optimal isodose on the tumour volume while sparing normal tissues. The benefits are threefold: patient cure, organ preservation and cost-efficiency. The efficacy and tolerance of radiotherapy were demonstrated by randomized trials in many different types of cancer (including breast, prostate and rectum) with a high level of scientific evidence. Such achievements, of major importance for the quality of life of patients, have been fostered during the past decade by linear accelerators with computer-assisted technology. More recently, these developments were augmented by proton and particle beam radiotherapy, usually combined with surgery and medical treatment in a multidisciplinary and personalized strategy against cancer. This article reviews the timeline of 100 years of radiotherapy with a focus on breakthroughs in the physics of radiotherapy and technology during the past two decades, and the associated clinical benefits.

PET and Radiation Therapy Planning and Delivery for Prostate Cancer

PET imaging has become an integral component of the diagnosis and management of a substantial number of lymphatic and solid malignancies. One of the greatest dilemmas in prostate cancer remains the need for greater personalization of treatment recommendations based on the true extent of disease, so that patients with extraprostatic, micrometastatic disease can be identified early and managed accordingly. These sites currently remain under the level of detection with standard imaging and continue to confound clinicians. Novel PET tracers to complement anatomic data from CT and MR imaging can truly make a difference, and ongoing research holds the greatest promise.

Testing the new ICRU 62 ‘Planning Organ at Risk Volume’ concept for the rectum

BACKGROUND AND PURPOSE

To study the impact of the new ICRU 62 'Planning organ at Risk Volume' (PRV) concept on the relationship between rectum dose-volume histogram (DVH) data and toxicity.

PATIENTS AND METHODS

The acute gastro-intestinal (GI) RTOG toxicity in 127 prostate cancer patients prescribed a total dose of 70 Gy with conformal irradiation to either the prostate, the prostate and seminal vesicles or the whole pelvis (initial 50 Gy only) were analysed. DVHs were derived for the rectum only and for rectum extended with six PRV margin sets (narrow/intermediate/wide; anterior/anterior and posterior). The data was analysed using permutation tests, logistic regression and effective uniform dose (EUD) calculations.

RESULTS

Acute Grade 2 GI toxicity was seen in 22 of 127 cases (17%). Permutation tests showed that the difference between DVHs for patients with and without Grade 2 effects was significant, both for rectum only and rectum PRVs (P-value range: 0.02-0.04), with generally lower P-values for the PRVs. In the logistic regression, the fractional DVH variables (i.e. volumes) were significantly related to toxicity, with approximately 2-3 times as many significant dose levels for the PRVs as for rectum only. E.g. with wide anterior and posterior margins (16 and 11 mm, respectively) the relation was significant at 26 different dose levels (6-7, 13-14, 35-43, 60-71 and 73 Gy), compared to nine levels (38-40, 43-44 and 71-74 Gy) for rectum only. EUDs were significantly different for patients with and without Grade 2 effects both for rectum only and the PRVs (95% confidence interval for EUD increase with Grade 2 effects: 0.1-3.1 Gy).

CONCLUSIONS

All statistical methods applied indicated a small, but definite difference in DVH parameters between patients with versus those without Grade 2 effects. The difference was most pronounced when margins of 16 mm anterior and 11 mm posterior were applied.

An algorithm for systematic selection of beam directions for IMRT

Selection of the number of beams and their directions can be an important problem in radiation therapy, especially when a tumor surrounds a critical organ or is surrounded by multiple critical organs. Beam directions, in this sense, are chosen to not only avoid critical organs, but also to achieve better target dose uniformity. In intensity-modulated radiation therapy (IMRT), optimization of beam directions is further complicated due to the dependence of one beam direction on its corresponding beamlet intensities and the beamlet intensities of all other beam directions. The result is an excessively enlarged search space, even when the number of beams is small (two to three). Until now, only a handful of publications exist regarding beam direction optimization in IMRT. Here, we report a new systematic approach that determines a suitable number of "more optimal" beam directions without optimizing a complicated objective function or resorting to brute force. We start by assuming that beam directions chosen for an N-beam plan are candidates for beam directions in the search for an (N + 1)-beam plan. Knowing that beam directions in an N-beam plan are not always the best choices for the (N + 1)-beam plan, we introduce into the beam direction selection process an analysis of the beamlet weights of every beam direction set sampled. If the relative weights of any particular beam compared to other beams are insignificant and hence have no significant effect on the quality of the treatment plan, then we eliminate this beam from the plan. The algorithm terminates basically when the relative weights of the last beam compared to other beams are insignificant or the replacement of an eliminated beam does not improve the plan. This concept was applied to three two-dimensional phantoms and each plan was compared to a standard equally spaced IMRT plan in terms of dose distributions, dose-volume histograms, and objective function values. The results show improvements in both target dose uniformity and critical organ sparing often with a fewer number of beams than standard equally spaced beam plans.

Effet volume en radiothérapie [II]. Deuxième partie : volume et tolérance des tissus sains

The first part of our work has focused on the relationship between tumor volume and tumor control. Indeed, it is well known that the importance of irradiated volume could be a main parameter of radiation-induced complications. Numerous mathematical models have described the correlation between the irradiated volume and the risk of adverse effects. These models should predict the complication rate of each treatment planning. At the present time late effects have been the most studied. In this report we firstly propose a review of different mathematical models described for volume effect. Secondly, we will discuss whether these theoretical considerations can influence our view of radiation treatment planning optimization.

Estimation of parameters of dose-volume models and their confidence limits

Predictions of the normal-tissue complication probability (NTCP) for the ranking of treatment plans are based on fits of dose-volume models to clinical and/or experimental data. In the literature several different fit methods are used. In this work frequently used methods and techniques to fit NTCP models to dose response data for establishing dose-volume effects, are discussed. The techniques are tested for their usability with dose-volume data and NTCP models. Different methods to estimate the confidence intervals of the model parameters are part of this study. From a critical-volume (CV) model with biologically realistic parameters a primary dataset was generated, serving as the reference for this study and describable by the NTCP model. The CV model was fitted to this dataset. From the resulting parameters and the CV model, 1000 secondary datasets were generated by Monte Carlo simulation. All secondary datasets were fitted to obtain 1000 parameter sets of the CV model. Thus the 'real' spread in fit results due to statistical spreading in the data is obtained and has been compared with estimates of the confidence intervals obtained by different methods applied to the primary dataset. The confidence limits of the parameters of one dataset were estimated using the methods, employing the covariance matrix, the jackknife method and directly from the likelihood landscape. These results were compared with the spread of the parameters, obtained from the secondary parameter sets. For the estimation of confidence intervals on NTCP predictions, three methods were tested. Firstly, propagation of errors using the covariance matrix was used. Secondly, the meaning of the width of a bundle of curves that resulted from parameters that were within the one standard deviation region in the likelihood space was investigated. Thirdly, many parameter sets and their likelihood were used to create a likelihood-weighted probability distribution of the NTCP. It is concluded that for the type of dose response data used here, only a full likelihood analysis will produce reliable results. The often-used approximations, such as the usage of the covariance matrix, produce inconsistent confidence limits on both the parameter sets and the resulting NTCP values.

Elimination of importance factors for clinically accurate selection of beam orientations, beam weights and wedge angles in conformal radiation therapy

A method of simultaneously optimizing beam orientations, beam weights, and wedge angles for conformal radiotherapy is presented. This method removes the need for importance factors by optimizing one objective only, subject to a set of rigid constraints. This facilitates the production of inverse solutions which, without trial-and-error modification of importance factors, precisely satisfy the specified constraints. The algorithm minimizes an objective function which is based upon the single objective to be optimized, but which is forced to an artificially high value when the constraints are not met, so that only satisfactory solutions are allowed. Due to the complex nature of the objective function space, including multiple local minima separated by large regions of plateau, a random search technique equivalent to fast simulated annealing is used for producing inverse plans. To illustrate the novel features of the new algorithm, a simulation is first presented, for the case of a cylindrical phantom. The morphology of the objective function space is shown to be significantly different for the new algorithm, compared to that for a conventional quadratic objective function. Clinical cases for prostate and craniopharyngioma are then presented. For the prostate case, the objective is to reduce irradiated rectal volume. Three-field, four-field, and six-field optimizations, with or without orientation optimization, are shown to provide solutions which are consistent with previously reported plans and class solutions. For the craniopharyngioma case, which involves the use of a high-precision stereotactic conformal technique, the objective is to reduce the irradiated volume of normal brain. Practically feasible beam angles are produced which, compared to a standard plan, provide a small but worthwhile sparing of normal brain. The algorithm is thereby shown to be robust and suitable for clinical application.

Class solutions for conformal external beam prostate radiotherapy

PURPOSE

To determine a class solution coplanar plan from comparisons of three-field (3F), four-field (4F), and six-field (6F) plans in conformal non-intensity-modulated prostate radiotherapy.

METHODS AND MATERIALS

Doses to two clinical target volumes, prostate only (PO) and prostate plus seminal vesicles (PSV) were evaluated in each of 10 patients using a variety of 3F, 4F, and 6F plans with a planning target volume margin of 10 mm. All plans were prescribed to 64 and 74 Gy. The class solution plan for each of 3F, 4F, and 6F was chosen from a variety of symmetrical and asymmetrical field arrangements that had been previously assessed. The class solution plans, 3F (0, 90, 270 degrees ), 4F (35, 90, 270, 325 degrees ), and 6F (50/lat/25) were compared with reference plans: 3F (0, 120, 240 degrees ), 4F (0, 90, 180, 270 degrees ), and 6F (55, 90, 125, 235, 270, 305 degrees ). Rectal volumes irradiated to greater than 50% (V(50)), 80% (V(80)), and 90% (V(90)) of the prescribed dose, normal tissue complication probabilities (NTCP) for rectum, bladder, and femoral heads (FH), and tumor control probabilities (TCP) were assessed. FH tolerance was set at 52 Gy to 10% volume.

RESULTS

The field arrangement that gave the lowest irradiated rectal volume with acceptable bladder and FH doses was a 3F (0, 90, 270 degrees ) class solution plan. This plan gave a reduction in rectal V(80) of 1.2-12.4% for the PO group and 2.3-23.8% for the PSV group compared with the other plans. The reduction in rectal V(90) was 0.2-11.9% for the PO group and 1.5-23.3% for the PSV group using the 3F (0, 90, 270 degrees ) plan. This plan provided one of the lowest rectal NTCPs, but the difference was not significant when compared with the 4F class solution plan. When target volumes with 10-mm margins remain unchanged to 74 Gy, the irradiated rectal volumes for all plans were higher and rectal NTCPs can be trebled.

CONCLUSION

The use of appropriate beam arrangements can provide a class solution plan using only 3 fields compared with 4 or 6 fields for the parameters considered. Both 3F (0, 90, 270 degrees ) and 4F (35, 90, 270, 325 degrees ) plans can be used as a class solution plan. Other practical issues that may influence the choice of class solution include delivery time with smaller number of fields, ease of verification, the use of 10-mm multileaf collimation vs. conformal blocks, and field shape fitting limitations when using dynamic wedges.

Radiobiological indices that consider volume: a review

Understanding and predicting the impact of any radiotherapy treatment is critical if patients are to receive treatment with a high likelihood of eliminating the tumour and low likelihood of complications. One of the major contributing factors in determining these effects is the volume treated. This review assesses the current use and accuracy of a series of models which consider volume, building on a previous review which investigated the impact of fractionation particularly with respect to the linear quadratic model. Volume is particularly important in assessing the overall effect with respect to destroying the clonogenic cells and preventing damage to the normal tissues. Dose volume histograms are one of the simplest and most useful forms of representing volume information, however it is difficult to correlate plans based only on DVHs. For this reason various reduction schemes have been introduced and tumour control probability and normal tissues complication probability models adjusted to use this information. Many of these models have proved quite useful in the clinic although they are limited by the available radiobiological data.

Comparison of two algorithms for determining beam weights and wedge filters

This article compares two algorithms for determining beam weights and wedge filters for conformal treatment planning. One algorithm, which is based on dose-gradient analysis, provides analytic formulas for determining beam weights, wedge angles, and collimator angles (i.e., wedge orientations) so that the dose distribution is homogeneous in the target volume. The second algorithm is based on the concept of the super-omni wedge (i.e., the arrangement of two pairs of orthogonal nominal wedged beams), numerically optimize beam weights, wedge angles, and collimator angles so that the dose requirements to targets and organs at risk are satisfied to the best. Three clinical cases were tested. For the first case, both algorithms resulted in comparable homogeneous dose distributions in the target volume. For the second case, the second algorithm resulted in much lower doses to the eyes plus a better homogeneous dose distribution in the target volume. For the third case, only the second algorithm was applicable, and the treatment plan it developed met the prescribed requirements. The results show that the first algorithm is better in terms of feasibility, whereas the second is better in terms of applicability and the quality of treatment plans.

Partially wedged beams improve radiotherapy treatment of urinary bladder cancer

BACKGROUND AND PURPOSE

Partially wedged beams (PWBs) having wedge in one part of the field only, can be shaped using dynamic jaw intensity modulation. The possible clinical benefit of PWBs was tested in treatment plans for muscle-infiltrating bladder cancer.

MATERIAL AND METHODS

Three-dimensional treatment plans for 25 bladder cancer patients were analyzed. The originally prescribed standard conformal four-field box technique, which includes the use of lateral ordinary wedge beams, was compared to a modified conformal treatment using customized lateral PWBs. In these modified treatment plans, only the anterior parts of the two lateral beams had a wedge. To analyze the potential clinical benefit of treatment with PWBs, treatment plans were scored and compared using both physical parameters and biological dose response models. One tumour control probability model and two normal tissue complication probability (NTCP) models were applied. Different parameters for normal tissue radiation tolerance presented in the literature were used.

RESULTS

By PWBs the dose homogeneity throughout the target volume was improved for all patients, reducing the average relative standard deviation of the target dose distribution from 2.3 to 1.8%. A consistent reduction in the maximum doses to surrounding normal tissue volumes was also found. The most notable improvement was demonstrated in the rectum where the volume receiving more than the prescribed tumour dose was halved. Treatment with PWBs would permit a target dose escalation of 2-6 Gy in several of the patients analyzed, without increasing the overall risk for complications. The number of patients suitable for dose escalation ranged from 3 to 15, depending on whether support from all or only one of the five applied NTCP model/parameter combinations were required in each case to recommend dose escalation.

CONCLUSION

PWBs represent a simple dose conformation tool that may allow radiation dose escalation in the treatment of muscle-infiltrating urinary bladder tumours.

Planning, computer optimization, and dosimetric verification of a segmented irradiation technique for prostate cancer

PURPOSE

To develop and verify a multisegment technique for prostate irradiation that results in better sparing of the rectal wall compared to a conventional three-field technique, for patients with a concave-shaped planning target volume (PTV) overlapping the rectal wall.

METHODS AND MATERIALS

Five patients have been selected with various degrees of overlap between PTV and rectal wall. The planned dose to the ICRU reference point is 78 Gy. The new technique consists of five beams, each having an open segment covering the entire PTV and several smaller segments in which the rectum is shielded. Segment weights are computer-optimized using an algorithm based on simulated annealing. The score function to be minimized consists of dose-volume constraints for PTV, rectal wall, and femoral heads. The resulting dose distribution is verified for each patient by using point measurements and line scans made with an ionization chamber in a water tank and by using film in a cylindrical polystyrene phantom.

RESULTS

The final number of segments in the five-field technique ranges from 7 to 9 after optimization. Compared to the standard three-field technique, the maximum dose to the rectal wall decreases by approximately 3 Gy for patients with a large overlap and 1 Gy for patients with no overlap, resulting in a reduction of the normal tissue complication probability (NTCP) by a factor of 1.3 and 1.2, respectively. The mean dose to the PTV is the same for the two techniques, but the dose distribution is slightly less homogeneous with the five-field technique (Average standard deviation of five patients is 1.1 Gy and 1.7 Gy for the three-field and five-field technique, respectively). Ionization chamber measurements show that in the PTV, the calculated dose is in general within 1% of the measured dose. Outside the PTV, systematic dose deviations of up to 3% exist. Film measurements show that for the complete treatment, the position of the isodose lines in sagittal and coronal planes is calculated fairly accurately, the maximum distance between measured and calculated isodoses being 4 mm.

CONCLUSIONS

We developed a relatively simple multisegment "step-and-shoot" technique that can be delivered within an acceptable time frame at the treatment machine (Extra time needed is approximately 3 minutes). The technique results in better sparing of the rectal wall compared to the conventional three-field technique. The technique can be planned and optimized relatively easily using automated procedures and a predefined score function. Dose calculation is accurate and can be verified for each patient individually.

Evaluation of the optimal co-planar field arrangement for use in the boost phase of dose escalated conformal radiotherapy for localized prostate cancer

The aim of this study was to determine the optimal co-planar beam arrangement from a variety of three-field (3F), four-field (4F) and six-field (6F) plans for the boost phase of a dose escalated conformal radiotherapy schedule. Three selected plans (3F 0 degrees, 90 degrees, 270 degrees plan, 4F 45 degrees, 90 degrees, 270 degrees, 315 degrees plan and 6F 40 degrees, 90 degrees, 115 degrees, 245 degrees, 270 degrees, 320 degrees plan) were compared with reference plans (3F 0 degrees, 120 degrees, 240 degrees plan, 4F 0 degrees, 90 degrees, 180 degrees, 270 degrees plan, 6F 55 degrees, 90 degrees, 125 degrees, 235 degrees, 270 degrees, 305 degrees plan and 6F 50 degrees, 90 degrees, 130 degrees, 230 degrees, 270 degrees, 310 degrees plan) in 10 patients. Doses of 64 Gy and 74 Gy were prescribed to the isocentre using 6 MV photons. The boost planning target volume comprised the prostate gland alone without a margin. Plans were compared by means of rectal volumes irradiated to >50% (V50), >80% (V80) and >90% (V90) of the prescribed dose. Irradiated volumes were also measured for the bladder (V90) and the femoral heads (V70). All optimal 3F, 4F and 6F plans gave lower irradiated rectal V80 and V90 levels than their corresponding reference plan. The 3F (0 degrees, 90 degrees, 270 degrees) plan consistently provided lower irradiated rectal levels at V50 to V90, with acceptable bladder and femoral head doses compared with the other plans in the study. When the 6F (50 degrees, 90 degrees, 130 degrees, 230 degrees, 270 degrees, 310 degrees) plan used at our institution for the boost phase was compared with the 3F (0 degrees, 90 degrees, 270 degrees) plan, the rectal V50 was reduced from 20.8+/-5.2%, to 12.6+/-5.1%, the rectal V80 was reduced from 8.7+/-2.9% to 6.5+/-3.1% and the rectal V90 was reduced from 5.5+/-2.1% to 3.9+/-2.0% (all p<0.001). The bladder V90 and the femoral heads V70 levels were equivalent. For the boost phase when escalating the dose from 64 Gy to 74 Gy, the co-planar plan that allowed optimal rectal sparing was a 3F beam arrangement using gantry angles of 0 degrees, 90 degrees and 270 degrees. This 3F plan provided improved rectal sparing compared with the 6F (50 degrees, 90 degrees, 130 degrees, 230 degrees, 270 degrees, 310 degrees) beam arrangement currently used at our institution, with equivalent and acceptable bladder and femoral head doses.

Radiotherapy of prostate cancer with or without intensity modulated beams: a planning comparison

PURPOSE

To evaluate whether intensity modulated radiotherapy (IMRT) by static segmented beams allows the dose to the main portion of the prostate target to escalate while keeping the maximal dose at the anterior rectal wall at 72 Gy. The value of such IMRT plans was analyzed by comparison with non-IMRT plans using the same beam incidences.

METHODS AND MATERIALS

We performed a planning study on the CT data of 32 consecutive patients with localized adenocarcinoma of the prostate. Three fields in the transverse plane with gantry angles of 0 degrees, 116 degrees, and 244 degrees were isocentered at the center of gravity of the target volume (prostate and seminal vesicles). The geometry of the beams was determined by beam's eye view autocontouring of the target volume with a margin of 1.5 cm. In study 1, the beam weights were determined by a human planner (3D-man) or by computer optimization using a biological objective function with (3D-optim-lim) or without (3D-optim-unlim) a physical term to limit target dose inhomogeneity. In study 2, the 3 beam incidences mentioned above were used and in-field uniform segments were added to allow IMRT. Plans with (IMRT-lim) or without (IMRT-unlim) constraints on target dose inhomogeneity were compared. In the IMRT-lim plan, target dose inhomogeneity was constrained between 15% and 20%. After optimization, plans in both studies were normalized to a maximal rectal dose of 72 Gy. Biological (tumor control probability [TCP], normal tissue complication probability [NTCP]) and physical indices for tumor control and normal tissue complication probabilities were computed, as well as the probability of the uncomplicated local control (P+).

RESULTS

The IMRT-lim plan was superior to all other plans concerning TCP (p < 0.0001). The IMRT-unlim plan had the worst TCP. Within the 3D plans, the 3D-optim-unlim had the best TCP, which was significantly different from the 3D-optim-lim plan (p = 0.0003). For rectal NTCP, both IMRT plans were superior to all other plans (p < 0.0001). The IMRT-unlim plan was significantly better than the IMRT-lim plan (p < 0.0001). Again, 3D-optim-unlim was superior to the other 3D plans (p < 0. 0007). Physical endpoints for target showed the mean minimal target dose to be the lowest in the IMRT-unlim plan, caused by a large target dose inhomogeneity (TDI). Medial target dose, 90th percentile, and maximal target dose were significantly higher in both IMRT plans. Physical endpoints for the rectum showed the IMRT-unlim plan to be superior compared to all other plans. There was a strong correlation between the 65th percentile (Rp65) and rectal NTCP (correlation coefficient > or =89%). For bladder, maximal bladder dose was significantly higher in the IMRT-unlim plan compared to all other plans (p < or = 0.0001).P+ was significantly higher in both IMRT-plans than in all other plans. The 3D-optim-unlim plan was significantly better than the two other 3D plans (p < 0.0001).

CONCLUSION

IMRT significantly increases the ratio of TCP over NTCP of the rectum in the treatment of prostate cancer. However, constraints for TDI are needed, because a high degree of TDI reduced minimal target dose. IMRT improved uncomplicated local control probability. In our department, IMRT by static segmented beams is planned and delivered in a cost-effective way. IMRT-lim has replaced non-modulated conformal radiotherapy as the standard treatment for prostate cancer.

An evaluation of three-field coplanar plans for conformal radiotherapy of prostate cancer

BACKGROUND AND PURPOSE

A series of coplanar three-field configurations for two different clinical treatment volumes, prostate only (PO) and prostate plus seminal vesicles (PSV) were studied to determine the optimal three-field plan arrangement for prostate radiotherapy.

MATERIALS AND METHODS

A variety of conformal three-field 6 MV plans prescribed to both 64 and 74 Gy were created for PO and PSV volumes in each of ten patients. For description, the orientation of each sequential beam was named in a clockwise fashion. Plans included series with arrangements of 0 degrees, 60-150 degrees, 210-300 degrees; 0 degrees, 90 degrees, 225-255 degrees; 90 degrees, 210-240 degrees, 300-330 degrees and a four-field (4F) box plan for comparison. Six-hundred and eighty plans were compared using the rectal volume irradiated to greater than 50% (V(50)), 80% (V(80)), and 90% (V(90)) of the prescribed dose, normal tissue complications (NTCP) for rectum, bladder, and femoral heads (FH), and tumour control probabilities (TCP). FH tolerance was set at 52 Gy to 10% volume.

RESULTS

In comparing the 34 different three-field configurations for each of the PO and PSV groups, the greatest rectal sparing was achieved by a three-field plan with gantry angles of 0 degrees, 90 degrees, 270 degrees (PO: rectal V(80)=22.8+/-5.5% (1S.D.), V(90)=18.4+/-5.7%, and PSV: rectal V(80)=41.9+/-5.8%, V(90)=35.5+/-5.9%). This also improved on the 4F-box plan (PO: rectal V(80)=26.0+/-5.8%, V(90)=21.4+/-5.2%, P<0.001; and PSV: rectal V(80)=47.3+/-5.5%, V(90)=41.6+/-5.1%, P<0.001). The worst rectal sparing was seen with the 0 degrees, 120 degrees, 240 degrees plan (PO: rectal V(80)=35.2+/-8.0%, V(90)=30.3+/-7.1%, P<0.001; and PSV: rectal V(80)=65.7+/-9.0%, V(90)=58.8+/-8.8%, P<0.001). In the PO group, the increase in predicted rectal NTCP with dose escalation from 64 to 74 Gy was 3.3% using the 0 degrees, 90 degrees, 270 degrees plan, 4.7% with the 4F-box plan, and 6.9% with the 0 degrees, 120 degrees, 240 degrees plan. In the PSV group, dose escalation increased the predicted rectal NTCP by 7.9, 10.1 and 15.7% for the 0 degrees, 90 degrees, 270 degrees plan, 4F-box plan, and 0 degrees, 120 degrees, 240 degrees plan, respectively.

CONCLUSIONS

For both PO and PSV volumes, the three-field plan which afforded the greatest rectal sparing with acceptable bladder and femoral head doses was the 0 degrees, 90 degrees, 270 degrees plan. This plan also improved on the 4F-box. The increase in predicted rectal NTCP when escalating dose from 64 to 74 Gy was smaller using this plan compared to either the three-field 0 degrees, 120 degrees, 240 degrees plan or the 4F-box plan.

Optimization of inverse treatment planning using a fuzzy weight function

A fuzzy approach has been applied to inverse treatment planning optimization in radiation therapy. The proposed inverse-planning algorithm optimizes both the intensity-modulated beam (IMB) and the normal tissue prescription. In the IMB optimization, we developed a fast-monotonic-descent (FMD) method that has the property of fast and monotonic convergence to the minimum for a constrained quadratic objective function. In addition, a fuzzy weight function is employed to express the vague knowledge about the importance of matching the calculated dose to the prescribed dose in the normal tissue. Then, a validity function is established to optimize the normal tissue prescription. The performance of this new fuzzy prescription algorithm has been compared to that based on hard prescription methods for two treatment geometries. The FMD method presented here both provides a full-analytical solution to the optimization of intensity-modulated beams, and guarantees fast and monotonic convergence to the minimum. It has been shown that the fuzzy inverse planning technique is capable of achieving an optimal balance between the objective of matching the calculated dose to the prescribed dose for the target volume and the objective of minimizing the normal tissue dose.

Optimization of coplanar six-field techniques for conformal radiotherapy of the prostate

PURPOSE

To determine the optimal coplanar treatment technique for six-field conformal radiotherapy of prostate only (PO) or prostate plus seminal vesicles (PSV).

METHODS AND MATERIALS

A series of 6-MV six-field coplanar treatment plans were created for PO and PSV volumes in 10 patients prescribed to both 64 and 74 Gy. All plans consisted of laterally-symmetric anterior oblique, lateral, and posterior oblique fields. The posterior oblique fields were varied through 20-45 degrees relative to the lateral fields, and for each of these angles, the anterior oblique fields were varied through 25-65 degrees relative to lateral. The plans were compared by means of rectal volumes irradiated to 80% or more of the prescribed dose (V80); normal tissue complication probability (NTCP) for rectum, bladder, and femoral heads; and tumor control probability (TCP). Femoral head tolerance was designated as 52 Gy to no more than 10% volume.

RESULTS

For the PO group, anterior oblique fields at 50 degrees from lateral and posterior oblique fields at 25 degrees from lateral produced the lowest V80, together with femoral head doses which were appropriate for most patients (V80 = 24.4+/-5.3% [1 SD]). Compared to a commonly-used six-field (reference) plan with both anterior and posterior oblique fields at 35 degrees from lateral (V80 = 26.3+/-5.9%), this represented an improvement (p = 0.001). For the PSV group, the optimal anterior and posterior oblique fields were at 65 degrees and 30 degrees from lateral, respectively (V80 = 47.5+/-6.3%). Relative to the reference plan (V80 = 49.4+/-5.6%), this was a marginal improvement (p = 0.07).

CONCLUSION

The optimized six-field plans provide increased rectal sparing at both standard and escalated doses. Moreover, the gain in TCP resulting from dose escalation can be achieved with a smaller increase in rectal NTCP using the optimized six-field plans.

Monte Carlo-based inverse treatment planning

A Monte Carlo based inverse treatment planning system (MCI) has been developed which combines arguably the most accurate dose calculation method (Monte Carlo particle transport) with a 'guaranteed' optimization method (simulated annealing). A distribution of photons is specified in the tumour volume; they are transported using an adjoint calculation method to outside the patient surface to build up an intensity distribution. This intensity distribution is used as the initial input into an optimization algorithm. The dose distribution from each beam element from a number of fields is pre-calculated using Monte Carlo transport. Simulated annealing optimization is then used to find the weighting of each beam element, to yield the optimal dose distribution for the given criteria and constraints. MCI plans have been generated in various theoretical phantoms and patient geometries. These plans show conformation of the dose to the target volume and avoidance of critical structures. To verify the code, an experiment was performed on an anthropomorphic phantom.

A comparison of coplanar four-field techniques for conformal radiotherapy of the prostate

BACKGROUND AND PURPOSE

Conformal radiotherapy of the prostate is an increasingly common technique, but the optimal choice of beam configuration remains unclear. This study systematically compares a number of coplanar treatment plans for four-field irradiation of two different clinical treatment volumes: prostate only (PO) and the prostate plus seminal vesicles (PSV).

MATERIALS AND METHODS

A variety of four-field coplanar treatment plans were created for PO and PSV volumes in each of ten patients. Plans included a four-field 'box' plan, a symmetric plan having bilateral anterior and posterior oblique fields, a plan with anterior oblique and lateral fields, a series of asymmetric plans, and a three-field plan having anterior and bilateral fields for comparison. Doses of 64 and 74 Gy were prescribed to the isocentre. Plans were compared using the volume of rectum irradiated to greater than 50% (V50), 80% (V80) and 90% (V90) of the prescribed dose. Tumour control probabilities (TCP) and normal tissue complication probabilities (NTCP) for the rectum, bladder and femoral heads were also evaluated. Femoral head dose was limited such that less than 10% of each femoral head received 70% of the prescribed dose.

RESULTS

For the PO group, the optimal plan consisted of anterior oblique and lateral fields (Rectal V80 = 23.8+/-5.0% (1 SD)), while the box technique (V80 = 26.0+/-5.8%) was less advantageous in terms of rectal sparing (P = 0.001). Similar results were obtained for the PSV group (Rectal V80 = 43.9+/-5.0% and 47.3+/-5.5% for the two plan types, respectively, P = 0.001). The three-field plan was comparable to the optimal four-field plan but gave higher superficial body dose. With dose escalation from 64 to 74 Gy, the mean TCP for the optimal plan rose from 52.0+/-2.8% to 74.1+/-2.0%. Meanwhile, rectal NTCP for the optimal plan rose by 3.5% (PO) or 8.4% (PSV), compared to 4.7% (PO) or 10.1% (PSV) for the box plan.

CONCLUSIONS

For PO volumes, a plan with gantry angles of 35 degrees, 90 degrees, 270 degrees and 325 degrees offers a high level of rectal sparing and acceptable dose to the femoral heads for all patients, while for PSV volumes, the corresponding plan has gantry angles of 20 degrees, 90 degrees , 270 degrees and 340 degrees. Using these plans, the gain in TCP resulting from dose escalation can be achieved with a smaller increase in anticipated rectal NTCP.

Constrained customization of non-coplanar beam orientations in radiotherapy of brain tumours

A methodology for the constrained customization of non-coplanar beam orientations in radiotherapy treatment planning has been developed and tested on a cohort of five patients with tumours of the brain. The methodology employed a combination of single and multibeam cost functions to produce customized beam orientations. The single-beam cost function was used to reduce the search space for the multibeam cost function, which was minimized using a fast simulated annealing algorithm. The scheme aims to produce well-spaced, customized beam orientations for each patient that produce low dose to organs at risk (OARs). The customized plans were compared with standard plans containing the number and orientation of beams chosen by a human planner. The beam orientation constraint-customized plans employed the same number of treatment beams as the standard plan but with beam orientations chosen by the constrained-customization scheme. Improvements from beam orientation constraint-customization were studied in isolation by customizing the beam weights of both plans using a dose-based downhill simplex algorithm. The results show that beam orientation constraint-customization reduced the maximum dose to the orbits by an average of 18.8 (+/-3.8, ISD)% and to the optic nerves by 11.4 (+/-4.8, ISD)% with no degradation of the planning target volume (PTV) dose distribution. The mean doses, averaged over the patient cohort, were reduced by 4.2 (+/-1.1, ISD)% and 12.4 (+/-3.1, ISD)% for the orbits and optic nerves respectively. In conclusion, the beam orientation constraint-customization can reduce the dose to OARs, for few-beam treatment plans, when compared with standard treatment plans developed by a human planner.

Intensity modulation methods for proton radiotherapy

The characteristic Bragg peak of protons or heavy ions provides a good localization of dose in three dimensions. Through their ability to deliver laterally and distally shaped homogenous fields, protons have been shown to be a precise and practical method for delivering highly conformal radiotherapy. However, in an analogous manner to intensity modulation for photons, protons can be used to construct dose distributions through the application of many individually inhomogeneous fields, but with the localization of dose in the Bragg peak providing the possibility of modulating intensity within each field in two or three dimensions. We describe four different methods of intensity modulation for protons and describe how these have been implemented in an existing proton planning system. As a preliminary evaluation of the efficacy of these methods, each has been applied to an example case using a variety of field combinations. Dose-volume histogram analysis of the resulting dose distributions shows that when large numbers of fields are used, all techniques exhibit both good target homogeneity and sparing of neighbouring critical structures, with little difference between the four techniques being discerned. As the number of fields is decreased, however, only a full 3D modulation of individual Bragg peaks can preserve both target coverage and sparing of normal tissues. We conclude that the 3D method provides the greatest flexibility for constructing conformal doses in challenging situations, but that when large numbers of beam ports are available, little advantage may be gained from the additional modulation of intensity in depth.

A case study comparing the relative benefit of optimizing beam weights, wedge angles, beam orientations and tomotherapy in stereotactic radiotherapy of the brain

A treatment-planning case study has been performed on a patient with a medium-sized, convex brain tumour. The study involved the application of advanced treatment-plan optimization techniques to improve on the dose distribution of the 'standard plan' used to treat the patient. The standard plan was created according to conventional protocol at the Royal Marsden NHS Trust, and consisted of a three-field (one open and two wedged) non-coplanar arrangement, with field shaping to the beam's-eye view of the planning target volume (PTV). Three optimized treatment plans were created corresponding to (i) the optimization of the beam weights and wedge angles of the standard plan, (ii) the optimization of the beam orientations, beam weights and wedge angles of the standard plan, and (iii) a full fluence tomotherapy optimization of 1 cm wide (at isocentre), 270 degree arcs. (i) and (ii) were created on the VOXELPLAN research 3D treatment-planning system, using in-house developed optimization algorithms, and (iii) was created on the PEACOCK tomotherapy planning system. The downhill-simplex optimization algorithm is used, in conjunction with 'threshold-dose' cost-function terms enabling the algorithm to optimize specific regions of the dose-volume histogram (DVH) curve. The 'beam-cost plot' tool is presented as a visual aid to the selection of beneficial beam directions. The methods and pitfalls in the transfer of plans and patient data between the two planning systems are discussed. Each optimization approach was evaluated, relative to the standard plan, on the basis of DVH and dose statistics in the PTV and organs at risk (OARs). All three optimization approaches were able to improve on the dose distribution of the standard plan. The magnitude of the improvement was greater for the optimized beam-orientation and tomotherapy plans (up to 15% and 30% for the maximum and mean OAR doses). A smaller improvement was observed in the beam-weight and wedge-angle optimized plan (up to 5% and 10% in the maximum and mean OAR doses). In the tomotherapy plan, difficulty was encountered achieving an acceptable homogeneity of dose in the PTV. This was improved by treating the gross tumour volume (GTV) and (PTV - GTV) regions as separate targets in the inverse planning, with the latter region prescribed a slightly higher dose to reduce edge under-dosing. In conclusion, for the medium-sized convex tumour studied, the tomotherapy dose distribution showed a significant improvement on the standard plan, but no significant improvement over a conventional three-field plan where the beam orientations, beam weights and wedge angles had been optimized.

Improvements in prostate radiotherapy from the customization of beam directions

A methodology for optimizing the beam directions in radiotherapy treatment planning has been developed and tested on a cohort of twelve prostate patients. An optimization algorithm employing a an objective cost function was used, based on beam's-eye-view volumetrics but also employing a simple dose model and biological considerations for organs-at-risk (OARs). The cost function embodies information about the volume of OARs in a single field and their position relative to the planning target volume (PTV). The proximity of the PTV to the surface of the patient is also included. Within the algorithm "importance factor" were used to model the clinical importance of different organs-at-risk so that all organs-at-risk were included in a single objective score. "Gantry-angle-windows" were introduced to restrict the available beam directions. The methodology was applied to twelve prostate patients to determine the optimum beam directions for three-field direction plans. Orientation-optimized and standard treatment plans were compared via measures of tumor control probability (TCP) and normal tissue complication probability (NTCP). Standard plans had fixed beam directions whereas orientation-optimized plans contained beam directions chosen by the algorithm. The beam-weights of both the orientation-optimized and standard plans were optimized using a dose-based simulated annealing algorithm to allow the improvements by optimizing the beam directions to be studied in isolation. The results of the comparison show that optimization of the beam directions yielded better plans, in terms of TCP and NTCP, than the standard plans. When the dose to the isocenter was scaled to produce a rectal NTCP of 1%, the average TCP of the orientation-optimized plans was (5.7 +/- 1.4)% greater than that for the standard plans. In conclusion, the customization of beam directions in the treatment planning of prostate patients using and objective cost function and allowed gantry-angle-windows produces superior three-field direction plans compared to standard treatment plans.

3D conformal intensity-modulated radiotherapy planning: interactive optimization by constrained matrix inversion

BACKGROUND AND PURPOSE

This paper presents a method for interactive optimization of 3D conformal intensity-modulated radiotherapy plans employing a quadratic objective that also contains dose limitations in the organs at risk. This objective function is minimized by constrained matrix inversion (CMI) that follows the same approach as the gradient technique using matrix notation.

MATERIALS AND METHODS

Sherouse's GRATIS radiotherapy design system is used to determine the outlines of the target volume and the organs at risk and to input beam segments which are given by the beam segmentation technique. This technique defines the beam incidences and the beam segmentation. The weights of the segments are then calculated using a quadratic objective function and CMI. The objective function to be minimized consists of two components based on the planning target volume (PTV) and the organ at risk (OAR) with an importance factor w associated with the OAR.

RESULTS

Optimization is tested for concave targets in the head and neck region wrapping around the spinal cord. For a predefined w-value, segment weights are optimized within a few seconds on a DEC Alpha 3000. In practice, 5-10 w-values have to be tested, making optimization a less than 5 min procedure. This optimization procedure predicts the possibility of target dose escalation for a tumour in the lower neck to 120-150 Gy without exceeding the spinal cord tolerance, whereas human planners could not increase the dose above 65-80 Gy.

CONCLUSIONS

Treatment plans optimized using a quadratic objective function and the CMI algorithm are superior to those which are generated by human planners. The optimization algorithm is very fast and allows interactive use. Quadratic optimization by CMI is routinely used by clinicians at the Division of Radiotherapy, U.Z.-Gent.

A gradient inverse planning algorithm with dose-volume constraints

An inverse planning algorithm for determining the intensity-modulated beams that will most closely generate a desired dose distribution is presented. The algorithm is three-dimensional and does not explicitly depend on beam energies and modalities. It allows a single prescription dose or a window of acceptable doses to be specified for the target, with additional constraints to account for under- or over-dosing. For the protection of organs at risk, it provides maximum-dose and dose-volume constraints. The latter apply to the entire volume of the organ exposed to the corresponding dose levels. Several levels of each type of constraint, with varying penalty weights, may be specified for each organ. The objective function that serves as the measure of the goodness of the solution is of the least-squares type and is minimized using conjugate gradient methods. Typical clinical cases involving 40,000 points and 4000 rays to be determined require about 10 min of CPU time on a DEC AlphaStation. Results are presented for two clinical sites, prostate and lung. The optimization algorithm yielded plans that featured higher target dose homogeneity, compared with the human planner's plan, while selectively sparing more of the normal organs at the desired dose regions.

Multiobjective decision theory for computational optimization in radiation therapy

Machine-guided iterative optimization in radiation oncology requires ordinal or cardinal ranking of competing treatment plans. When the clinical objectives are multifaceted and incommensurable, the ranking formalism must take into account the decision maker's tradeoff strategies in a multidimensional decision space. To capture the decision processes in treatment planning, a multiobjective decision-theoretic scheme is formulated. Ranking among a group of candidate plans is based on a generalized distance metric. A dynamic metric weighting function is defined based on the state energy of the decision system, which is assumed to undergo thermodynamic cooling with iteration time. The decision maker is required to specify a baseline ranking of the objectives, which is taken to be the ground state of the decision system. This decision-theoretic formalism was applied to idealized cases in stereotactic radiosurgery and prostatic implantation, using the genetic algorithm as the optimization engine. The optimization pathways and the outcome at limited horizons indicated that the combined scheme of decision-theoretic steering and iterative optimization was robust and produced treatment plans consistent with the user's expectation. The effect of treatment uncertainties was simulated using imperfect objectives; however, certain recurring plans could be identified as optimized baseline solutions. Overall, the present formalism provides a realistic alternative to complete utility assessment or human-guided exploration of the efficient solution set.

Methods for transferring patient and plan data between radiotherapy treatment planning systems

The effectiveness of conformal radiotherapy can ultimately only be assessed by the use of clinical trials. As large multicentre clinical trials become more widespread, methods of transferring patient and plan data between radiotherapy treatment planning systems become increasingly important. In this paper, the general strategy for the transfer of data is discussed, and also illustrated with reference to two specific systems: TARGET 2 (GE Medical Systems) and VOXELPLAN (DKFZ-Heidelberg). The transfer method involves using a computer program to translate the data formats used by each of the two systems for CT scans, patient outlines, plan information and block descriptions. This paper does not address the question of transferring beam data between systems: beam data must first be entered separately into both machines. The physical concepts encountered when transferring plans are described, with specific reference to the two planning systems used. Differences in the strategies used by the two planning systems for definition of irregular field shapes are compared. The dose calculations used by the two systems are also briefly evaluated. Isodoses produced by VOXELPLAN around a circular target volume are found to be up to 3 mm different in location to those produced by TARGET 2, owing to the use of a smooth field shape contour as opposed to a stepped field shape which closely models the leaves of a multileaf collimator. In general, dose distributions generated by both systems are comparable, but some differences are found in the presence of large tissue inhomogeneities. It is concluded that the transfer of patient and plan data between two different treatment planning systems is feasible, provided that any differences in field shape definition methods or dose calculation methods between the two systems are understood.

Selection of coplanar or noncoplanar beams using three-dimensional optimization based on maximum beam separation and minimized nontarget irradiation

PURPOSE

The design of an appropriate set of multiple fixed fields to achieve a steep dose gradient at the tumor edge, with minimal normal tissue exposure, is a very difficult problem, since a virtually infinite number of possible beam orientations exists. In practice we have selected beams in an iterative and often time-consuming process. This work proposes an optimization method, based on geometric and dose elements, to effectively arrive at a set of beam orientations.

METHODS AND MATERIALS

Beams are selected by minimizing a goal function including an angle function (beam separation for steep dose gradient at target edge) and a length function (related to normal tissue dose volume histogram). The relative importance of these two factors may be adjusted depending on the clinic situation. The model is flexible and can include case specific practical anatomic and physical considerations.

RESULTS

In extremely simple situations, the goal function yields results consistent with well-known analytical solutions. When applied to more complex clinical situations, it provides clinically reasonable solutions similar to those empirically developed by the clinician. The optimization process takes approximately 25 min on a UNIX workstation.

CONCLUSION

The optimization scheme provides a practical means for rapidly designing multiple field coplanar or noncoplanar treatments. It overcomes limitations in human three-dimensional visualization such as trying to visualize beam directions and keeping track of the hinge angle between beams while accounting for anatomic/machine constraints. In practice, it has been used as a starting point for physicians to make modifications, based on their clinical judgment.

Predicting late rectal complications following prostate conformal radiotherapy using biologically effective doses and normalized dose-surface histograms

A model to predict the late normal tissue complication probability (NTCP) of the rectum following conformal therapy is described. The model evaluates the biological consequence of inhomogeneities in the physical dose by computing dose histograms of the biologically effective dose to the surface of the rectum for a given fractionation scheme. A method of normalizing the surface area of the rectum is employed so that the predicted NTCP is independent of the differing cross-sectional size of sections of the rectum, ensuring the NTCP is dependent only on the dose delivered to sensitive rectal tissues. The model has been used to assess severe late rectal complications and the milder RTOG grades 2 and 3 reactions. This model was found to predict severe toxicity levels of 1.7 +/- 0.6% for an accelerated treatment of 50 Gy in 16 fractions commonly employed at this centre. This result lies between the severe toxicities predicted for 60 and 62 Gy delivered in 2 Gy fractions. The model predicts that the average NTCP for severe late effects for nine prostate patients becomes greater than 5% with a fractionation scheme of 70 Gy in 35 fractions, for the four fields treatment. The effects of not treating all fields at each therapy session on rectal toxicity were also investigated. Biologically effective dose-surface histograms show that the dose to the lower surface of the rectum is increased by not treating all fields at each therapy session, but the predicted differences in rectal NTCP are negligible.

The optimisation of wedge filters in radiotherapy of the prostate

A treatment plan optimisation algorithm has been applied to 12 patients with early prostate cancer in order to determine the optimum beam-weights and wedge angles for a standard conformal three-field treatment technique. The optimisation algorithm was based on fast-stimulated-annealing using a cost function designed to achieve a uniform dose in the planning-target-volume (PTV) and to minimise the integral doses to the organs-at-risk. The algorithm has been applied to standard conformal three-field plans created by an experienced human planner, and run in three PLAN MODES: (1) where the wedge angles were fixed by the human planner and only the beam-weights were optimised; (2) where both the wedge angles and beam-weights were optimised; and (3) where both the wedge angles and beam-weights were optimised and a non-uniform dose was prescribed to the PTV. In the latter PLAN MODE, a uniform 100% dose was prescribed to all of the PTV except for that region that overlaps with the rectum where a lower (e.g., 90%) dose was prescribed. The resulting optimised plans have been compared with those of the human planner who found beam-weights by conventional forward planning techniques. Plans were compared on the basis of dose statistics, normal-tissue-complication-probability (NTCP) and tumour-control-probability (TCP). The results of the comparison showed that all three PLAN MODES produced plans with slightly higher TCP for the same rectal NTCP, than the human planner. The best results were observed for PLAN MODE 3, where an average increase in TCP of 0.73% (+/- 0.20, 95% confidence interval) was predicted by the biological models. This increase arises from a beneficial dose gradient which is produced across the tumour. Although the TCP gain is small it comes with no increase in treatment complexity, and could translate into increased cures given the large numbers of patients being referred. A study of the beam-weights and wedge angles chosen by the optimisation algorithm revealed significant inter-patient variability the implications of which are examined. Probably the most significant benefit of the algorithm is the time saved (about a factor of 10) in computing optimised beam-weights and wedge angles for this simple plan.

Comparison of treatment techniques for conformal radiotherapy of the prostate using dose-volume histograms and normal tissue complication probabilities

The aim of this study was to evaluate the relative merits of the coplanar field arrangements most frequently used for conformal radiotherapy of the prostate using dose-volume histograms and normal tissue complication probabilities (NTCPs). Twelve patients with early prostate cancer underwent a planning CT scan of the pelvis. Isocentric plans for each patient were devised using three, four, six and eight conformal fields and beam-weights optimised using fast simulated annealing to give a dose homogeneity across the planning target volume of +/- 5% or better while minimising irradiation of the relevant organs at risk. The plans were then evaluated using dose-volume histograms of the organs at risk (bladder, rectum and both femoral heads) and the Lyman model of normal tissue complication probability for the rectum. Analysis of dose-volume histogram data averaged over the 12 patients indicates an advantage for six (p = 0.002) and eight (p = 0.0001) fields with respect to the percentage volume of the femoral heads receiving > 50% of the prescribed dose compared with three fields. There was a similar advantage for six (p = 0.0007) and eight (p = 0.0001) fields compared with four fields. Ranking of the treatment techniques indicates that the four-field technique is the worst with respect to femoral head irradiation but the best with respect to reducing rectal irradiation. A higher dose can be prescribed to the isocentre with the four-field technique for a 5% rectal NTCP. The six-field technique led to sparing of the bladder when the different treatment techniques were ranked using bladder dose-volume histogram data. We conclude that none of the techniques studied consistently proved to be superior when applied to this sample of patients with prostate cancer with respect to sparing all the organs at risk. The absolute differences between techniques are small and would be very difficult to detect with respect to clinically relevant endpoints.