Geometric and dosimetric evaluation of the differences between rigid and deformable registration to assess interfraction motion during pelvic radiotherapy

Improving hybrid image and structure-based deformable image registration for large internal deformations

Objective.Deformable image registration (DIR) is a widely used technique in radiotherapy. Complex deformations, resulting from large anatomical changes, are a regular challenge. DIR algorithms generally seek a balance between capturing large deformations and preserving a smooth deformation vector field (DVF). We propose a novel structure-based term that can enhance the registration efficacy while ensuring a smooth DVF.Approach.The proposed novel similarity metric for controlling structures was introduced as a new term into a commercially available algorithm. Its performance was compared to the original algorithm using a dataset of 46 patients who received pelvic re-irradiation, many of which exhibited complex deformations.Main results.The mean Dice Similarity Coefficient (DSC) under the improved algorithm was 0.96, 0.94, 0.76, and 0.91 for bladder, rectum, colon, and bone respectively, compared to 0.69, 0.89, 0.62, and 0.88 for the original algorithm. The improvement was more pronounced for complex deformations.Significance.With this work, we have demonstrated that the proposed term is able to improve registration accuracy for complex cases while maintaining realistic deformations.

Evaluation of a cone-beam computed tomography system calibrated for accurate radiotherapy dose calculation

Background and purpose

Dose calculation on cone-beam computed tomography (CBCT) images has been less accurate than on computed tomography (CT) images due to lower image quality and discrepancies in CT numbers for CBCT. As increasing interest arises in offline and online re-planning, dose calculation accuracy was evaluated for a novel CBCT imager integrated into a ring gantry treatment machine.

Materials and methods

The new CBCT system allowed fast image acquisition (5.9 s) by using new hardware, including a large-size flat panel detector, and incorporated image-processing algorithms with iterative reconstruction techniques, leading to accurate CT numbers allowing dose calculation. In this study, CBCT- and CT-based dose calculations were compared based on three anthropomorphic phantoms, after CBCT-to-mass-density calibration was performed. Six plans were created on the CT scans covering various target locations and complexities, followed by CBCT to CT registrations, copying of contours, and re-calculation of the plans on the CBCT scans. Dose-volume histogram metrics for target volumes and organs-at-risk (OARs) were evaluated, and global gamma analyses were performed.

Results

Target coverage differences were consistently below 1.2 %, demonstrating the agreement between CT and re-calculated CBCT dose distributions. Differences in Dmean for OARs were below 0.5 Gy for all plans, except for three OARs, which were below 0.8 Gy (<1.1 %). All plans had a 3 %/1mm gamma pass rate > 97 %.

Conclusions

This study demonstrated comparable results between dose calculations performed on CBCT and CT acquisitions. The new CBCT system with enhanced image quality and CT number accuracy opens possibilities for off-line and on-line re-planning.

Joint radial trajectory correction for accelerated T2 * mapping on an MR-Linac

BACKGROUND

T2 * mapping can characterize tumor hypoxia, which may be associated with resistance to therapy. Acquiring T2 * maps during MR-guided radiotherapy could inform treatment adaptation by, for example, escalating the dose to resistant sub-volumes.

PURPOSE

The purpose of this work is to demonstrate the feasibility of the accelerated T2 * mapping technique using model-based image reconstruction with integrated trajectory auto-correction (TrACR) for MR-guided radiotherapy on an MR-Linear accelerator (MR-Linac).

MATERIALS AND METHODS

The proposed method was validated in a numerical phantom, where two T2 * mapping approaches (sequential and joint) were compared for different noise levels (0,0.1,0.5,1) and gradient delays ([1, -1] and [1, -2] in units of dwell time for x- and y-axis, respectively). Fully sampled k-space was retrospectively undersampled using two different undersampling patterns. Root mean square errors (RMSEs) were calculated between reconstructed T2 * maps and ground truth. In vivo data was acquired twice weekly in one prostate and one head and neck cancer patient undergoing treatment on a 1.5 T MR-Linac. Data were retrospectively undersampled and T2 * maps reconstructed, with and without trajectory corrections were compared.

RESULTS

Numerical simulations demonstrated that, for all noise levels, T2 * maps reconstructed with a joint approach demonstrated less error compared to an uncorrected and sequential approach. For a noise level of 0.1, uniform undersampling and gradient delay [1, -1] (in units of dwell time for x- and y-axis, respectively), RMSEs for sequential and joint approaches were 13.01 and 9.32 ms, respectively, which reduced to 10.92 and 5.89 ms for a gradient delay of [1, 2]. Similarly, for alternate undersampling and gradient delay [1, -1], RMSEs for sequential and joint approaches were 9.80 and 8.90 ms, respectively, which reduced to 9.10 and 5.40 ms for gradient delay [1, 2]. For in vivo data, T2 * maps reconstructed with our proposed approach resulted in less artifacts and improved visual appearance compared to the uncorrected approach. For both prostate and head and neck cancer patients, T2 * maps reconstructed from different treatment fractions showed changes within the planning target volume (PTV).

CONCLUSION

Using the proposed approach, a retrospective data-driven gradient delay correction can be performed, which is particularly relevant for hybrid devices, where full information on the machine configuration is not available for image reconstruction. T2 * maps were acquired in under 5 min and can be integrated into MR-guided radiotherapy treatment workflows, which minimizes patient burden and leaves time for additional imaging for online adaptive radiotherapy on an MR-Linac.

Evaluation of DIR algorithm performance in real patients for radiotherapy treatments: A systematic review of operator-dependent strategies

PURPOSE

The performance of deformable medical image registration (DIR) algorithms has become a major concern.

METHODS

We aimed to obtain updated information on DIR algorithm performance quantification through a literature review of articles published between 2010 and 2022. We focused only on studies using operator-based methods to treat real patients. The PubMed, Google Scholar and Embase databases were searched following PRISMA guidelines.

RESULTS

One hundred and seven articles were identified. The mean number of patients and registrations per publication was 20 and 63, respectively. We found 23 different geometric metrics appearing at least twice, and the dosimetric impact of DIR was quantified in 32 articles. Forty-eight different at-risk organs were described, and target volumes were studied in 43 publications. Prostate, head-and-neck and thoracic locations represented more than ¾ of the studied locations. We summarized the type of DIR and the images used, and other key elements. Intra/interobserver variability, threshold values and the correlation between metrics were also discussed.

CONCLUSIONS

This literature review covers the past decade and should facilitate the implementation of DIR algorithms in clinical practice by providing practical and pertinent information to quantify their performance on real patients.

[Material Investigation for the Development of Non-rigid Phantoms for CT-MRI Image Registration]

PURPOSE

In radiotherapy, deformable image registration (DIR) has been frequently used in different imaging examinations in recent years. However, no phantom has been established for quality assurance for DIR. In order to develop a non-rigid phantom for accuracy control between CT and MRI images, we investigated the suitability of 3D printing materials and gel materials in this study.

METHODS

We measured CT values, T₁ values, T₂ values, and the proton densities of 31 3D printer materials-purchased from three manufacturers-and one gel material. The dice coefficient after DIR was calculated for the CT-MRI images using a prototype phantom made of a gel material compatible with CT-MRI.

RESULTS

The CT number of the 3D printing materials ranged from -6.8 to 146.4 HU. On MRI, T₁ values were not measurable in most cases, whereas T₂ values were not measurable in all cases; proton density (PD) ranged from 2.51% to 4.9%. The gel material had a CT number of 111.16 HU, T₁ value of 813.65 ms, and T₂ value of 27.19 ms. The prototype phantom was flexible, and the usefulness of DIR with CT and MRI images was demonstrated using this phantom.

CONCLUSION

The CT number and T₁ and T₂ values of the gel material are close to those of the human body and may therefore be developed as a DIR verification phantom between CT and MRI. These findings may contribute to the development of non-rigid phantoms for DIR in the future.

An endovaginal MRI array with a forward-looking coil for advanced gynecological cancer brachytherapy procedures: Design and initial results

PURPOSE

To develop an endovaginal MRI array that provides signal enhancement forward into the posterior parametrium and sideways into the vaginal wall, accelerating multiple-contrast detection of residual tumors that survive external beam radiation. The array's enclosure should form an obturator for cervical cancer brachytherapy, allowing integration with MRI-guided catheter placement, CT, and interstitial radiation dose delivery.

METHODS

The endovaginal array consisted of forward-looking and sideways-looking components. The forward-looking element imaged the cervix and posterior endometrium, and the sideways-looking elements imaged the vaginal wall. Electromagnetic simulation was performed to optimize the geometry of a forward-looking coil placed on a conductive-metallic substrate, extending the forward penetration above the coil's tip. Thereafter, an endovaginal array with one forward-looking coil and four sideways-looking elements was constructed and tested at 1.5 Tesla in saline and gel phantoms, and three sexually mature swine. Each coil's tuning, matching, and decoupling were optimized theoretically, implemented with electronic circuits, and validated with network-analyzer measurements. The array enclosure emulates a conventional brachytherapy obturator, allowing use of the internal imaging array together with tandem coils and interstitial catheters, as well as use of the enclosure alone during CT and radiation delivery. To evaluate the receive magnetic field ( B 1 - ) spatial profile, the endovaginal array's specific absorption-rate (SAR) distribution was simulated inside a gel ASTM phantom to determine extreme heating locations in advance of a heating test. Heating tests were then performed during high SAR imaging in a gel phantom at the predetermined locations, testing compliance with MRI safety standards. To assess array imaging performance, signal-to-noise-ratios (SNR) were calculated in a saline phantom and in vivo. Swine images were acquired with the endovaginal array combined with the scanner's body and spine arrays.

RESULTS

Simulated B 1 - profiles for the forward-looking lobe pattern, obtained while varying several geometric parameters, disclosed that a forward-looking coil placed on a metal-backed substrate could double the effective forward penetration from approximately 25 to ∼40 mm. An endovaginal array, enclosed in an obturator enclosure was then constructed, with all coils tuned, matched, and decoupled. The ASTM gel-phantom SAR test showed that peak local SAR was 1.2 W/kg in the forward-looking coil and 0.3 W/kg in the sideways-looking elements, well within ASTM/FDA/IEC guidelines. A 15-min 4 W/kg average SAR imaging experiment resulted in less than 2^o C temperature increase, also within ASTM/FDA/IEC heating limits. In a saline phantom, the forward-looking coil and sideways-looking array's SNR was four to eight times, over a 20-30 mm field-of-view (FOV), and five to eight times, over a 15-25 mm FOV, relative to the spine array's SNR, respectively. In three sexually mature swine, the forward-looking coil provided a 5 + 0.2 SNR enhancement factor within the cervix and posterior endometrium, and the sideways-looking array provided a 4 + 0.2 SNR gain factor in the vaginal wall, relative to the Siemens spine array, demonstrating that the array could significantly reduce imaging time.

CONCLUSIONS

Higher SNR gynecological imaging is supported by forward-looking and sideways-looking coils. A forward-looking endovaginal coil for cervix and parametrium imaging was built with optimized metal backing. Array placement within an obturator enhanced integration with the brachytherapy procedure and accelerated imaging for detecting postexternal-beam residual tumors.

Interfraction dose deviation and catheter position in cervical interstitial and intracavitary image guided HDR brachytherapy

Interstitial and intracavitary gynecological HDR brachytherapy involve precise, localized delivery to targets with high dose gradients, sparing adjacent organs at risk (OAR). Due to the proximity of the rectum, bowel and bladder to the target, deviations in the applicator or catheter with respect to patient anatomy can significantly increase dose to OAR. The magnitude and direction of applicator and catheter migration at each fraction was assessed for template interstitial and tandem and ring (T&R) cohorts. The cohort included twelve gynecological patients with intact cervical lesions treated with external beam and brachytherapy. Pre-treatment CT images were registered to the simulation CT with respect to the target. Treatment catheter positions transformed into the planning CT coordinate system to evaluate localized catheter displacement and dose distributions calculated at each fraction. Dose was evaluated on the planning CT with planning contours and dwell locations at treatment position. Absolute deviation, depth and deflection angle for all patients were 4.6 ± 4.2 mm, -1.4 ± 4.0 mm, and 3.1 ± 2.3° respectively (n = 516 catheter positions for all treatment fractions and patients, mean ± SD). Absolute catheter deviation and deflection magnitude for interstitial treatments increased overall with each subsequent fraction with an overall increase of catheter retraction at each fraction (p < 0.005, n = 492 catheters, Kruskal-Wallis). A target EQD2 D90 reduction of 10 ± 10% and 7.7 ± 8.7% of the planned dose for interstitial and T&R cohorts respectively. There was an overall increase in bladder and rectal doses at each fraction. Catheter tracking in interstitial and intracavitary gynecological treatments with CT imaging revealed significant changes in catheter positioning with respect to the target volume. Overall deviations increased in magnitude with each subsequent fraction in the interstitial treatments. This caused patient dosimetry deviations, including target dose reduction and adjacent OAR doses changes.

Impact of interfractional target motion in locally advanced cervical cancer patients treated with spot scanning proton therapy using an internal target volume strategy

Background and purpose

The more localized dose deposition of proton therapy (PT) compared to photon therapy might allow a reduction in treatment-related side effects but induces additional challenges to address. The aim of this study was to evaluate the impact of interfractional motion on the target and organs at risk (OARs) in cervical cancer patients treated with spot scanning PT using an internal target volume (ITV) strategy.

Methods and materials

For ten locally advanced cervical cancer patients, empty and full bladder planning computed tomography (pCT) as well as 25 daily cone beam CTs (CBCTs) were available. The Clinical Target Volume (CTV), the High Risk CTV (CTV_HR) (gross tumor volume and whole cervix), the non-involved uterus as well as the OARs (bowel, bladder and rectum) were contoured on the daily CBCTs and transferred to the pCT through rigid bony match. Using synthetic CTs derived from pCTs, four-beam spot scanning PT plans were generated to target the patient-specific ITV with 45 Gy(RBE) in 25 fractions. This structure was defined based on pre-treatment MRI and CT to anticipate potential target motion throughout the treatment. D98% of the targets and V40Gy(RBE) of the OARs were extracted from the daily anatomies, accumulated and analyzed. In addition, the impact of bladder volume deviations from planning values on target and bowel dose was investigated.

Results

The ITV strategy ensured a total accumulated dose >42.75 Gy(RBE) to the CTV_HR for all ten patients. Two patients with large bladder-related uterus motion had accumulated dose to the non-involved uterus of 35.7 Gy(RBE) and 41.1 Gy(RBE). Variations in bowel V40Gy(RBE) were found to be correlated (Pearson r = −0.55; p-value <0.0001) with changes in bladder volume during treatment.

Conclusion

The ITV concept ensured adequate dose to the CTV_HR, but was insufficient for the non-involved uterus of patients subject to large target interfractional motion. CBCT monitoring and occasional replanning is recommended along the same lines as with photon radiotherapy in cervical cancer.

Dosimetric impact of organ at risk daily variation during prostate stereotactic ablative radiotherapy

OBJECTIVE

Prostate stereotactic ablative radiotherapy (SABR) delivers large doses using a fast dose rate. This amplifies the effect geometric uncertainties have on normal tissue dose. The aim of this study was to determine whether the treatment dose-volume histogram (DVH) agrees with the planned dose to organs at risk (OAR).

METHODS

41 low-intermediate risk prostate cancer patients were treated with SABR using a linac based technique. Dose prescribed was 35 Gy in five fractions delivered on alternate days, planned using volumetric modulated arc therapy (VMAT) with 10X flattening filter free (FFF). On treatment, prostate was matched to fiducial markers on cone beam CT (CBCT). OAR were retrospectively delineated on 205 pre-treatment CBCT images. Daily CBCT contours were overlaid on the planning CT for dosimetric analysis. Verification plan used to evaluate the daily DVH for each structure. The daily doses received by OAR were recorded using the D%.

RESULTS

The median rectum and bladder volumes at planning were 67.1 cm3 (interquartile range 56.4-78.2) and 164.4 cm3 (interquartile range 120.3-213.4) respectively. There was no statistically significant difference in median rectal volume at each of the five treatment scans compared to the planning scan (p = 0.99). This was also the case for median bladder volume (p = 0.79). The median dose received by rectum and bladder at each fraction was higher than planned, at the majority of dose levels. For rectum the increase ranged from 0.78-1.64Gy and for bladder 0.14-1.07Gy. The percentage of patients failing for rectum D35% < 18 Gy (p = 0.016), D10% < 28 Gy (p = 0.004), D5% < 32 Gy (p = 0.0001), D1% < 35 Gy (p = 0.0001) and bladder D1% < 35 Gy (p = 0.001) at treatment were all statistically significant.

CONCLUSION

In this cohort of prostate SABR patients, we estimate the OAR treatment DVH was higher than planned. This was due to rectal and bladder organ variation.

ADVANCES IN KNOWLEDGE

OAR variation in prostate SABR using a FFF technique, may cause the treatment DVH to be higher than planned.