On-line re-optimization of prostate IMRT plans for adaptive radiation therapy

Dosimetric impact of variable air cavity within PTV for rectum cancer

PURPOSE

The aim of this study is to determine the impact of rectal air volume changes on treatment plan quality, and subsequently inform daily cone-beam computed tomography (CBCT) evaluation constraints, in terms of acceptable rectal air volume during treatment.

METHODS

Twelve rectal cancer patients who exhibited rectal air within the PTV on their planning CT were selected. A study was conducted to evaluate the deterioration in plan quality due to expanding air volume. For each case, the air cavity volume was isotropically expanded in three dimensions using predefined margins of 3, 5, 7, and 10 mm, while deforming bladder and rectum contours. A constraint was applied to the bony anatomy to restrict the deformation. Treatment plans were then generated for all twelve patients by recalculating the reference plan with the expanded air cavity volume.

RESULTS

As the air cavity expanded, the maximum relative change in D98% coverage, compared to the reference plan, decreased by 10.8% ± 3.5%, while the D2% increased by 3.5% ± 0.9%. The positioning of the air cavity notably influenced the D98% variability with the 3 mm expansion. D98% coverage falls below 95% when the air cavity volume exceeds 17 cm3. On average, D2% coverage increased by 0.5% with each expansion. At the largest expansion, extensive coverage of 102% and 105% isodoses was observed compared to the reference plan.

CONCLUSION

Air cavity volumes above 17 cm3 can potentially degrade the high-dose PTV coverage while increasing the regions covered by the 102% and 105% isodoses. Clinical CBCT guidelines were deduced, recommending a maximum threshold of 3.2 cm in diameter in any direction.

The geometric and dosimetric accuracy of kilovoltage cone beam computed tomography images for adaptive treatment: a systematic review

OBJECTIVES

To provide an overview and meta-analysis of different techniques adopted to accomplish kVCBCT for dose calculation and automated segmentation.

METHODS

A systematic review and meta-analysis were performed on eligible studies demonstrating kVCBCT-based dose calculation and automated contouring of different tumor features. Meta-analysis of the performance was accomplished on the reported γ analysis and dice similarity coefficient (DSC) score of both collected results as three subgroups (head and neck, chest, and abdomen).

RESULTS

After the literature scrutinization (n = 1008), 52 papers were recognized for the systematic review. Nine studies of dosimtric studies and eleven studies of geometric analysis were suitable for inclusion in meta-analysis. Using kVCBCT for treatment replanning depends on a method used. Deformable Image Registration (DIR) methods yielded small dosimetric error (≤2%), γ pass rate (≥90%) and DSC (≥0.8). Hounsfield Unit (HU) override and calibration curve-based methods also achieved satisfactory yielded small dosimetric error (≤2%) and γ pass rate ((≥90%), but they are prone to error due to their sensitivity to a vendor-specific variation in kVCBCT image quality.

CONCLUSIONS

Large cohorts of patients ought to be undertaken to validate methods achieving low levels of dosimetric and geometric errors. Quality guidelines should be established when reporting on kVCBCT, which include agreed metrics for reporting on the quality of corrected kVCBCT and defines protocols of new site-specific standardized imaging used when obtaining kVCBCT images for adaptive radiotherapy.

ADVANCES IN KNOWLEDGE

This review gives useful knowledge about methods making kVCBCT feasible for kVCBCT-based adaptive radiotherapy, simplifying patient pathway and reducing concomitant imaging dose to the patient.

Online adaptive planning methods for intensity-modulated radiotherapy

Online adaptive radiation therapy aims at adapting a patient's treatment plan to their current anatomy to account for inter-fraction variations before daily treatment delivery. As this process needs to be accomplished while the patient is immobilized on the treatment couch, it requires time-efficient adaptive planning methods to generate a quality daily treatment plan rapidly. The conventional planning methods do not meet the time requirement of online adaptive radiation therapy because they often involve excessive human intervention, significantly prolonging the planning phase. This article reviews the planning strategies employed by current commercial online adaptive radiation therapy systems, research on online adaptive planning, and artificial intelligence's potential application to online adaptive planning.

Dosimetric impact of deep learning-based CT auto-segmentation on radiation therapy treatment planning for prostate cancer

Background

The evaluation of automatic segmentation algorithms is commonly performed using geometric metrics. An analysis based on dosimetric parameters might be more relevant in clinical practice but is often lacking in the literature. The aim of this study was to investigate the impact of state-of-the-art 3D U-Net-generated organ delineations on dose optimization in radiation therapy (RT) for prostate cancer patients.

Methods

A database of 69 computed tomography images with prostate, bladder, and rectum delineations was used for single-label 3D U-Net training with dice similarity coefficient (DSC)-based loss. Volumetric modulated arc therapy (VMAT) plans have been generated for both manual and automatic segmentations with the same optimization settings. These were chosen to give consistent plans when applying perturbations to the manual segmentations. Contours were evaluated in terms of DSC, average and 95% Hausdorff distance (HD). Dose distributions were evaluated with the manual segmentation as reference using dose volume histogram (DVH) parameters and a 3%/3 mm gamma-criterion with 10% dose cut-off. A Pearson correlation coefficient between DSC and dosimetric metrics, i.e. gamma index and DVH parameters, has been calculated.

Results

3D U-Net-based segmentation achieved a DSC of 0.87 (0.03) for prostate, 0.97 (0.01) for bladder and 0.89 (0.04) for rectum. The mean and 95% HD were below 1.6 (0.4) and below 5 (4) mm, respectively. The DVH parameters, V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{60/65/70\,{\mathrm{Gy}}}$$\end{document}60/65/70Gy for the bladder and V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{50/65/70\,{\mathrm{Gy}}}$$\end{document}50/65/70Gy for the rectum, showed agreement between dose distributions within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \, 5\%$$\end{document}±5% and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm \,2\%$$\end{document}±2%, respectively. The D\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{98/2\%}$$\end{document}98/2% and V\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{95\%}$$\end{document}95%, for prostate and its 3 mm expansion (surrogate clinical target volume) showed agreement with the reference dose distribution within 2% and 3 Gy with the exception of one case. The average gamma pass-rate was 85%. The comparison between geometric and dosimetric metrics showed no strong statistically significant correlation.

Conclusions

The 3D U-Net developed for this work achieved state-of-the-art geometrical performance. Analysis based on clinically relevant DVH parameters of VMAT plans demonstrated neither excessive dose increase to OARs nor substantial under/over-dosage of the target in all but one case. Yet the gamma analysis indicated several cases with low pass rates. The study highlighted the importance of adding dosimetric analysis to the standard geometric evaluation.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

Clinical Experience in Prostate Ultrahypofractionated Radiation Therapy With an Online Adaptive Method

PURPOSE

This study aimed to describe the feasibility of the online adaptive radiation therapy (oART) method developed at the Hospital Quirónsalud Barcelona for prostate cancer, using a standard C-arm linear accelerator (linac) and without the support of artificial intelligence.

METHODS AND MATERIALS

The first 18 patients treated at the Hospital Quirónsalud Barcelona with the developed oART method were included. An ultrahypofractionated radiation therapy scheme consisting of 7 × 6.1 Gy was used. Patients were treated on 2 conventional Varian C-arm linacs. For each patient, a reference plan based on a planning computed tomography (pCT) scan was generated using the Eclipse system. On each treatment session, the pCT scan was rigidly registered with the daily cone beam computed tomography (CT) scan. The pCT-based target (prostate) and organs at risk were mapped onto the cone beam CT images and manually adapted to take into account the anatomy of the day. The reference plan was then copied to the cone beam CT scan, and a full reoptimization was done for the current anatomy (adapted plan). For each treatment session, the unaltered reference plan was recomputed on the daily cone beam CT scan by mimicking the soft-tissue alignment performed per our standard procedure (nonadapted plan). Over the 126 adapted sessions from the 18 patients, a dosimetric comparison of adapted against nonadapted plans was done.

RESULTS

A significant difference in the target coverage was found between the adapted and nonadapted plans (97.1 vs 90.4; P < .001) in favor of adapting. Without online adaptation, the optimal coverage of the prostate was not attained in 35% of fractions. Adapting allows for the improvement of the target coverage with compliance of all organ-at-risk dose constraints in all treatment fractions.

CONCLUSIONS

The oART technique described in this study is technically feasible with a C-arm linac. To our knowledge, this is the first clinical experience with oART for prostate cancer including full replanning and delivered with a C-arm linac without artificial intelligence capability.

Prospects for daily online adaptive radiotherapy via ethos for prostate cancer patients without nodal involvement using unedited CBCT auto-segmentation

Purpose

 Implementing new online adaptive radiation therapy technologies is challenging because extra clinical resources are required particularly expert contour review. Here, we provide the first assessment of Varian's Ethos™ adaptive platform for prostate cancer using no manual edits after auto‐segmentation to minimize this impact on clinical efficiency.

Methods

Twenty‐five prostate patients previously treated at our clinic were re‐planned using an Ethos™ emulator. Clinical target volumes (CTV) included intact prostate and proximal seminal vesicles. The following clinical margins were used: 3 mm posterior, 5 mm left/right/anterior, and 7 mm superior/inferior. Adapted plans were calculated for 10 fractions per patient using Ethos's auto‐segmentation and auto‐planning workflow without manual contouring edits. Doses and auto‐segmented structures were exported to our clinical treatment planning system where contours were modified as needed for all 250 CTVs and organs‐at‐risk. Dose metrics from adapted plans were compared to unadapted plans to evaluate CTV and OAR dose changes.

Results

Overall 96% of fractions required auto‐segmentation edits, although corrections were generally minor (<10% of the volume for 70% of CTVs, 88% of bladders, and 90% of rectums). However, for one patient the auto‐segmented CTV failed to include the superior portion of prostate that extended into the bladder at all 10 fractions resulting in under‐contouring of the CTV by 31.3% ± 6.7%. For the 24 patients with minor auto‐segmentation corrections, adaptation improved CTV D98% by 2.9% ± 5.3%. For non‐adapted fractions where bladder or rectum V90% exceeded clinical thresholds, adaptation reduced them by 13.1% ± 1.0% and 6.5% ± 7.3%, respectively.

Conclusion

 For most patients, Ethos's online adaptive radiation therapy workflow improved CTV D98% and reduced normal tissue dose when structures would otherwise exceed clinical thresholds, even without time‐consuming manual edits. However, for one in 25 patients, large contour edits were required and thus scrutiny of the daily auto‐segmentation is necessary and not all patients will be good candidates for adaptation.

Synthetic CT generation from CBCT images via unsupervised deep learning

Adaptive-radiation-therapy (ART) is applied to account for anatomical variations observed over the treatment course. Daily or weekly cone-beam computed tomography (CBCT) is commonly used in clinic for patient positioning, but CBCT's inaccuracy in Hounsfield units (HU) prevents its application to dose calculation and treatment planning. Adaptive re-planning can be performed by deformably registering planning CT (pCT) to CBCT. However, scattering artifacts and noise in CBCT decrease the accuracy of deformable registration and induce uncertainty in treatment plan. Hence, generating from CBCT a synthetic CT (sCT) that has the same anatomical structure as CBCT but accurate HU values is desirable for ART. We proposed an unsupervised style-transfer-based approach to generate sCT based on CBCT and pCT. Unsupervised learning was desired because exactly matched CBCT and CT are rarely available, even when they are taken a few minutes apart. In the proposed model, CBCT and pCT are two inputs that provide anatomical structure and accurate HU information, respectively. The training objective function is designed to simultaneously minimize (1) contextual loss between sCT and CBCT to maintain the content and structure of CBCT in sCT and (2) style loss between sCT and pCT to achieve pCT-like image quality in sCT. We used CBCT and pCT images of 114 patients to train and validate the designed model, and another 29 independent patient cases to test the model's effectiveness. We quantitatively compared the resulting sCT with the original CBCT using the deformed same-day pCT as reference. Structure-similarity-index, peak-signal-to-noise-ratio, and mean-absolute-error in HU of sCT were 0.9723, 33.68, and 28.52, respectively, while those of CBCT were 0.9182, 29.67, and 49.90, respectively. We have demonstrated the effectiveness of the proposed model in using CBCT and pCT to synthesize CT-quality images. This model may permit using CBCT for advanced applications such as adaptive treatment planning.

Measurements of human tolerance to horizontal rotation within an MRI scanner: Towards gantry-free radiation therapy

INTRODUCTION

Recent advances in image guidance and adaptive radiotherapy could enable gantry-free radiotherapy using patient rotation. Gantry-free radiotherapy could substantially reduce the cost of radiotherapy systems and facilities. MRI guidance complements a gantry-free approach because of its ability to visualise soft tissue deformation during rotation. A potential barrier to gantry-free radiotherapy is patient acceptability, especially when combined with MRI. This study investigates human experiences of horizontal rotation within an MRI scanner.

METHODS

Ten healthy human participants and nine participants previously treated with radiotherapy were rotated within an MRI scanner. Participants' anxiety and motion sickness was assessed before being rotated in 45-degree increments and paused, representing a multi-field intensity-modulated radiotherapy treatment. An MR image was acquired at each 45-degree angle. Following imaging, anxiety and motion sickness were re-assessed, followed by a comfort questionnaire and exit interview. The significance of the differences in anxiety and motion sickness pre- versus post-imaging was assessed using Wilcoxon signed-rank tests. Content analysis was performed on exit interview transcripts.

RESULTS

Eight of ten healthy and eight of nine patient participants completed the imaging session. Mean anxiety scores before and after imaging were 7.9/100 and 11.8/100, respectively (P = 0.26), and mean motion sickness scores were 5.3/100 and 13.7/100, respectively (P = 0.02). Most participants indicated likely acceptance of rotation if MRI were to be used in a hypothetical treatment. Physical discomfort was reported to be the biggest concern.

CONCLUSIONS

Horizontal rotation within an MRI scanner was acceptable for most (17/19) participants.

Initial Evaluation of a Novel Cone-Beam CT-Based Semi-Automated Online Adaptive Radiotherapy System for Head and Neck Cancer Treatment - A Timing and Automation Quality Study

Introduction

A novel on-line adaptive radiotherapy (ART) system based on O-ring linear accelerator (LINAC) and cone-beam CT (CBCT) was evaluated for treatment and management of head & neck (H&N) cancer in an emulated environment accessed via remote desktop connection. In this on-line ART system, organs-at-risk (OARs) and target contours and radiotherapy (RT) plans are semi-automatically generated based on the patient CBCT, expediting a typically hours-long RT planning session to under half an hour. In this paper, we describe our initial experiences with the system and explore optimization strategies to expedite the process further.

Methods

We retroactively studied five patients with head and neck cancers, treated 16-35 fractions to 50-70 Gys. For each patient, on-line ART was simulated with one planning CT and three daily CBCT images taken beginning, middle, and end of treatment (tx). Key OAR (mandible, parotids, and spinal cord) and target (planning target volume (PTV) = clinical target volume (CTV) + 3 mm margin) contours were auto-generated and adjusted as needed by therapist/dosimetrist and attending physician, respectively. Duration of OAR contouring, target contouring, and plan review was recorded. Key OAR auto-contours were qualitatively rated from 1 (unacceptable) - 5 (perfect OAR delineation), and then quantitatively compared to human-adjusted “ground truth” contours via dice similarity coefficient (DSC) and 95-percentile Hausdorff distance (HD95%). Once contours were approved, adapted RT plans were auto-generated for physician review. Simulated doses to OARs and targets from the adapted plan were compared to that from the original (un-adapted) plan.

Results

Median on-line ART planning duration in the remote emulated environment was 19 min 34 sec (range: 13 min 10 sec - 31 min 20 sec). Automated key OAR quality was satisfactory overall (98% scored ≥3; 82% ≥4), though mandible was rated lower than others (p < 0.05). Most key OARs and all targets were within 2 mm margin of human-adjusted contours, but a few parotid and spinal cord contours deviated up to 5 mm. Anatomical changes over tx course further increased auto-contour error (p < 0.05, ΔHD95% = 0.77 mm comparing start and end of tx). Further optimizing auto-contoured OAR and target quality could reduce the on-line treatment planning duration by ~5 min and ~4.5 min, respectively. Dosimetrically, adapted plan spared OARs at a rate much greater than random chance compared to the original plan (χ² = 22.3, p << 0.001), while maintaining similar therapeutic dose to treatment target CTV (χ² = 1.14, p > 0.05). In addition, a general decrease in accumulated OAR dose was observed with adaptation. Unsupervised adapted plans where contours were auto-generated without human review still spared OAR at a greater rate than the original plans, suggesting benefits of adaptation can be maintained even with some leniency in contour accuracy.

Conclusion

Feasibility of a novel, semi-automated on-line ART system for various head and neck (H&N) cancer sites was demonstrated in terms of treatment duration, dosimetric benefits, and automated contour accuracy in a remote emulator environment. Adaptive planning duration was clinically viable at 19 min and 34 sec, but further improvements in automated contour accuracy and performance improvements of plan auto-generation may reduce adaptive planning duration by up to 10 minutes.

Comparison of radiographer interobserver image registration variability using cone beam CT and MR for cervix radiotherapy

OBJECTIVES

The aim of this study was to assess the consistency of therapy radiographers performing image registration using cone beam computed tomography (CBCT)-CT, magnetic resonance (MR)-CT, and MR-MR image guidance for cervix cancer radiotherapy and to assess that MR-based image guidance is not inferior to CBCT standard practice.

METHODS

10 patients receiving cervix radiation therapy underwent daily CBCT guidance and magnetic resonance (MR) imaging weekly during treatment. Offline registration of each MR image, and corresponding CBCT, to planning CT was performed by five radiographers. MR images were also registered to the earliest MR interobserver variation was assessed using modified Bland-Altman analysis with clinically acceptable 95% limits of agreement (LoA) defined as ±5.0 mm.

RESULTS

30 CBCT-CT, 30 MR-CT and 20 MR-MR registrations were performed by each observer. Registration variations between CBCT-CT and MR-CT were minor and both strategies resulted in 95% LoA over the clinical threshold in the anteroposterior direction (CBCT-CT ±5.8 mm, MR-CT ±5.4 mm). MR-MR registrations achieved a significantly improved 95% LoA in the anteroposterior direction (±4.3 mm). All strategies demonstrated similar results in lateral and longitudinal directions.

CONCLUSION

The magnitude of interobserver variations between CBCT-CT and MR-CT were similar, confirming that MR-CT radiotherapy workflows are comparable to CBCT-CT image-guided radiotherapy. Our results suggest MR-MR radiotherapy workflows may be a superior registration strategy.

ADVANCES IN KNOWLEDGE

This is the first publication quantifying interobserver registration of multimodality image registration strategies for cervix radical radiotherapy patients.

Pelvic multi-organ segmentation on cone-beam CT for prostate adaptive radiotherapy

BACKGROUND AND PURPOSE

The purpose of this study is to develop a deep learning-based approach to simultaneously segment five pelvic organs including prostate, bladder, rectum, left and right femoral heads on cone-beam CT (CBCT), as required elements for prostate adaptive radiotherapy planning.

MATERIALS AND METHODS

We propose to utilize both CBCT and CBCT-based synthetic MRI (sMRI) for the segmentation of soft tissue and bony structures, as they provide complementary information for pelvic organ segmentation. CBCT images have superior bony structure contrast and sMRIs have superior soft tissue contrast. Prior to segmentation, sMRI was generated using a cycle-consistent adversarial networks (CycleGAN), which was trained using paired CBCT-MR images. To combine the advantages of both CBCT and sMRI, we developed a cross-modality attention pyramid network with late feature fusion. Our method processes CBCT and sMRI inputs separately to extract CBCT-specific and sMRI-specific features prior to combining them in a late-fusion network for final segmentation. The network was trained and tested using 100 patients' datasets, with each dataset including the CBCT and manual physician contours. For comparison, we trained another two networks with different network inputs and architectures. The segmentation results were compared to manual contours for evaluations.

RESULTS

For the proposed method, dice similarity coefficients and mean surface distances between the segmentation results and the ground truth were 0.96 ± 0.03, 0.65 ± 0.67 mm; 0.91 ± 0.08, 0.93 ± 0.96 mm; 0.93 ± 0.04, 0.72 ± 0.61 mm; 0.95 ± 0.05, 1.05 ± 1.40 mm; and 0.95 ± 0.05, 1.08 ± 1.48 mm for bladder, prostate, rectum, left and right femoral heads, respectively. As compared to the other two competing methods, our method has shown superior performance in terms of the segmentation accuracy.

CONCLUSION

We developed a deep learning-based segmentation method to rapidly and accurately segment five pelvic organs simultaneously from daily CBCTs. The proposed method could be used in the clinic to support rapid target and organs-at-risk contouring for prostate adaptive radiation therapy.

Comparison between Conventional IMRT Planning and a Novel Real-Time Adaptive Planning Strategy in Hypofractionated Regimes for Prostate Cancer: A Proof-of-Concept Planning Study

In prostate cancer external beam radiation therapy (EBRT), intra-fraction prostate drifts may compromise the treatment efficacy by underdosing the target and/or overdosing the organs at risk. In this study, a recently developed real-time adaptive planning strategy for intensity-modulated radiation therapy (IMRT) for prostate cancer was evaluated in hypofractionated regimes against traditional treatment planning based on a treatment volume margin expansion. The proposed workflow makes use of a “library of plans” corresponding to possible intra-fraction prostate positions. During delivery, at each beam end, the plan prepared for the position of the prostate closest to the current one is selected and the corresponding beam delivered. This adaptive planning strategy was compared with the traditional approach on a clinical prostate cancer case where different prostate shift magnitudes were considered. Five, six and fifteen fraction hypofractionated schemes were considered for each of these scenarios. When shifts larger than the treatment margin were present, using the traditional approach the seminal vesicles were underdosed by 3–4% of the prescribed dose. The adaptive approach instead allowed for correct target dose coverage and lowered the dose on the rectum for each dosimetric endpoint on average by 3–4% in all the fractionation schemes. Standard intensity-modulated radiation therapy planning did not always guarantee a correct dose distribution on the seminal vesicles and the rectum. The adaptive planning strategy proposed resulted insensitive to the intra-fraction prostate drifts, produced a dose distribution in agreement with the dosimetric requirements in every case analysed and significantly lowered the dose on the rectum.

ESTRO ACROP consensus guideline on the use of image guided radiation therapy for localized prostate cancer

Use of image-guided radiation therapy (IGRT) helps to account for daily prostate position changes during radiation therapy for prostate cancer. However, guidelines for the use of IGRT are scarce. An ESTRO panel consisting of leading radiation oncologists and medical physicists was assembled to review the literature and formulate a consensus guideline of methods and procedure for IGRT in prostate cases. Advanced methods and procedures are also described which the committee judged relevant to further improve clinical practice. Moreover, ranges for margins for the three most popular IGRT scenarios have been suggested as examples.

Properties of the anisotropy of dose contributions: A planning study on prostate cases

Purpose

To characterize the static properties of the anisotropy of dose contributions for different treatment techniques on real patient data (prostate cases). From this, we aim to define a class of treatment techniques with invariant anisotropy distribution carrying information of target coverage and organ‐at‐risk (OAR) sparing. The anisotropy presumably is a helpful quantity for plan adaptation problems.

Methods

The anisotropy field is analyzed for different intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT) techniques for a total of ten planning CTs of prostate cases. Primary irradiation directions ranged from 5 to 15. The uniqueness of anisotropy was explored: In particular, the anisotropy distribution inside the planning treatment volume (PTV) and in its vicinity was investigated. Furthermore, deviations of the anisotropy under beam rotations were explored by direct plan comparison as an indicating the susceptibility of each planned technique to changes in the geometric plan configuration. In addition, plan comparisons enabled the categorization of treatment techniques in terms of their anisotropy distribution.

Results

The anisotropy profile inside the PTV and in the transition between OAR and PTV is independent of the treatment technique as long as a sufficient number of beams contribute to the dose distribution. Techniques with multiple beams constitute a class of almost identical and technique‐independent anisotropy distribution. For this class of techniques, substructures of the anisotropy are particularly pronounced in the PTV, thus offering good options for applying adaptation rules. Additionally, the techniques forming the mentioned class fortunately allow a better OAR sparing at constant PTV coverage. Besides the characterization of the distribution, a pairwise plan comparison reveals each technique's susceptibility to deviations which decreases for an increasing number of primary irradiation directions.

Conclusions

Techniques using many irradiation directions form a class of almost identical anisotropy distributions which are assumed to provide a basis for improved adaptation procedures. Encouragingly, these techniques deliver quite invariant anisotropy distributions with respect to rotations correlated with good plan qualities than techniques using few gantry angles. The following will be the next steps toward anisotropy‐based adaptation: first, the quantification of anisotropy regarding organ deformations; and second, establishing the interrelation between the anisotropy and beam shaping.

[Image-guided radiotherapy contribution and patient setup for anorectal cancer treatment]

Intensity-modulated radiation therapy is recommended in anal squamous cell carcinoma treatment and is increasingly used in rectal cancer. It adapts the dose to target volumes, with a high doses gradient. Intensity-modulated radiation therapy allows to reduce toxicity to critical normal structures and to consider dose-escalation studies or systemic treatment intensification. Image-guided radiation therapy is a warrant of quality for intensity-modulated radiation therapy, especially for successful delivery of the dose as planned. There is no recommended international or national anorectal cancer image-guided radiation therapy protocol currently available. Dose-escalation trials or expert opinions about intensity-modulated/image-guided radiation therapy good practice guidelines recommend daily volumetric imaging throughout the treatment or during the five first fractions and weekly thereafter as a minimum. Image-guided radiation therapy allows to reduce margins related to patient setup errors. Internal margin, related to the internal organ motion, needs to be adapted according to short- or long-course radiotherapy, gender, rectal location; it can be higher than current recommended planning target volume margins, particularly in the upper and anterior part of mesorectum, which has the most significant movement. Image-guided radiation therapy based on volumetric imaging allows to take target volume shrinkage into account and to develop adaptive strategies, in particular for mesorectum shrinkage during rectal cancer treatment. Lastly, the emergence of new image-guided radiation therapy technologies including MRI (which plays a major role in pelvic tumours assessment and diagnosis) opens up interesting perspectives for adaptive radiotherapy, taking into account both organs' movements and tumour shrinkage.

Volumetric-based image guidance is superior to marker-based alignments for stereotactic body radiotherapy of prostate cancer

Purposes

The aim of this study was to evaluate a dual marker‐based and soft‐tissue based image guidance for inter‐fractional corrections in stereotactic body radiotherapy (SBRT) of prostate cancer.

Methods/Materials

We reviewed 18 patients treated with SBRT for prostate cancer. An endorectal balloon was inserted at simulation and each treatment. Planning margins were 3 mm/0 mm posteriorly. Prior to each treatment, a dual image guidance protocol was applied to align three makers using stereoscopic x ray images and then to the soft tissue using kilo‐voltage cone beam CT (kV‐CBCT). After treatment, prostate (CTV), rectal wall, and bladder were delineated on each kV‐CBCT, and delivered dose was recalculated. Dosimetric endpoints were analyzed, including V_{36.25 Gy} for prostate, and D_{0.03 cc} for bladder and rectal wall.

Results

Following initial marker alignment, additional translational shifts were applied to 22 of 84 fractions after kV‐CBCT. Among the 22 fractions, ten fractions exceeded 3 mm shifts in any direction, including one in the left‐right direction, four in the superior‐inferior direction, and five in the anterior‐posterior direction. With and without the additional kV‐CBCT shifts, the average V_{36.25 Gy} of the prostate for the 22 fractions was 97.6 ± 2.6% with the kV x ray image alone, and was 98.1 ± 2.4% after applying the additional kV‐CBCT shifts. The improvement was borderline statistical significance using Wilcoxon signed‐rank test (P = 0.007). D_{0.03 cc} was 45.8 ± 6.3 Gy vs. 45.1 ± 4.9 Gy for the rectal wall; and 49.5 ± 8.6 Gy vs. 49.3 ± 7.9 Gy for the bladder before and after applying kV‐CBCT shifts.

Conclusions

Marker‐based alignment alone is not sufficient. Additional adjustments are needed for some patients based kV‐CBCT.

Magnetic Resonance Imaging-Guided Adaptive Radiation Therapy: A "Game Changer" for Prostate Treatment?

Radiation therapy to the prostate involves increasingly sophisticated delivery techniques and changing fractionation schedules. With a low estimated α/β ratio, a larger dose per fraction would be beneficial, with moderate fractionation schedules rapidly becoming a standard of care. The integration of a magnetic resonance imaging (MRI) scanner and linear accelerator allows for accurate soft tissue tracking with the capacity to replan for the anatomy of the day. Extreme hypofractionation schedules become a possibility using the potentially automated steps of autosegmentation, MRI-only workflow, and real-time adaptive planning. The present report reviews the steps involved in hypofractionated adaptive MRI-guided prostate radiation therapy and addresses the challenges for implementation.

Application programming in C# environment with recorded user software interactions and its application in autopilot of VMAT/IMRT treatment planning

An autopilot scheme of volumetric‐modulated arc therapy (VMAT)/intensity‐modulated radiation therapy (IMRT) planning with the guidance of prior knowledge is established with recorded interactions between a planner and a commercial treatment planning system (TPS). Microsoft (MS) Visual Studio Coded UI is applied to record some common planner‐TPS interactions as subroutines. The TPS used in this study is a Windows‐based Eclipse system. The interactions of our application program with Eclipse TPS are realized through a series of subroutines obtained by prerecording the mouse clicks or keyboard strokes of a planner in operating the TPS. A strategy to autopilot Eclipse VMAT/IMRT plan selection process is developed as a specific example of the proposed “scripting” method. The autopiloted planning is navigated by a decision function constructed with a reference plan that has the same prescription and similar anatomy with the case at hand. The calculation proceeds by alternating between the Eclipse optimization and the outer‐loop optimization independent of the Eclipse. In the C# program, the dosimetric characteristics of a reference treatment plan are used to assess and modify the Eclipse planning parameters and to guide the search for a clinically sensible treatment plan. The approach is applied to plan a head and neck (HN) VMAT case and a prostate IMRT case. Our study demonstrated the feasibility of application programming method in C# environment with recorded interactions of planner‐TPS. The process mimics a planner's planning process and automatically provides clinically sensible treatment plans that would otherwise require a large amount of manual trial and error of a planner. The proposed technique enables us to harness a commercial TPS by application programming via the use of recorded human computer interactions and provides an effective tool to greatly facilitate the treatment planning process.

PACS number(s): 87.55.D‐, 87.55.kd, 87.55.de

Anisotropy of dose contributions-An instrument to upgrade real time IMRT and VMAT adaptation?

PURPOSE

To suggest a definition of dose deposition anisotropy for the purpose of ad hoc adaptation of intensity modulated arc therapy (IMRT) and volumetric arc therapy (VMAT), particularly in the vicinity of important organs at risk (OAR), also for large deformations.

METHODS

Beam's-eye-view (BEV) based fluence warping is a standard adaptation method with disadvantages for strongly varying OAR shapes. 2-Step-adaptation overcomes these difficulties by a deeper analysis of the 3D properties of adaptation processes, but requires separate arcs for every OAR to spare, which makes it impractical for cases with multiple OARs. The authors aim to extend the 2-Step method to arbitrary intensity modulated plan by analyzing the anisotropy of dose contributions. Anisotropy was defined as a second term of Fourier transformation of gantry angle dependent dose contributions. For a cylindrical planning target volume (PTV) surrounding an OAR of varying diameter, the anisotropy and the dose-normalized anisotropy were analyzed for several scenarios of optimized fluence distributions. 2-Step adaptation to decreasing and increasing OAR diameter was performed, and compared to a usual fluence based adaptation method. For two clinical cases, prostate and neck, the VMAT was generated and the behavior of anisotropy was qualitatively explored for deformed organs at risk.

RESULTS

Dose contribution anisotropy in the PTV peaks around nearby OARs. The thickness of the "anisotropy wall" around OAR increases for increasing OAR radius, as also does the width of 2-Step dose saturating fluence peak adjacent to the OAR [K. Bratengeier et al., "A comparison between 2-Step IMRT and conventional IMRT planning," Radiother. Oncol. 84, 298-306 (2007)]. Different optimized beam fluence profiles resulted in comparable radial dependence of normalized anisotropy. As predicted, even for patient cases, anisotropy was inflated even more than increasing diameters of OAR.

CONCLUSIONS

For cylindrically symmetric cases, the dose distribution anisotropy defined in the present work implicitly contains adaptation-relevant information about 3D relationships between PTV and OAR and degree of OAR sparing. For more complex realistic cases, it shows the predicted behavior qualitatively. The authors claim to have found a first component for advancing a 2-Step adaptation to a universal adaptation algorithm based on the BEV projection of the dose anisotropy. Further planning studies to explore the potential of anisotropy for adaptation algorithms using phantoms and clinical cases of differing complexity will follow.

Dosimetric impact of different bladder and rectum filling during prostate cancer radiotherapy

Background

The aim of this study was to analyze the influence of volumetric changes of bladder and rectum filling on the 3D dose distribution in prostate cancer radiotherapy.

Methods

A total of 314 cone-beam CT (CBCT) image data sets from 19 patients were enrolled in this study. For each CBCT, the bladder and rectum were contoured and volume sizes were normalized to those on their original CT. The daily delivered dose was recalculated on the CBCT images and the doses to bladder and rectum were investigated. Linear regression analysis was performed to identify the mean dose change of the volume change using SPSS 19.

Results

The data show that the variances of the normalized volume of the bladder and the rectum are 0.13–0.58 and 0.12–0.50 respectively. The variances of V_70Gy, V_60Gy, V_50Gy, V_40Gy and V_30Gy of bladder are bigger than those of rectum for 17 patients. The linear regression analysis indicates a 10 % increase in bladder volume will cause a 5.6 % (±4.9 %) reduction in mean dose (p <0.05).

Conclusions

The bladder’s volume change is more significant than that of the rectum for the prostate cancer patient. The rectum volume variations are not significant except for air bubbles, which change the shape and the position of the rectum. The bladder volume variations may cause dose changes proportionately. Monitoring the bladder’s volume before fractional treatment delivery will be crucial for accurate dose delivery.

Adaptive radiotherapy strategies for pelvic tumors - a systematic review of clinical implementations

Introdution: Variation in shape, position and treatment response of both tumor and organs at risk are major challenges for accurate dose delivery in radiotherapy. Adaptive radiotherapy (ART) has been proposed to customize the treatment to these motion/response patterns of the individual patients, but increases workload and thereby challenges clinical implementation. This paper reviews strategies and workflows for clinical and in silico implemented ART for prostate, bladder, gynecological (gyne) and ano-rectal cancers.

MATERIAL AND METHODS

Initial identification of papers was based on searches in PubMed. For each tumor site, the identified papers were screened independently by two researches for selection of studies describing all processes of an ART workflow: treatment monitoring and evaluation, decision and execution of adaptations. Both brachytherapy and external beam studies were eligible for review.

RESULTS

The review consisted of 43 clinical studies and 51 in silico studies. For prostate, 1219 patients were treated with offline re-planning, mainly to adapt prostate motion relative to bony anatomy. For gyne 1155 patients were treated with online brachytherapy re-planning while 25 ano-rectal cancer patients were treated with offline re-planning, all to account for tumor regression detected by magnetic resonance imaging (MRI)/computed tomography (CT). For bladder and gyne, 161 and 64 patients, respectively, were treated with library-based online plan selection to account for target volume and shape variations. The studies reported sparing of rectum (prostate and bladder cancer), bladder (ano-rectal cancer) and bowel cavity (gyne and bladder cancer) as compared to non-ART.

CONCLUSION

Implementations of ART were dominated by offline re-planning and online brachytherapy re-planning strategies, although recently online plan selection workflows have increased with the availability of cone-beam CT. Advantageous dosimetric and outcome patterns using ART was documented by the studies of this review. Despite this, clinical implementations were scarce due to challenges in target/organ re-contouring and suboptimal patient selection in the ART workflows.

A deformable head and neck phantom with in-vivo dosimetry for adaptive radiotherapy quality assurance

PURPOSE

Patients' interfractional anatomic changes can compromise the initial treatment plan quality. To overcome this issue, adaptive radiotherapy (ART) has been introduced. Deformable image registration (DIR) is an important tool for ART and several deformable phantoms have been built to evaluate the algorithms' accuracy. However, there is a lack of deformable phantoms that can also provide dosimetric information to verify the accuracy of the whole ART process. The goal of this work is to design and construct a deformable head and neck (HN) ART quality assurance (QA) phantom with in vivo dosimetry.

METHODS

An axial slice of a HN patient is taken as a model for the phantom construction. Six anatomic materials are considered, with HU numbers similar to a real patient. A filled balloon inside the phantom tissue is inserted to simulate tumor. Deflation of the balloon simulates tumor shrinkage. Nonradiopaque surface markers, which do not influence DIR algorithms, provide the deformation ground truth. Fixed and movable holders are built in the phantom to hold a diode for dosimetric measurements.

RESULTS

The measured deformations at the surface marker positions can be compared with deformations calculated by a DIR algorithm to evaluate its accuracy. In this study, the authors selected a Demons algorithm as a DIR algorithm example for demonstration purposes. The average error magnitude is 2.1 mm. The point dose measurements from the in vivo diode dosimeters show a good agreement with the calculated doses from the treatment planning system with a maximum difference of 3.1% of prescription dose, when the treatment plans are delivered to the phantom with original or deformed geometry.

CONCLUSIONS

In this study, the authors have presented the functionality of this deformable HN phantom for testing the accuracy of DIR algorithms and verifying the ART dosimetric accuracy. The authors' experiments demonstrate the feasibility of this phantom serving as an end-to-end ART QA phantom.

Atlas-guided prostate intensity modulated radiation therapy (IMRT) planning

An atlas-based IMRT planning technique for prostate cancer was developed and evaluated. A multi-dose atlas was built based on the anatomy patterns of the patients, more specifically, the percent distance to the prostate and the concaveness angle formed by the seminal vesicles relative to the anterior-posterior axis. A 70-case dataset was classified using a k-medoids clustering analysis to recognize anatomy pattern variations in the dataset. The best classification, defined by the number of classes or medoids, was determined by the largest value of the average silhouette width. Reference plans from each class formed a multi-dose atlas. The atlas-guided planning (AGP) technique started with matching the new case anatomy pattern to one of the reference cases in the atlas; then a deformable registration between the atlas and new case anatomies transferred the dose from the atlas to the new case to guide inverse planning with full automation. 20 additional clinical cases were re-planned to evaluate the AGP technique. Dosimetric properties between AGP and clinical plans were evaluated. The classification analysis determined that the 5-case atlas would best represent anatomy patterns for the patient cohort. AGP took approximately 1 min on average (corresponding to 70 iterations of optimization) for all cases. When dosimetric parameters were compared, the differences between AGP and clinical plans were less than 3.5%, albeit some statistical significances observed: homogeneity index (p  >  0.05), conformity index (p  <  0.01), bladder gEUD (p  <  0.01), and rectum gEUD (p  =  0.02). Atlas-guided treatment planning is feasible and efficient. Atlas predicted dose can effectively guide the optimizer to achieve plan quality comparable to that of clinical plans.

Adaptive optimization by 6 DOF robotic couch in prostate volumetric IMRT treatment: rototranslational shift and dosimetric consequences

The purpose of this study was to investigate the magnitude and dosimetric relevance of translational and rotational shifts on IGRT prostate volumetric‐modulated arc therapy (VMAT) using Protura six degrees of freedom (DOF) Robotic Patient Positioning System. Patients with cT3aN0M0 prostate cancer, treated with VMAT simultaneous integrated boost (VMAT‐SIB), were enrolled. PTV2 was obtained adding 0.7 cm margin to seminal vesicles base (CTV2), while PTV1 adding to prostate (CTV1) 0.7 cm margin in all directions, except 1.2 cm, as caudal margin. A daily CBCT was acquired before dose delivery. The translational and rotational displacements were corrected through Protura Robotic Couch, collected and applied to the simulation CT to obtain a translated CT (tCT) and a rototranslated CT (rtCT) on which we recalculated the initial treatment plan (TP). We analyzed the correlation between dosimetric coverage, organs at risk (OAR) sparing, and translational or rotational displacements. The dosimetric impact of a rototranslational correction was calculated. From October 2012 to September 2013, a total of 263 CBCT scans from 12 patients were collected. Translational shifts were <5mm in 81% of patients and the rotational shifts were <2∘ in 93% of patient scans. The dosimetric analysis was performed on 172 CBCT scans and calculating 344 VMAT‐TP. Two significant linear correlations were observed between yaw and the V20 femoral heads and between pitch rotation and V50 rectum (p<0.001); rototranslational correction seems to impact more on PTV2 than on PTV1, especially when margins are reduced. Rotational errors are of dosimetric significance in sparing OAR and in target coverage. This is relevant for femoral heads and rectum because of major distance from isocenter, and for seminal vesicles because of irregular shape. No correlation was observed between translational and rotational errors. A study considering the intrafractional error and the deformable registration is ongoing.

PACS number: 87.55.de

Adaptive and robust radiation therapy in the presence of drift

Combining adaptive and robust optimization in radiation therapy has the potential to mitigate the negative effects of both intrafraction and interfraction uncertainty over a fractionated treatment course. A previously developed adaptive and robust radiation therapy (ARRT) method for lung cancer was demonstrated to be effective when the sequence of breathing patterns was well-behaved. In this paper, we examine the applicability of the ARRT method to less well-behaved breathing patterns. We develop a novel method to generate sequences of probability mass functions that represent different types of drift in the underlying breathing pattern. Computational results derived from applying the ARRT method to these sequences demonstrate that the ARRT method is effective for a much broader class of breathing patterns than previously demonstrated.

Accelerated gradient-based free form deformable registration for online adaptive radiotherapy

The registration of planning fan-beam computed tomography (FBCT) and daily cone-beam CT (CBCT) is a crucial step in adaptive radiation therapy. The current intensity-based registration algorithms, such as Demons, may fail when they are used to register FBCT and CBCT, because the CT numbers in CBCT cannot exactly correspond to the electron densities. In this paper, we investigated the effects of CBCT intensity inaccuracy on the registration accuracy and developed an accurate gradient-based free form deformation algorithm (GFFD). GFFD distinguishes itself from other free form deformable registration algorithms by (a) measuring the similarity using the 3D gradient vector fields to avoid the effect of inconsistent intensities between the two modalities; (b) accommodating image sampling anisotropy using the local polynomial approximation-intersection of confidence intervals (LPA-ICI) algorithm to ensure a smooth and continuous displacement field; and (c) introducing a 'bi-directional' force along with an adaptive force strength adjustment to accelerate the convergence process. It is expected that such a strategy can decrease the effect of the inconsistent intensities between the two modalities, thus improving the registration accuracy and robustness. Moreover, for clinical application, the algorithm was implemented by graphics processing units (GPU) through OpenCL framework. The registration time of the GFFD algorithm for each set of CT data ranges from 8 to 13 s. The applications of on-line adaptive image-guided radiation therapy, including auto-propagation of contours, aperture-optimization and dose volume histogram (DVH) in the course of radiation therapy were also studied by in-house-developed software.

A DVH-guided IMRT optimization algorithm for automatic treatment planning and adaptive radiotherapy replanning

PURPOSE

To develop a novel algorithm that incorporates prior treatment knowledge into intensity modulated radiation therapy optimization to facilitate automatic treatment planning and adaptive radiotherapy (ART) replanning.

METHODS

The algorithm automatically creates a treatment plan guided by the DVH curves of a reference plan that contains information on the clinician-approved dose-volume trade-offs among different targets/organs and among different portions of a DVH curve for an organ. In ART, the reference plan is the initial plan for the same patient, while for automatic treatment planning the reference plan is selected from a library of clinically approved and delivered plans of previously treated patients with similar medical conditions and geometry. The proposed algorithm employs a voxel-based optimization model and navigates the large voxel-based Pareto surface. The voxel weights are iteratively adjusted to approach a plan that is similar to the reference plan in terms of the DVHs. If the reference plan is feasible but not Pareto optimal, the algorithm generates a Pareto optimal plan with the DVHs better than the reference ones. If the reference plan is too restricting for the new geometry, the algorithm generates a Pareto plan with DVHs close to the reference ones. In both cases, the new plans have similar DVH trade-offs as the reference plans.

RESULTS

The algorithm was tested using three patient cases and found to be able to automatically adjust the voxel-weighting factors in order to generate a Pareto plan with similar DVH trade-offs as the reference plan. The algorithm has also been implemented on a GPU for high efficiency.

CONCLUSIONS

A novel prior-knowledge-based optimization algorithm has been developed that automatically adjust the voxel weights and generate a clinical optimal plan at high efficiency. It is found that the new algorithm can significantly improve the plan quality and planning efficiency in ART replanning and automatic treatment planning.

Metal artifact correction for x-ray computed tomography using kV and selective MV imaging

PURPOSE

The overall goal of this work is to improve the computed tomography (CT) image quality for patients with metal implants or fillings by completing the missing kilovoltage (kV) projection data with selectively acquired megavoltage (MV) data that do not suffer from photon starvation. When both of these imaging systems, which are available on current radiotherapy devices, are used, metal streak artifacts are avoided, and the soft-tissue contrast is restored, even for regions in which the kV data cannot contribute any information.

METHODS

Three image-reconstruction methods, including two filtered back-projection (FBP)-based analytic methods and one iterative method, for combining kV and MV projection data from the two on-board imaging systems of a radiotherapy device are presented in this work. The analytic reconstruction methods modify the MV data based on the information in the projection or image domains and then patch the data onto the kV projections for a FBP reconstruction. In the iterative reconstruction, the authors used dual-energy (DE) penalized weighted least-squares (PWLS) methods to simultaneously combine the kV/MV data and perform the reconstruction.

RESULTS

The authors compared kV/MV reconstructions to kV-only reconstructions using a dental phantom with fillings and a hip-implant numerical phantom. Simulation results indicated that dual-energy sinogram patch FBP and the modified dual-energy PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in the kV projections. The root-mean-square errors of soft-tissue patterns obtained using combined kV/MV data are 10-15 Hounsfield units smaller than those of the kV-only images, and the structural similarity index measure also indicates a 5%-10% improvement in the image quality. The added dose from the MV scan is much less than the dose from the kV scan if a high efficiency MV detector is assumed.

CONCLUSIONS

The authors have shown that it is possible to improve the image quality of kV CTs for patients with metal implants or fillings by completing the missing kV projection data with selectively acquired MV data that do not suffer from photon starvation. Numerical simulations demonstrated that dual-energy sinogram patch FBP and a modified kV/MV PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in kV projections. Combined kV/MV images may permit the improved delineation of structures of interest in CT images for patients with metal implants or fillings.

Investigation into image quality and dose for different patient geometries with multiple cone-beam CT systems

PURPOSE

To provide quantitative and qualitative image quality metrics and imaging dose for modern Varian On-board Imager (OBI) (ver. 1.5) and Elekta X-ray Volume Imager (XVI) (ver. 4.5R) cone-beam computed tomography (CBCT) systems in a clinical adaptive radiation therapy environment by accounting for varying patient thickness.

METHODS

Image quality measurements were acquired with Catphan 504 phantom (nominal diameter and with additional 10 cm thickness) for OBI and XVI systems and compared to planning CT (pCT) (GE LightSpeed). Various clinical protocols were analyzed for the OBI and XVI systems and analyzed using image quality metrics, including spatial resolution, low contrast detectability, uniformity, and HU sensitivity. Imaging dose measurements were acquired in Wellhofer Scanditronix i'mRT phantom at nominal phantom diameter and with additional 4 cm phantom diameter using GafChromic XRQA2 film. Calibration curves were generated using previously published in-air Air Kerma calibration method.

RESULTS

The OBI system full trajectory scans exhibited very little dependence on phantom thickness for accurate HU calculation, while half-trajectory scans with full-fan filter exhibited dependence of HU calculation on phantom thickness. The contrast-to-noise ratio (CNR) for the OBI scans decreased with additional phantom thickness. The uniformity of Head protocol scan was most significantly affected with additional phantom thickness. The spatial resolution and CNR compared favorably with pCT, while the uniformity of the OBI system was slightly inferior to pCT. The OBI scan protocol dose levels for nominal phantom thickness at the central portion of the phantom were 2.61, 0.72, and 1.88 cGy, and for additional phantom thickness were 1.95, 0.48, and 1.52 cGy, for the Pelvis, Thorax, and Spotlight protocols, respectively. The XVI system scans exhibited dependence on phantom thickness for accurate HU calculation regardless of trajectory. The CNR for the XVI scans decreased with additional phantom thickness. The uniformity of the XVI scans was significantly dependent on the selection of the proper FOV setting for all phantom geometries. The spatial resolution, CNR, and uniformity for XVI were lower than values measured for pCT. The XVI scan protocol dose levels at the central portion of the phantom for nominal phantom thickness were 2.14, 2.15, and 0.33 cGy, and for additional phantom thickness were 1.56, 1.68, and 0.21 cGy, for the Pelvis M20, Chest M20, and Prostate Seed S10 scan protocols, respectively.

CONCLUSIONS

The OBI system offered comparable spatial resolution and CNR results to the results for pCT. Full trajectory scans with the OBI system need little-to-no correction for HU calculation based on HU stability with changing phantom thickness. The XVI system offered lower spatial resolution and CNR results than pCT. In addition, the HU calculation for all scan protocols was dependent on the phantom thickness. The uniformity for each CBCT system was inferior to that of pCT for each phantom geometry. The dose for each system and scan protocol in the interior of the phantom tended to decrease by approximately 25% with 4 cm additional phantom thickness.

Reconstructing cone-beam CT with spatially varying qualities for adaptive radiotherapy: a proof-of-principle study

With the aim of maximally reducing imaging dose while meeting requirements for adaptive radiation therapy (ART), we propose in this paper a new cone beam CT (CBCT) acquisition and reconstruction method that delivers images with a low noise level inside a region of interest (ROI) and a relatively high noise level outside the ROI. The acquired projection images include two groups: densely sampled projections at a low exposure with a large field of view (FOV) and sparsely sampled projections at a high exposure with a small FOV corresponding to the ROI. A new algorithm combining the conventional filtered back-projection algorithm and the tight-frame iterative reconstruction algorithm is also designed to reconstruct the CBCT based on these projection data. We have validated our method on a simulated head-and-neck (HN) patient case, a semi-real experiment conducted on a HN cancer patient under a full-fan scan mode, as well as a Catphan phantom under a half-fan scan mode. Relative root-mean-square errors (RRMSEs) of less than 3% for the entire image and ~1% within the ROI compared to the ground truth have been observed. These numbers demonstrate the ability of our proposed method to reconstruct high-quality images inside the ROI. As for the part outside ROI, although the images are relatively noisy, it can still provide sufficient information for radiation dose calculations in ART. Dose distributions calculated on our CBCT image and on a standard CBCT image are in agreement, with a mean relative difference of 0.082% inside the ROI and 0.038% outside the ROI. Compared with the standard clinical CBCT scheme, an imaging dose reduction of approximately 3-6 times inside the ROI was achieved, as well as an 8 times outside the ROI. Regarding computational efficiency, it takes 1-3 min to reconstruct a CBCT image depending on the number of projections used. These results indicate that the proposed method has the potential for application in ART.

Uncertainties encountered in implementation of adaptive planning with in vivo dosimeters

Under a previously approved institutional review board protocol for prostate cancer patients, implanted metal-oxide semiconductor field-effect transistor dosimeters (dose verification system, Sicel Technologies) were used for measurement of the in vivo delivered daily dose. This dosimetric information provided the ability to adapt the plan if the measured doses did not match the dose expected from the planning system. Due to the inherent uncertainty in the dosimeters, the decision to adapt the treatment plan was made only if readings differed by more than 7 % for three consecutive days. To validate this method, we acquired daily cone beam computed tomography images for five patients, and the dose delivered to the dosimeters was calculated by use of (1) an automated procedure (MIM Maestro, MIM Software) and (2) the treatment planning system (XIO, Elekta). 72 % of the doses calculated automatically fell within 1 % of the doses calculated in the planning system, and 99 % agreed within 2 %. When compared to the calculated dose, 53 % of the in vivo measurements fell within 3 % of the calculated dose, 80 % fell within 5 %, and 9.8 % were greater than 7 %, but never on three consecutive days. The measured doses agreed reasonably well with the calculated doses, supporting the decision to adapt the plan only if there were discrepancies of more than 7 % over three consecutive days. Even with the inherent uncertainty in the dosimeters, this adaptive planning method can detect delivery inaccuracies that would not otherwise be caught with the use of only daily image guidance or other dose calculation surrogates.

Applications of linac-mounted kilovoltage Cone-beam Computed Tomography in modern radiation therapy: A review

The use of Cone-beam Computed Tomography (CBCT) in radiotherapy is increasing due to the widespread implementation of kilovoltage systems on the currently available linear accelerators. Cone beam CT acts as an effective Image-Guided Radiotherapy (IGRT) tool for the verification of patient position. It also opens up the possibility of real-time re-optimization of treatment plans for Adaptive Radiotherapy (ART). This paper reviews the most prominent applications of CBCT (linac-mounted) in radiation therapy, focusing on CBCT-based planning and dose calculation studies. This is followed by a concise review of the main issues associated with CBCT, such as imaging artifacts, dose and image quality. It explores how medical physicists and oncologists can best apply CBCT for therapeutic applications.

Artifact-resistant motion estimation with a patient-specific artifact model for motion-compensated cone-beam CT

PURPOSE

In image-guided radiation therapy (IGRT) valuable information for patient positioning, dose verification, and adaptive treatment planning is provided by an additional kV imaging unit. However, due to the limited gantry rotation speed during treatment the typical acquisition time is quite long. Tomographic images of the thorax suffer from motion blurring or, if a gated 4D reconstruction is performed, from significant streak artifacts. Our purpose is to provide a method that reliably estimates respiratory motion in presence of severe artifacts. The estimated motion vector fields are then used for motion-compensated image reconstruction to provide high quality respiratory-correlated 4D volumes for on-board cone-beam CT (CBCT) scans.

METHODS

The proposed motion estimation method consists of a model that explicitly addresses image artifacts because in presence of severe artifacts state-of-the-art registration methods tend to register artifacts rather than anatomy. Our artifact model, e.g., generates streak artifacts very similar to those included in the gated 4D CBCT images, but it does not include respiratory motion. In combination with a registration strategy, the model gives an error estimate that is used to compensate the corresponding errors of the motion vector fields that are estimated from the gated 4D CBCT images. The algorithm is tested in combination with a cyclic registration approach using temporal constraints and with a standard 3D-3D registration approach. A qualitative and quantitative evaluation of the motion-compensated results was performed using simulated rawdata created on basis of clinical CT data. Further evaluation includes patient data which were scanned with an on-board CBCT system.

RESULTS

The model-based motion estimation method is nearly insensitive to image artifacts of gated 4D reconstructions as they are caused by angular undersampling. The motion is accurately estimated and our motion-compensated image reconstruction algorithm can correct for it. Motion artifacts of 3D standard reconstruction are significantly reduced, while almost no new artifacts are introduced.

CONCLUSIONS

Using the artifact model allows to accurately estimate and compensate for patient motion, even if the initial reconstructions are of very low image quality. Using our approach together with a cyclic registration algorithm yields a combination which shows almost no sensitivity to sparse-view artifacts and thus ensures both high spatial and high temporal resolution.

Target volume coverage and dose to organs at risk in prostate cancer patients. Dose calculation on daily cone-beam CT data sets

PURPOSE

On the basis of correct Hounsfield unit to electron density calibration, cone-beam computed tomography (CBCT) data provide the opportunity for retrospective dose recalculation in the patient. Therefore, the consequences of translational positioning corrections and of morphological changes in the patient anatomy can be quantified for prostate cancer patients.

MATERIALS AND METHODS

The organs at risk were newly contoured on the CBCT data sets of 7 patients so as to evaluate the actual applied dose. The daily dose to the planning target volume (PTV) was recalculated with and without the translation data, which result from the real patient repositioning.

RESULTS

A CBCT-based dose recalculation with uncertainties less than 3 % is possible. The deviations between the planning CT and the CBCT without the translational positioning correction vector show an average dose difference of - 8 % inside the PTV. An inverse proportional relation between the mean bladder dose and the actual volume of the bladder could be established. The daily applied dose to the rectum is about 1-54 % higher than predicted by the planning CT.

CONCLUSION

A dose calculation based on CBCT data is possible. The daily positioning correction of the patient is necessary to avoid an underdosage in the PTV. The new contouring of the organs at risk - the bladder and rectum - allows a better appraisal to be made of the total applied dose to these organs.

Automatic treatment plan re-optimization for adaptive radiotherapy guided with the initial plan DVHs

Adaptive radiation therapy (ART) can reduce normal tissue toxicity and/or improve tumor control through treatment adaptations based on the current patient anatomy. Developing an efficient and effective re-planning algorithm is an important step toward the clinical realization of ART. For the re-planning process, manual trial-and-error approach to fine-tune planning parameters is time-consuming and is usually considered unpractical, especially for online ART. It is desirable to automate this step to yield a plan of acceptable quality with minimal interventions. In ART, prior information in the original plan is available, such as dose-volume histogram (DVH), which can be employed to facilitate the automatic re-planning process. The goal of this work is to develop an automatic re-planning algorithm to generate a plan with similar, or possibly better, DVH curves compared with the clinically delivered original plan. Specifically, our algorithm iterates the following two loops. An inner loop is the traditional fluence map optimization, in which we optimize a quadratic objective function penalizing the deviation of the dose received by each voxel from its prescribed or threshold dose with a set of fixed voxel weighting factors. In outer loop, the voxel weighting factors in the objective function are adjusted according to the deviation of the current DVH curves from those in the original plan. The process is repeated until the DVH curves are acceptable or maximum iteration step is reached. The whole algorithm is implemented on GPU for high efficiency. The feasibility of our algorithm has been demonstrated with three head-and-neck cancer IMRT cases, each having an initial planning CT scan and another treatment CT scan acquired in the middle of treatment course. Compared with the DVH curves in the original plan, the DVH curves in the resulting plan using our algorithm with 30 iterations are better for almost all structures. The re-optimization process takes about 30 s using our in-house optimization engine.