Improved efficiency of multi-criteria IMPT treatment planning using iterative resampling of randomly placed pencil beams

This study investigates whether 'pencil beam resampling', i.e. iterative selection and weight optimization of randomly placed pencil beams (PBs), reduces optimization time and improves plan quality for multi-criteria optimization in intensity-modulated proton therapy, compared with traditional modes in which PBs are distributed over a regular grid. Resampling consisted of repeatedly performing: (1) random selection of candidate PBs from a very fine grid, (2) inverse multi-criteria optimization, and (3) exclusion of low-weight PBs. The newly selected candidate PBs were added to the PBs in the existing solution, causing the solution to improve with each iteration. Resampling and traditional regular grid planning were implemented into our in-house developed multi-criteria treatment planning system 'Erasmus iCycle'. The system optimizes objectives successively according to their priorities as defined in the so-called 'wish-list'. For five head-and-neck cancer patients and two PB widths (3 and 6 mm sigma at 230 MeV), treatment plans were generated using: (1) resampling, (2) anisotropic regular grids and (3) isotropic regular grids, while using varying sample sizes (resampling) or grid spacings (regular grid). We assessed differences in optimization time (for comparable plan quality) and in plan quality parameters (for comparable optimization time). Resampling reduced optimization time by a factor of 2.8 and 5.6 on average (7.8 and 17.0 at maximum) compared with the use of anisotropic and isotropic grids, respectively. Doses to organs-at-risk were generally reduced when using resampling, with median dose reductions ranging from 0.0 to 3.0 Gy (maximum: 14.3 Gy, relative: 0%-42%) compared with anisotropic grids and from -0.3 to 2.6 Gy (maximum: 11.4 Gy, relative: -4%-19%) compared with isotropic grids. Resampling was especially effective when using thin PBs (3 mm sigma). Resampling plans contained on average fewer PBs, energy layers and protons than anisotropic grid plans and more energy layers and protons than isotropic grid plans. In conclusion, resampling resulted in improved plan quality and in considerable optimization time reduction compared with traditional regular grid planning.

Accurate IMRT fluence verification for prostate cancer patients using 'in-vivo' measured EPID images and in-room acquired kilovoltage cone-beam CT scans

Background

To investigate for prostate cancer patients the comparison of ‘in-vivo’ measured portal dose images (PDIs) with predictions based on a kilovoltage cone-beam CT scan (CBCT), acquired during the same treatment fraction, as an alternative for pre-treatment verification. For evaluation purposes, predictions were also performed using the patients’ planning CTs (pCT).

Methods

To get reliable CBCT electron densities for PDI predictions, Hounsfield units from the pCT were mapped onto the CBCT, while accounting for non-rigidity in patient anatomy in an approximate way. PDI prediction accuracy was first validated for an anatomical phantom, using IMRT treatment plans of ten prostate cancer patients. Clinical performance was studied using data acquired for 50 prostate cancer patients. For each patient, 4–5 CBCTs were available, resulting in a total of 1413 evaluated images. Measured and predicted PDIs were compared using γ-analyses with 3% global dose difference and 3 mm distance to agreement as reference criteria. Moreover, the pass rate for automated PDI comparison was assessed. To quantify improvements in IMRT fluence verification accuracy results from multiple fractions were combined by generating a γ-image with values halfway the minimum and median γ values, pixel by pixel.

Results

For patients, CBCT-based PDI predictions showed a high agreement with measurements, with an average percentage of rejected pixels of 1.41% only. In spite of possible intra-fraction motion and anatomy changes, this was only slightly larger than for phantom measurements (0.86%). For pCT-based predictions, the agreement deteriorated (average percentage of rejected pixels 2.98%), due to an enhanced impact of anatomy variations. For predictions based on CBCT, combination of the first 2 fractions yielded gamma results in close agreement with pre-treatment analyses (average percentage of rejected pixels 0.63% versus 0.35%, percentage of rejected beams 0.6% versus 0%). For the pCT-based approach, only combination of the first 5 fractions resulted in acceptable agreement with pre-treatment results.

Conclusion

In-room acquired CBCT scans can be used for high accuracy IMRT fluence verification based on in-vivo measured EPID images. Combination of γ results for the first 2 fractions can largely compensate for small accuracy reductions, with respect to pre-treatment verification, related to intra-fraction motion and anatomy changes.

Integrated multicriterial optimization of beam angles and intensity profiles for coplanar and noncoplanar head and neck IMRT and implications for VMAT

PURPOSE

To quantify improved salivary gland sparing for head and neck cancer patients using intensity-modulated radiotherapy (IMRT) plans based on integrated computerized optimization of beam orientations and intensity profiles. To assess if optimized nonzero couch angles also improve VMAT plans.

METHODS

Our in-house developed algorithm iCycle was used for automated generation of multicriterial optimized plans with optimized beam orientations and intensity profiles, and plans with optimized profiles for preselected beam arrangements. For 20 patients, five IMRT plans, based on one "wish-list," were compared: (i) and (ii) seven- and nine-beam equiangular coplanar plans (iCycle(7equi), iCycle(9equi)), (iii) and (iv) nine-beam plans with optimized coplanar and noncoplanar beam orientations (iCycle(copl), iCycle(noncopl)), and (v) a nine-beam coplanar plan with optimized gantry angles and one optimized couch rotation (iCycle(couch)). VMAT plans without and with this optimized couch rotation were evaluated.

RESULTS

iCycle(noncopl) resulted in the best salivary gland sparing, while iCycle(couch) yielded similar results for 18 patients. For iCycle(7equi), submandibular gland NTCP values were on average 5% higher. iCycle(9equi) performed better than iCycle(7equi). iCycle(copl) showed further improvement. Application of the optimized couch angle from iCycle(couch) also improved NTCP values in VMAT plans.

CONCLUSIONS

iCycle allows objective comparison of competing planning strategies. Integrated optimization of beam profiles and angles can significantly improve normal tissue sparing, yielding optimal results for iCycle(noncopl).

On the beam direction search space in computerized non-coplanar beam angle optimization for IMRT-prostate SBRT

In a recent paper, we have published a new algorithm, designated 'iCycle', for fully automated multi-criterial optimization of beam angles and intensity profiles. In this study, we have used this algorithm to investigate the relationship between plan quality and the extent of the beam direction search space, i.e. the set of candidate beam directions that may be selected for generating an optimal plan. For a group of ten prostate cancer patients, optimal IMRT plans were made for stereotactic body radiation therapy (SBRT), mimicking high dose rate brachytherapy dosimetry. Plans were generated for five different beam direction input sets: a coplanar (CP) set and four non-coplanar (NCP) sets. For CP treatments, the search space consisted of 72 orientations (5° separations). The NCP CyberKnife (CK) space contained all directions available in the robotic CK treatment unit. The fully non-coplanar (F-NCP) set facilitated the highest possible degree of freedom in selecting optimal directions. CK(+) and CK(++) were subsets of F-NCP to investigate some aspects of the CK space. For each input set, plans were generated with up to 30 selected beam directions. Generated plans were clinically acceptable, according to an assessment of our clinicians. Convergence in plan quality occurred only after around 20 included beams. For individual patients, variations in PTV dose delivery between the five generated plans were minimal, as aimed for (average spread in V(95): 0.4%). This allowed plan comparisons based on organ at risk (OAR) doses, with the rectum considered most important. Plans generated with the NCP search spaces had improved OAR sparing compared to the CP search space, especially for the rectum. OAR sparing was best with the F-NCP, with reductions in rectum D(Mean), V(40Gy), V(60Gy) and D(2%) compared to CP of 25%, 35%, 37% and 8%, respectively. Reduced rectum sparing with the CK search space compared to F-NCP could be largely compensated by expanding CK with beams with relatively large direction components along the superior-inferior axis (CK(++)). Addition of posterior beams (CK(++) → F-NCP) did not lead to further improvements in OAR sparing. Plans with 25 beams clearly performed better than 11-beam plans. For CP plans, an increase from 11 to 25 involved beams resulted in reductions in rectum D(Mean), V(40Gy), V(60Gy) and D(2%) of 39%, 57%, 64% and 13%, respectively.

Toward fully automated multicriterial plan generation: a prospective clinical study

PURPOSE

To prospectively compare plans generated with iCycle, an in-house-developed algorithm for fully automated multicriterial intensity modulated radiation therapy (IMRT) beam profile and beam orientation optimization, with plans manually generated by dosimetrists using the clinical treatment planning system.

METHODS AND MATERIALS

For 20 randomly selected head-and-neck cancer patients with various tumor locations (of whom 13 received sequential boost treatments), we offered the treating physician the choice between an automatically generated iCycle plan and a manually optimized plan using standard clinical procedures. Although iCycle used a fixed "wish list" with hard constraints and prioritized objectives, the dosimetrists manually selected the beam configuration and fine tuned the constraints and objectives for each IMRT plan. Dosimetrists were not informed in advance whether a competing iCycle plan was made. The 2 plans were simultaneously presented to the physician, who then selected the plan to be used for treatment. For the patient group, differences in planning target volume coverage and sparing of critical tissues were quantified.

RESULTS

In 32 of 33 plan comparisons, the physician selected the iCycle plan for treatment. This highly consistent preference for the automatically generated plans was mainly caused by the improved sparing for the large majority of critical structures. With iCycle, the normal tissue complication probabilities for the parotid and submandibular glands were reduced by 2.4% ± 4.9% (maximum, 18.5%, P=.001) and 6.5% ± 8.3% (maximum, 27%, P=.005), respectively. The reduction in the mean oral cavity dose was 2.8 ± 2.8 Gy (maximum, 8.1 Gy, P=.005). For the swallowing muscles, the esophagus and larynx, the mean dose reduction was 3.3 ± 1.1 Gy (maximum, 9.2 Gy, P<.001). For 15 of the 20 patients, target coverage was also improved.

CONCLUSIONS

In 97% of cases, automatically generated plans were selected for treatment because of the superior quality. Apart from the improved plan quality, automatic plan generation is economically attractive because of the reduced workload.

iCycle: Integrated, multicriterial beam angle, and profile optimization for generation of coplanar and noncoplanar IMRT plans

PURPOSE

To introduce iCycle, a novel algorithm for integrated, multicriterial optimization of beam angles, and intensity modulated radiotherapy (IMRT) profiles.

METHODS

A multicriterial plan optimization with iCycle is based on a prescription called wish-list, containing hard constraints and objectives with ascribed priorities. Priorities are ordinal parameters used for relative importance ranking of the objectives. The higher an objective priority is, the higher the probability that the corresponding objective will be met. Beam directions are selected from an input set of candidate directions. Input sets can be restricted, e.g., to allow only generation of coplanar plans, or to avoid collisions between patient/couch and the gantry in a noncoplanar setup. Obtaining clinically feasible calculation times was an important design criterium for development of iCycle. This could be realized by sequentially adding beams to the treatment plan in an iterative procedure. Each iteration loop starts with selection of the optimal direction to be added. Then, a Pareto-optimal IMRT plan is generated for the (fixed) beam setup that includes all so far selected directions, using a previously published algorithm for multicriterial optimization of fluence profiles for a fixed beam arrangement Breedveld et al. [Phys. Med. Biol. 54, 7199-7209 (2009)]. To select the next direction, each not yet selected candidate direction is temporarily added to the plan and an optimization problem, derived from the Lagrangian obtained from the just performed optimization for establishing the Pareto-optimal plan, is solved. For each patient, a single one-beam, two-beam, three-beam, etc. Pareto-optimal plan is generated until addition of beams does no longer result in significant plan quality improvement. Plan generation with iCycle is fully automated.

RESULTS

Performance and characteristics of iCycle are demonstrated by generating plans for a maxillary sinus case, a cervical cancer patient, and a liver patient treated with SBRT. Plans generated with beam angle optimization did better meet the clinical goals than equiangular or manually selected configurations. For the maxillary sinus and liver cases, significant improvements for noncoplanar setups were seen. The cervix case showed that also in IMRT with coplanar setups, beam angle optimization with iCycle may improve plan quality. Computation times for coplanar plans were around 1-2 h and for noncoplanar plans 4-7 h, depending on the number of beams and the complexity of the site.

CONCLUSIONS

Integrated beam angle and profile optimization with iCycle may result in significant improvements in treatment plan quality. Due to automation, the plan generation workload is minimal. Clinical application has started.

Treatment precision of image-guided liver SBRT using implanted fiducial markers depends on marker-tumour distance

The purpose of this study is to assess the accuracy of day-to-day predictions of liver tumour position using implanted gold markers as surrogates and to compare the method with alternative set-up strategies, i.e. no correction, vertebrae and 3D diaphragm-based set-up. Twenty patients undergoing stereotactic body radiation therapy (SBRT) with abdominal compression for primary or metastatic liver cancer were analysed. We determined the day-to-day correlation between gold marker and tumour positions in contrast-enhanced CT scans acquired at treatment preparation and before each treatment session. The influence of marker-tumour distance on the accuracy of prediction was estimated by introducing a method extension of the set-up error paradigm. The distance between gold markers and the centre of the tumour varied between 5 and 96 mm. Marker-guidance was superior to guiding treatment using other surrogates, although both the random and systematic components of the prediction error SD depended on the tumour-marker distance. For a marker-tumour distance of 4 cm, we observed σ = 1.3 mm and Σ = 1.6 mm. The 3D position of the diaphragm dome was the second best predictor. In conclusion, the tumour position can be predicted accurately using implanted markers, but marker-guided set-up accuracy decreases with increasing distance between implanted markers and the tumour.

Intensity modulated radiation therapy planning for patients with a metal hip prosthesis based on class solutions

PURPOSE

With the aging of the population, an increasing number of patients with metallic hip implants are referred for radiotherapy treatment. Class solutions for intensity modulated radiation therapy (IMRT) treatment planning are generally not applicable for these patients due to the required avoidance of dose delivery through prostheses. In this work a new approach for IMRT planning is presented, allowing the use of a default beam setup.

METHODS AND MATERIALS

For IMRT planning, Monaco (Elekta; CMS Software, Maryland Heights, MO) was used. In addition to the target and organs at risk, so-called prosthesis avoidance volumes (PAVs) were delineated in the beam's eye view projection for beams in which the prosthesis was partially in front of the target. By putting strict constraints on these virtual organs at risk, entrance dose delivery through a prosthesis is avoided while exit dose delivery is allowed. In this way, uncertainties in the dose delivery to the target and organs at risk, as derived by the treatment planning system, are largely minimized. To show the advantages of this IMRT-PAV technique, for 2 prostate cancer patients, 1 with bilateral and the other with unilateral metallic hip prostheses, obtained IMRT plans were compared with conventional IMRT plans using a prosthesis-avoiding beam setup.

RESULTS

For both IMRT techniques a similar planning target volume coverage was achieved, but with the IMRT-PAV technique the mean doses to the bladder and the rectum were reduced by up to 25%. While the IMRT-PAV technique required more time for delineation, the time for treatment planning reduced because the default beam setup could be applied. The number of segments needed for dose delivery was comparable for both techniques.

CONCLUSIONS

With the new IMRT-PAV technique IMRT class solutions can safely be applied for cancer patients with metallic hip prostheses, generally yielding a reduced dose delivery to organs at risk or improved target coverage.

IMRT for image-guided single vocal cord irradiation

PURPOSE

We have been developing an image-guided single vocal cord irradiation technique to treat patients with stage T1a glottic carcinoma. In the present study, we compared the dose coverage to the affected vocal cord and the dose delivered to the organs at risk using conventional, intensity-modulated radiotherapy (IMRT) coplanar, and IMRT non-coplanar techniques.

METHODS AND MATERIALS

For 10 patients, conventional treatment plans using two laterally opposed wedged 6-MV photon beams were calculated in XiO (Elekta-CMS treatment planning system). An in-house IMRT/beam angle optimization algorithm was used to obtain the coplanar and non-coplanar optimized beam angles. Using these angles, the IMRT plans were generated in Monaco (IMRT treatment planning system, Elekta-CMS) with the implemented Monte Carlo dose calculation algorithm. The organs at risk included the contralateral vocal cord, arytenoids, swallowing muscles, carotid arteries, and spinal cord. The prescription dose was 66 Gy in 33 fractions.

RESULTS

For the conventional plans and coplanar and non-coplanar IMRT plans, the population-averaged mean dose ± standard deviation to the planning target volume was 67 ± 1 Gy. The contralateral vocal cord dose was reduced from 66 ± 1 Gy in the conventional plans to 39 ± 8 Gy and 36 ± 6 Gy in the coplanar and non-coplanar IMRT plans, respectively. IMRT consistently reduced the doses to the other organs at risk.

CONCLUSIONS

Single vocal cord irradiation with IMRT resulted in good target coverage and provided significant sparing of the critical structures. This has the potential to improve the quality-of-life outcomes after RT and maintain the same local control rates.

Increasing treatment accuracy for cervical cancer patients using correlations between bladder-filling change and cervix-uterus displacements: proof of principle

PURPOSE

To investigate application of pre-treatment established correlations between bladder-filling changes and cervix-uterus displacements in adaptive therapy.

MATERIALS AND METHODS

Thirteen cervical cancer patients participated in this prospective study. Pre-treatment, and after delivery of 40 Gy, a full bladder CT-scan was acquired, followed by voiding the bladder and acquisition of 4 other 3D scans in a 1h period with a naturally filling bladder (variable bladder filling CT-scans, VBF-scans). For the pre-treatment VBF-scans, linear correlations between bladder volume change and displacements of the tip of the uterus (ToU) and the center of mass (CoM) of markers implanted in the fornices of the vagina relative to the full bladder planning scan were established. Prediction accuracy of these correlation models was assessed by comparison with actual displacements in CT-scans, both pre-treatment and after 40 Gy. Inter-fraction ToU and marker-CoM displacements were derived from the established correlations and twice-weekly performed in-room bladder volume measurements, using a 3D ultrasound scanner.

RESULTS

Target displacement in VBF-scans ranged from up to 65 mm in a single direction to almost 0mm, depending on the patient. For pre-treatment VBF-scans, the linear correlation models predicted the mean 3D position change for the ToU of 26.1 mm±10.8 with a residual of only 2.2 mm±1.7. For the marker-CoM, the 8.4 mm±5.3 mean positioning error was predicted with a residual of 0.9 mm±0.7. After 40Gy, the mean ToU displacement was 26.8 mm±15.8, while prediction based on the pre-treatment correlation models yielded a mean residual error of 9.0 mm±3.7. Target positioning errors in the fractioned treatments were very large, especially for the ToU (-18.5mm±11.2 for systematic errors in SI-direction).

CONCLUSIONS

Pre-treatment acquired VBF-scans may be used to substantially enhance treatment precision of cervical cancer patients. Application in adaptive therapy is promising and warrants further investigation. For highly conformal (IMRT) treatments, the use of a full bladder drinking protocol results in unacceptably large systematic set-up errors.

Comparison of macroscopic pathology measurements with magnetic resonance imaging and assessment of microscopic pathology extension for colorectal liver metastases

PURPOSE

To compare pathology macroscopic tumor dimensions with magnetic resonance imaging (MRI) measurements and to establish the microscopic tumor extension of colorectal liver metastases.

METHODS AND MATERIALS

In a prospective pilot study we included patients with colorectal liver metastases planned for surgery and eligible for MRI. A liver MRI was performed within 48 hours before surgery. Directly after surgery, an MRI of the specimen was acquired to measure the degree of tumor shrinkage. The specimen was fixed in formalin for 48 hours, and another MRI was performed to assess the specimen/tumor shrinkage. All MRI sequences were imported into our radiotherapy treatment planning system, where the tumor and the specimen were delineated. For the macroscopic pathology analyses, photographs of the sliced specimens were used to delineate and reconstruct the tumor and the specimen volumes. Microscopic pathology analyses were conducted to assess the infiltration depth of tumor cell nests.

RESULTS

Between February 2009 and January 2010 we included 13 patients for analysis with 21 colorectal liver metastases. Specimen and tumor shrinkage after resection and fixation was negligible. The best tumor volume correlations between MRI and pathology were found for T1-weighted (w) echo gradient sequence (r(s) = 0.99, slope = 1.06), and the T2-w fast spin echo (FSE) single-shot sequence (r(s) = 0.99, slope = 1.08), followed by the T2-w FSE fat saturation sequence (r(s) = 0.99, slope = 1.23), and the T1-w gadolinium-enhanced sequence (r(s) = 0.98, slope = 1.24). We observed 39 tumor cell nests beyond the tumor border in 12 metastases. Microscopic extension was found between 0.2 and 10 mm from the main tumor, with 90% of the cases within 6 mm.

CONCLUSIONS

MRI tumor dimensions showed a good agreement with the macroscopic pathology suggesting that MRI can be used for accurate tumor delineation. However, microscopic extensions found beyond the tumor border indicate that caution is needed in selecting appropriate tumor margins.

Three-dimensional dose addition of external beam radiotherapy and brachytherapy for oropharyngeal patients using nonrigid registration

PURPOSE

To develop and evaluate a method for adding dose distributions of combined external beam radiotherapy (EBRT) and brachytherapy (BT) for oropharyngeal patients.

METHODS AND MATERIALS

Two computed tomography (CT) scans were used for 5 patients: the EBRT CT, used for EBRT planning, and the BT CT, acquired after catheter implantation. For each scan, the salivary glands and the chewing and swallowing muscles were contoured, and a dose distribution was calculated. A nonrigid transformation was obtained by registering the organs' surfaces. Then the BT dose distribution was mapped onto the EBRT dose distribution by applying the transformation obtained. To account for differences in fractionation, the physical doses were converted to equivalent dose in 2 Gy (EQD(2)), and the total dose was found by adding dose voxel by voxel. The robustness of the dose addition was investigated by varying delineations and input parameters of the registration method and by varying the α/β parameter for EQD(2). The effect of the perturbations was quantified using dose-volume histograms (DVH) and gamma analyses (distance-to-agreement/dose-difference = 1 mm/1 Gy).

RESULTS

The variations in input parameters and delineations caused only small perturbations in the DVH of the added dose distributions. For most organs the gamma index was low, and it was moderately elevated for organs lying in areas with a steep gradient (median gamma index ≤ 2.3 for constrictor muscles, ≤ 0.7 for all other organs).

CONCLUSIONS

The presented method allows adding dose distributions of combined EBRT and BT for oropharyngeal patients. In general, the method is reliable and robust with respect to uncertainties in organ delineation, perturbations in input parameters of the method, and α/β values.

Implications of artefacts reduction in the planning CT originating from implanted fiducial markers

The efficacy of metal artefact reduction (MAR) software to suppress artefacts in reconstructed computed tomography (CT) images originating from small metal objects, like tumor markers and surgical clips, was evaluated. In addition, possible implications of using digital reconstructed radiographs (DRRs), based on the MAR CT images, for setup verification were analyzed. A phantom and 15 patients with different tumor sites and implanted markers were imaged with a multislice CT scanner. The raw image data was reconstructed both with the clinically used filtered-backprojection (FBP) and with the MAR software. Using the MAR software, improvements in image quality were often observed in CT slices with markers or clips. Especially when several markers were located near to each other, fewer streak artefacts were observed than with the FBP algorithm. In addition, the shape and size of markers could be identified more accurately, reducing the contoured marker volumes by a factor of 2. For the phantom study, the CT numbers measured near to the markers corresponded more closely to the expected values. However, the MAR images were slightly more smoothed compared with the images reconstructed with FBP. For 8 prostate cancer patients in this study, the interobserver variation in 3D marker definition was similar (<0.4 mm) when using DRRs based on either FBP or MAR CT scans. Automatic marker matches also showed a similar success rate. However, differences in automatic match results up to 1 mm, caused by differences in the marker definition, were observed, which turned out to be (borderline) statistically significant (p = 0.06) for 2 patients. In conclusion, the MAR software might improve image quality by suppressing metal artefacts, probably allowing for a more reliable delineation of structures. When implanted markers or clips are used for setup verification, the accuracy may slightly be improved as well, which is relevant when using very tight clinical target volume (CTV) to planning target volume (PTV) margins for planning.

The equivalence of multi-criteria methods for radiotherapy plan optimization

Several methods can be used to achieve multi-criteria optimization of radiation therapy treatment planning, which strive for Pareto-optimality. The property of the solution being Pareto optimal is desired, because it guarantees that no criteria can be improved without deteriorating another criteria. The most widely used methods are the weighted-sum method, in which the different treatment objectives are weighted, and constrained optimization methods, in which treatment goals are set and the algorithm has to find the best plan fulfilling these goals. The constrained method used in this paper, the 2p element of c (2-phase element-constraint) method is based on the element-constraint method, which generates Pareto-optimal solutions. Both approaches are uniquely related to each other. In this paper, we will show that it is possible to switch from the constrained method to the weighted-sum method by using the Lagrange multipliers from the constrained optimization problem, and vice versa by setting the appropriate constraints. In general, the theory presented in this paper can be useful in cases where a new situation is slightly different from the original situation, e.g. in online treatment planning, with deformations of the volumes of interest, or in automated treatment planning, where changes to the automated plan have to be made. An example of the latter is given where the planner is not satisfied with the result from the constrained method and wishes to decrease the dose in a structure. By using the Lagrange multipliers, a weighted-sum optimization problem is constructed, which generates a Pareto-optimal solution in the neighbourhood of the original plan, but fulfills the new treatment objectives.

Automated non-coplanar beam direction optimization improves IMRT in SBRT of liver metastasis

PURPOSE

To investigate whether automatically optimized coplanar, or non-coplanar beam setups improve intensity modulated radiotherapy (IMRT) treatment plans for stereotactic body radiotherapy (SBRT) of liver tumors, compared to a reference equi-angular IMRT plan.

METHODS

For a group of 13 liver patients, an in-house developed beam selection algorithm (Cycle) was used for generation of 3D-CRT plans with either optimized coplanar-, or non-coplanar beam setups. These 10 field, coplanar and non-coplanar setups, and an 11 field, equi-angular coplanar reference setup were then used as input for generation of IMRT plans. For all plans, the PTV dose was maximized in an iterative procedure by increasing the prescribed PTV dose in small steps until further increase was prevented by constraint violation(s).

RESULTS

For optimized non-coplanar setups, D(PTV, max) increased by on average 30% (range 8-64%) compared to the corresponding reference IMRT plan. Similar increases were observed for D(PTV, 99%) and gEUD(a). For optimized coplanar setups, mean PTV dose increases were only approximately 4%. After re-scaling all plans to the clinically applied dose, optimized non-coplanar configurations resulted in the best sparing of organs at risk (healthy liver, spinal cord, bowel).

CONCLUSION

Compared to an equi-angular beam setup, computer optimized non-coplanar setups do result in substantial improvements in IMRT plans for SBRT of liver tumors.

PTV dose prescription strategies for SBRT of metastatic liver tumours

PURPOSE

Recently we have demonstrated that our in-house developed algorithm for automated plan generation for fully non-coplanar SBRT of liver patients (designated Cycle) yields plans that are superior to conventionally generated plans of experienced dosimetrists. Here we use Cycle in the comparison of plans with prescription isodoses of 65% or 80% of the isocentre dose.

METHODS

Plans were generated using CT-data of 15 previously treated patients. For each patient, both for the 65%- and the 80% strategy, Cycle was used to generate a plan with the maximum isocentre dose, D(isoc), while strictly obeying a set of hard constraints for the organs at risk (OAR). Plans for the two strategies were compared using D(isoc), D(PTV,99%) (the minimum dose delivered to 99% of the PTV), and the generalised equivalent uniform dose, gEUD(PTV)(a), for several values of the parameter a. Moreover, for the OARs, the distance to the constraint values was analysed.

RESULTS

The 65% strategy resulted in treatment plans with a higher D(isoc) (average 17.6%, range 7.6-31.1%) than the 80% strategy, at the cost of a somewhat lower D(PTV,99%) (average -2.0%, range -9.6% to 9.3%). On average, voxels with a dose in the 65% strategy, lower than the minimum PTV dose in the 80% strategy, were within 0.2cm from the PTV surface. For a-10, the 65% strategy yielded on average a significantly (P<0.01) higher gEUD(PTV)(a) than the 80% strategy, whereas for highly negative a-values the 80% approach was slightly better, although not significantly. Large variations between patients were observed. Generally, for the OAR the approach to the constraint levels was similar for the two strategies.

CONCLUSION

On average, PTV dose delivery is superior with the 65% strategy. However, apart from the isocentre dose, for each applied PTV dose parameter at least one patient would have been better off with the 80% dose prescription strategy.

A novel approach to multi-criteria inverse planning for IMRT

Treatment plan optimization is a multi-criteria process. Optimizing solely on one objective or on a sum of a priori weighted objectives may result in inferior treatment plans. Manually adjusting weights or constraints in a trial and error procedure is time consuming. In this paper we introduce a novel multi-criteria optimization approach to automatically optimize treatment constraints (dose-volume and maximum-dose). The algorithm tries to meet these constraints as well as possible, but in the case of conflicts it relaxes lower priority constraints so that higher priority constraints can be met. Afterwards, all constraints are tightened, starting with the highest priority constraints. Applied constraint priority lists can be used as class solutions for patients with similar tumour types. The presented algorithm does iteratively apply an underlying algorithm for beam profile optimization, based on a quadratic objective function with voxel-dependent importance factors. These voxel-dependent importance factors are automatically adjusted to reduce dose-volume and maximum-dose constraint violations.