Magnetic resonance imaging radiomic features stability in brain metastases: Impact of image preprocessing, image-, and feature-level harmonization.

Background and purpose: Magnetic resonance imaging (MRI) scans are highly sensitive to acquisition and reconstruction parameters which affect feature stability and model generalizability in radiomic research. This work aims to investigate the effect of image pre-processing and harmonization methods on the stability of brain MRI radiomic features and the prediction performance of radiomic models in patients with brain metastases (BMs). Materials and methods: Two T1 contrast enhanced brain MRI data-sets were used in this study. The first contained 25 BMs patients with scans at two different time points and was used for features stability analysis. The effect of gray level discretization (GLD), intensity normalization (Z-score, Nyul, WhiteStripe, and in house-developed method named N-Peaks), and ComBat harmonization on features stability was investigated and features with intraclass correlation coefficient >0.8 were considered as stable. The second data-set containing 64 BMs patients was used for a classification task to investigate the informativeness of stable features and the effects of harmonization methods on radiomic model performance. Results: Applying fixed bin number (FBN) GLD, resulted in higher number of stable features compare to fixed bin size (FBS) discretization (10 ± 5.5 % higher). `Harmonization in feature domain improved the stability for non-normalized and normalized images with Z-score and WhiteStripe methods. For the classification task, keeping the stable features resulted in good performance only for normalized images with N-Peaks along with FBS discretization. Conclusions: To develop a robust MRI based radiomic model we recommend using an intensity normalization method based on a reference tissue (e.g N-Peaks) and then using FBS discretization.

Validation of echo planar imaging based diffusion-weighted magnetic resonance imaging on a 0.35 T MR-Linac.

Background and Purpose: The feasibility of acquiring diffusion-weighted imaging (DWI) images on an MR-Linac for quantitative response assessment during radiotherapy was explored. DWI data obtained with a Spin Echo Echo Planar Imaging sequence adapted for a 0.35 T MR-Linac were examined and compared with DWI data from a conventional 3 T scanner. Materials and Methods: Apparent diffusion coefficient (ADC) measurements and a distortion correction technique were investigated using DWI-calibrated phantoms and in the brains of seven volunteers. All DWI utilized two phase-encoding directions for distortion correction and off-resonance field estimation. ADC maps in the brain were analyzed for automatically segmented normal tissues. Results: Phantom ADC measurements on the MR-Linac were within a 3 % margin of those recorded by the 3 T scanner. The maximum distortion observed in the phantom was 2.0 mm prior to correction and 1.1 mm post-correction on the MR-Linac, compared to 6.0 mm before correction and 3.6 mm after correction at 3 T. In vivo, the average ADC values for gray and white matter exhibited variations of 14 % and 4 %, respectively, for different selections of b-values on the MR-Linac. Distortions in brain images before correction, estimated through the off-resonance field, reached 2.7 mm on the MR-Linac and 12 mm at 3 T. Conclusion: Accurate ADC measurements are achievable on a 0.35 T MR-Linac, both in phantom and in vivo. The selection of b-values significantly influences ADC values in vivo. DWI on the MR-Linac demonstrated lower distortion levels, with a maximum distortion reduced to 1.1 mm after correction.

Patient-specific quality assurance strategies for synthetic computed tomography in magnetic resonance-only radiotherapy of the abdomen.

Background and purpose: The superior tissue contrast of magnetic resonance (MR) compared to computed tomography (CT) led to an increasing interest towards MR-only radiotherapy. For the latter, the dose calculation should be performed on a synthetic CT (sCT). Patient-specific quality assurance (PSQA) methods have not been established yet and this study aimed to assess several software-based solutions. Materials and methods: A retrospective study was performed on 20 patients treated at an MR-Linac, which were selected to evenly cover four subcategories: (i) standard, (ii) air pockets, (iii) lung and (iv) implant cases. The neural network (NN) CycleGAN was adopted to generate a reference sCT, which was then compared to four PSQA methods: (A) water override of body, (B) five tissue classes with bulk densities, (C) sCT generated by a separate NN (pix2pix) and (D) deformed CT. Results: The evaluation of the dose endpoints demonstrated that while all methods A-D provided statistically equivalent results (p = 0.05) within the 2% level for the standard cases (i), only the methods C-D guaranteed the same result over the whole cohort. The bulk densities override was shown to be a valuable method in absence of lung tissue within the beam path. Conclusion: The observations of this study suggested that the use of an additional sCT generated by a separate NN was an appropriate tool to perform PSQA of a sCT in an MR-only workflow at an MR-Linac. The time and dose endpoints requirements were respected, namely within 10 min and 2%.

Impact of the gradient in gantry-table rotation on dynamic trajectory radiotherapy plan quality.

BACKGROUND: To improve organ at risk (OAR) sparing, dynamic trajectory radiotherapy (DTRT) extends VMAT by dynamic table and collimator rotation during beam-on. However, comprehensive investigations regarding the impact of the gantry-table (GT) rotation gradient on the DTRT plan quality have not been conducted. PURPOSE: To investigate the impact of a user-defined GT rotation gradient on plan quality of DTRT plans in terms of dosimetric plan quality, dosimetric robustness, deliverability, and delivery time. METHODS: The dynamic trajectories of DTRT are described by GT and gantry-collimator paths. The GT path is determined by minimizing the overlap of OARs with planning target volume (PTV). This approach is extended to consider a GT rotation gradient by means of a maximum gradient of the path ( G m a x ${G}_{max}$ ) between two adjacent control points ( G = | Δ table angle / Δ gantry angle | $G = | \Delta {{\mathrm{table\ angle}}/\Delta {\mathrm{gantry\ angle}}} |$ ) and maximum absolute change of G ( Δ G m a x ${{\Delta}}{G}_{max}$ ). Four DTRT plans are created with different maximum G&∆G: G m a x ${G}_{max}$ & Δ G m a x ${{\Delta}}{G}_{max}$  = 0.5&0.125 (DTRT-1), 1&0.125 (DTRT-2), 3&0.125 (DTRT-3) and 3&1|(DTRT-4), including 3-4 dynamic trajectories, for three clinically motivated cases in the head and neck and brain region (A, B, and C). A reference VMAT plan for each case is created. For all plans, plan quality is assessed and compared. Dosimetric plan quality is evaluated by target coverage, conformity, and OAR sparing. Dosimetric robustness is evaluated against systematic and random patient-setup uncertainties between ± 3 mm $ \pm 3\ {\mathrm{mm}}$ in the lateral, longitudinal, and vertical directions, and machine uncertainties between ± 4 ∘ $ \pm 4^\circ \ $ in the dynamically rotating machine components (gantry, table, collimator rotation). Delivery time is recorded. Deliverability and delivery accuracy on a TrueBeam are assessed by logfile analysis for all plans and additionally verified by film measurements for one case. All dose calculations are Monte Carlo based. RESULTS: The extension of the DTRT planning process with user-defined G m a x & Δ G m a x ${G}_{max}\& {{\Delta}}{G}_{max}$ to investigate the impact of the GT rotation gradient on plan quality is successfully demonstrated. With increasing G m a x & Δ G m a x ${G}_{max}\& {{\Delta}}{G}_{max}$ , slight (case C, D m e a n , p a r o t i d l . ${D}_{mean,\ parotid\ l.}$ : up to|-1|Gy) and substantial (case A, D 0.03 c m 3 , o p t i c n e r v e r . ${D}_{0.03c{m}^3,\ optic\ nerve\ r.}$ : up to -9.3 Gy, case|B, D m e a n , b r a i n $\ {D}_{mean,\ brain}$ : up to -4.7|Gy) improvements in OAR sparing are observed compared to VMAT, while maintaining similar target coverage. All plans are delivered on the TrueBeam. Expected and actual machine position values recorded in the logfiles deviated by <0.2° for gantry, table and collimator rotation. The film measurements agreed by >96% (2%|global/2 mm Gamma passing rate) with the dose calculation. With increasing G m a x & Δ G m a x ${G}_{max}\& {{\Delta}}{G}_{max}$ , delivery time is prolonged by <2 min/trajectory (DTRT-4) compared to VMAT and DTRT-1. The DTRT plans for case A and B and the VMAT plan for case C plan reveal the best dosimetric robustness for the considered uncertainties. CONCLUSION: The impact of the GT rotation gradient on DTRT plan quality is comprehensively investigated for three cases in the head and neck and brain region. Increasing freedom in this gradient improves dosimetric plan quality at the cost of increased delivery time for the investigated cases. No clear dependency of GT rotation gradient on dosimetric robustness is observed.

Synthetic computed tomography for low-field magnetic resonance-guided radiotherapy in the abdomen.

Background and purpose: The requirement of computed tomography (CT) for radiotherapy planning may be bypassed by synthetic CT (sCT) generated from magnetic resonance (MR), which has recently led to the clinical introduction of MR-only radiotherapy for specific sites. Further developments are required for abdominal sCT, mostly due to the presence of mobile air pockets affecting the dose calculation. In this study we aimed to overcome this limitation for abdominal sCT at a low field (0.35 T) hybrid MR-Linac. Materials and methods: A retrospective analysis was conducted enrolling 168 patients corresponding to 215 MR-CT pairs. After the exclusion criteria, 152 volumetric images were used to train the cycle-consistent generative adversarial network (CycleGAN) and 34 to test the sCT. Image similarity metrics and dose recalculation analysis were performed. Results: The generated sCT faithfully reproduced the original CT and the location of the air pockets agreed with the MR scan. The dose calculation did not require manual bulk density overrides and the mean deviations of the dose-volume histogram dosimetric points were within 1 % of the CT, without any outlier above 2 %. The mean gamma passing rates were above 99 % for the 2 %/ 2 mm analysis and no cases below 95 % were observed. Conclusions: This study presented the implementation of CycleGAN to perform sCT generation in the abdominal region for a low field hybrid MR-Linac. The sCT was shown to correctly allocate the electron density for the mobile air pockets and the dosimetric analysis demonstrated the potential for future implementation of MR-only radiotherapy in the abdomen.