Cancer Informatics for Cancer Centers: Sharing Ideas on How to Build an Artificial Intelligence-Ready Informatics Ecosystem for Radiation Oncology.

In August 2022, the Cancer Informatics for Cancer Centers brought together cancer informatics leaders for its biannual symposium, Precision Medicine Applications in Radiation Oncology, co-chaired by Quynh-Thu Le, MD (Stanford University), and Walter J. Curran, MD (GenesisCare). Over the course of 3 days, presenters discussed a range of topics relevant to radiation oncology and the cancer informatics community more broadly, including biomarker development, decision support algorithms, novel imaging tools, theranostics, and artificial intelligence (AI) for the radiotherapy workflow. Since the symposium, there has been an impressive shift in the promise and potential for integration of AI in clinical care, accelerated in large part by major advances in generative AI. AI is now poised more than ever to revolutionize cancer care. Radiation oncology is a field that uses and generates a large amount of digital data and is therefore likely to be one of the first fields to be transformed by AI. As experts in the collection, management, and analysis of these data, the informatics community will take a leading role in ensuring that radiation oncology is prepared to take full advantage of these technological advances. In this report, we provide highlights from the symposium, which took place in Santa Barbara, California, from August 29 to 31, 2022. We discuss lessons learned from the symposium for data acquisition, management, representation, and sharing, and put these themes into context to prepare radiation oncology for the successful and safe integration of AI and informatics technologies.

Volumetric Modulated Arc Therapy Enabled Total Body Irradiation (VMAT-TBI): Six-year Clinical Experience and Treatment Outcomes.

Total body irradiation is an important part of the conditioning regimens frequently used to prepare patients for allogeneic hematopoietic stem cell transplantation (SCT). Volumetric-modulated arc therapy enabled total body irradiation (VMAT-TBI), an alternative to conventional TBI (cTBI), is a novel radiotherapy treatment technique that has been implemented and investigated in our institution. The purpose of this study is to (1) report our six-year clinical experience in terms of treatment planning strategy and delivery time and (2) evaluate the clinical outcomes and toxicities in our cohort of patients treated with VMAT-TBI. This is a retrospective single center study. Forty-four patients at our institution received VMAT-TBI and chemotherapy conditioning followed by allogeneic SCT between 2014 and 2020. Thirty-two patients (73%) received standard-dose TBI (12-13.2 Gy in 6-8 fractions twice daily), whereas 12 (27%) received low-dose TBI (2-4 Gy in one fraction). Treatment planning, delivery, and treatment outcome data including overall survival (OS), relapse-free survival (RFS), and toxicities were analyzed. The developed VMAT-TBI planning strategy consistently generated plans satisfying our dose constraints, with planning target volume coverage >90%, mean lung dose ∼50% to 75% of prescription dose, and minimal hotspots in critical organs. Most of the treatment deliveries were <100 minutes (range 33-147, mean 72). The median follow-up was 26 months. At the last follow-up, 34 of 44 (77%) of patients were alive, with 1- and 2-year OS of 90% and 79% and RFS of 88% and 71%, respectively. The most common grade 3+ toxicities observed were mucositis (31 patients [71%]) and nephrotoxicity (6 patients [13%]), both of which were deemed multifactorial in cause. Four patients (9%) in standard-dose cohort developed grade 3+ pneumonitis, with 3 cases in the setting of documented respiratory infection and only 1 (2%) deemed likely related to radiation alone. VMAT-TBI provides a safe alternative to cTBI. The dose modulation capability of VMAT-TBI may lead to new treatment strategies, such as simultaneous boost and further critical organ sparing, for better malignant cell eradication, immune suppression, and lower toxicities.

Semi-automatic sigmoid colon segmentation in CT for radiation therapy treatment planning via an iterative 2.5-D deep learning approach.

Automatic sigmoid colon segmentation in CT for radiotherapy treatment planning is challenging due to complex organ shape, close distances to other organs, and large variations in size, shape, and filling status. The patient bowel is often not evacuated, and CT contrast enhancement is not used, which further increase problem difficulty. Deep learning (DL) has demonstrated its power in many segmentation problems. However, standard 2-D approaches cannot handle the sigmoid segmentation problem due to incomplete geometry information and 3-D approaches often encounters the challenge of a limited training data size. Motivated by human's behavior that segments the sigmoid slice by slice while considering connectivity between adjacent slices, we proposed an iterative 2.5-D DL approach to solve this problem. We constructed a network that took an axial CT slice, the sigmoid mask in this slice, and an adjacent CT slice to segment as input and output the predicted mask on the adjacent slice. We also considered other organ masks as prior information. We trained the iterative network with 50 patient cases using five-fold cross validation. The trained network was repeatedly applied to generate masks slice by slice. The method achieved average Dice similarity coefficients of 0.82 0.06 and 0.88 0.02 in 10 test cases without and with using prior information.

On the robustness of deep learning-based lung-nodule classification for CT images with respect to image noise.

Robustness is an important aspect when evaluating a method of medical image analysis. In this study, we investigated the robustness of a deep learning (DL)-based lung-nodule classification model for CT images with respect to noise perturbations. A deep neural network (DNN) was established to classify 3D CT images of lung nodules into malignant or benign groups. The established DNN was able to predict malignancy rate of lung nodules based on CT images, achieving the area under the curve of 0.91 for the testing dataset in a tenfold cross validation as compared to radiologists' prediction. We then evaluated its robustness against noise perturbations. We added to the input CT images noise signals generated randomly or via an optimization scheme using a realistic noise model based on a noise power spectrum for a given mAs level, and monitored the DNN's output. The results showed that the CT noise was able to affect the prediction results of the established DNN model. With random noise perturbations at 100 mAs, DNN's predictions for 11.2% of training data and 17.4% of testing data were successfully altered by at least once. The percentage increased to 23.4% and 34.3%, respectively, for optimization-based perturbations. We further evaluated robustness of models with different architectures, parameters, number of output labels, etc, and robustness concern was found in these models to different degrees. To improve model robustness, we empirically proposed an adaptive training scheme. It fine-tuned the DNN model by including perturbations in the training dataset that successfully altered the DNN's perturbations. The adaptive scheme was repeatedly performed to gradually improve DNN's robustness. The numbers of perturbations at 100 mAs affecting DNN's predictions were reduced to 10.8% for training and 21.1% for testing by the adaptive training scheme after two iterations. Our study illustrated that robustness may potentially be a concern for an exemplary DL-based lung-nodule classification model for CT images, indicating the needs for evaluating and ensuring model robustness when developing similar models. The proposed adaptive training scheme may be able to improve model robustness.

Operating a treatment planning system using a deep-reinforcement learning-based virtual treatment planner for prostate cancer intensity-modulated radiation therapy treatment planning.

PURPOSE: In the treatment planning process of intensity-modulated radiation therapy (IMRT), a human planner operates the treatment planning system (TPS) to adjust treatment planning parameters, for example, dose volume histogram (DVH) constraints' locations and weights, to achieve a satisfactory plan for each patient. This process is usually time-consuming, and the plan quality depends on planer's experience and available planning time. In this study, we proposed to model the behaviors of human planners in treatment planning by a deep reinforcement learning (DRL)-based virtual treatment planner network (VTPN), such that it can operate the TPS in a human-like manner for treatment planning. METHODS AND MATERIALS: Using prostate cancer IMRT as an example, we established the VTPN using a deep neural network developed. We considered an in-house optimization engine with a weighted quadratic objective function. Virtual treatment planner network was designed to observe an intermediate plan DVHs and decide the action to improve the plan by changing weights and threshold dose in the objective function. We trained the VTPN in an end-to-end DRL process in 10 patient cases. A plan score was used to measure plan quality. We demonstrated the feasibility and effectiveness of the trained VTPN in another 64 patient cases. RESULTS: Virtual treatment planner network was trained to spontaneously learn how to adjust treatment planning parameters to generate high-quality treatment plans. In the 64 testing cases, with initialized parameters, quality score was 4.97 (±2.02), with 9.0 being the highest possible score. Using VTPN to perform treatment planning improved quality score to 8.44 (±0.48). CONCLUSIONS: To our knowledge, this was the first time that intelligent treatment planning behaviors of human planner in external beam IMRT are autonomously encoded in an artificial intelligence system. The trained VTPN is capable of behaving in a human-like way to produce high-quality plans.

An introduction to deep learning in medical physics: advantages, potential, and challenges.

As one of the most popular approaches in artificial intelligence, deep learning (DL) has attracted a lot of attention in the medical physics field over the past few years. The goals of this topical review article are twofold. First, we will provide an overview of the method to medical physics researchers interested in DL to help them start the endeavor. Second, we will give in-depth discussions on the DL technology to make researchers aware of its potential challenges and possible solutions. As such, we divide the article into two major parts. The first part introduces general concepts and principles of DL and summarizes major research resources, such as computational tools and databases. The second part discusses challenges faced by DL, present available methods to mitigate some of these challenges, as well as our recommendations.

Beam modeling and beam model commissioning for Monte Carlo dose calculation-based radiation therapy treatment planning: Report of AAPM Task Group 157.

Dose calculation plays an important role in the accuracy of radiotherapy treatment planning and beam delivery. The Monte Carlo (MC) method is capable of achieving the highest accuracy in radiotherapy dose calculation and has been implemented in many commercial systems for radiotherapy treatment planning. The objective of this task group was to assist clinical physicists with the potentially complex task of acceptance testing and commissioning MC-based treatment planning systems (TPS) for photon and electron beam dose calculations. This report provides an overview on the general approach of clinical implementation and testing of MC-based TPS with a specific focus on models of clinical photon and electron beams. Different types of beam models are described including those that utilize MC simulation of the treatment head and those that rely on analytical methods and measurements. The trade-off between accuracy and efficiency in the various source-modeling approaches is discussed together with guidelines for acceptance testing of MC-based TPS from the clinical standpoint. Specific recommendations are given on methods and practical procedures to commission clinical beam models for MC-based TPS.

Intelligent inverse treatment planning via deep reinforcement learning, a proof-of-principle study in high dose-rate brachytherapy for cervical cancer.

Inverse treatment planning in radiation therapy is formulated as solving optimization problems. The objective function and constraints consist of multiple terms designed for different clinical and practical considerations. Weighting factors of these terms are needed to define the optimization problem. While a treatment planning optimization engine can solve the optimization problem with given weights, adjusting the weights to yield a high-quality plan is typically performed by a human planner. Yet the weight-tuning task is labor intensive, time consuming, and it critically affects the final plan quality. An automatic weight-tuning approach is strongly desired. The procedure of weight adjustment to improve the plan quality is essentially a decision-making problem. Motivated by the tremendous success in deep learning for decision making with human-level intelligence, we propose a novel framework to adjust the weights in a human-like manner. This study used inverse treatment planning in high-dose-rate brachytherapy (HDRBT) for cervical cancer as an example. We developed a weight-tuning policy network (WTPN) that observes dose volume histograms of a plan and outputs an action to adjust organ weighting factors, similar to the behaviors of a human planner. We trained the WTPN via end-to-end deep reinforcement learning. Experience replay was performed with the epsilon greedy algorithm. After training was completed, we applied the trained WTPN to guide treatment planning of five testing patient cases. It was found that the trained WTPN successfully learnt the treatment planning goals and was able to guide the weight tuning process. On average, the quality score of plans generated under the WTPN's guidance was improved by ~8.5% compared to the initial plan with arbitrarily set weights, and by 10.7% compared to the plans generated by human planners. To our knowledge, this was the first time that a tool was developed to adjust organ weights for the treatment planning optimization problem in a human-like fashion based on intelligence learnt from a training process, which was different from existing strategies based on pre-defined rules. The study demonstrated potential feasibility to develop intelligent treatment planning approaches via deep reinforcement learning.

Intelligent Parameter Tuning in Optimization-Based Iterative CT Reconstruction via Deep Reinforcement Learning.

A number of image-processing problems can be formulated as optimization problems. The objective function typically contains several terms specifically designed for different purposes. Parameters in front of these terms are used to control the relative importance among them. It is of critical importance to adjust these parameters, as quality of the solution depends on their values. Tuning parameters are a relatively straight forward task for a human, as one can intuitively determine the direction of parameter adjustment based on the solution quality. Yet manual parameter tuning is not only tedious in many cases, but also becomes impractical when a number of parameters exist in a problem. Aiming at solving this problem, this paper proposes an approach that employs deep reinforcement learning to train a system that can automatically adjust parameters in a human-like manner. We demonstrate our idea in an example problem of optimization-based iterative computed tomography (CT) reconstruction with a pixel-wise total-variation regularization term. We set up a parameter-tuning policy network (PTPN), which maps a CT image patch to an output that specifies the direction and amplitude by which the parameter at the patch center is adjusted. We train the PTPN via an end-to-end reinforcement learning procedure. We demonstrate that under the guidance of the trained PTPN, reconstructed CT images attain quality similar or better than those reconstructed with manually tuned parameters.

MR-CBCT image-guided system for radiotherapy of orthotopic rat prostate tumors.

Multi-modality image-guided radiotherapy is the standard of care in contemporary cancer management; however, it is not common in preclinical settings due to both hardware and software limitations. Soft tissue lesions, such as orthotopic prostate tumors, are difficult to identify using cone beam computed tomography (CBCT) imaging alone. In this study, we characterized a research magnetic resonance (MR) scanner for preclinical studies and created a protocol for combined MR-CBCT image-guided small animal radiotherapy. Two in-house dual-modality, MR and CBCT compatible, phantoms were designed and manufactured using 3D printing technology. The phantoms were used for quality assurance tests and to facilitate end-to-end testing for combined preclinical MR and CBCT based treatment planning. MR and CBCT images of the phantoms were acquired utilizing a Varian 4.7 T scanner and XRad-225Cx irradiator, respectively. The geometry distortion was assessed by comparing MR images to phantom blueprints and CBCT. The corrected MR scans were co-registered with CBCT and subsequently used for treatment planning. The fidelity of 3D printed phantoms compared to the blueprint design yielded favorable agreement as verified with the CBCT measurements. The geometric distortion, which varied between -5% and 11% throughout the scanning volume, was substantially reduced to within 0.4% after correction. The distortion free MR images were co-registered with the corresponding CBCT images and imported into a commercial treatment planning software SmART Plan. The planning target volume (PTV) was on average 19% smaller when contoured on the corrected MR-CBCT images relative to raw images without distortion correction. An MR-CBCT based preclinical workflow was successfully designed and implemented for small animal radiotherapy. Combined MR-CBCT image-guided radiotherapy for preclinical research potentially delivers enhanced relevance to human radiotherapy for various disease sites. This novel protocol is wide-ranging and not limited to the orthotopic prostate tumor study presented in the study.

Real-Time Whole-Brain Radiation Therapy: A Single-Institution Experience.

PURPOSE: To demonstrate the feasibility of a real-time whole-brain radiation therapy (WBRT) workflow, taking advantage of contemporary radiation therapy capabilities and seeking to optimize clinical workflow for WBRT. METHODS AND MATERIALS: We developed a method incorporating the linear accelerator's on-board imaging system for patient simulation, used cone-beam computed tomography (CBCT) data for treatment planning, and delivered the first fraction of prescribed therapy, all during the patient's initial appointment. Simulation was performed in the linear accelerator vault. An acquired CBCT data set was used for scripted treatment planning protocol, providing inversely planned, automated treatment plan generation. The osseous boundaries of the brain were auto-contoured to create a target volume. Two parallel-opposed beams using field-in-field intensity modulate radiation therapy covered this target to the user-defined inferior level (C1 or C2). The method was commissioned using an anthropomorphic head phantom and verified using 100 clinically treated patients. RESULTS: Whole-brain target heterogeneity was within 95%-107% of the prescription dose, and target coverage compared favorably to standard, manually created 3-dimensional plans. For the commissioning CBCT datasets, the secondary monitor unit verification and independent 3-dimensional dose distribution comparison for computed and delivered doses were within 2% agreement relative to the scripted auto-plans. On average, time needed to complete the entire process was 35.1 ± 10.3 minutes from CBCT start to last beam delivered. CONCLUSIONS: The real-time WBRT workflow using integrated on-site imaging, planning, quality assurance, and delivery was tested and deemed clinically feasible. The design necessitates a synchronized team consisting of physician, physicist, dosimetrist, and therapists. This work serves as a proof of concept of real-time planning and delivery for other treatment sites.

Prototype volumetric ultrasound tomography image guidance system for prone stereotactic partial breast irradiation: proof-of-concept.

Accurate dose delivery in stereotactic partial breast irradiation (S-PBI) is challenging because of the target position uncertainty caused by breast deformation, the target volume changes caused by lumpectomy cavity shrinkage, and the target delineation uncertainty on simulation computed tomography (CT) images caused by poor soft tissue contrast. We have developed a volumetric ultrasound tomography (UST) image guidance system for prone position S-PBI. The system is composed of a novel 3D printed rotation water tank, a patient-specific resin breast immobilization cup, and a 1D array ultrasound transducer. Coronal 2D US images were acquired in 5° increments over a 360° range, and planes were acquired every 2 mm in elevation. A super-compounding technique was used to reconstruct the image volume. The image quality of UST was evaluated with a BB-1 breast phantom and BioZorb surgical marker, and the results revealed that UST offered better soft tissue contrast than CT and similar image quality to MR. In the evaluated plane, the size and location of five embedded objects were measured and compared to MR, which is considered as the ground truth. Objects' diameters and the distances between objects in UST differ by approximately 1 to 2 mm from those in MR, which showed that UST offers the image quality required for S-PBI. In future work we will develop a robotic system that will be ultimately implemented in the clinic.

Full Monte Carlo-Based Biologic Treatment Plan Optimization System for Intensity Modulated Carbon Ion Therapy on Graphics Processing Unit.

PURPOSE: One of the major benefits of carbon ion therapy is enhanced biological effectiveness at the Bragg peak region. For intensity modulated carbon ion therapy (IMCT), it is desirable to use Monte Carlo (MC) methods to compute the properties of each pencil beam spot for treatment planning, because of their accuracy in modeling physics processes and estimating biological effects. We previously developed goCMC, a graphics processing unit (GPU)-oriented MC engine for carbon ion therapy. The purpose of the present study was to build a biological treatment plan optimization system using goCMC. METHODS AND MATERIALS: The repair-misrepair-fixation model was implemented to compute the spatial distribution of linear-quadratic model parameters for each spot. A treatment plan optimization module was developed to minimize the difference between the prescribed and actual biological effect. We used a gradient-based algorithm to solve the optimization problem. The system was embedded in the Varian Eclipse treatment planning system under a client-server architecture to achieve a user-friendly planning environment. We tested the system with a 1-dimensional homogeneous water case and 3 3-dimensional patient cases. RESULTS: Our system generated treatment plans with biological spread-out Bragg peaks covering the targeted regions and sparing critical structures. Using 4 NVidia GTX 1080 GPUs, the total computation time, including spot simulation, optimization, and final dose calculation, was 0.6 hour for the prostate case (8282 spots), 0.2 hour for the pancreas case (3795 spots), and 0.3 hour for the brain case (6724 spots). The computation time was dominated by MC spot simulation. CONCLUSIONS: We built a biological treatment plan optimization system for IMCT that performs simulations using a fast MC engine, goCMC. To the best of our knowledge, this is the first time that full MC-based IMCT inverse planning has been achieved in a clinically viable time frame.

A deep convolutional neural network-based automatic delineation strategy for multiple brain metastases stereotactic radiosurgery.

Accurate and automatic brain metastases target delineation is a key step for efficient and effective stereotactic radiosurgery (SRS) treatment planning. In this work, we developed a deep learning convolutional neural network (CNN) algorithm for segmenting brain metastases on contrast-enhanced T1-weighted magnetic resonance imaging (MRI) datasets. We integrated the CNN-based algorithm into an automatic brain metastases segmentation workflow and validated on both Multimodal Brain Tumor Image Segmentation challenge (BRATS) data and clinical patients' data. Validation on BRATS data yielded average DICE coefficients (DCs) of 0.75±0.07 in the tumor core and 0.81±0.04 in the enhancing tumor, which outperformed most techniques in the 2015 BRATS challenge. Segmentation results of patient cases showed an average of DCs 0.67±0.03 and achieved an area under the receiver operating characteristic curve of 0.98±0.01. The developed automatic segmentation strategy surpasses current benchmark levels and offers a promising tool for SRS treatment planning for multiple brain metastases.

Data correlation based noise level estimation for cone beam projection data.

BACKGROUND: In regularized iterative reconstruction algorithms, the selection of regularization parameter depends on the noise level of cone beam projection data. OBJECTIVE: Our aim is to propose an algorithm to estimate the noise level of cone beam projection data. METHODS: We first derived the data correlation of cone beam projection data in the Fourier domain, based on which, the signal and the noise were decoupled. Then the noise was extracted and averaged for estimation. An adaptive regularization parameter selection strategy was introduced based on the estimated noise level. Simulation and real data studies were conducted for performance validation. RESULTS: There exists an approximately zero-energy double-wedge area in the 3D Fourier domain of cone beam projection data. As for the noise level estimation results, the averaged relative errors of the proposed algorithm in the analytical/MC/spotlight-mode simulation experiments were 0.8%, 0.14% and 0.24%, respectively, and outperformed the homogeneous area based as well as the transformation based algorithms. Real studies indicated that the estimated noise levels were inversely proportional to the exposure levels, i.e., the slopes in the log-log plot were -1.0197 and -1.049 with respect to the short-scan and half-fan modes. The introduced regularization parameter selection strategy could deliver promising reconstructed image qualities. CONCLUSIONS: Based on the data correlation of cone beam projection data in Fourier domain, the proposed algorithm could estimate the noise level of cone beam projection data accurately and robustly. The estimated noise level could be used to adaptively select the regularization parameter.

Multistage stereotactic radiosurgery for large cerebral arteriovenous malformations using the Gamma Knife platform.

PURPOSE: Radiosurgery is an established technique to treat cerebral arteriovenous malformations (AVMs). Obliteration of larger AVMs (> 10-15 cm3 or diameter > 3 cm) in a single session is challenging with current radiosurgery platforms due to toxicity. We present a novel technique of multistage stereotactic radiosurgery (SRS) for large intracranial arteriovenous malformations (AVM) using the Gamma Knife system. MATERIALS/METHODS: Eighteen patients with large (> 10-15 cm3 or diameter > 3 cm) AVMs, which were previously treated using a staged SRS technique on the Cyberknife platform, were retrospectively selected for this study. The AVMs were contoured and divided into 3-8 subtargets to be treated sequentially in a staged approach at half to 4 week intervals. The prescription dose ranged from 15 Gy to 20 Gy, depending on the subtarget number, volume, and location. Gamma Knife plans using multiple collimator settings were generated and optimized. The coordinates of each shot from the initial plan covering the total AVM target were extracted based on their relative positions within the frame system. The shots were regrouped based on their location with respect to the subtarget contours to generate subplans for each stage. The delivery time of each shot for a subtarget was decay corrected with 60 Co for staging the treatment course to generate the same dose distribution as that planned for the total AVM target. Conformality indices and dose-volume analysis were performed to evaluate treatment plans. RESULTS: With the shot redistribution technique, the composite dose for the multistaged treatment of multiple subtargets is equivalent to the initial plan for total AVM target. Gamma Knife plans resulted in an average PTV coverage of 96.3 ± 0.9% and a PITV of 1.23 ± 0.1. The resulting Conformality indices, V12Gy and R50 dose spillage values were 0.76 ± 0.05, 3.4 ± 1.8, and 3.1 ± 0.5 respectively. CONCLUSION: The Gamma Knife system can deliver a multistaged conformal dose to treat large AVMs when correcting for translational setup errors of each shot at each staged treatment.

A novel geometry-dosimetry label fusion method in multi-atlas segmentation for radiotherapy: a proof-of-concept study.

Multi-atlas segmentation (MAS) has been widely used to automate the delineation of organs at risk (OARs) for radiotherapy. Label fusion is a crucial step in MAS to cope with the segmentation variabilities among multiple atlases. However, most existing label fusion methods do not consider the potential dosimetric impact of the segmentation result. In this proof-of-concept study, we propose a novel geometry-dosimetry label fusion method for MAS-based OAR auto-contouring, which evaluates the segmentation performance in terms of both geometric accuracy and the dosimetric impact of the segmentation accuracy on the resulting treatment plan. Differently from the original selective and iterative method for performance level estimation (SIMPLE), we evaluated and rejected the atlases based on both Dice similarity coefficient and the predicted error of the dosimetric endpoints. The dosimetric error was predicted using our previously developed geometry-dosimetry model. We tested our method in MAS-based rectum auto-contouring on 20 prostate cancer patients. The accuracy in the rectum sub-volume close to the planning tumor volume (PTV), which was found to be a dosimetric sensitive region of the rectum, was greatly improved. The mean absolute distance between the obtained contour and the physician-drawn contour in the rectum sub-volume 2 mm away from PTV was reduced from 3.96 mm to 3.36 mm on average for the 20 patients, with the maximum decrease found to be from 9.22 mm to 3.75 mm. We also compared the dosimetric endpoints predicted for the obtained contours with those predicted for the physician-drawn contours. Our method led to smaller dosimetric endpoint errors than the SIMPLE method in 15 patients, comparable errors in 2 patients, and slightly larger errors in 3 patients. These results indicated the efficacy of our method in terms of considering both geometric accuracy and dosimetric impact during label fusion. Our algorithm can be applied to different tumor sites and radiation treatments, given a specifically trained geometry-dosimetry model.

New concept on an integrated interior magnetic resonance imaging and medical linear accelerator system for radiation therapy.

Image guidance plays a critical role in radiotherapy. Currently, cone-beam computed tomography (CBCT) is routinely used in clinics for this purpose. While this modality can provide an attenuation image for therapeutic planning, low soft-tissue contrast affects the delineation of anatomical and pathological features. Efforts have recently been devoted to several MRI linear accelerator (LINAC) projects that lead to the successful combination of a full diagnostic MRI scanner with a radiotherapy machine. We present a new concept for the development of the MRI-LINAC system. Instead of combining a full MRI scanner with the LINAC platform, we propose using an interior MRI (iMRI) approach to image a specific region of interest (RoI) containing the radiation treatment target. While the conventional CBCT component still delivers a global image of the patient's anatomy, the iMRI offers local imaging of high soft-tissue contrast for tumor delineation. We describe a top-level system design for the integration of an iMRI component into an existing LINAC platform. We performed numerical analyses of the magnetic field for the iMRI to show potentially acceptable field properties in a spherical RoI with a diameter of 15 cm. This field could be shielded to a sufficiently low level around the LINAC region to avoid electromagnetic interference. Furthermore, we investigate the dosimetric impacts of this integration on the radiotherapy beam.

Accelerated Monte Carlo simulation on the chemical stage in water radiolysis using GPU.

The accurate simulation of water radiolysis is an important step to understand the mechanisms of radiobiology and quantitatively test some hypotheses regarding radiobiological effects. However, the simulation of water radiolysis is highly time consuming, taking hours or even days to be completed by a conventional CPU processor. This time limitation hinders cell-level simulations for a number of research studies. We recently initiated efforts to develop gMicroMC, a GPU-based fast microscopic MC simulation package for water radiolysis. The first step of this project focused on accelerating the simulation of the chemical stage, the most time consuming stage in the entire water radiolysis process. A GPU-friendly parallelization strategy was designed to address the highly correlated many-body simulation problem caused by the mutual competitive chemical reactions between the radiolytic molecules. Two cases were tested, using a 750 keV electron and a 5 MeV proton incident in pure water, respectively. The time-dependent yields of all the radiolytic species during the chemical stage were used to evaluate the accuracy of the simulation. The relative differences between our simulation and the Geant4-DNA simulation were on average 5.3% and 4.4% for the two cases. Our package, executed on an Nvidia Titan black GPU card, successfully completed the chemical stage simulation of the two cases within 599.2 s and 489.0 s. As compared with Geant4-DNA that was executed on an Intel i7-5500U CPU processor and needed 28.6 h and 26.8 h for the two cases using a single CPU core, our package achieved a speed-up factor of 171.1-197.2.