101
|
Yang Y, Cao M, Gao Y, Kamrava M, Kalbasi A, Lamb J, Agazaryan N, Sheng K, Low D, Hu P. Longitudinal Diffusion MRI for Early Assessment of Treatment Response in Sarcoma Patients After Preoperative Radiation Therapy. Int J Radiat Oncol Biol Phys 2017. [DOI: 10.1016/j.ijrobp.2017.06.2389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
102
|
Nguyen D, O'Connor D, Ruan D, Sheng K. Deterministic direct aperture optimization using multiphase piecewise constant segmentation. Med Phys 2017; 44:5596-5609. [PMID: 28834556 DOI: 10.1002/mp.12529] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2016] [Revised: 06/07/2017] [Accepted: 08/11/2017] [Indexed: 11/11/2022] Open
Abstract
PURPOSE Direct aperture optimization (DAO) attempts to incorporate machine constraints in the inverse optimization to eliminate the post-processing steps in fluence map optimization (FMO) that degrade plan quality. Current commercial DAO methods utilize a stochastic or greedy approach to search a small aperture solution space. In this study, we propose a novel deterministic direct aperture optimization that integrates the segmentation of fluence map in the optimization problem using the multiphase piecewise constant Mumford-Shah formulation. METHODS The Mumford-Shah based direct aperture optimization problem was formulated to include an L2-norm dose fidelity term to penalize differences between the projected dose and the prescribed dose, an anisotropic total variation term to promote piecewise continuity in the fluence maps, and the multiphase piecewise constant Mumford-Shah function to partition the fluence into pairwise discrete segments. A proximal-class, first-order primal-dual solver was implemented to solve the large scale optimization problem, and an alternating module strategy was implemented to update fluence and delivery segments. Three patients of varying complexity-one glioblastoma multiforme (GBM) patient, one lung (LNG) patient, and one bilateral head and neck (H&N) patient with 3 PTVs-were selected to test the new DAO method. For each patient, 20 non-coplanar beams were first selected using column generation, followed by the Mumford-Shah based DAO (DAOMS ). For comparison, a popular and successful approach to DAO known as simulated annealing-a stochastic approach-was replicated. The simulated annealing DAO (DAOSA ) plans were then created using the same beam angles and maximum number of segments per beam. PTV coverage, PTV homogeneity D95D5, and OAR sparing were assessed for each plan. In addition, high dose spillage, defined as the 50% isodose volume divided by the tumor volume, as well as conformity, defined as the van't Riet conformation number, were evaluated. RESULTS DAOMS achieved essentially the same OAR doses compared with the DAOSA plans for the GBM case. The average difference of OAR Dmax and Dmean between the two plans were within 0.05% of the plan prescription dose. The lung case showed slightly improved critical structure sparing using the DAOMS approach, where the average OAR Dmax and Dmean were reduced by 3.67% and 1.08%, respectively, of the prescription dose. The DAOMS plan substantially improved OAR dose sparing for the H&N patient, where the average OAR Dmax and Dmean were reduced by over 10% of the prescription dose. The DAOMS and DAOSA plans were comparable for the GBM and LNG PTV coverage, while the DAOMS plan substantially improved the H&N PTV coverage, increasing D99 by 6.98% of the prescription dose. For the GBM and LNG patients, the DAOMS and DAOSA plans had comparable high dose spillage but slightly worse conformity with the DAOMS approach. For the H&N plan, DAOMS was considerably superior in high dose spillage and conformity to the DAOSA . The deterministic approach is able to solve the DAO problem substantially faster than the simulated annealing approach, with a 9.5- to 40-fold decrease in total solve time, depending on the patient case. CONCLUSIONS A novel deterministic direct aperture optimization formulation was developed and evaluated. It combines fluence map optimization and the multiphase piecewise constant Mumford-Shah segmentation into a unified framework, and the resulting optimization problem can be solved efficiently. Compared to the widely and commercially used simulated annealing DAO approach, it showed comparable dosimetry behavior for simple plans, and substantially improved OAR sparing, PTV coverage, PTV homogeneity, high dose spillage, and conformity for the more complex head and neck plan.
Collapse
|
103
|
Johnstone CD, Lindsay P, Graves EE, Wong E, Perez JR, Poirier Y, Ben-Bouchta Y, Kanesalingam T, Chen H, Rubinstein AE, Sheng K, Bazalova-Carter M. Multi-institutional MicroCT image comparison of image-guided small animal irradiators. Phys Med Biol 2017; 62:5760-5776. [PMID: 28574405 DOI: 10.1088/1361-6560/aa76b4] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
To recommend imaging protocols and establish tolerance levels for microCT image quality assurance (QA) performed on conformal image-guided small animal irradiators. A fully automated QA software SAPA (small animal phantom analyzer) for image analysis of the commercial Shelley micro-CT MCTP 610 phantom was developed, in which quantitative analyses of CT number linearity, signal-to-noise ratio (SNR), uniformity and noise, geometric accuracy, spatial resolution by means of modulation transfer function (MTF), and CT contrast were performed. Phantom microCT scans from eleven institutions acquired with four image-guided small animal irradiator units (including the commercial PXi X-RAD SmART and Xstrahl SARRP systems) with varying parameters used for routine small animal imaging were analyzed. Multi-institutional data sets were compared using SAPA, based on which tolerance levels for each QA test were established and imaging protocols for QA were recommended. By analyzing microCT data from 11 institutions, we established image QA tolerance levels for all image quality tests. CT number linearity set to R 2 > 0.990 was acceptable in microCT data acquired at all but three institutions. Acceptable SNR > 36 and noise levels <55 HU were obtained at five of the eleven institutions, where failing scans were acquired with current-exposure time of less than 120 mAs. Acceptable spatial resolution (>1.5 lp mm-1 for MTF = 0.2) was obtained at all but four institutions due to their large image voxel size used (>0.275 mm). Ten of the eleven institutions passed the set QA tolerance for geometric accuracy (<1.5%) and nine of the eleven institutions passed the QA tolerance for contrast (>2000 HU for 30 mgI ml-1). We recommend performing imaging QA with 70 kVp, 1.5 mA, 120 s imaging time, 0.20 mm voxel size, and a frame rate of 5 fps for the PXi X-RAD SmART. For the Xstrahl SARRP, we recommend using 60 kVp, 1.0 mA, 240 s imaging time, 0.20 mm voxel size, and 6 fps. These imaging protocols should result in high quality images that pass the set tolerance levels on all systems. Average SAPA computation time for complete QA analysis for a 0.20 mm voxel, 400 slice Shelley phantom microCT data set was less than 20 s. We present image quality assurance recommendations for image-guided small animal radiotherapy systems that can aid researchers in maintaining high image quality, allowing for spatially precise conformal dose delivery to small animals.
Collapse
|
104
|
Sheng K. EP-1521: Non-coplanar beam orientation and fluence map optimization based on group sparsity. Radiother Oncol 2017. [DOI: 10.1016/s0167-8140(17)31956-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
105
|
Nguyen D, Thomas D, Cao M, O'Connor D, Lamb J, Sheng K. Computerized triplet beam orientation optimization for MRI-guided Co-60 radiotherapy. Med Phys 2017; 43:5667. [PMID: 27782726 DOI: 10.1118/1.4963212] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Magnetic resonance imaging (MRI)-guided Co-60 provides daily and intrafractional MRI soft tissue imaging for improved target and critical organ tracking. To increase delivery efficiency, the system uses three Co-60 sources at 120° apart, allowing up to 600 cGy combined dose rate at isocenter. Despite the potential tripling in output, creating a delivery plan that uses all three sources is considerably unintuitive. Here, the authors computerize the triplet orientation optimization using column generation, an approach that was demonstrated effective in integrated beam orientation and fluence optimization for noncoplanar therapies. To achieve a better plan quality without increasing the treatment time, the authors then solve a fluence map optimization (FMO) problem while regularizing the fluence maps to reduce the number of deliverable MLC segments. METHODS Three patients-one prostate, one lung, and one head and neck boost plan (H&NBoost)-were evaluated in this study. For each patient, the beamlet doses were calculated using Monte Carlo, under a 0.35 T magnetic field, for 180 equally spaced coplanar beams grouped into 60 triplets. The beamlet size is 1.05 × 0.5 cm determined by the MLC leaf thickness and step size. The triplets were selected using the column generation algorithm. The FMO problem was formulated using an L2-norm dose fidelity term and an L1-norm anisotropic total variation regularization term, which allows controlling the number of MLC segments, and hence the treatment time, with minimal degradation to the dose. The authors' Fluence Regularization and Optimized Selection of Triplets (FROST) plans were compared against the clinical treatment plans (CLNs) produced by an experienced dosimetrist. PTV homogeneity, max dose, mean dose, D95, D98, and D99 were evaluated. OAR max and mean doses, as well as R50, defined as the ratio of the 50% isodose volume over the planning target volume were investigated. RESULTS The mean PTV D95, D98, and D99 differ by +0.04%, +0.07%, and +0.25% of the prescription dose between planning methods. The mean PTV homogeneity was virtually same with values at 0.8788 (FROST) and 0.8812 (CLN). R50 decreased by 0.67 comparing FROST to CLN. On average, FROST reduced Dmax and Dmean of OARs by 7.30% and 6.08% of the prescription dose, respectively. The manual CLN planning processes required numerous trial and error runs. The FROST plans on the other hand required minimal human intervention. CONCLUSIONS Efficient delivery of MRI-guided Co-60 therapy needs the output of multiple sources yet suffers from unintuitive and laborious manual beam selection processes. Computerized triplet orientation optimization improves both planning efficiency and plan dosimetry. The novel fluence map regularization provides additional controls over the number of MLC segments and treatment time.
Collapse
|
106
|
Han F, Zhou Z, Cao M, Yang Y, Sheng K, Hu P. Respiratory motion-resolved, self-gated 4D-MRI using rotating cartesian k-space (ROCK). Med Phys 2017; 44:1359-1368. [PMID: 28133752 DOI: 10.1002/mp.12139] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2016] [Revised: 12/20/2016] [Accepted: 01/23/2017] [Indexed: 12/26/2022] Open
Abstract
PURPOSE To propose and validate a respiratory motion resolved, self-gated (SG) 4D-MRI technique to assess patient-specific breathing motion of abdominal organs for radiation treatment planning. METHODS The proposed 4D-MRI technique was based on the balanced steady-state free-precession (bSSFP) technique and 3D k-space encoding. A novel rotating cartesian k-space (ROCK) reordering method was designed which incorporates repeatedly sampled k-space centerline as the SG motion surrogate and allows for retrospective k-space data binning into different respiratory positions based on the amplitude of the surrogate. The multiple respiratory-resolved 3D k-space data were subsequently reconstructed using a joint parallel imaging and compressed sensing method with spatial and temporal regularization. The proposed 4D-MRI technique was validated using a custom-made dynamic motion phantom and was tested in six healthy volunteers, in whom quantitative diaphragm and kidney motion measurements based on 4D-MRI images were compared with those based on 2D-CINE images. RESULTS The 5-minute 4D-MRI scan offers high-quality volumetric images in 1.2 × 1.2 × 1.6 mm3 and eight respiratory positions, with good soft-tissue contrast. In phantom experiments with triangular motion waveform, the motion amplitude measurements based on 4D-MRI were 11.89% smaller than the ground truth, whereas a -12.5% difference was expected due to data binning effects. In healthy volunteers, the difference between the measurements based on 4D-MRI and the ones based on 2D-CINE were 6.2 ± 4.5% for the diaphragm, 8.2 ± 4.9% and 8.9 ± 5.1% for the right and left kidney. CONCLUSION The proposed 4D-MRI technique could provide high-resolution, high-quality, respiratory motion-resolved 4D images with good soft-tissue contrast and are free of the "stitching" artifacts usually seen on 4D-CT and 4D-MRI based on resorting 2D-CINE. It could be used to visualize and quantify abdominal organ motion for MRI-based radiation treatment planning.
Collapse
|
107
|
Nguyen D, Lyu Q, Ruan D, O'Connor D, Low DA, Sheng K. A comprehensive formulation for volumetric modulated arc therapy planning. Med Phys 2017; 43:4263. [PMID: 27370141 DOI: 10.1118/1.4953832] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
PURPOSE Volumetric modulated arc therapy (VMAT) is a widely employed radiation therapy technique, showing comparable dosimetry to static beam intensity modulated radiation therapy (IMRT) with reduced monitor units and treatment time. However, the current VMAT optimization has various greedy heuristics employed for an empirical solution, which jeopardizes plan consistency and quality. The authors introduce a novel direct aperture optimization method for VMAT to overcome these limitations. METHODS The comprehensive VMAT (comVMAT) planning was formulated as an optimization problem with an L2-norm fidelity term to penalize the difference between the optimized dose and the prescribed dose, as well as an anisotropic total variation term to promote piecewise continuity in the fluence maps, preparing it for direct aperture optimization. A level set function was used to describe the aperture shapes and the difference between aperture shapes at adjacent angles was penalized to control MLC motion range. A proximal-class optimization solver was adopted to solve the large scale optimization problem, and an alternating optimization strategy was implemented to solve the fluence intensity and aperture shapes simultaneously. Single arc comVMAT plans, utilizing 180 beams with 2° angular resolution, were generated for a glioblastoma multiforme case, a lung (LNG) case, and two head and neck cases-one with three PTVs (H&N3PTV) and one with foue PTVs (H&N4PTV)-to test the efficacy. The plans were optimized using an alternating optimization strategy. The plans were compared against the clinical VMAT (clnVMAT) plans utilizing two overlapping coplanar arcs for treatment. RESULTS The optimization of the comVMAT plans had converged within 600 iterations of the block minimization algorithm. comVMAT plans were able to consistently reduce the dose to all organs-at-risk (OARs) as compared to the clnVMAT plans. On average, comVMAT plans reduced the max and mean OAR dose by 6.59% and 7.45%, respectively, of the prescription dose. Reductions in max dose and mean dose were as high as 14.5 Gy in the LNG case and 15.3 Gy in the H&N3PTV case. PTV coverages measured by D95, D98, and D99 were within 0.25% of the prescription dose. By comprehensively optimizing all beams, the comVMAT optimizer gained the freedom to allow some selected beams to deliver higher intensities, yielding a dose distribution that resembles a static beam IMRT plan with beam orientation optimization. CONCLUSIONS The novel nongreedy VMAT approach simultaneously optimizes all beams in an arc and then directly generates deliverable apertures. The single arc VMAT approach thus fully utilizes the digital Linac's capability in dose rate and gantry rotation speed modulation. In practice, the new single VMAT algorithm generates plans superior to existing VMAT algorithms utilizing two arcs.
Collapse
|
108
|
Shaverdian N, Yang Y, Hu P, Hart S, Sheng K, Lamb J, Cao M, Agazaryan N, Thomas D, Steinberg M, Low DA, Lee P. Feasibility evaluation of diffusion-weighted imaging using an integrated MRI-radiotherapy system for response assessment to neoadjuvant therapy in rectal cancer. Br J Radiol 2017; 90:20160739. [PMID: 28079398 DOI: 10.1259/bjr.20160739] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
OBJECTIVE To evaluate the feasibility of on-board diffusion-weighted imaging (DWI) with an integrated low-field MRI radiotherapy system to assess responses to neoadjuvant chemoradiation (NAC) in rectal cancer. METHODS A spin echo-based planar imaging diffusion sequence on a 0.35-T MRI radiotherapy system was acquired over the course of NAC. The apparent diffusion coefficients (ADCs) from the tumour regions of interest (ROIs) were calculated. A functional diffusion map (fDM) was created showing a pixelwise ADC analysis of the ROI over the course of treatment. Surgical pathology was correlated with ADC data. RESULTS Consecutive patients treated on a 0.35-T MRI radiotherapy system were evaluated. Patient A had the worst pathological response to NAC with a tumour regression score of 1 and was the only patient with a negative slope in the change of ADC values over the entire course of NAC, and during both the first and second half of NAC. The fDM from the first half of NAC for Patient A showed discrete dark areas in the tumour ROI, reflecting subregions with decreasing ADC values during NAC. Patient C had the most favourable pathological response to NAC with a Grade 3 response and was the only patient who had an increase in the slope in the change of ADC values from the first to the second half of NAC. CONCLUSION DWI using a low-field MRI radiotherapy system for evaluating the responses to NAC is feasible. Advances in knowledge: ADC values obtained using a 0.35-T MRI radiotherapy system over the course of NAC for rectal cancer correlate with pathological responses.
Collapse
|
109
|
Tran A, Zhang J, Woods K, Yu V, Nguyen D, Gustafson G, Rosen L, Sheng K. Treatment planning comparison of IMPT, VMAT and 4π radiotherapy for prostate cases. Radiat Oncol 2017; 12:10. [PMID: 28077128 PMCID: PMC5225526 DOI: 10.1186/s13014-016-0761-0] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2016] [Accepted: 12/29/2016] [Indexed: 11/10/2022] Open
Abstract
Background Intensity-modulated proton therapy (IMPT), non-coplanar 4π intensity-modulated radiation therapy (IMRT), and volumetric-modulated arc therapy (VMAT) represent the most advanced treatment methods based on heavy ion and X-rays, respectively. Here we compare their performance for prostate cancer treatment. Methods Ten prostate patients were planned using IMPT with robustness optimization, VMAT, and 4π to an initial dose of 54 Gy to a clinical target volume (CTV) that encompassed the prostate and seminal vesicles, then a boost prescription dose of 25.2 Gy to the prostate for a total dose of 79.2 Gy. The IMPT plans utilized two coplanar, oblique scanning beams 10° posterior of the lateral beam positions. Range uncertainties were taken into consideration in the IMPT plans. VMAT plans used two full, coplanar arcs to ensure sufficient PTV coverage. 4π plans were created by inversely selecting and optimizing 30 beams from 1162 candidate non-coplanar beams using a greedy column generation algorithm. CTV doses, bladder and rectum dose volumes (V40, V45, V60, V65, V70, V75, and V80), R100, R50, R10, and CTV homogeneity index (D95/D5) were evaluated. Results Compared to IMPT, 4π resulted in lower anterior rectal wall mean dose as well as lower rectum V40, V45, V60, V65, V70, and V75. Due to the opposing beam arrangement, IMPT resulted in significantly (p < 0.05) greater femoral head doses. However, IMPT plans had significantly lower bladder, rectum, and anterior rectal wall max dose. IMPT doses were also significantly more homogeneous than 4π and VMAT doses. Conclusion Compared to the VMAT and 4π plans, IMPT treatment plans are superior in CTV homogeneity and maximum point organ-at-risk (OAR) doses with the exception of femur heads. IMPT is inferior in rectum and bladder volumes receiving intermediate to high doses, particularly to the 4π plans, but significantly reduced low dose spillage and integral dose, which are correlated to secondary cancer for patients with expected long survival. The dosimetric benefits of 4π plans over VMAT are consistent with the previous publication.
Collapse
|
110
|
Han Y, Sheng K, Su M, Yang N, Wan D. Local mild hypothermia therapy as an augmentation strategy for minimally invasive surgery of hypertensive intracerebral hemorrhage: a meta-analysis of randomized clinical trials. Neuropsychiatr Dis Treat 2017; 13:41-49. [PMID: 28096671 PMCID: PMC5207467 DOI: 10.2147/ndt.s123501] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/02/2022] Open
Abstract
BACKGROUND Previous studies reported that the mild hypothermia therapy (MHT) could significantly improve the clinical outcomes for patients with hypertensive intracerebral hemorrhage (HICH). Therefore, this meta-analysis was conducted to systematically assess whether the addition of local MHT (LMHT) could significantly improve the efficacy of minimally invasive surgery (MIS) in treating HICH. METHODS Randomized clinical trials on the combined application of MIS and LMHT (MIS+LMHT) vs MIS alone for treating HICH were searched up to September 2016 in databases. Response rate and mortality rate were the primary outcomes, and the neurologic function and Barthel index were the secondary outcomes. Side effects were also analyzed. RESULTS Totally, 28 studies composed of 2,325 patients were included to compare the efficacy of MIS+LMHT to MIS alone. The therapeutic effects of MIS+LMHT were significantly better than MIS alone. The pooled odds ratio of response rate and mortality rate was 2.68 (95% confidence interval [CI]=2.22-3.24) and 0.43 (95% CI=0.32-0.57), respectively. In addition, the MIS+LMHT led to a significantly better improvement in the neurologic function and activities of daily living. The incidence of pneumonia was similar between the two treatment methods. CONCLUSION These results indicated that compared to MIS alone, the MIS+LMHT could be more effective for the acute treatment of patients with HICH. This treatment modality should be further explored and optimized.
Collapse
|
111
|
Hasse K, Neylon J, Sheng K, Santhanam AP. Systematic feasibility analysis of a quantitative elasticity estimation for breast anatomy using supine/prone patient postures. Med Phys 2016; 43:1299-1311. [PMID: 26936715 DOI: 10.1118/1.4941745] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
PURPOSE Breast elastography is a critical tool for improving the targeted radiotherapy treatment of breast tumors. Current breast radiotherapy imaging protocols only involve prone and supine CT scans. There is a lack of knowledge on the quantitative accuracy with which breast elasticity can be systematically measured using only prone and supine CT datasets. The purpose of this paper is to describe a quantitative elasticity estimation technique for breast anatomy using only these supine/prone patient postures. Using biomechanical, high-resolution breast geometry obtained from CT scans, a systematic assessment was performed in order to determine the feasibility of this methodology for clinically relevant elasticity distributions. METHODS A model-guided inverse analysis approach is presented in this paper. A graphics processing unit (GPU)-based linear elastic biomechanical model was employed as a forward model for the inverse analysis with the breast geometry in a prone position. The elasticity estimation was performed using a gradient-based iterative optimization scheme and a fast-simulated annealing (FSA) algorithm. Numerical studies were conducted to systematically analyze the feasibility of elasticity estimation. For simulating gravity-induced breast deformation, the breast geometry was anchored at its base, resembling the chest-wall/breast tissue interface. Ground-truth elasticity distributions were assigned to the model, representing tumor presence within breast tissue. Model geometry resolution was varied to estimate its influence on convergence of the system. A priori information was approximated and utilized to record the effect on time and accuracy of convergence. The role of the FSA process was also recorded. A novel error metric that combined elasticity and displacement error was used to quantify the systematic feasibility study. For the authors' purposes, convergence was set to be obtained when each voxel of tissue was within 1 mm of ground-truth deformation. RESULTS The authors' analyses showed that a ∼97% model convergence was systematically observed with no-a priori information. Varying the model geometry resolution showed no significant accuracy improvements. The GPU-based forward model enabled the inverse analysis to be completed within 10-70 min. Using a priori information about the underlying anatomy, the computation time decreased by as much as 50%, while accuracy improved from 96.81% to 98.26%. The use of FSA was observed to allow the iterative estimation methodology to converge more precisely. CONCLUSIONS By utilizing a forward iterative approach to solve the inverse elasticity problem, this work indicates the feasibility and potential of the fast reconstruction of breast tissue elasticity using supine/prone patient postures.
Collapse
|
112
|
Yang Y, Cao M, Sheng K, Gao Y, Chen A, Kamrava M, Lee P, Agazaryan N, Lamb J, Thomas D, Low D, Hu P. Longitudinal diffusion MRI for treatment response assessment: Preliminary experience using an MRI-guided tri-cobalt 60 radiotherapy system. Med Phys 2016; 43:1369-73. [PMID: 26936721 DOI: 10.1118/1.4942381] [Citation(s) in RCA: 79] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
PURPOSE To demonstrate the preliminary feasibility of a longitudinal diffusion magnetic resonance imaging (MRI) strategy for assessing patient response to radiotherapy at 0.35 T using an MRI-guided radiotherapy system (ViewRay). METHODS Six patients (three head and neck cancer, three sarcoma) who underwent fractionated radiotherapy were enrolled in this study. A 2D multislice spin echo single-shot echo planar imaging diffusion pulse sequence was implemented on the ViewRay system and tested in phantom studies. The same pulse sequence was used to acquire longitudinal diffusion data (every 2-5 fractions) on the six patients throughout the entire course of radiotherapy. The reproducibility of the apparent diffusion coefficient (ADC) measurements was assessed using reference regions and the temporal variations of the tumor ADC values were evaluated. RESULTS In diffusion phantom studies, the ADC values measured on the ViewRay system matched well with reference ADC values with <5% error for a range of ground truth diffusion coefficients of 0.4-1.1 × 10(-3) mm(2)/s. The remote reference regions (i.e., brainstem in head and neck patients) had consistent ADC values throughout the therapy for all three head and neck patients, indicating acceptable reproducibility of the diffusion imaging sequence. The tumor ADC values changed throughout therapy, with the change differing between patients, ranging from a 40% drop in ADC within the first week of therapy to gradually increasing throughout therapy. For larger tumors, intratumoral heterogeneity was observed. For one sarcoma patient, postradiotherapy biopsy showed less than 10% necrosis score, which correlated with the observed 40% decrease in ADC from the fifth fraction to the eighth treatment fraction. CONCLUSIONS This pilot study demonstrated that longitudinal diffusion MRI is feasible using the 0.35 T ViewRay MRI. Larger patient cohort studies are warranted to correlate the longitudinal diffusion measurements to patient outcomes. Such an approach may enable response-guided adaptive radiotherapy.
Collapse
|
113
|
Nguyen D, Ruan D, O'Connor D, Woods K, Low DA, Boucher S, Sheng K. A novel software and conceptual design of the hardware platform for intensity modulated radiation therapy. Med Phys 2016; 43:917-29. [PMID: 26843252 DOI: 10.1118/1.4940353] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
PURPOSE To deliver high quality intensity modulated radiotherapy (IMRT) using a novel generalized sparse orthogonal collimators (SOCs), the authors introduce a novel direct aperture optimization (DAO) approach based on discrete rectangular representation. METHODS A total of seven patients-two glioblastoma multiforme, three head & neck (including one with three prescription doses), and two lung-were included. 20 noncoplanar beams were selected using a column generation and pricing optimization method. The SOC is a generalized conventional orthogonal collimators with N leaves in each collimator bank, where N = 1, 2, or 4. SOC degenerates to conventional jaws when N = 1. For SOC-based IMRT, rectangular aperture optimization (RAO) was performed to optimize the fluence maps using rectangular representation, producing fluence maps that can be directly converted into a set of deliverable rectangular apertures. In order to optimize the dose distribution and minimize the number of apertures used, the overall objective was formulated to incorporate an L2 penalty reflecting the difference between the prescription and the projected doses, and an L1 sparsity regularization term to encourage a low number of nonzero rectangular basis coefficients. The optimization problem was solved using the Chambolle-Pock algorithm, a first-order primal-dual algorithm. Performance of RAO was compared to conventional two-step IMRT optimization including fluence map optimization and direct stratification for multileaf collimator (MLC) segmentation (DMS) using the same number of segments. For the RAO plans, segment travel time for SOC delivery was evaluated for the N = 1, N = 2, and N = 4 SOC designs to characterize the improvement in delivery efficiency as a function of N. RESULTS Comparable PTV dose homogeneity and coverage were observed between the RAO and the DMS plans. The RAO plans were slightly superior to the DMS plans in sparing critical structures. On average, the maximum and mean critical organ doses were reduced by 1.94% and 1.44% of the prescription dose. The average number of delivery segments was 12.68 segments per beam for both the RAO and DMS plans. The N = 2 and N = 4 SOC designs were, on average, 1.56 and 1.80 times more efficient than the N = 1 SOC design to deliver. The mean aperture size produced by the RAO plans was 3.9 times larger than that of the DMS plans. CONCLUSIONS The DAO and dose domain optimization approach enabled high quality IMRT plans using a low-complexity collimator setup. The dosimetric quality is comparable or slightly superior to conventional MLC-based IMRT plans using the same number of delivery segments. The SOC IMRT delivery efficiency can be significantly improved by increasing the leaf numbers, but the number is still significantly lower than the number of leaves in a typical MLC.
Collapse
|
114
|
Yang Y, Cao M, Kaprealian T, Sheng K, Gao Y, Han F, Gomez C, Santhanam A, Tenn S, Agazaryan N, Low DA, Hu P. Accuracy of UTE-MRI-based patient setup for brain cancer radiation therapy. Med Phys 2016; 43:262. [PMID: 26745919 DOI: 10.1118/1.4938266] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
PURPOSE Radiation therapy simulations solely based on MRI have advantages compared to CT-based approaches. One feature readily available from computed tomography (CT) that would need to be reproduced with MR is the ability to compute digitally reconstructed radiographs (DRRs) for comparison against on-board radiographs commonly used for patient positioning. In this study, the authors generate MR-based bone images using a single ultrashort echo time (UTE) pulse sequence and quantify their 3D and 2D image registration accuracy to CT and radiographic images for treatments in the cranium. METHODS Seven brain cancer patients were scanned at 1.5 T using a radial UTE sequence. The sequence acquired two images at two different echo times. The two images were processed using an in-house software to generate the UTE bone images. The resultant bone images were rigidly registered to simulation CT data and the registration error was determined using manually annotated landmarks as references. DRRs were created based on UTE-MRI and registered to simulated on-board images (OBIs) and actual clinical 2D oblique images from ExacTrac™. RESULTS UTE-MRI resulted in well visualized cranial, facial, and vertebral bones that quantitatively matched the bones in the CT images with geometric measurement errors of less than 1 mm. The registration error between DRRs generated from 3D UTE-MRI and the simulated 2D OBIs or the clinical oblique x-ray images was also less than 1 mm for all patients. CONCLUSIONS UTE-MRI-based DRRs appear to be promising for daily patient setup of brain cancer radiotherapy with kV on-board imaging.
Collapse
|
115
|
Nguyen D, Thomas D, Cao M, O'Connor D, Lamb J, Sheng K. Automated Triplet Beam Orientation Optimization for Magnetic Resonance Imaging–Guided Co-60 Radiation Therapy. Int J Radiat Oncol Biol Phys 2016. [DOI: 10.1016/j.ijrobp.2016.06.155] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
116
|
Yu V, Tran A, Nguyen D, Woods K, Kaprealian T, Chin R, Low D, Sheng K. Significant Cord and Esophagus Dose Reduction by 4π Non-Coplanar Spine Stereotactic Body Radiation Therapy and Stereotactic Radiosurgery. Int J Radiat Oncol Biol Phys 2016. [DOI: 10.1016/j.ijrobp.2016.06.2246] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
117
|
Yu VY, Tran A, Nguyen D, Cao M, Ruan D, Low DA, Sheng K. The development and verification of a highly accurate collision prediction model for automated noncoplanar plan delivery. Med Phys 2016; 42:6457-67. [PMID: 26520735 DOI: 10.1118/1.4932631] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
PURPOSE Significant dosimetric benefits had been previously demonstrated in highly noncoplanar treatment plans. In this study, the authors developed and verified an individualized collision model for the purpose of delivering highly noncoplanar radiotherapy and tested the feasibility of total delivery automation with Varian TrueBeam developer mode. METHODS A hand-held 3D scanner was used to capture the surfaces of an anthropomorphic phantom and a human subject, which were positioned with a computer-aided design model of a TrueBeam machine to create a detailed virtual geometrical collision model. The collision model included gantry, collimator, and couch motion degrees of freedom. The accuracy of the 3D scanner was validated by scanning a rigid cubical phantom with known dimensions. The collision model was then validated by generating 300 linear accelerator orientations corresponding to 300 gantry-to-couch and gantry-to-phantom distances, and comparing the corresponding distance measurements to their corresponding models. The linear accelerator orientations reflected uniformly sampled noncoplanar beam angles to the head, lung, and prostate. The distance discrepancies between measurements on the physical and virtual systems were used to estimate treatment-site-specific safety buffer distances with 0.1%, 0.01%, and 0.001% probability of collision between the gantry and couch or phantom. Plans containing 20 noncoplanar beams to the brain, lung, and prostate optimized via an in-house noncoplanar radiotherapy platform were converted into XML script for automated delivery and the entire delivery was recorded and timed to demonstrate the feasibility of automated delivery. RESULTS The 3D scanner measured the dimension of the 14 cm cubic phantom within 0.5 mm. The maximal absolute discrepancy between machine and model measurements for gantry-to-couch and gantry-to-phantom was 0.95 and 2.97 cm, respectively. The reduced accuracy of gantry-to-phantom measurements was attributed to phantom setup errors due to the slightly deformable and flexible phantom extremities. The estimated site-specific safety buffer distance with 0.001% probability of collision for (gantry-to-couch, gantry-to-phantom) was (1.23 cm, 3.35 cm), (1.01 cm, 3.99 cm), and (2.19 cm, 5.73 cm) for treatment to the head, lung, and prostate, respectively. Automated delivery to all three treatment sites was completed in 15 min and collision free using a digital Linac. CONCLUSIONS An individualized collision prediction model for the purpose of noncoplanar beam delivery was developed and verified. With the model, the study has demonstrated the feasibility of predicting deliverable beams for an individual patient and then guiding fully automated noncoplanar treatment delivery. This work motivates development of clinical workflows and quality assurance procedures to allow more extensive use and automation of noncoplanar beam geometries.
Collapse
|
118
|
Sheng K, Fang W, Zhu Y, Shuai G, Zou D, Su M, Han Y, Cheng O. Different Alterations of Cerebral Regional Homogeneity in Early-Onset and Late-Onset Parkinson's Disease. Front Aging Neurosci 2016; 8:165. [PMID: 27462265 PMCID: PMC4940400 DOI: 10.3389/fnagi.2016.00165] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2016] [Accepted: 06/22/2016] [Indexed: 11/13/2022] Open
Abstract
UNLABELLED HIGHLIGHTS Eighteen EOPD, 21 LOPD and 37 age-matched normal control subjects participated in the resting state fMRI scans.Age at onset of PD modulates the distribution of cerebral regional homogeneity during resting state.Disproportionate putamen alterations are more prominent in PD patients with a younger age of onset. OBJECTIVE Early-onset Parkinson's disease (EOPD) is distinct from late-onset PD (LOPD) as it relates to the clinical profile and response to medication. The objective of current paper is to investigate whether characteristics of spontaneous brain activity in the resting state are associated with the age of disease onset. METHODS We assessed the correlation between neural activity and age-at-onset in a sample of 39 PD patients (18 EOPD and 21 LOPD) and 37 age-matched normal control subjects. Regional homogeneity (ReHo) approaches were employed using ANOVA with two factors: PD and age. RESULTS In the comparisons between LOPD and EOPD, EOPD revealed lower ReHo values in the right putamen and higher ReHo values in the left superior frontal gyrus. Compared with age-matched control subjects, EOPD exhibited lower ReHo values in the right putamen and higher ReHo values in the left inferior temporal gyrus; However, LOPD showed lower ReHo values in the right putamen and left insula. The ReHo values were negatively correlated with the UPDRS total scores in the right putamen in LOPD, but a correlation between the ReHo value and UPDRS score was not detected in EOPD. CONCLUSIONS Our findings support the notion that age at onset is associated with the distribution of cerebral regional homogeneity in the resting state and suggest that disproportionate putamen alterations are more prominent in patients with a younger age of onset.
Collapse
|
119
|
Yu V, Ruan D, Nguyen D, Kaprealian T, Chin R, Sheng K. SU-F-R-17: Advancing Glioblastoma Multiforme (GBM) Recurrence Detection with MRI Image Texture Feature Extraction and Machine Learning. Med Phys 2016. [DOI: 10.1118/1.4955789] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
120
|
Zhang J, Nguyen D, Woods K, Tran A, Li X, Ding X, Kabolizadeh P, Guerrero T, Sheng K. SU-F-T-186: A Treatment Planning Study of Normal Tissue Sparing with Robustness Optimized IMPT, 4Pi IMRT, and VMAT for Head and Neck Cases. Med Phys 2016. [DOI: 10.1118/1.4956323] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
121
|
Woods K, Karunamuni R, Tran A, Yu V, Nguyen D, Hattangadi-Gluth J, Sheng K. TH-EF-BRB-01: BEST IN PHYSICS (THERAPY): Dosimetric Comparison of 4π and Clinical IMRT for Cortex-Sparing High-Grade Glioma Treatment. Med Phys 2016. [DOI: 10.1118/1.4958247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
122
|
O'Connor D, Voronenko Y, Nguyen D, Yin W, Sheng K. TH-EF-BRB-05: 4pi Non-Coplanar IMRT Beam Angle Selection by Convex Optimization with Group Sparsity Penalty. Med Phys 2016. [DOI: 10.1118/1.4958251] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
123
|
Han F, Zhou Z, Yang Y, Sheng K, Hu P. SU-F-J-158: Respiratory Motion Resolved, Self-Gated 4D-MRI Using Rotating Cartesian K-Space Sampling. Med Phys 2016. [DOI: 10.1118/1.4956066] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
124
|
Nguyen D, Thomas D, Cao M, O'Connor D, Lamb J, Sheng K. TH-AB-BRA-02: Automated Triplet Beam Orientation Optimization for MRI-Guided Co-60 Radiotherapy. Med Phys 2016. [DOI: 10.1118/1.4958053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
125
|
Yang Y, Cao M, Kamrava M, Low D, Sheng K, Lamb J, Agazaryan N, Thomas D, Hu P. WE-FG-202-11: Longitudinal Diffusion MRI for Treatment Assessment of Sarcoma Patients with Pre-Operative Radiation Therapy. Med Phys 2016. [DOI: 10.1118/1.4957923] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|