Generation of synthetic CT data using patient specific daily MR image data and image registration

A body mass index-based method for "MR-only" abdominal MR-guided adaptive radiotherapy

PURPOSE

Dose calculation for MR-guided radiotherapy (MRgRT) at the 0.35 T MR-Linac is currently based on deformation of planning CTs (defCT) acquired for each patient. We present a simple and robust bulk density overwrite synthetic CT (sCT) method for abdominal treatments in order to streamline clinical workflows.

METHOD

Fifty-six abdominal patient treatment plans were retrospectively evaluated. All patients had been treated at the MR-Linac using MR datasets for treatment planning and plan adaption and defCT for dose calculation. Bulk density CTs (4M-sCT) were generated from MR images with four material compartments (bone, lung, air, soft tissue). The relative electron densities (RED) for bone and lung were extracted from contoured CT structure average REDs. For soft tissue, a correlation between BMI and RED was evaluated. Dose was recalculated on 4M-sCT and compared to dose distributions on defCTs assessing dose differences in the PTV and organs at risk (OAR).

RESULTS

Mean RED of bone was 1.17 ± 0.02, mean RED of lung 0.17 ± 0.05. The correlation between BMI and RED for soft tissue was statistically significant (p < 0.01). PTV dose differences between 4M-sCT and defCT were Dmean: -0.4 ± 1.0%, D1%: -0.3 ± 1.1% and D95%: -0.5 ± 1.0%. OARs showed D2%: -0.3 ± 1.9% and Dmean: -0.1 ± 1.4% differences. Local 3D gamma index pass rates (2%/2mm) between dose calculated using 4M-sCT and defCT were 96.8 ± 2.6% (range 89.9-99.6%).

CONCLUSION

The presented method for sCT generation enables precise dose calculation for MR-only abdominal MRgRT.

SC-GAN: Structure-completion generative adversarial network for synthetic CT generation from MR images with truncated anatomy

Creating synthetic CT (sCT) from magnetic resonance (MR) images enables MR-based treatment planning in radiation therapy. However, the MR images used for MR-guided adaptive planning are often truncated in the boundary regions due to the limited field of view and the need for sequence optimization. Consequently, the sCT generated from these truncated MR images lacks complete anatomic information, leading to dose calculation error for MR-based adaptive planning. We propose a novel structure-completion generative adversarial network (SC-GAN) to generate sCT with full anatomic details from the truncated MR images. To enable anatomy compensation, we expand input channels of the CT generator by including a body mask and introduce a truncation loss between sCT and real CT. The body mask for each patient was automatically created from the simulation CT scans and transformed to daily MR images by rigid registration as another input for our SC-GAN in addition to the MR images. The truncation loss was constructed by implementing either an auto-segmentor or an edge detector to penalize the difference in body outlines between sCT and real CT. The experimental results show that our SC-GAN achieved much improved accuracy of sCT generation in both truncated and untruncated regions compared to the original cycleGAN and conditional GAN methods.

A systematic literature review: deep learning techniques for synthetic medical image generation and their applications in radiotherapy

The aim of this systematic review is to determine whether Deep Learning (DL) algorithms can provide a clinically feasible alternative to classic algorithms for synthetic Computer Tomography (sCT). The following categories are presented in this study: ∙ MR-based treatment planning and synthetic CT generation techniques. ∙ Generation of synthetic CT images based on Cone Beam CT images. ∙ Low-dose CT to High-dose CT generation. ∙ Attenuation correction for PET images. To perform appropriate database searches, we reviewed journal articles published between January 2018 and June 2023. Current methodology, study strategies, and results with relevant clinical applications were analyzed as we outlined the state-of-the-art of deep learning based approaches to inter-modality and intra-modality image synthesis. This was accomplished by contrasting the provided methodologies with traditional research approaches. The key contributions of each category were highlighted, specific challenges were identified, and accomplishments were summarized. As a final step, the statistics of all the cited works from various aspects were analyzed, which revealed that DL-based sCTs have achieved considerable popularity, while also showing the potential of this technology. In order to assess the clinical readiness of the presented methods, we examined the current status of DL-based sCT generation.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Region of interest focused MRI to synthetic CT translation using regression and segmentation multi-task network

Objective. In MR-only clinical workflow, replacing CT with MR image is of advantage for workflow efficiency and reduces radiation to the patient. An important step required to eliminate CT scan from the workflow is to generate the information provided by CT via an MR image. In this work, we aim to demonstrate a method to generate accurate synthetic CT (sCT) from an MR image to suit the radiation therapy (RT) treatment planning workflow. We show the feasibility of the method and make way for a broader clinical evaluation.Approach. We present a machine learning method for sCT generation from zero-echo-time (ZTE) MRI aimed at structural and quantitative accuracies of the image, with a particular focus on the accurate bone density value prediction. The misestimation of bone density in the radiation path could lead to unintended dose delivery to the target volume and results in suboptimal treatment outcome. We propose a loss function that favors a spatially sparse bone region in the image. We harness the ability of the multi-task network to produce correlated outputs as a framework to enable localization of region of interest (RoI) via segmentation, emphasize regression of values within RoI and still retain the overall accuracy via global regression. The network is optimized by a composite loss function that combines a dedicated loss from each task.Main results. We have included 54 brain patient images in this study and tested the sCT images against reference CT on a subset of 20 cases. A pilot dose evaluation was performed on 9 of the 20 test cases to demonstrate the viability of the generated sCT in RT planning. The average quantitative metrics produced by the proposed method over the test set were-(a) mean absolute error (MAE) of 70 ± 8.6 HU; (b) peak signal-to-noise ratio (PSNR) of 29.4 ± 2.8 dB; structural similarity metric (SSIM) of 0.95 ± 0.02; and (d) Dice coefficient of the body region of 0.984 ± 0.Significance. We demonstrate that the proposed method generates sCT images that resemble visual characteristics of a real CT image and has a quantitative accuracy that suits RT dose planning application. We compare the dose calculation from the proposed sCT and the real CT in a radiation therapy treatment planning setup and show that sCT based planning falls within 0.5% target dose error. The method presented here with an initial dose evaluation makes an encouraging precursor to a broader clinical evaluation of sCT based RT planning on different anatomical regions.

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%.

Global democratisation of proton radiotherapy

Proton radiotherapy is an advanced treatment option compared with conventional x-ray treatment, delivering much lower doses of radiation to healthy tissues surrounding the tumour. However, proton therapy is currently not widely available. In this Review, we summarise the evolution of proton therapy to date, together with the benefits to patients and society. These developments have led to an exponential growth in the number of hospitals using proton radiotherapy worldwide. However, the gap between the number of patients who should be treated with proton radiotherapy and those who have access to it remains large. We summarise the ongoing research and development that is contributing to closing this gap, including the improvement of treatment efficiency and efficacy, and advances in fixed-beam treatments that do not require an enormously large, heavy, and costly gantry. The ultimate goal of decreasing the size of proton therapy machines to fit into standard treatment rooms appears to be within reach, and we discuss future research and development opportunities to achieve this goal.

Evaluation of MRI-only based online adaptive radiotherapy of abdominal region on MR-linac

PURPOSE

A hybrid magnetic resonance linear accelerator (MRL) can perform magnetic resonance imaging (MRI) with high soft-tissue contrast to be used for online adaptive radiotherapy (oART). To obtain electron densities needed for the oART dose calculation, a computed tomography (CT) is often deformably registered to MRI. Our aim was to evaluate an MRI-only based synthetic CT (sCT) generation as an alternative to the deformed CT (dCT)-based oART in the abdominal region.

METHODS

The study data consisted of 57 patients who were treated on a 0.35 T MRL system mainly for abdominal tumors. Simulation MRI-CT pairs of 43 patients were used for training and validation of a prototype convolutional neural network sCT-generation algorithm, based on HighRes3DNet, for the abdominal region. For remaining test patients, sCT images were produced from simulation MRIs and daily MRIs. The dCT-based plans were re-calculated on sCT with identical calculation parameters. The sCT and dCT were compared in terms of geometric agreement and calculated dose.

RESULTS

The mean and one standard deviation of the geometric agreement metrics over dCT-sCT-pairs were: mean error of 8 ± 10 HU, mean absolute error of 49 ± 10 HU, and Dice similarity coefficient of 55 ± 12%, 60 ± 5%, and 82 ± 15% for bone, fat, and lung tissues, respectively. The dose differences between the sCT and dCT-based dose for planning target volumes were 0.5 ± 0.9%, 0.6 ± 0.8%, and 0.5 ± 0.8% at D_2% , D_50% , and D_98% in physical dose and 0.8 ± 1.4%, 0.8 ± 1.2%, and 0.6 ± 1.1% in biologically effective dose (BED). For organs-at-risk, the dose differences of all evaluated dose-volume histogram points were within [-4.5%, 7.8%] and [-1.1 Gy, 3.5 Gy] in both physical dose and BED.

CONCLUSIONS

The geometric agreement metrics were within typically reported values and most average relative dose differences were within 1%. Thus, an MRI-only sCT-based approach is a promising alternative to the current clinical practice of the abdominal oART on MRL.

Deep Learning for Synthetic CT from Bone MRI in the Head and Neck

BACKGROUND AND PURPOSE

Bone MR imaging techniques enable visualization of cortical bone without the need for ionizing radiation. Automated conversion of bone MR imaging to synthetic CT is highly desirable for downstream image processing and eventual clinical adoption. Given the complex anatomy and pathology of the head and neck, deep learning models are ideally suited for learning such mapping.

MATERIALS AND METHODS

This was a retrospective study of 39 pediatric and adult patients with bone MR imaging and CT examinations of the head and neck. For each patient, MR imaging and CT data sets were spatially coregistered using multiple-point affine transformation. Paired MR imaging and CT slices were generated for model training, using 4-fold cross-validation. We trained 3 different encoder-decoder models: Light_U-Net (2 million parameters) and VGG-16 U-Net (29 million parameters) without and with transfer learning. Loss functions included mean absolute error, mean squared error, and a weighted average. Performance metrics included Pearson R, mean absolute error, mean squared error, bone precision, and bone recall. We investigated model generalizability by training and validating across different conditions.

RESULTS

The Light_U-Net architecture quantitatively outperformed VGG-16 models. Mean absolute error loss resulted in higher bone precision, while mean squared error yielded higher bone recall. Performance metrics decreased when using training data captured only in a different environment but increased when local training data were augmented with those from different hospitals, vendors, or MR imaging techniques.

CONCLUSIONS

We have optimized a robust deep learning model for conversion of bone MR imaging to synthetic CT, which shows good performance and generalizability when trained on different hospitals, vendors, and MR imaging techniques. This approach shows promise for facilitating downstream image processing and adoption into clinical practice.

Adaptive Proton Therapy of Pediatric Head and Neck Cases Using MRI-Based Synthetic CTs: Initial Experience of the Prospective KiAPT Study

BACKGROUND AND PURPOSE

Interfractional anatomical changes might affect the outcome of proton therapy (PT). We aimed to prospectively evaluate the role of Magnetic Resonance Imaging (MRI) based adaptive PT for children with tumors of the head and neck and base of skull.

METHODS

MRI verification images were acquired at half of the treatment course. A synthetic computed tomography (CT) image was created using this MRI and a deformable image registration (DIR) to the reference MRI. The methodology was verified with in-silico phantoms and validated using a clinical case with a shrinking cystic hygroma on the basis of dosimetric quantities of contoured structures. The dose distributions on the verification X-ray CT and on the synthetic CT were compared with a gamma-index test using global 2 mm/2% criteria.

RESULTS

Regarding the clinical validation case, the gamma-index pass rate was 98.3%. Eleven patients were included in the clinical study. The most common diagnosis was rhabdomyosarcoma (73%). Craniofacial tumor site was predominant in 64% of patients, followed by base of skull (18%). For one individual case the synthetic CT showed an increase in the median D2 and Dmax dose on the spinal cord from 20.5 GyRBE to 24.8 GyRBE and 14.7 GyRBE to 25.1 GyRBE, respectively. Otherwise, doses received by OARs remained relatively stable. Similarly, the target volume coverage seen by D95% and V95% remained unchanged.

CONCLUSIONS

The method of transferring anatomical changes from MRIs to a synthetic CTs was successfully implemented and validated with simple, commonly available tools. In the frame of our early results on a small cohort, no clinical relevant deterioration for neither PTV coverage nor an increased dose burden to OARs occurred. However, the study will be continued to identify a pediatric patient cohort, which benefits from adaptive treatment planning.

Development of an anthropomorphic multimodality pelvic phantom for quantitative evaluation of a deep-learning-based synthetic computed tomography generation technique

PURPOSE

The objective of this study was to fabricate an anthropomorphic multimodality pelvic phantom to evaluate a deep-learning-based synthetic computed tomography (CT) algorithm for magnetic resonance (MR)-only radiotherapy.

METHODS

Polyurethane-based and silicone-based materials with various silicone oil concentrations were scanned using 0.35 T MR and CT scanner to determine the tissue surrogate. Five tissue surrogates were determined by comparing the organ intensity with patient CT and MR images. Patient-specific organ modeling for three-dimensional printing was performed by manually delineating the structures of interest. The phantom was finally fabricated by casting materials for each structure. For the quantitative evaluation, the mean and standard deviations were measured within the regions of interest on the MR, simulation CT (CT_sim ), and synthetic CT (CT_syn ) images. Intensity-modulated radiation therapy plans were generated to assess the impact of different electron density assignments on plan quality using CT_sim and CT_syn . The dose calculation accuracy was investigated in terms of gamma analysis and dose-volume histogram parameters.

RESULTS

For the prostate site, the mean MR intensities for the patient and phantom were 78.1 ± 13.8 and 86.5 ± 19.3, respectively. The mean intensity of the synthetic image was 30.9 Hounsfield unit (HU), which was comparable to that of the real CT phantom image. The original and synthetic CT intensities of the fat tissue in the phantom were -105.8 ± 4.9 HU and -107.8 ± 7.8 HU, respectively. For the target volume, the difference in D_95% was 0.32 Gy using CT_syn with respect to CT_sim values. The V_65Gy values for the bladder in the plans using CT_sim and CT_syn were 0.31% and 0.15%, respectively.

CONCLUSION

This work demonstrated that the anthropomorphic phantom was physiologically and geometrically similar to the patient organs and was employed to quantitatively evaluate the deep-learning-based synthetic CT algorithm.

Generation of Synthetic-Pseudo MR Images from Real CT Images

This study aimed to generate synthetic MR images from real CT images. CT# mean and standard deviation of a moving window across every pixel in the reconstructed CT images were mapped to their corresponding tissue-mimicking types. Identification of the tissue enabled remapping it to its corresponding intrinsic parameters: T1, T2, and proton density (ρ). Lastly, synthetic weighted MR images of a selected slice were generated by simulating a spin-echo sequence using the intrinsic parameters and proper contrast parameters (TE and TR). Experiments were performed on a 3D multimodality abdominal phantom and on human knees at different TE and TR parameters to confirm the clinical effectiveness of the approach. Results demonstrated the validity of the approach of generating synthetic MR images at different weightings using only CT images and the three predefined mapping functions. The slope of the fitting line and percentage root-mean-square difference (PRD) between real and synthetic image vector representations were (0.73, 10%), (0.9, 18%), and (0.2, 8.7%) for T1-, T2-, and ρ-weighted images of the phantom, respectively. The slope and PRD for human knee images, on average, were 0.89% and 18.8%, respectively. The generated MR images provide valuable guidance for physicians with regard to deciding whether acquiring real MR images is crucial.

Future technological developments in proton therapy - A predicted technological breakthrough

In the current spectrum of cancer treatments, despite high costs, a lack of robust evidence based on clinical outcomes or technical and radiobiological uncertainties, particle therapy and in particular proton therapy (PT) is rapidly growing. Despite proton therapy being more than fifty years old (first proposed by Wilson in 1946) and more than 220,000 patients having been treated with in 2020, many technological challenges remain and numerous new technical developments that must be integrated into existing systems. This article presents an overview of on-going technical developments and innovations that we felt were most important today, as well as those that have the potential to significantly shape the future of proton therapy. Indeed, efforts have been done continuously to improve the efficiency of a PT system, in terms of cost, technology and delivery technics, and a number of different developments pursued in the accelerator field will first be presented. Significant developments are also underway in terms of transport and spatial resolution achievable with pencil beam scanning, or conformation of the dose to the target: we will therefore discuss beam focusing and collimation issues which are important parameters for the development of these techniques, as well as proton arc therapy. State of the art and alternative approaches to adaptive PT and the future of adaptive PT will finally be reviewed. Through these overviews, we will finally see how advances in these different areas will allow the potential for robust dose shaping in proton therapy to be maximised, probably foreshadowing a future era of maturity for the PT technique.

Applying MRI Intensity Normalization on Non-Bone Tissues to Facilitate Pseudo-CT Synthesis from MRI

This study aimed to facilitate pseudo-CT synthesis from MRI by normalizing MRI intensity of the same tissue type to a similar intensity level. MRI intensity normalization was conducted through dividing MRI by a shading map, which is a smoothed ratio image between MRI and a three-intensity mask. Regarding pseudo-CT synthesis from MRI, a conversion model based on a three-layer convolutional neural network was trained and validated. Before MRI intensity normalization, the mean value ± standard deviation of fat tissue in 0.35 T chest MRI was 297 ± 73 (coefficient of variation (CV) = 24.58%), which was 533 ± 91 (CV = 17.07%) in 1.5 T abdominal MRI. The corresponding results were 149 ± 32 (CV = 21.48%) and 148 ± 28 (CV = 18.92%) after intensity normalization. With regards to pseudo-CT synthesis from MRI, the differences in mean values between pseudo-CT and real CT were 3, 15, and 12 HU for soft tissue, fat, and lung/air in 0.35 T chest imaging, respectively, while the corresponding results were 3, 14, and 15 HU in 1.5 T abdominal imaging. Overall, the proposed workflow is reliable in pseudo-CT synthesis from MRI and is more practicable in clinical routine practice compared with deep learning methods, which demand a high level of resources for building a conversion model.

Advanced Magnetic Resonance Imaging of the Skull Base

Magnetic resonance (MR) imaging is a crucial tool for evaluation of the skull base, enabling characterization of complex anatomy by utilizing multiple image contrasts. Recent technical MR advances have greatly enhanced radiologists' capability to diagnose skull base pathology and help direct management. In this paper, we will summarize cutting-edge clinical and emerging research MR techniques for the skull base, including high-resolution, phase-contrast, diffusion, perfusion, vascular, zero echo-time, elastography, spectroscopy, chemical exchange saturation transfer, PET/MR, ultra-high-field, and 3D visualization. For each imaging technique, we provide a high-level summary of underlying technical principles accompanied by relevant literature review and clinical imaging examples.

Zero-TE MRI: Potential Applications in the Oral Cavity and Oropharynx

ABSTRACT

Zero-echo time (ZTE) magnetic resonance imaging (MRI) is the newest in a family of MRI pulse sequences that involve ultrafast sequence readouts, permitting visualization of short-T2 tissues such as cortical bone. Inherent sequence properties enable rapid, high-resolution, quiet, and artifact-resistant imaging. ZTE can be performed as part of a "one-stop-shop" MRI examination for comprehensive evaluation of head and neck pathology. As a potential alternative to computed tomography for bone imaging, this approach could help reduce patient exposure to ionizing radiation and improve radiology resource utilization. Because ZTE is not yet widely used clinically, it is important to understand the technical limitations and pitfalls for diagnosis. Imaging cases are presented to demonstrate potential applications of ZTE for imaging of oral cavity, oropharynx, and jaw anatomy and pathology in adult and pediatric patients. Emerging studies indicate promise for future clinical implementation based on synthetic computed tomography image generation, 3D printing, and interventional applications.

MR-guided proton therapy: a review and a preview

Background

The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics.

Methods

Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches.

Conclusions

Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.

Online daily adaptive proton therapy

It is recognized that the use of a single plan calculated on an image acquired some time before the treatment is generally insufficient to accurately represent the daily dose to the target and to the organs at risk. This is particularly true for protons, due to the physical finite range. Although this characteristic enables the generation of steep dose gradients, which is essential for highly conformal radiotherapy, it also tightens the dependency of the delivered dose to the range accuracy. In particular, the use of an outdated patient anatomy is one of the most significant sources of range inaccuracy, thus affecting the quality of the planned dose distribution. A plan should be ideally adapted as soon as anatomical variations occur, ideally online. In this review, we describe in detail the different steps of the adaptive workflow and discuss the challenges and corresponding state-of-the art developments in particular for an online adaptive strategy.

Towards a generalised development of synthetic CT images and assessment of their dosimetric accuracy

Background: The interest in generating "synthetic computed tomography (CT) images" from magnetic resonance (MR) images has been increasing over the past years due to advances in MR guidance for radiotherapy. A variety of methods for synthetic CT creation have been developed, from simple bulk density assignment to complex machine learning algorithms.Material and methods: In this study, we present a general method to determine simplistic synthetic CTs and evaluate them according to their dosimetric accuracy. It separates the requirements on the MR image and the associated calculation effort to generate a synthetic CT. To evaluate the significance of the dosimetric accuracy under realistic conditions, clinically common uncertainties including position shifts and Hounsfield lookup table (HLUT) errors were simulated. To illustrate our approach, we first translated CT images from a test set of six pelvic cancer patients to relative electron density (ED) via a clinical HLUT. For each patient, seven simplified ED images (simED) were generated at different levels of complexity, ranging from one to four tissue classes. Then, dose distributions optimised on the reference ED image and the simEDs were compared to each other in terms of gamma pass rates (2 mm/2% criteria) and dose volume metrics.Results: For our test set, best results were obtained for simEDs with four tissue classes representing fat, soft tissue, air, and bone. For this simED, gamma pass rates of 99.95% (range: 99.72-100%) were achieved. The decrease in accuracy from ED simplification was smaller in this case than the influence of the uncertainty scenarios on the reference image, both for gamma pass rates and dose volume metrics.Conclusions: The presented workflow helps to determine the required complexity of synthetic CTs with respect to their dosimetric accuracy. The investigated cases showed potential simplifications, based on which the synthetic CT generation could be faster and more reproducible.

Automated air region delineation on MRI for synthetic CT creation

Automatically and accurately separating air from other low signal regions (especially bone, liver, etc) in an MRI is difficult because these tissues produce similar MR intensities, resulting in errors in synthetic CT generation for MRI-based radiation therapy planning. This work aims to develop a technique to accurately and automatically determine air-regions for MR-guided adaptive radiation therapy. CT and MRI scans (T2-weighted) of phantoms with fabricated air-cavities and abdominal cancer patients were used to establish an MR intensity threshold for air delineation. From the phantom data, air/tissue boundaries in MRI were identified by CT-MRI registration. A formula relating the MRI intensities of air and surrounding materials was established to auto-threshold air-regions. The air-regions were further refined by using quantitative image texture features. A naive Bayesian classifier was trained using the extracted features with a leave-one-out cross validation technique to differentiate air from non-air voxels. The multi-step air auto-segmentation method was tested against the manually segmented air-regions. The dosimetry impacts of the air-segmentation methods were studied. Air-regions in the abdomen can be segmented on MRI within 1 mm accuracy using a multi-step auto-segmentation method as compared to manually delineated contours. The air delineation based on the MR threshold formula was improved using the MRI texture differences between air and tissues, as judged by the area under the receiver operating characteristic curve of 81% when two texture features (autocorrelation and contrast) were used. The performance increased to 82% with using three features (autocorrelation, sum-variance, and contrast). Dosimetric analysis showed no significant difference between the auto-segmentation and manual MR air delineation on commonly used dose volume parameters. The proposed techniques consisting of intensity-based auto-thresholding and image texture-based voxel classification can automatically and accurately segment air-regions on MRI, allowing synthetic CT to be generated quickly and precisely for MR-guided adaptive radiation therapy.

On the accuracy of bulk synthetic CT for MR-guided online adaptive radiotherapy

PURPOSE

MR-guided radiotherapy (MRgRT) relies on the daily assignment of a relative electron density (RED) map to allow the fraction specific dose calculation. One approach to assign the RED map consists of segmenting the daily magnetic resonance image into five different density levels and assigning a RED bulk value to each level to generate a synthetic CT (sCT). The aim of this study is to evaluate the dose calculation accuracy of this approach for applications in MRgRT.

METHODS

A planning CT (pCT) was acquired for 26 patients with abdominal and pelvic lesions and segmented in five levels similar to an online approach: air, lung, fat, soft tissue and bone. For each patient, the median RED value was calculated for fat, soft tissue and bone. Two sCTs were generated assigning different bulk values to the segmented levels on pCT: The sCT_ICRU uses the RED values recommended by ICRU46, and the sCT_tailor uses the median patient-specific RED values. The same treatment plan was calculated on two the sCTs and the pCT. The dose calculation accuracy was investigated in terms of gamma analysis and dose volume histogram parameters.

RESULTS

Good agreement was found between dose calculated on sCTs and pCT (gamma passing rate 1%/1 mm equal to 91.2% ± 6.9% for sCT_ICRU and 93.7% ± 5.3% b or sCT_tailor). The mean difference in estimating V95 (PTV) was equal to 0.2% using sCT_tailor and 1.2% using sCT_ICRU, respect to pCT values CONCLUSIONS: The bulk sCT guarantees a high level of dose calculation accuracy also in presence of magnetic field, making this approach suitable to MRgRT. This accuracy can be improved by using patient-specific RED values.

Multiatlas Fusion with a Hybrid CT Number Correction Technique for Subject-Specific Pseudo-CT Estimation in the Context of MRI-Only Radiation Therapy

OBJECTIVE

To propose a hybrid multiatlas fusion and correction approach to estimate a pseudo-computed tomography (pCT) image from T2-weighted brain magnetic resonance (MR) images in the context of MRI-only radiotherapy.

MATERIALS AND METHODS

A set of eleven pairs of T2-weighted MR and CT brain images was included. Using leave-one-out cross-validation, atlas MR images were registered to the target MRI with multimetric, multiresolution deformable registration. The subsequent deformations were applied to the atlas CT images, producing uncorrected pCT images. Afterward, a three-dimensional hybrid CT number correction technique was used. This technique uses information about MR intensity, spatial location, and tissue label from segmented MR images with the fuzzy c-means algorithm and combines them in a weighted fashion to correct Hounsfield unit values of the uncorrected pCT images. The corrected pCT images were then fused into a final pCT image.

RESULTS

The proposed hybrid approach proved to be performant in correcting Hounsfield unit values in terms of qualitative and quantitative measures. Average correlation was 0.92 and 0.91 for the proposed approach by taking the mean and the median, respectively, compared with 0.86 for the uncorrected unfused version. Average values of dice similarity coefficient for bone were 0.68 and 0.72 for the fused corrected pCT images by taking the mean and the median, respectively, compared with 0.65 for the uncorrected unfused version indicating a significant bone estimation improvement.

CONCLUSION

A hybrid fusion and correction method is presented to estimate a pCT image from T2-weighted brain MR images.

Synthetic 4D-CT of the thorax for treatment plan adaptation on MR-guided radiotherapy systems

MR-guided radiotherapy treatment planning utilises the high soft-tissue contrast of MRI to reduce uncertainty in delineation of the target and organs at risk. Replacing 4D-CT with MRI-derived synthetic 4D-CT would support treatment plan adaptation on hybrid MR-guided radiotherapy systems for inter- and intrafractional differences in anatomy and respiration, whilst mitigating the risk of CT to MRI registration errors. Three methods were devised to calculate synthetic 4D and midposition (time-weighted mean position of the respiratory cycle) CT from 4D-T1w and Dixon MRI. The first approach employed intensity-based segmentation of Dixon MRI for bulk-density assignment (sCTD). The second step added spine density information using an atlas of CT and Dixon MRI (sCTDS). The third iteration used a polynomial function relating Hounsfield units and normalised T1w image intensity to account for variable lung density (sCTDSL). Motion information in 4D-T1w MRI was applied to generate synthetic CT in midposition and in twenty respiratory phases. For six lung cancer patients, synthetic 4D-CT was validated against 4D-CT in midposition by comparison of Hounsfield units and dose-volume metrics. Dosimetric differences found by comparing sCTD,DS,DSL and CT were evaluated using a Wilcoxon signed-rank test (p   =  0.05). Compared to sCTD and sCTDS, planning on sCTDSL significantly reduced absolute dosimetric differences in the planning target volume metrics to less than 98 cGy (1.7% of the prescribed dose) on average. When comparing sCTDSL and CT, average radiodensity differences were within 97 Hounsfield units and dosimetric differences were significant only for the planning target volume D99% metric. All methods produced clinically acceptable results for the organs at risk in accordance with the UK SABR consensus guidelines and the LungTech EORTC phase II trial. The overall good agreement between sCTDSL and CT demonstrates the feasibility of employing synthetic 4D-CT for plan adaptation on hybrid MR-guided radiotherapy systems.

Quality assurance measurements of geometric accuracy for magnetic resonance imaging-based radiotherapy treatment planning

BACKGROUND

Using magnetic resonance imaging (MRI) as the only imaging method for radiotherapy treatment planning (RTP) is becoming more common as MRI-only RTP solutions have evolved. The geometric accuracy of MR images is an essential factor of image quality when determining the suitability of MRI for RTP. The need is therefore clear for clinically feasible quality assurance (QA) methods for the geometric accuracy measurement.

MATERIALS AND METHODS

This work evaluates long-term stability of geometric accuracy and the validity of a 2D geometric accuracy QA method compared to a prototype 3D method and analysis software in routine QA. The long-term follow-up measurements were conducted on one of the 1.5 T scanners over a period of 19 months using both methods. Inter-scanner variability of geometric distortions was also evaluated in three 1.5 T and three 3 T MRI scanners from a single vendor by using the prototype 3D QA method.

RESULTS

The geometric accuracy of the magnetic resonance for radiotherapy (MR-RT) platform remained stable within 2 mm at distances of <250 mm from isocenter. All scanners achieved good geometric accuracy with mean geometric distortions of <1 mm at <150 mm and <2 mm at <250 mm from the isocenter. Both measurement methods provided relevant information about geometric distortions.

CONCLUSIONS

Geometric distortions are often considered a limitation of MRI-only RTP. Results indicate that geometric accuracy of modern scanners remain within acceptable limits by default even after many years of clinical use based on the 3D QA evaluation.

Evolutionary Machine Learning for Multi-Objective Class Solutions in Medical Deformable Image Registration

Current state-of-the-art medical deformable image registration (DIR) methods optimize a weighted sum of key objectives of interest. Having a pre-determined weight combination that leads to high-quality results for any instance of a specific DIR problem (i.e., a class solution) would facilitate clinical application of DIR. However, such a combination can vary widely for each instance and is currently often manually determined. A multi-objective optimization approach for DIR removes the need for manual tuning, providing a set of high-quality trade-off solutions. Here, we investigate machine learning for a multi-objective class solution, i.e., not a single weight combination, but a set thereof, that, when used on any instance of a specific DIR problem, approximates such a set of trade-off solutions. To this end, we employed a multi-objective evolutionary algorithm to learn sets of weight combinations for three breast DIR problems of increasing difficulty: 10 prone-prone cases, 4 prone-supine cases with limited deformations and 6 prone-supine cases with larger deformations and image artefacts. Clinically-acceptable results were obtained for the first two problems. Therefore, for DIR problems with limited deformations, a multi-objective class solution can be machine learned and used to compute straightforwardly multiple high-quality DIR outcomes, potentially leading to more efficient use of DIR in clinical practice.

Current state and future applications of radiological image guidance for particle therapy

In this review paper, we first give a short overview of radiological image guidance in photon radiotherapy, placing emphasis on the fact that linac based radiotherapy has outpaced particle therapy in the adoption of volumetric image guidance. While cone beam computed tomography (CBCT) has been an established technique in linac treatment rooms for almost two decades, the widespread adoption of volumetric image guidance in particle therapy, whether by means of CBCT or in-room CT imaging, is recent. This lag may be attributable to the bespoke nature and lower number of particle therapy installations, as well as the differences in geometry between those installations and linac treatment rooms. In addition, for particle therapy the so called shift invariance of the dose distribution rarely applies. An overview of the different volumetric image guidance solutions found at modern particle therapy facilities is provided, covering gantry, nozzle, C-arm, and couch-mounted CBCT as well different in-room CT configurations. A summary of the use of in-room volumetric imaging data beyond anatomy-based positioning is also presented as well as the necessary corrections to CBCT images for accurate water equivalent thickness calculation. Finally, the use of non-ionizing imaging modalities is discussed.

A phantom study to create synthetic CT from orthogonal twodimensional cine MRI and evaluate the effect of irregular breathing

An exciting innovation in radiotherapy is the use of real-time MRI for treatment adaptation. This study proposes an in-silico framework for the generation of 3D synthetic CT (sCT) from orthogonal interleaved 2D cine MRI data to overcome the lack of electron density information in MR images. The method uses pre-treatment data to build a patient breathing motion model. This model is then driven by surrogates extracted from cine MR images during the treatment. The effect of irregular breathing on the motion model is also evaluated by simulating different motion components related to uncorrelated diaphragm, chest and tumor motion. 3D sCT were successfully created for each of the 512 cine MRI pairs in the digital phantom study. The analysis showed that the diaphragm position was a good surrogate to rescale the 3D breathing motion for the current regular breathing phase. However, respiratory and tumor motion were correlated in only 59% of the phases, resulting in tumor position uncertainties of up to 3mm. The inclusion of additional chest and tumor motion information improved the accuracy for irregular changes in breathing pattern and enhanced the tumor position uncertainties to less than 1mm. This study successfully demonstrated a proof-ofprinciple for a digital phantom dataset based on patient parameters, providing a way to create real-time 3D electron density volumes and enhancing the need to account for irregular breathing pattern.

Teaching treatment planning for protons with educational open-source software: experience with FoCa and matRad

Open‐source, MATLAB‐based treatment planning systems FoCa and matRAD were used in a pilot project for training prospective medical physicists and postgraduate physics students in treatment planning and beam modeling techniques for proton therapy. In the four exercises designed, students learnt how proton pencil beams are modeled and how dose is calculated in three‐dimensional voxelized geometries, how pencil beam scanning plans (PBS) are constructed, the rationale behind the choice of spot spacing in patient plans, and the dosimetric differences between photon IMRT and proton PBS plans. Sixty students of two courses participated in the pilot project, with over 90% of satisfactory rating from student surveys. The pilot experience will certainly be continued.

An uncertainty metric to evaluate deformation vector fields for dose accumulation in radiotherapy

Background and purpose

In adaptive radiotherapy, deformable image registration (DIR) is used to propagate delineations of tumors and organs into a new therapy plan and to calculate the accumulated total dose. Many DIR accuracy metrics have been proposed. An alternative proposed here could be a local uncertainty (LU) metric for DIR results.

Materials and methods

The LU represented the uncertainty of each DIR position and was focused on deformation evaluation in uniformly-dense regions. Four cases demonstrated LU calculations: two head and neck cancer cases, a lung cancer case, and a prostate cancer case. Each underwent two CT examinations for radiotherapy planning.

Results

LU maps were calculated from each DIR of the clinical cases. Reduced fat regions had LUs of 4.6 ± 0.9 mm, 4.8 ± 1.0 mm, and 4.5 ± 0.7 mm, while the shrunken left parotid gland had a LU of 4.1 ± 0.8 mm and the shrunken lung tumor had a LU of 3.7 ± 0.7 mm. The bowels in the pelvic region had a LU of 10.2 ± 3.7 mm. LU histograms for the cases were similar and 99% of the voxels had a LU < 3 mm.

Conclusions

LU is a new uncertainty metric for DIR that was demonstrated for clinical cases. It had a tolerance of <3 mm.

A Novel method to generate on-board 4D MRI using prior 4D MRI and on-board kV projections from a conventional LINAC for target localization in liver SBRT

PURPOSE

On-board MRI can provide superb soft tissue contrast for improving liver SBRT localization. However, the availability of on-board MRI in clinics is extremely limited. On the contrary, on-board kV imaging systems are widely available on radiotherapy machines, but its capability to localize tumors in soft tissue is limited due to its poor soft tissue contrast. This study aims to explore the feasibility of using an on-board kV imaging system and patient prior knowledge to generate on-board four-dimensional (4D)-MRI for target localization in liver SBRT.

METHODS

Prior 4D MRI volumes were separated into end of expiration (EOE) phase (MRI_prior ) and all other phases. MRI_prior was used to generate a synthetic CT at EOE phase (sCT_prior ). On-board 4D MRI at each respiratory phase was considered a deformation of MRI_prior . The deformation field map (DFM) was estimated by matching DRRs of the deformed sCT_prior to on-board kV projections using a motion modeling and free-form deformation optimization algorithm. The on-board 4D MRI method was evaluated using both XCAT simulation and real patient data. The accuracy of the estimated on-board 4D MRI was quantitatively evaluated using Volume Percent Difference (VPD), Volume Dice Coefficient (VDC), and Center of Mass Shift (COMS). Effects of scan angle and number of projections were also evaluated.

RESULTS

In the XCAT study, VPD/VDC/COMS among all XCAT scenarios were 10.16 ± 1.31%/0.95 ± 0.01/0.88 ± 0.15 mm using orthogonal-view 30° scan angles with 102 projections. The on-board 4D MRI method was robust against the various scan angles and projection numbers evaluated. In the patient study, estimated on-board 4D MRI was generated successfully when compared to the "reference on-board 4D MRI" for the liver patient case.

CONCLUSIONS

A method was developed to generate on-board 4D MRI using prior 4D MRI and on-board limited kV projections. Preliminary results demonstrated the potential for MRI-based image guidance for liver SBRT using only a kV imaging system on a conventional LINAC.

Evaluation of Dosimetric Robustness of Carbon Ion Boost Therapy for Anal Carcinoma

Purpose

The radiation therapy treatment outcome of human papillomavirus-negative anal carcinoma may be improved by the biological effectiveness of carbon ions. However, abdominal tissue motion can compromise the precision of carbon ion therapy. This work aims to evaluate the dosimetric feasibility of carbon ion boost (CIB) therapy for anal carcinoma.

Materials and Methods

An algorithm to generate computed tomographies based on daily magnetic resonance imaging data and deformable image registration was developed. By means of this algorithm, fractional computed tomography data for 54 treatment fractions for 3 different patients with anal carcinoma were derived. The dose for a sequential CIB (CIB_seq) treatment plan was recalculated on the fractional computed tomography data and accumulated over the number of fractions. The resulting dose distributions were compared to standard intensity-modulated radiation therapy treatment with an integrated photon boost.

Results

For the investigated patient cases, similar dosimetric results for CIB_seq treatment and for intensity-modulated radiation therapy with an integrated photon boost were found. For CIB_seq treatment, bladder-filling variation had the strongest influence on the dose distribution. However, the detrimental effects on the mean target dose remained below 1 Gy (RBE) as compared to photon therapy.

Conclusion

This study shows the dosimetric feasibility of CIB therapy for anal carcinoma for the first time and gives reason for clinical exploitation of the enhanced biological effect of carbon ions for patients with human papillomavirus-negative anal cancer.