Feasibility of MR-only proton dose calculations for prostate cancer radiotherapy using a commercial pseudo-CT generation method

Generating synthetic computed tomography for radiotherapy: SynthRAD2023 challenge report

Radiation therapy plays a crucial role in cancer treatment, necessitating precise delivery of radiation to tumors while sparing healthy tissues over multiple days. Computed tomography (CT) is integral for treatment planning, offering electron density data crucial for accurate dose calculations. However, accurately representing patient anatomy is challenging, especially in adaptive radiotherapy, where CT is not acquired daily. Magnetic resonance imaging (MRI) provides superior soft-tissue contrast. Still, it lacks electron density information, while cone beam CT (CBCT) lacks direct electron density calibration and is mainly used for patient positioning. Adopting MRI-only or CBCT-based adaptive radiotherapy eliminates the need for CT planning but presents challenges. Synthetic CT (sCT) generation techniques aim to address these challenges by using image synthesis to bridge the gap between MRI, CBCT, and CT. The SynthRAD2023 challenge was organized to compare synthetic CT generation methods using multi-center ground truth data from 1080 patients, divided into two tasks: (1) MRI-to-CT and (2) CBCT-to-CT. The evaluation included image similarity and dose-based metrics from proton and photon plans. The challenge attracted significant participation, with 617 registrations and 22/17 valid submissions for tasks 1/2. Top-performing teams achieved high structural similarity indices (≥0.87/0.90) and gamma pass rates for photon (≥98.1%/99.0%) and proton (≥97.3%/97.0%) plans. However, no significant correlation was found between image similarity metrics and dose accuracy, emphasizing the need for dose evaluation when assessing the clinical applicability of sCT. SynthRAD2023 facilitated the investigation and benchmarking of sCT generation techniques, providing insights for developing MRI-only and CBCT-based adaptive radiotherapy. It showcased the growing capacity of deep learning to produce high-quality sCT, reducing reliance on conventional CT for treatment planning.

Evaluation of magnetic resonance imaging derived synthetic computed tomography for proton therapy planning in prostate cancer

Background and purpose

Magnetic resonance imaging (MRI)-only workflow is used in photon radiotherapy (RT) today, but not yet for protons. To bring MRI-only proton RT into clinical use, proton dose calculation on MRI-derived synthetic CT (sCT) must be validated. We evaluated proton dose calculation accuracy of prostate cancer proton plans using a commercially available sCT generator already validated for photon planning.

Materials and methods

The retrospective planning study included 10 prostate cancer patients who underwent MRI and planning CT (pCT) before RT. sCT were generated from the MRI with MRI Planner v2.3, and compared to pCT using structural mean absolute error (MAE). The pCT was used to create one-arc volumetric modulated arc therapy (VMAT) photon plan and two-field intensity modulated proton therapy (IMPT) proton plan. Each plan was recalculated on the sCT and compared to pCT doses. Dose volume histogram parameters, gamma analyses and range differences were evaluated.

Results

Median MAE for the body contour was 71 HU. Dose differences between pCT and sCT were small and similar for VMAT and IMPT plans. Median (range) gamma pass rates were lower for IMPT plans with 95.8 (89.3-98.7) % compared to VMAT plans with 99.4 (91.2-99.6) %. The proton range difference was 1.0 (interquartile range -0.1 - 0.2) mm deeper for sCT compared to the reference.

Conclusion

MRI-only IMPT planning for prostate cancer seems feasible in a clinical setting for the evaluated beam arrangement and sCT generator. More patients and evaluation of other beam arrangements are needed for a more general conclusion.

Validation of an MR-based multimodal method for molecular composition and proton stopping power ratio determination using ex vivo animal tissues and tissue-mimicking phantoms

Objective. Range uncertainty in proton therapy is an important factor limiting clinical effectiveness. Magnetic resonance imaging (MRI) can measure voxel-wise molecular composition and, when combined with kilovoltage CT (kVCT), accurately determine mean ionization potential (Im), electron density, and stopping power ratio (SPR). We aimed to develop a novel MR-based multimodal method to accurately determine SPR and molecular compositions. This method was evaluated in tissue-mimicking andex vivoporcine phantoms, and in a brain radiotherapy patient.Approach. Four tissue-mimicking phantoms with known compositions, two porcine tissue phantoms, and a brain cancer patient were imaged with kVCT and MRI. Three imaging-based values were determined: SPRCM(CT-based Multimodal), SPRMM(MR-based Multimodal), and SPRstoich(stoichiometric calibration). MRI was used to determine two tissue-specific quantities of the Bethe Bloch equation (Im, electron density) to compute SPRCMand SPRMM. Imaging-based SPRs were compared to measurements for phantoms in a proton beam using a multilayer ionization chamber (SPRMLIC).Main results. Root mean square errors relative to SPRMLICwere 0.0104(0.86%), 0.0046(0.45%), and 0.0142(1.31%) for SPRCM, SPRMM, and SPRstoich, respectively. The largest errors were in bony phantoms, while soft tissue and porcine tissue phantoms had <1% errors across all SPR values. Relative to known physical molecular compositions, imaging-determined compositions differed by approximately ≤10%. In the brain case, the largest differences between SPRstoichand SPRMMwere in bone and high lipids/fat tissue. The magnitudes and trends of these differences matched phantom results.Significance. Our MR-based multimodal method determined molecular compositions and SPR in various tissue-mimicking phantoms with high accuracy, as confirmed with proton beam measurements. This method also revealed significant SPR differences compared to stoichiometric kVCT-only calculation in a clinical case, with the largest differences in bone. These findings support that including MRI in proton therapy treatment planning can improve the accuracy of calculated SPR values and reduce range uncertainties.

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.

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.

A study to determine the impact of IMPT optimization techniques on prostate synthetic CT image sets dose comparison against CT image sets

Background

The objective of this study is to determine the impact of intensity modulated proton therapty (IMPT) optimization techniques on the proton dose comparison of commercially available magnetic resonance for calculating attenuation (MRCA T) images, a synthetic computed tomography CT (sCT) based on magnetic resonance imaging (MRI) scan against the CT images and find out the optimization technique which creates plans with the least dose differences against the regular CT image sets.

Material and methods

Regular CT data sets and sCT image sets were obtained for 10 prostate patients for the study. Six plans were created using six distinct IMPT optimization techniques including multi-field optimization (MFO), single field uniform dose (SFUD) optimization, and robust optimization (RO) in CT image sets. These plans were copied to MRCA T, sCT datasets and doses were computed. Doses from CT and MRCA T data sets were compared for each patient using 2D dose distribution display, dose volume histograms (DVH), homogeneity index (HI), conformation number (CN) and 3D gamma analysis. A two tailed t-test was conducted on HI and CN with 5% significance level with a null hypothesis for CT and sCT image sets.

Results

Analysis of ten CT and sCT image sets with different IMPT optimization techniques shows that a few of the techniques show significant differences between plans for a few evaluation parameters. Isodose lines, DVH, HI, CN and t-test analysis shows that robust optimizations with 2% range error incorporated results in plans, when re-computed in sCT image sets results in the least dose differences against CT plans compared to other optimization techniques. The second best optimization technique with the least dose differences was robust optimization with 5% range error.

Conclusion

This study affirmatively demonstrates the impact of IMPT optimization techniques on synthetic CT image sets dose comparison against CT images and determines the robust optimization with 2% range error as the optimization technique which gives the least dose difference when compared to CT plans.

An MRI sequence independent convolutional neural network for synthetic head CT generation in proton therapy

A magnetic resonance imaging (MRI) sequence independent deep learning technique was developed and validated to generate synthetic computed tomography (sCT) scans for MR guided proton therapy. 47 meningioma patients previously undergoing proton therapy based on pencil beam scanning were divided into training (33), validation (6), and test (8) cohorts. T₁, T₂, and contrast enhanced T₁ (T1CM) MRI sequences were used in combination with the planning CT (pCT) data to train a 3D U-Net architecture with ResNet-Blocks. A hyperparameter search was performed including two loss functions, two group sizes of normalisation, and depth of the network. Training outcome was compared between models trained for each individual MRI sequence and for all sequences combined. The performance was evaluated based on a metric and dosimetric analysis as well as spot difference maps. Furthermore, the influence of immobilisation masks that are not visible on MRIs was investigated. Based on the hyperparameter search, the final model was trained with fixed features per group for the group normalisation, six down-convolution steps, an input size of 128×192×192, and feature loss. For the test dataset for body/bone the mean absolute error (MAE) values were on average 79.8/216.3Houndsfield unit (HU) when trained using T1 images, 71.1/186.1HU for T2, and 82.9/236.4HU for T1CM. The structural similarity metric (SSIM) ranged from 0.95 to 0.98 for all sequences. The investigated dose parameters of the target structures agreed within 1% between original proton treatment plans and plans recalculated on sCTs. The spot difference maps had peaks at ±0.2cm and for 98% of all spots the difference was less than 1cm. A novel MRI sequence independent sCT generator was developed, which suggests that the training phase of neural networks can be disengaged from specific MRI acquisition protocols. In contrast to previous studies, the patient cohort consisted exclusively of actual proton therapy patients (i.e. "real-world data").

Intra-fractional per-beam adaptive workflow to mitigate the need for a rotating gantry during MRI-guided proton therapy

The integration of real-time magnetic resonance imaging (MRI)-guidance and proton therapy would potentially improve the proton dose steering capability by reducing daily uncertainties due to anatomical variations. The use of a fixed beamline coupled with an axial patient couch rotation would greatly simplify the proton delivery with MRI-guidance. Nonetheless, it is mandatory to assure that the plan quality is not deteriorated by the anatomical deformations due to patient rotation. In this work, an in-house tool allowing for intra-fractional per-beam adaptation of intensity-modulated proton plans (BeamAdapt) was implemented through features available in RayStation. A set of three MRIs was acquired for two healthy volunteers (V1, V2): (1) no rotation/static, (2) rotation to the right and (3) left. V1 was rotated by 15º, to simulate a clinical pediatric abdominal case and V2 by 45º, to simulate an extreme patient rotation case. For each volunteer, a total of four intensity-modulated pencil beam scanning plans were optimized on the static MRI using virtual abdominal targets and 2-3 posterior-oblique beams. Beam angles were defined according to the angulations on the rotated MRIs. With BeamAdapt, each original plan was first converted into separate plans with one beam per plan. In an iterative order, individual beam doses were non-rigidly deformed to the rotated anatomies and re-optimized accounting for the consequent deformations and the beam doses delivered so far. For evaluation, the final adapted dose distribution was propagated back to the static MRI. Planned and adapted dose distributions were compared by computing relative differences between dose-volume histogram (DVH) metrics. Absolute target dose differences were on average below 1% and mean dose organs-at-risk differences were below 3%. With BeamAdapt, not only intra-fractional per-beam proton plan adaptation coupled with axial patient rotation is possible but also the need for a rotating gantry during MRI-guidance might be mitigated.

Geometrical imaging accuracy and imaging and plan quality for prostate cancer on a 1.5T MRLinac in patients with a unilateral hip implant

PURPOSE

To assess the feasibility of prostate cancer radiotherapy for patients with a hip implant on an 1.5T MRI-Linac (MRL) in terms of geometrical image accuracy, image quality, and plan quality.

METHODS

Pretreatment MRI images on a 1.5T MRL and 3T MRI consisting of a T2-weighted 3D delineation scan and main magnetic field homogeneity (B0) scan were performed in 6 patients with a unilateral hip implant. System specific geometrical errors due to gradient non-linearity were determined for the MRL. Within the prostate and skin contour, B0 inhomogeneity, gradient non-linearity error and the total geometrical error (vector summation of the prior two) was determined. Image quality was determined by visually scoring the extent of implant-born image artifacts. A treatment planning study was performed on 5 patients to quantify the impact of the implant on plan quality, in which conventional MRL IMRT plans were created, as well as plans which avoid radiation through the left or right femur.

RESULTS

The total maximum geometrical error in the prostate was < 1 mm and the skin contour < 1.7 mm; in all cases the machine-specific gradient error was most dominant. The B0 error for the MRlinac MRI could partly be predicted based on the pre-treatment 3T scan. Image quality for all patients was sufficient at 1.5T MRL. Plan comparison showed that, even with avoidance of the hips, in all cases sufficient target coverage could be obtained with similar D1cc and D5cc to rectum and bladder, while V28Gy was slightly poorer in only the rectum for femur avoidance.

CONCLUSION

We showed that geometrical accuracy, image quality and plan quality for six prostate patients with a hip implant or hip fixation treated on a 1.5T MRL did not show relevant deterioration for the used image settings, which allowed safe treatment.

Feasibility of using a single MRI acquisition for fiducial marker localization and synthetic CT generation towards MRI-only prostate radiation therapy treatment planning

Purpose.To investigate the feasibility of using a single MRI acquisition for fiducial marker identification and synthetic CT (sCT) generation towards MRI-only treatment planning for prostate external beam radiation therapy (EBRT).Methods.Seven prostate cancer patients undergoing EBRT, each with three implanted gold fiducial markers, participated in this study. In addition to the planning CT scan, all patients were scanned on a 3 T MR scanner with a 3D double-echo gradient echo (GRE) sequence. Quantitative susceptibility mapping (QSM) was performed for marker localization. QSM-derived marker positions were compared to those from CT. The bulk density assignment technique for sCT generation was adopted. The magnitude GRE images were segmented into muscle, bone, fat, and air using a combination of unsupervised intensity-based classification of soft tissue and convolutional neural networks (CNN) for bone segmentation.Results.All implanted markers were visualized and accurately identified (average error: 0.7 ± 0.5 mm). QSM generated distinctive contrast for hemorrhage, calcifications, and gold fiducial markers. The estimated susceptibility/HU values on QSM/CT for gold and calcifications were 31.5 ± 2.9 ppm/1220 ± 100 HU and 14.6 ± 0.9 ppm/440 ± 100 HU, respectively. The intensity-based soft tissue classification resulted in an average Dice score of 0.97 ± 0.02; bone segmentation using CNN resulted in an average Dice score of 0.93 ± 0.03.Conclusion.This work indicates the feasibility of simultaneous fiducial marker identification and sCT generation using a single MRI acquisition. Future works includes evaluation of the proposed method in a large cohort of patients with optimized acquisition parameters as well as dosimetric evaluations.

Synthetic CT for single-fraction neoadjuvant partial breast irradiation on an MRI-linac

A synthetic computed tomography (sCT) is required for daily plan optimization on an MRI-linac. Yet, only limited information is available on the accuracy of dose calculations on sCT for breast radiotherapy. This work aimed to (1) evaluate dosimetric accuracy of treatment plans for single-fraction neoadjuvant partial breast irradiation (PBI) on a 1.5 T MRI-linac calculated on a) bulk-density sCT mimicking the current MRI-linac workflow and b) deep learning-generated sCT, and (2) investigate the number of bulk-density levels required. For ten breast cancer patients we created three bulk-density sCTs of increasing complexity from the planning-CT, using bulk-density for: (1) body, lungs, and GTV (sCT_BD1); (2) volumes for sCT_BD1plus chest wall and ipsilateral breast (sCT_BD2); (3) volumes for sCT_BD2plus ribs (sCT_BD3); and a deep learning-generated sCT (sCT_DL) from a 1.5 T MRI in supine position. Single-fraction neoadjuvant PBI treatment plans for a 1.5 T MRI-linac were optimized on each sCT and recalculated on the planning-CT. Image evaluation was performed by assessing mean absolute error (MAE) and mean error (ME) in Hounsfield Units (HU) between the sCTs and the planning-CT. Dosimetric evaluation was performed by assessing dose differences, gamma pass rates, and dose-volume histogram (DVH) differences. The following results were obtained (median across patients for sCT_BD1/sCT_BD2/sCT_BD3/sCT_DLrespectively): MAE inside the body contour was 106/104/104/75 HU and ME was 8/9/6/28 HU, mean dose difference in the PTV_GTVwas 0.15/0.00/0.00/-0.07 Gy, median gamma pass rate (2%/2 mm, 10% dose threshold) was 98.9/98.9/98.7/99.4%, and differences in DVH parameters were well below 2% for all structures except for the skin in the sCT_DL. Accurate dose calculations for single-fraction neoadjuvant PBI on an MRI-linac could be performed on both bulk-density and deep learning sCT, facilitating further implementation of MRI-guided radiotherapy for breast cancer. Balancing simplicity and accuracy, sCT_BD2showed the optimal number of bulk-density levels for a bulk-density approach.

An MR-only acquisition and artificial intelligence based image-processing protocol for photon and proton therapy using a low field MR

OBJECTIVE

Recent developments on synthetically generated CTs (sCT), hybrid MRI linacs and MR-only simulations underlined the clinical feasibility and acceptance of MR guided radiation therapy. However, considering clinical application of open and low field MR with a limited field of view can result in truncation of the patient's anatomy which further affects the MR to sCT conversion. In this study an acquisition protocol and subsequent MR image stitching is proposed to overcome the limited field of view restriction of open MR scanners, for MR-only photon and proton therapy.

MATERIAL AND METHODS

12 prostate cancer patients scanned with an open 0.35T scanner were included. To obtain the full body contour an enhanced imaging protocol including two repeated scans after bilateral table movement was introduced. All required structures (patient contour, target and organ at risk) were delineated on a post-processed combined transversal image set (stitched MRI). The postprocessed MR was converted into a sCT by a pretrained neural network generator. Inversely planned photon and proton plans (VMAT and SFUD) were designed using the sCT and recalculated for rigidly and deformably registered CT images and compared based on D2%, D50%, V70Gy for organs at risk and based on D2%, D50%, D98% for the CTV and PTV. The stitched MRI and the untruncated MRI were compared to the CT, and the maximum surface distance was calculated. The sCT was evaluated with respect to delineation accuracy by comparing on stitched MRI and sCT using the DICE coefficient for femoral bones and the whole body.

RESULTS

Maximum surface distance analysis revealed uncertainties in lateral direction of 1-3mm on average. DICE coefficient analysis confirms good performance of the sCT conversion, i.e. 92%, 93%, and 100% were obtained for femoral bone left and right and whole body. Dose comparison resulted in uncertainties below 1% between deformed CT and sCT and below 2% between rigidly registered CT and sCT in the CTV for photon and proton treatment plans.

DISCUSSION

A newly developed acquisition protocol for open MR scanners and subsequent Sct generation revealed good acceptance for photon and proton therapy. Moreover, this protocol tackles the restriction of the limited FOVs and expands the capacities towards MR guided proton therapy with horizontal beam lines.

Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas

In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD_2%, ΔD_98%, and ΔD_mean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD_2%≥0.3%, ΔD_98%≥15.9%, ΔD_mean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD_2%≤0.4%, ΔD_98%≤5.6%, ΔD_mean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.

Clinical experience and workflow challenges with magnetic resonance-only radiation therapy simulation and planning for prostate cancer

BACKGROUND AND PURPOSE

Magnetic Resonance (MR)-only planning has been implemented clinically for radiotherapy of prostate cancer. However, fewer studies exist regarding the overall success rate of MR-only workflows. We report on successes and challenges of implementing MR-only workflows for prostate.

MATERIALS AND METHODS

A total of 585 patients with prostate cancer underwent an MR-only simulation and planning between 06/2016-06/2018. MR simulation included images for contouring, synthetic-CT generation and fiducial identification. Workflow interruptions occurred that required a backup CT, a re-simulation or an update to our current quality assurance (QA) process. The challenges were prospectively evaluated and classified into syn-CT generation, motion/artifacts in the MRs, fiducial QA and bowel preparation guidelines.

RESULTS

MR-only simulation was successful in 544 (93.2 %) patients. . In seventeen patients (2.9%), reconstruction of synthetic-CT failed due to patient size, femur angulation, or failure to determine the body contour. Twenty-four patients (4.1%) underwent a repeat/backup CT scan because of artifacts on the MR such as image blur due to patient motion or biopsy/surgical artifacts that hampered identification of the implanted fiducial markers. In patients requiring large coverage due to nodal involvement, inhomogeneity artifacts were resolved by using a two-stack acquisition and adaptive inhomogeneity correction. Bowel preparation guidelines were modified to address frequent rectum/gas issues due to longer MR scan time.

CONCLUSIONS

MR-only simulation has been successfully implemented for a majority of patients in the clinic. However, MR-CT or CT-only pathway may still be needed for patients where MR-only solution fails or patients with MR contraindications.

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.

Dosimetric evaluation of synthetic CT generated with GANs for MRI-only proton therapy treatment planning of brain tumors

PURPOSE

The purpose of this study was to address the dosimetric accuracy of synthetic computed tomography (sCT) images of patients with brain tumor generated using a modified generative adversarial network (GAN) method, for their use in magnetic resonance imaging (MRI)-only treatment planning for proton therapy.

METHODS

Dose volume histogram (DVH) analysis was performed on CT and sCT images of patients with brain tumor for plans generated for intensity-modulated proton therapy (IMPT). All plans were robustly optimized using a commercially available treatment planning system (RayStation, from RaySearch Laboratories) and standard robust parameters reported in the literature. The IMPT plan was then used to compute the dose on CT and sCT images for dosimetric comparison, using RayStation analytical (pencil beam) dose algorithm. We used a second, independent Monte Carlo dose calculation engine to recompute the dose on both CT and sCT images to ensure a proper analysis of the dosimetric accuracy of the sCT images.

RESULTS

The results extracted from RayStation showed excellent agreement for most DVH metrics computed on the CT and sCT for the nominal case, with a mean absolute difference below 0.5% (0.3 Gy) of the prescription dose for the clinical target volume (CTV) and below 2% (1.2 Gy) for the organs at risk (OARs) considered. This demonstrates a high dosimetric accuracy for the generated sCT images, especially in the target volume. The metrics obtained from the Monte Carlo doses mostly agreed with the values extracted from RayStation for the nominal and worst-case scenarios (mean difference below 3%).

CONCLUSIONS

This work demonstrated the feasibility of using sCT generated with a GAN-based deep learning method for MRI-only treatment planning of patients with brain tumor in intensity-modulated proton therapy.

Status and innovations in pre-treatment CT imaging for proton therapy

Pre-treatment CT imaging is a topic of growing importance in particle therapy. Improvements in the accuracy of stopping-power prediction are demanded to allow for a dose conformality that is not inferior to state-of-the-art image-guided photon therapy. Although range uncertainty has been kept practically constant over the last decades, recent technological and methodological developments, like the clinical application of dual-energy CT, have been introduced or arise at least on the horizon to improve the accuracy and precision of range prediction. This review gives an overview of the current status, summarizes the innovations in dual-energy CT and its potential impact on the field as well as potential alternative technologies for stopping-power prediction.

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.

Dosimetric evaluation of MR-derived synthetic-CTs for MR-only proton treatment planning

PURPOSE

To evaluate proton dose calculation accuracy of optimized pencil beam scanning (PBS) plans on MR-derived synthetic-CTs for prostate patients.

MATERIAL AND METHODS

Ten patient datasets with both a CT and an MRI were planned with opposed lateral proton beams optimized to single field uniform dose under an IRB-approved study. The proton plans were created on CT datasets generated by a commercial synthetic CT-based software called MRCAT (MR for Calculating ATtenuation) routinely used in our clinic for photon-based MR-only planning. A standard prescription of 79.2 Gy (RBE) and 68.4 Gy (RBE) was used for intact prostate and prostate bed cases, respectively. Proton plans were first generated and optimized using the MRCAT synthetic-CT (syn-CT), and then recalculated on the planning CT rigidly aligned with the syn-CT (aligned-CT) and a deformed planning CT (deformed-CT), which was deformed to match outer contour between syn-CT and aligned-CT. The same beam arrangement, total MUs, MUs/spot, spot positions were used to recalculate dose on deformed-CT and aligned-CT without renormalization. DVH analysis was performed on aligned-CT, deformed-CT, and syn-CT to compare D_98%, V_100%, V_95% for PTV, PTVeval, and GTV as well as V_70Gy, V_50Gy for OARs.

RESULTS

The relative percentage dose difference between syn-CT and deformed-CT, were (0.17 ± 0.33 %) for PTVeval D_98% and (0.07 ± 0.1 %) for CTV D_98%. Rectum V_70Gy, V_50Gy, and Bladder V_70Gy were (2.76 ± 4.01 %), (11.6 ± 11.2 %), and (3.41 ± 2.86 %), respectively for the syn-CT, and (3.23 ± 3.63 %), (11.3 ± 8.18 %), and (3.29 ± 2.76 %), respectively for the deformed-CT, and (1.37 ± 1.84 %), (8.48 ± 6.67 %), and (4.91 ± 3.65 %), respectively for aligned-CT.

CONCLUSION

Dosimetric analysis shows that MR-only proton planning is feasible using syn-CT based on current clinical margins that account for a range uncertainty.

MRI-based IMPT planning for prostate cancer

PURPOSE

Treatment planning for proton therapy requires the relative proton stopping power ratio (RSP) information of the patient for accurate dose calculations. RSP are conventionally obtained after mapping of the Hounsfield units (HU) from a calibrated patient computed tomography (CT). One or multiple CT are needed for a given treatment which represents additional, undesired dose to the patient. For prostate cancer, magnetic resonance imaging (MRI) scans are the gold standard for segmentation while offering dose-less imaging. We here quantify the clinical applicability of converted MR images as a substitute for intensity modulated proton therapy (IMPT) treatment of the prostate.

METHODS

MRCAT (Magnetic Resonance for Calculating ATtenuation) is a Philips-developed technology which produces a synthetic CT image consisting of five HU from a specific set of MRI acquisitions. MRCAT and original planning CT data sets were obtained for ten patients. An IMPT plan was generated on the MRCAT for each patient. Plans were produced such that they fulfill the prostate protocol in use at Massachusetts General Hospital (MGH). The plans were then recomputed onto the nominal planning CT for each patient. Robustness analyses (±5 mm setup shifts and ±3.5 % range uncertainties) were also performed.

RESULTS

Comparison of MRCAT plans and their recomputation onto the planning CT plan showed excellent agreement. Likewise, dose perturbations due to setup shifts and range uncertainties were well within clinical acceptance demonstrating the clinical viability of the approach.

CONCLUSIONS

This work demonstrate the clinical acceptability of substituting MR converted RSP images instead of CT for IMPT planning of prostate cancer. This further translates into higher contouring accuracy along with lesser imaging dose.

MRI-guided proton therapy planning: accounting for an inline MRI fringe field

MRI-guided proton therapy is being pursued for its promise to provide a more conformal, accurate proton therapy. However, the presence of the magnetic field imposes a challenge for the beam delivery as protons are deflected due to the Lorenz force. In this study, the impact of realistic inline MRI fringe field on IMPT plan delivery is investigated for a water phantom, liver tumor and prostate cancer differing in target volume, shape, and field configuration using Monte Carlo simulations. A method to correct for the shift of the beam spot positions in the presence of the inline magnetic field is presented. Results show that when not accounting for the effect of the magnetic field on the pencil beam delivery, the spot positions are substantially shifted and the quality of delivered plans is significantly deteriorated leading to dose inhomogeneities and creation of hot and cold spots. However, by correcting the pencil beam delivery, the dose quality of the IMPT plans is restored to a high degree. Nevertheless, adaptation of beam delivery alone is not robust regarding different treatment sites. By fully accounting during plan optimization for the dose distortions caused by the fringe and imaging fields, highly conformal IMPT plans are achieved. These results demonstrate proton pencil beam scanning and treatment planning can be adapted for precise delivery of state-of-the-art IMPT plans in MR-guided proton therapy in the presence of an inline MRI fringe field.

Evaluation of a deep learning-based pelvic synthetic CT generation technique for MRI-based prostate proton treatment planning

The purpose of this work is to validate the application of a deep learning-based method for pelvic synthetic CT (sCT) generation that can be used for prostate proton beam therapy treatment planning. We propose to integrate dense block minimization into 3D cycle-consistent generative adversarial networks (cycleGAN) framework to effectively learn the nonlinear mapping between MRI and CT pairs. A cohort of 17 patients with co-registered CT and MR pairs were used to test the deep learning-based sCT generation method by leave-one-out cross-validation. Image quality between the sCT and CT images, gamma analysis passing rate, dose-volume metrics, distal range displacement, and the individual pencil beam Bragg peak shift between sCT- and CT-based proton plans were evaluated. The average mean absolute error (MAE) was 51.32  ±  16.91 HU. The relative differences of the statistics of the PTV dose-volume histogram (DVH) metrics in between sCT and CT were generally less than 1%. Mean values of dose difference, absolute dose difference (in percent of the prescribed dose) were  -0.07%  ±  0.07% and 0.23%  ±  0.08%. Mean gamma analysis pass rate of 1 mm/1%, 2 mm/2%, 3 mm/3% criteria with 10% dose threshold were 92.39%  ±  5.97%, 97.95%  ±  2.95% and 98.97%  ±  1.62% respectively. The median, mean and standard deviation of absolute maximum range differences were 0.09 cm and 0.23  ±  0.25 cm. The median and mean Bragg peak shifts among the 17 patients were 0.09 cm and 0.18  ±  0.07 cm. The image similarity, dosimetric and distal range agreement between sCT and original CT suggests the feasibility of further development of an MRI-only workflow for prostate proton radiotherapy.

Evaluation of proton and photon dose distributions recalculated on 2D and 3D Unet-generated pseudoCTs from T1-weighted MR head scans

Introduction: The recent developments of magnetic resonance (MR) based adaptive strategies for photon and, potentially for proton therapy, require a fast and reliable conversion of MR images to X-ray computed tomography (CT) values. CT values are needed for photon and proton dose calculation. The improvement of conversion results employing a 3D deep learning approach is evaluated. Material and methods: A database of 89 T1-weighted MR head scans with about 100 slices each, including rigidly registered CTs, was created. Twenty-eight validation patients were randomly sampled, and four patients were selected for application. The remaining patients were used to train a 2D and a 3D U-shaped convolutional neural network (Unet). A stack size of 32 slices was used for 3D training. For all application cases, volumetric modulated arc therapy photon and single-field uniform dose pencil-beam scanning proton plans at four different gantry angles were optimized for a generic target on the CT and recalculated on 2D and 3D Unet-based pseudoCTs. Mean (absolute) error (MAE/ME) and a gradient sharpness estimate were used to quantify the image quality. Three-dimensional gamma and dose difference analyses were performed for photon (gamma criteria: 1%, 1 mm) and proton dose distributions (gamma criteria: 2%, 2 mm). Range (80% fall off) differences for beam's eye view profiles were evaluated for protons. Results: Training 36 h for 1000 epochs in 3D (6 h for 200 epochs in 2D) yielded a maximum MAE of 147 HU (135 HU) for the application patients. Except for one patient gamma pass rates for photon and proton dose distributions were above 96% for both Unets. Slice discontinuities were reduced for 3D training at the cost of sharpness. Conclusions: Image analysis revealed a slight advantage of 2D Unets compared to 3D Unets. Similar dose calculation performance was reached for the 2D and 3D network.

MRI-based treatment planning for proton radiotherapy: dosimetric validation of a deep learning-based liver synthetic CT generation method

Magnetic resonance imaging (MRI) has been widely used in combination with computed tomography (CT) radiation therapy because MRI improves the accuracy and reliability of target delineation due to its superior soft tissue contrast over CT. The MRI-only treatment process is currently an active field of research since it could eliminate systematic MR-CT co-registration errors, reduce medical cost, avoid diagnostic radiation exposure, and simplify clinical workflow. The purpose of this work is to validate the application of a deep learning-based method for abdominal synthetic CT (sCT) generation by image evaluation and dosimetric assessment in a commercial proton pencil beam treatment planning system (TPS). This study proposes to integrate dense block into a 3D cycle-consistent generative adversarial networks (cycle GAN) framework in an effort to effectively learn the nonlinear mapping between MRI and CT pairs. A cohort of 21 patients with co-registered CT and MR pairs were used to test the deep learning-based sCT image quality by leave-one-out cross validation. The CT image quality, dosimetric accuracy and the distal range fidelity were rigorously checked, using side-by-side comparison against the corresponding original CT images. The average mean absolute error (MAE) was 72.87  ±  18.16 HU. The relative differences of the statistics of the PTV dose volume histogram (DVH) metrics between sCT and CT were generally less than 1%. Mean 3D gamma analysis passing rate of 1 mm/1%, 2 mm/2%, 3 mm/3% criteria with 10% dose threshold were 90.76%  ±  5.94%, 96.98%  ±  2.93% and 99.37%  ±  0.99%, respectively. The median, mean and standard deviation of absolute maximum range differences were 0.170 cm, 0.186 cm and 0.155 cm. The image similarity, dosimetric and distal range agreement between sCT and original CT suggests the feasibility of further development of an MRI-only workflow for liver proton radiotherapy.

MRI-based treatment planning for liver stereotactic body radiotherapy: validation of a deep learning-based synthetic CT generation method

OBJECTIVE

The purpose of this work is to develop and validate a learning-based method to derive electron density from routine anatomical MRI for potential MRI-based SBRT treatment planning.

METHODS

We proposed to integrate dense block into cycle generative adversarial network (GAN) to effectively capture the relationship between the CT and MRI for CT synthesis. A cohort of 21 patients with co-registered CT and MR pairs were used to evaluate our proposed method by the leave-one-out cross-validation. Mean absolute error, peak signal-to-noise ratio and normalized cross-correlation were used to quantify the imaging differences between the synthetic CT (sCT) and CT. The accuracy of Hounsfield unit (HU) values in sCT for dose calculation was evaluated by comparing the dose distribution in sCT-based and CT-based treatment planning. Clinically relevant dose-volume histogram metrics were then extracted from the sCT-based and CT-based plans for quantitative comparison.

RESULTS

The mean absolute error, peak signal-to-noise ratio and normalized cross-correlation of the sCT were 72.87 ± 18.16 HU, 22.65 ± 3.63 dB and 0.92 ± 0.04, respectively. No significant differences were observed in the majority of the planning target volume and organ at risk dose-volume histogram metrics ( p  >  0.05). The average pass rate of γ analysis was over 99% with 1%/1 mm acceptance criteria on the coronal plane that intersects with isocenter.

CONCLUSION

The image similarity and dosimetric agreement between sCT and original CT warrant further development of an MRI-only workflow for liver stereotactic body radiation therapy.

ADVANCES IN KNOWLEDGE

This work is the first deep-learning-based approach to generating abdominal sCT through dense-cycle-GAN. This method can successfully generate the small bony structures such as the rib bones and is able to predict the HU values for dose calculation with comparable accuracy to reference CT images.

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.

Feasibility of MRI-only photon and proton dose calculations for pediatric patients with abdominal tumors

The purpose of this study was to develop a method enabling synthetic computed tomography (sCT) generation of the whole abdomen using magnetic resonance imaging (MRI) scans of pediatric patients with abdominal tumors. The proposed method relies on an automatic atlas-based segmentation of bone and lungs followed by an MRI intensity to synthetic Hounsfield unit conversion. Separate conversion algorithms were used for bone, lungs and soft-tissue. Rigidly registered CT and T2-weighted MR images of 30 patients in treatment position and with the same field of view were used for the evaluation of the atlas and the conversion algorithms. The dose calculation accuracy of the generated sCTs was verified for volumetric modulated arc therapy (VMAT) and pencil beam scanning (PBS). VMAT and PBS plans were robust optimized on an internal target volume (ITV) against a patient set-up uncertainty of 5 mm. Average differences between CT and sCT dose calculations for the ITV V _95% were 0.5% (min 0.0%; max 5.0%) and 0.0% (min  -0.1%; max 0.1%) for VMAT and PBS dose distributions, respectively. Average differences for the mean dose to the organs at risk were  <0.2% (min  -0.6%; max 1.2%) and  <0.2% (min  -2.0%; max 2.6%) for VMAT and PBS dose distributions, respectively. Results show that MRI-only photon and proton dose calculations are feasible for children with abdominal tumors.

MRI-guidance for motion management in external beam radiotherapy: current status and future challenges

High precision conformal radiotherapy requires sophisticated imaging techniques to aid in target localisation for planning and treatment, particularly when organ motion due to respiration is involved. X-ray based imaging is a well-established standard for radiotherapy treatments. Over the last few years, the ability of magnetic resonance imaging (MRI) to provide radiation-free images with high-resolution and superb soft tissue contrast has highlighted the potential of this imaging modality for radiotherapy treatment planning and motion management. In addition, these advantageous properties motivated several recent developments towards combined MRI radiation therapy treatment units, enabling in-room MRI-guidance and treatment adaptation. The aim of this review is to provide an overview of the state-of-the-art in MRI-based image guidance for organ motion management in external beam radiotherapy. Methodological aspects of MRI for organ motion management are reviewed and their application in treatment planning, in-room guidance and adaptive radiotherapy described. Finally, a roadmap for an optimal use of MRI-guidance is highlighted and future challenges are discussed.

Two-dimensional noise reconstruction in proton computed tomography using distance-driven filtered back-projection of simulated projections

We present a formalism for two-dimensional (2D) noise reconstruction in proton computed tomography (pCT). This is necessary for the application of fluence modulated pCT (FMpCT) since it permits image noise prescription and the corresponding proton fuence optimization. We aimed at extending previously published formalisms to account for the impact of multiple Coulomb scattering (MCS) on projection noise, and the use of filtered back projection (FBP) reconstruction along curved paths with distance driven binning (DDB). 2D noise reconstruction for a beam of protons with parallel initial momentum vectors, and for projections binned both at the rear tracker and with DDB, was established. Monte Carlo (MC) simulations of pCT scans of a water cylinder were employed to generate pCT projections and to calculate their noise for use in 2D noise reconstruction. These were compared to results from an analytical model accounting for MCS for rear tracker binning as well as against the previously published central pixel model which ignores MCS. Image noise reconstructed with the formalism for rear tracker binning and DDB were compared to MC results using annular regions of interest (ROIs). Agreement better than 8% was obtained between the noise of projections calculated with MC simulation and our model. Noise from annular ROIs agreed with our noise reconstructions for rear tracker binning and DDB. The central pixel model ignoring MCS underestimated projection and thus image noise by up to 40% towards the object's edge. The use of DDB decreased the image noise towards the object's edge when compared to rear tracker binning and yielded more uniform noise throughout the image. MCS should not be neglected when predicting image noise for pixels away from the center of an object in a pCT scan due to the increasing influence of the gradient of the object's hull closer to the edges.

Magnetic Resonance Imaging only Workflow for Radiotherapy Simulation and Planning in Prostate Cancer

Magnetic resonance imaging (MRI) is often combined with computed tomography (CT) in prostate radiotherapy to optimise delineation of the target and organs-at-risk (OAR) while maintaining accurate dose calculation. Such a dual-modality workflow requires two separate imaging sessions, and it has some fundamental and logistical drawbacks. Due to the availability of new MRI hardware and software solutions, CT examinations can be omitted for prostate radiotherapy simulations. All information for treatment planning, including electron density maps and bony anatomy, can nowadays be obtained with MRI. Such an MRI-only simulation workflow reduces delineation ambiguities, eases planning logistics, and improves patient comfort; however, careful validation of the complete MRI-only workflow is warranted. The first institutes are now adopting this MRI-only workflow for prostate radiotherapy. In this article, we will review technology and workflow requirements for an MRI-only prostate simulation workflow.

Dose evaluation of fast synthetic-CT generation using a generative adversarial network for general pelvis MR-only radiotherapy

To enable magnetic resonance (MR)-only radiotherapy and facilitate modelling of radiation attenuation in humans, synthetic CT (sCT) images need to be generated. Considering the application of MR-guided radiotherapy and online adaptive replanning, sCT generation should occur within minutes. This work aims at assessing whether an existing deep learning network can rapidly generate sCT images for accurate MR-based dose calculations in the entire pelvis. A study was conducted on data of 91 patients with prostate (59), rectal (18) and cervical (14) cancer who underwent external beam radiotherapy acquiring both CT and MRI for patients' simulation. Dixon reconstructed water, fat and in-phase images obtained from a conventional dual gradient-recalled echo sequence were used to generate sCT images. A conditional generative adversarial network (cGAN) was trained in a paired fashion on 2D transverse slices of 32 prostate cancer patients. The trained network was tested on the remaining patients to generate sCT images. For 30 patients in the test set, dose recalculations of the clinical plan were performed on sCT images. Dose distributions were evaluated comparing voxel-based dose differences, gamma and dose-volume histogram (DVH) analysis. The sCT generation required 5.6 s and 21 s for a single patient volume on a GPU and CPU, respectively. On average, sCT images resulted in a higher dose to the target of maximum 0.3%. The average gamma pass rates using the 3%, 3 mm and 2%, 2 mm criteria were above 97 and 91%, respectively, for all volumes of interests considered. All DVH points calculated on sCT differed less than  ±2.5% from the corresponding points on CT. Results suggest that accurate MR-based dose calculation using sCT images generated with a cGAN trained on prostate cancer patients is feasible for the entire pelvis. The sCT generation was sufficiently fast for integration in an MR-guided radiotherapy workflow.

MR-Only Brain Radiation Therapy: Dosimetric Evaluation of Synthetic CTs Generated by a Dilated Convolutional Neural Network

PURPOSE

This work aims to facilitate a fast magnetic resonance (MR)-only workflow for radiation therapy of intracranial tumors. Here, we evaluate whether synthetic computed tomography (sCT) images generated with a dilated convolutional neural network (CNN) enable accurate MR-based dose calculations in the brain.

METHODS AND MATERIALS

We conducted a retrospective study of 52 patients with brain tumors who underwent both computed tomography (CT) and MR imaging for radiation therapy treatment planning. To generate the sCTs, a T1-weighted gradient echo MR sequence was selected from the clinical protocol for multiple types of brain tumors. sCTs were created for all 52 patients with a dilated CNN using 2-fold cross validation; in each fold, 26 patients were used for training and the remaining 26 patients were used for evaluation. For each patient, the clinical CT-based treatment plan was recalculated on sCT. We calculated dose differences and gamma pass rates between CT- and sCT-based plans inside body and planning target volume. Geometric fidelity of the sCT and differences in beam depth and equivalent path length were assessed between both treatment plans.

RESULTS

sCT generation took 1 minute per patient. Over the patient population, the mean absolute error of the sCT within the intersection of body contours was 67 ± 11 HU (±1 standard deviation [SD], range: 51-117 HU), and the mean error was 13 ± 9 HU (±1 SD, range: -2 to 38 HU). Dosimetric analysis showed mean deviations of 0.00% ± 0.02% (±1 SD, range: -0.05 to 0.03) for dose within the body contours and -0.13% ± 0.39% (±1 SD, range: -1.43 to 0.80) inside the planning target volume. Mean γ_1mm/1% was 98.8% ± 2.2% for doses >50% of the prescribed dose.

CONCLUSIONS

The presented dilated CNN generated sCTs from conventional MR images without adding scan time to the acquisition. Dosimetric evaluation suggests that dose calculations performed on the sCTs are accurate and can therefore be used for MR-only intracranial radiation therapy treatment planning.

Image guidance in proton therapy for lung cancer

Proton therapy is a promising but challenging treatment modality for the management of lung cancer. The technical challenges are due to respiratory motion, low dose tolerance of adjacent normal tissue and tissue density heterogeneity. Different imaging modalities are applied at various steps of lung proton therapy to provide information on target definition, target motion, proton range, patient setup and treatment outcome assessment. Imaging data is used to guide treatment design, treatment delivery, and treatment adaptation to ensure the treatment goal is achieved. This review article will summarize and compare various imaging techniques that can be used in every step of lung proton therapy to address these challenges.