Skull phantom-based methodology to validate MRI co-registration accuracy for Gamma Knife radiosurgery

PURPOSE

Target localization, for stereotactic radiosurgery (SRS) treatment with Gamma Knife, has become increasingly reliant on the co-registration between the planning MRI and the stereotactic cone-beam computed tomography (CBCT). Validating image registration between modalities would be particularly beneficial when considering the emergence of novel functional and metabolic MRI pulse sequences for target delineation. This study aimed to develop a phantom-based methodology to quantitatively compare the co-registration accuracy of the standard clinical imaging protocol to a representative MRI sequence that was likely to fail co-registration. The comparative methodology presented in this study may serve as a useful tool to evaluate the clinical translatability of novel MRI sequences.

METHODS

A realistic human skull phantom with fiducial marker columns was designed and manufactured to fit into a typical MRI head coil and the Gamma Knife patient positioning system. A series of "optimized" 3D MRI sequences-T₁ -weighted Dixon, T₁ -weighted fast field echo (FFE), and T₂ -weighted fluid-attenuated inversion recovery (FLAIR)-were acquired and co-registered to the CBCT. The same sequences were "compromised" by reconstructing without geometric distortion correction and re-collecting with lower signal-to-noise-ratio (SNR) to simulate a novel MRI sequence with poor co-registration accuracy. Image similarity metrics-structural similarity (SSIM) index, mean squared error (MSE), and peak SNR (PSNR)-were used to quantitatively compare the co-registration of the optimized and compromised MR images.

RESULTS

The ground truth fiducial positions were compared to positions measured from each optimized image volume revealing a maximum median geometric uncertainty of 0.39 mm (LR), 0.92 mm (AP), and 0.13 mm (SI) between the CT and CBCT, 0.60 mm (LR), 0.36 mm (AP), and 0.07 mm (SI) between the CT and T₁ -weighted Dixon, 0.42 mm (LR), 0.23 mm (AP), and 0.08 mm (SI) between the CT and T₁ -weighted FFE, and 0.45 mm (LR), 0.19 mm (AP), and 1.04 mm (SI) between the CT and T₂ -weighted FLAIR. Qualitatively, pairs of optimized and compromised image slices were compared using a fusion image where separable colors were used to differentiate between images. Quantitatively, MSE was the most predictive and SSIM the second most predictive metric for evaluating co-registration similarity. A clinically relevant threshold of MSE, SSIM, and/or PSNR may be defined beyond which point an MRI sequence should be rejected for target delineation based on its dissimilarity to an optimized sequence co-registration. All dissimilarity thresholds calculated using correlation coefficients with in-plane geometric uncertainty would need to be defined on a sequence-by-sequence basis and validated with patient data.

CONCLUSION

This study utilized a realistic skull phantom and image similarity metrics to develop a methodology capable of quantitatively assessing whether a modern research-based MRI sequence can be co-registered to the Gamma Knife CBCT with equal or less than equal accuracy when compared to a clinically accepted protocol.

Chemical exchange saturation transfer MRI in central nervous system tumours on a 1.5 T MR-Linac

PURPOSE

To describe the implementation and initial results of using Chemical Exchange Saturation Transfer (CEST) for monitoring patients with central nervous system (CNS) tumours treated using a 1.5 tesla MR-guided radiotherapy system.

METHODS

CNS patients were treated with up to 30 fractions (total dose up to 60 Gy) using a 1.5 T Elekta Unity MR-Linac. CEST scans were obtained in 54 subjects at one or more time points during treatment. CEST metrics, including the amide magnetization transfer ratio (MTR_Amide), nuclear Overhauser effect (NOE) MTR (MTR_NOE) and asymmetry, were quantified in phantoms and CNS patients. The signal was investigated between tumour and white matter, across time, and across disease categories including high- and low-grade tumours.

RESULTS

The gross tumour volume (GTV) exhibited lower MTR_Amide and MTR_NOE and higher asymmetry compared to contralateral normal appearing white matter. Signal changes in the GTV during fractionated radiotherapy were observed. There were differences between high- and low-grade tumours, with higher CEST asymmetry associated with higher grade disease.

CONCLUSION

CEST MRI using a 1.5 T MR-Linac was demonstrated to be feasible for in vivo imaging of CNS tumours. CEST images showed tumour/white-matter contrast, temporal CEST signal changes, and associations with tumour grade. These results show promise for the eventual goal of using metabolic imaging to inform the design of adaptive radiotherapy protocols.

A strategy to prevent a temperature-induced MRI artifact in warm liquid phantoms due to convection currents

MRI phantom studies often fail to mimic the temperature of the human body, which can negatively impact accuracy. An artifact induced by increasing temperature in liquid phantoms was observed, presenting a significant challenge to temperature-controlled experiments. In this study we characterize and provide a solution to eliminate this temperature-induced MRI artifact. Low concentration (0.5-2.5 mM) agar phantoms were prepared. Utilizing a temperature-controlled phantom holder, T₁ - and T₂ -weighted structural images were acquired at 7 T along with quantitative B₀ , B₁ , T₁ , T₂ and ADC maps at both 25 and 37°C. Additionally, computer simulations were conducted to demonstrate the fluid flow and thermal flux patterns in water to provide an insight into the origins of the artifact. Evidence from computer simulation and quantitative MRI strongly suggest the artifact was caused by heat transfer in the form of natural convection leading to structured patterns of signal loss in MR images. The artifact was present up to agar concentrations of 1.5 mM (T₁ = 3068 ± 16 ms, T₂ = 1052 ± 20 ms, ADC = 2.29 ± 0.36 × 10^-3 mm² /s at 25°C; T₁ = 3928 ± 44 ms, T₂ = 1122 ± 24 ms, ADC = 2.64 ± 0.49 × 10^-3 mm² /s at 37°C), above which point increased sample viscosity no longer allows for convection currents, thereby eliminating the artifact. The methodology described in this work simplifies quantitative MR acquisition of liquid phantoms at physiological temperature by suppressing convection currents with relatively small changes to intrinsic MR parameters (T₁ increased by 1.4% and T₂ decreased by 17% for 1.5 mM agar at 25°C).