Comparison of skin dose in IMRT and VMAT with TrueBeam and Halcyon linear accelerator for whole breast irradiation.

With the increasing use of flattening filter free (FFF) beams, it is important to evaluate the impact on the skin dose and target coverage of breast cancer treatments. This study aimed to compare skin doses of treatments using FFF and flattening filter (FF) beams for breast cancer. The study established treatment plans for left breast of an anthropomorphic phantom using Halcyon's 6-MV FFF beam and TrueBeam's 6-MV FF beam. Volumetric modulated arc therapy (VMAT) with varying numbers of arcs and intensity modulated radiation therapy (IMRT) were employed, and skin doses were measured at five points using Gafchromic EBT3 film. Each measurement was repeated three times, and averaged to reduce uncertainty. All plans were compared in terms of plan quality to ensure homogeneous target coverage. The study found that when using VMAT with two, four, and six arcs, in-field doses were 19%, 15%, and 6% higher, respectively, when using Halcyon compared to TrueBeam. Additionally, when using two arcs for VMAT, in-field doses were 10% and 15% higher compared to four and six arcs when using Halcyon. Finally, in-field dose from Halcyon using IMRT was about 1% higher than when using TrueBeam. Our research confirmed that when treating breast cancer with FFF beams, skin dose is higher than with traditional FF beams. Moreover, number of arcs used in VMAT treatment with FFF beams affects skin dose to the patient. To maintain a skin dose similar to that of FF beams when using Halcyon, it may be worth considering increasing the number of arcs.

Development of open access tool for automatic use factor calculation using DICOM-RT patient data.

Our study recalculated the use factor of linear accelerators (LINACs) by using an in-house program based on Digital Imaging and Communications in Medicine radiation therapy (DICOM-RT). We considered the impact of advancements and changes in treatment trends, including modality, technology, and radiation dose, on the use factor, which is one of the shielding parameters. In accordance with the methodology described in the NCRP 151 report, we computed the use factor for four linear accelerators (LINACs) across three hospitals. We analyzed the results based on the treatment techniques and treatment sites for three-dimensional conformal radiation therapy (3D-CRT) and intensity modulated radiation therapy or volumetric modulated arc therapy. Our findings revealed that the use factors obtained at 45° and 90° were 14.8% and 13.5% higher than those of the NCRP 151 report. In treatment rooms with a high 3D-CRT ratio, the use factor at a specific angle differed by up to 14.6% relative to the NCRP 151 report value. Our results showed a large difference in the use factor for specific sites such as the breast and spine, so it is recommended that each institution recalculate the use factor using patient's data.

Surgery vs. radiotherapy for locally advanced hypopharyngeal cancer in the contemporary era: A population-based study.

OBJECTIVES: To compare overall survival (OS) in locally advanced hypopharyngeal cancer treated with surgery or definitive chemoradiotherapy in the contemporary era. METHODS: From 2010 to 2015, data for patients diagnosed with hypopharyngeal cancer (T2-T4aM0) and treated with total pharyngectomy with lymph node dissection (surgery group) or definitive radiotherapy and chemotherapy (chemoradiotherapy group) was retrieved from the SEER database. Multivariate analyses were performed in each subgroup divided according to T category (T2-3 or T4a). RESULTS: The number of patients in the surgery and chemoradiotherapy groups was 209 and 648, respectively. Among them, the number of T4a patients was 111 and 126 in each group. Three-year OS rate in the surgery and chemoradiotherapy groups was 37.9% and 44.1%, respectively (P = 0.178). The 3-year OS rate for the T2-3 patients was 46.5% and 48.7% (P = 0.598), and the 3-year OS rate for the T4a patients was 29.9% and 26.1% in the surgery and chemoradiotherapy groups, respectively (P = 0.439). On multivariate analysis, the chemoradiotherapy group was not inferior to the surgery group in T2-T4a patients (Hazard ratio [HR] for the chemoradiotherapy group 0.889, 95% confidence interval [CI] 0.699-1.129, P = 0.334), in T2-3 patients (HR 0.932, 95% CI 0.699-1.297, P = 0.675), and in T4a patients (HR 0.880, 95% CI 0.617-1.256, P = 0.481). CONCLUSIONS: Chemoradiotherapy for locally advanced hypophagyngeal cancer showed a comparable OS rate to surgery. For patients with T4a category cancer with high possibility of preserving the laryngopharyngeal function, chemoradiotherapy may be a promising alternative treatment.

Analysis of Neutron Production in Passively Scattered Ion-Beam Therapy.

A new treatment facility for heavy ion therapy since 2010 was constructed. In the broad beam, a range shifter, ridge filter and multi leaf collimator (MLC) for the generation of the spread-out Bragg peak is used. In this case, secondary neutrons produced by the interactions of the ion field with beam-modifying devices (e.g. double-scattering system, beam shaping collimators and range compensators) are very important for patient safety. Therefore, these components must be carefully examined in the context of secondary neutron yield and associated secondary cancer risk. In this article, Monte Carlo simulation has been carried out with the FLUktuierende KAskade particle transport code, the fluence and distribution of neutron generation and the neutron dose equivalent from the broad beam components are compared using carbon and proton beams. As a result, it is confirmed that the yield of neutron production using a carbon beam from all components of the broad beam was higher than using a proton beam. The ambient dose by neutrons per heavy ion and proton ion from the MLC surface was 0.12-0.18 and 0.0067-0.0087 pSv, respectively, which shows that heavy ions generate more neutrons than protons. However, ambient dose per treatment 2 Gy, which means physical dose during treatment by ion beam, is higher than carbon beam because proton therapy needs more beam flux to make 2-Gy prescription dose. Therefore, the neutron production from the MLC, which is closed to the patient, is a very important parameter for patient safety.

Treatment outcomes after adjuvant radiotherapy following surgery for patients with stage I endometrial cancer.

PURPOSE: The purpose of this study is to evaluate the treatment outcomes of adjuvant radiotherapy using vaginal brachytherapy (VB) with a lower dose per fraction and/or external beam radiotherapy (EBRT) following surgery for patients with stage I endometrial carcinoma. MATERIALS AND METHODS: The subjects were 43 patients with the International Federation of Gynecology and Obstetrics (FIGO) stage I endometrial cancer who underwent adjuvant radiotherapy following surgery between March 2000 and April 2014. Of these, 25 received postoperative VB alone, while 18 received postoperative EBRT to the whole pelvis; 3 of these were treated with EBRT plus VB. The median EBRT dose was 50.0 Gy (45.0-50.4 Gy) and the VB dose was 24 Gy in 6 fractions. Tumor dose was prescribed at a depth of 5 mm from the cylinder surface and delivered twice per week. RESULTS: The median follow-up period for all patients was 57 months (range, 9 to 188 months). Five-year disease-free survival (DFS) and overall survival (OS) for all patients were 92.5% and 95.3%, respectively. Adjuvant radiotherapy was performed according to risk factors and stage IB, grade 3 and lymphovascular invasion were observed more frequently in the EBRT group. Five-year DFS for EBRT and VB alone were 88.1% and 96.0%, respectively (p = 0.42), and 5-year OS for EBRT and VB alone were 94.4% and 96%, respectively (p = 0.38). There was no locoregional recurrence in any patient. Two patients who received EBRT and 1 patient who received VB alone developed distant metastatic disease. Two patients who received EBRT had severe complications, one each of grade 3 gastrointestinal complication and pelvic bone insufficiency fracture. CONCLUSION: Adjuvant radiotherapy achieved high DFS and OS with acceptable toxicity in stage I endometrial cancer. VB (with a lower dose per fraction) may be a viable option for selected patients with early-stage endometrial cancer following surgery.

Two-year Follow-up of Volumetric-modulated Arc Therapy for Treating Internal Mammary Nodes in Locally Advanced Breast Cancer.

BACKGROUND/AIM: We evaluated dosimetric and clinical results of volumetric-modulated arc therapy (VMAT) treating internal mammary nodes (IMNs) in patients with locally advanced breast cancer. PATIENTS AND METHODS: Between July 2012 and December 2013, treatment outcomes of 31 patients with stage III breast cancer treated with postoperative radiotherapy including the IMN field using VMAT were analyzed. RESULTS: The median follow-up was 25.2 months. VMAT showed good homogeneity and improved conformity. The incidence of lung toxicity of grade 2 or more was 3.2%, and cardiac toxicity was not observed. No patient developed local recurrence. CONCLUSION: VMAT appears safe for the treatment of locally advanced breast cancer including the IMN field.

Automatic deformable surface registration for medical applications by radial basis function-based robust point-matching.

Deformable surface mesh registration is a useful technique for various medical applications, such as intra-operative treatment guidance and intra- or inter-patient study. In this paper, we propose an automatic deformable mesh registration technique. The proposed method iteratively deforms a source mesh to a target mesh without manual feature extraction. Each iteration of the registration consists of two steps, automatic correspondence finding using robust point-matching (RPM) and local deformation using a radial basis function (RBF). The proposed RBF-based RPM algorithm solves the interlocking problems of correspondence and deformation using a deterministic annealing framework with fuzzy correspondence and RBF interpolation. Simulation tests showed promising results, with the average deviations decreasing by factors of 21.2 and 11.9, respectively. In the human model test, the average deviation decreased from 1.72±1.88mm to 0.57±0.66mm. We demonstrate the effectiveness of the proposed method by presenting some medical applications.

Prognostic analysis of uterine cervical cancer treated with postoperative radiotherapy: importance of positive or close parametrial resection margin.

PURPOSE: To analyze prognostic factors for locoregional recurrence (LRR), distant metastasis (DM), and overall survival (OS) in cervical cancer patients who underwent radical hysterectomy followed by postoperative radiotherapy (PORT) in a single institute. MATERIALS AND METHODS: Clinicopathologic data of 135 patients with clinical stage IA2 to IIA2 cervical cancer treated with PORT from 2001 to 2012 were reviewed, retrospectively. Postoperative parametrial resection margin (PRM) and vaginal resection margin (VRM) were investigated separately. The median treatment dosage of external beam radiotherapy (EBRT) to the whole pelvis was 50.4 Gy in 1.8 Gy/fraction. High-dose-rate vaginal brachytherapy after EBRT was given to patients with positive or close VRMs. Concurrent platinum-based chemoradiotherapy (CCRT) was administered to 73 patients with positive resection margin, lymph node (LN) metastasis, or direct extension of parametrium. Kaplan-Meier method and log-rank test were used for analyzing LRR, DM, and OS; Cox regression was applied to analyze prognostic factors. RESULTS: The 5-year disease-free survival was 79% and 5-year OS was 91%. In univariate analysis, positive or close PRM, LN metastasis, direct extension of parametrium, lymphovascular invasion, histology of adenocarcinoma, and chemotherapy were related with more DM and poor OS. In multivariate analysis, PRM and LN metastasis remained independent prognostic factors for OS. CONCLUSION: PORT after radical hysterectomy in uterine cervical cancer showed excellent OS in this study. Positive or close PRM after radical hysterectomy in uterine cervical cancer correlates with poor prognosis even with CCRT. Therefore, additional treatments to improve local control such as radiation boosting need to be considered.

Beam's-eye-view dosimetrics (BEVD) guided rotational station parameter optimized radiation therapy (SPORT) planning based on reweighted total-variation minimization.

Conventional VMAT optimizes aperture shapes and weights at uniformly sampled stations, which is a generalization of the concept of a control point. Recently, rotational station parameter optimized radiation therapy (SPORT) has been proposed to improve the plan quality by inserting beams to the regions that demand additional intensity modulations, thus formulating nonuniform beam sampling. This work presents a new rotational SPORT planning strategy based on reweighted total-variation (TV) minimization (min.), using beam's-eye-view dosimetrics (BEVD) guided beam selection. The convex programming based reweighted TV min. assures the simplified fluence-map, which facilitates single-aperture selection at each station for single-arc delivery. For the rotational arc treatment planning and non-uniform beam angle setting, the mathematical model needs to be modified by additional penalty term describing the fluence-map similarity and by determination of appropriate angular weighting factors. The proposed algorithm with additional penalty term is capable of achieving more efficient and deliverable plans adaptive to the conventional VMAT and SPORT planning schemes by reducing the dose delivery time about 5 to 10 s in three clinical cases (one prostate and two head-and-neck (HN) cases with a single and multiple targets). The BEVD guided beam selection provides effective and yet easy calculating methodology to select angles for denser, non-uniform angular sampling in SPORT planning. Our BEVD guided SPORT treatment schemes improve the dose sparing to femoral heads in the prostate and brainstem, parotid glands and oral cavity in the two HN cases, where the mean dose reduction of those organs ranges from 0.5 to 2.5 Gy. Also, it increases the conformation number assessing the dose conformity to the target from 0.84, 0.75 and 0.74 to 0.86, 0.79 and 0.80 in the prostate and two HN cases, while preserving the delivery efficiency, relative to conventional single-arc VMAT plans.

Accuracy of surface registration compared to conventional volumetric registration in patient positioning for head-and-neck radiotherapy: a simulation study using patient data.

PURPOSE: 3D optical surface imaging has been applied to patient positioning in radiation therapy (RT). The optical patient positioning system is advantageous over conventional method using cone-beam computed tomography (CBCT) in that it is radiation free, frameless, and is capable of real-time monitoring. While the conventional radiographic method uses volumetric registration, the optical system uses surface matching for patient alignment. The relative accuracy of these two methods has not yet been sufficiently investigated. This study aims to investigate the theoretical accuracy of the surface registration based on a simulation study using patient data. METHODS: This study compares the relative accuracy of surface and volumetric registration in head-and-neck RT. The authors examined 26 patient data sets, each consisting of planning CT data acquired before treatment and patient setup CBCT data acquired at the time of treatment. As input data of surface registration, patient's skin surfaces were created by contouring patient skin from planning CT and treatment CBCT. Surface registration was performed using the iterative closest points algorithm by point-plane closest, which minimizes the normal distance between source points and target surfaces. Six degrees of freedom (three translations and three rotations) were used in both surface and volumetric registrations and the results were compared. The accuracy of each method was estimated by digital phantom tests. RESULTS: Based on the results of 26 patients, the authors found that the average and maximum root-mean-square translation deviation between the surface and volumetric registrations were 2.7 and 5.2 mm, respectively. The residual error of the surface registration was calculated to have an average of 0.9 mm and a maximum of 1.7 mm. CONCLUSIONS: Surface registration may lead to results different from those of the conventional volumetric registration. Only limited accuracy can be achieved for patient positioning with an approach based solely on surface information.

Motion management within two respiratory-gating windows: feasibility study of dual quasi-breath-hold technique in gated medical procedures.

A dual quasi-breath-hold (DQBH) technique is proposed for respiratory motion management (a hybrid technique combining breathing-guidance with breath-hold task in the middle). The aim of this study is to test a hypothesis that the DQBH biofeedback system improves both the capability of motion management and delivery efficiency. Fifteen healthy human subjects were recruited for two respiratory motion measurements (free breathing and DQBH biofeedback breathing for 15 min). In this study, the DQBH biofeedback system utilized the abdominal position obtained using an real-time position management (RPM) system (Varian Medical Systems, Palo Alto, USA) to audio-visually guide a human subject for 4 s breath-hold at EOI and 90% EOE (EOE90%) to improve delivery efficiency. We investigated the residual respiratory motion and the delivery efficiency (duty-cycle) of abdominal displacement within the gating window. The improvement of the abdominal motion reproducibility was evaluated in terms of cycle-to-cycle displacement variability, respiratory period and baseline drift. The DQBH biofeedback system improved the abdominal motion management capability compared to that with free breathing. With a phase based gating (mean ± std: 55 ± 5%), the averaged root mean square error (RMSE) of the abdominal displacement in the dual-gating windows decreased from 2.26 mm of free breathing to 1.16 mm of DQBH biofeedback (p-value = 0.007). The averaged RMSE of abdominal displacement over the entire respiratory cycles reduced from 2.23 mm of free breathing to 1.39 mm of DQBH biofeedback breathing in the dual-gating windows (p-value = 0.028). The averaged baseline drift dropped from 0.9 mm min(-1) with free breathing to 0.09 mm min(-1) with DQBH biofeedback (p-value = 0.048). The averaged duty-cycle with an 1 mm width of displacement bound increased from 15% of free breathing to 26% of DQBH biofeedback (p-value = 0.003). The study demonstrated that the DQBH biofeedback system has the potential to significantly reduce the residual respiratory motion with the improved duty cycle during the respiratory gating procedure.

Negative margin technique - a novel planning strategy to improve dose conformation in SBRT using dynamic conformal arc delivery.

The purpose of this study was to introduce a planning strategy for dynamic conformal arc therapy (DCAT), named negative margin technique (NMT), and evaluate its dosimetric gain in lung stereotactic body radiation therapy (SBRT). In DCAT, the field aperture is continuously conformed to the planning target volume (PTV) with an aperture margin (AM) to compensate for the penumbra effect with gantry rotation. It is a common belief the AM should be positive (or at least 'zero'). However, the radial penumbra width becomes significantly wider because of continuously overlapped beams in arc delivery. Therefore, we hypothesize if the 'negative margin' is applied in the radial direction, it would improve the PTV dose conformation while reducing normal tissue dose. For verification, trial plans were made using the NMT and compared with 'zero margin (ZM)' plans for five lung SBRT cases representing different situations depending on the location of the PTV and organs at risk. All plans met 95% PTV coverage with the prescription dose and spared the spinal cord below the tolerance. Two conventional conformation indices (the ratio of prescription isodose volume to the PTV (CI100) and the ratio of 50% prescription isodose volume to the PTV (CI50)) and a modified conformation index were investigated. The maximum dose at 2 cm from the PTV (Dmax-2cm) and the percent of lung volume receiving 20 Gy (V20) were also evaluated. Another planning simulation was performed with a total of ten randomly selected lung SBRT cases to mimic actual practice. In this simulation, optimization with ZM was first performed and further optimization using the NMT was processed for cases that could not meet a goal of CI100 = 1.2 with the ZM optimization. In all cases, both the CI100 and CI50 values were significantly reduced (overall, 9.4% ± 4.1% and 5.9%± 3.1% for CI100 and CI50, respectively). The modified conformation index values also showed similar improvement (overall, 10.1% ± 5.7% increase). Reduction of Dmax-2cm was also observed in all cases (4.5% ± 2.2%). V20 values decreased in all cases but one (5.7% ± 3.9%, excluding the increased case). In the random group simulation, it was possible to achieve the goal with just one NMT trial for five out of six cases that did not meet the goal in the ZM optimization. Interestingly, however, one case needed as many as six iterations to get the CI100 = 1.2 goal. The NMT turned out to be an effective planning strategy that could bring significant improvement of dose conformation. The NMT can be easily implemented in most clinics with no prerequisite.

Improving IMRT delivery efficiency with reweighted L1-minimization for inverse planning.

PURPOSE: This study presents an improved technique to further simplify the fluence-map in intensity modulated radiation therapy (IMRT) inverse planning, thereby reducing plan complexity and improving delivery efficiency, while maintaining the plan quality. METHODS: First-order total-variation (TV) minimization (min.) based on L1-norm has been proposed to reduce the complexity of fluence-map in IMRT by generating sparse fluence-map variations. However, with stronger dose sparing to the critical structures, the inevitable increase in the fluence-map complexity can lead to inefficient dose delivery. Theoretically, L0-min. is the ideal solution for the sparse signal recovery problem, yet practically intractable due to its nonconvexity of the objective function. As an alternative, the authors use the iteratively reweighted L1-min. technique to incorporate the benefits of the L0-norm into the tractability of L1-min. The weight multiplied to each element is inversely related to the magnitude of the corresponding element, which is iteratively updated by the reweighting process. The proposed penalizing process combined with TV min. further improves sparsity in the fluence-map variations, hence ultimately enhancing the delivery efficiency. To validate the proposed method, this work compares three treatment plans obtained from quadratic min. (generally used in clinic IMRT), conventional TV min., and our proposed reweighted TV min. techniques, implemented by a large-scale L1-solver (template for first-order conic solver), for five patient clinical data. Criteria such as conformation number (CN), modulation index (MI), and estimated treatment time are employed to assess the relationship between the plan quality and delivery efficiency. RESULTS: The proposed method yields simpler fluence-maps than the quadratic and conventional TV based techniques. To attain a given CN and dose sparing to the critical organs for 5 clinical cases, the proposed method reduces the number of segments by 10-15 and 30-35, relative to TV min. and quadratic min. based plans, while MIs decreases by about 20%-30% and 40%-60% over the plans by two existing techniques, respectively. With such conditions, the total treatment time of the plans obtained from our proposed method can be reduced by 12-30 s and 30-80 s mainly due to greatly shorter multileaf collimator (MLC) traveling time in IMRT step-and-shoot delivery. CONCLUSIONS: The reweighted L1-minimization technique provides a promising solution to simplify the fluence-map variations in IMRT inverse planning. It improves the delivery efficiency by reducing the entire segments and treatment time, while maintaining the plan quality in terms of target conformity and critical structure sparing.

Accuracy of volume measurement using 3D ultrasound and development of CT-3D US image fusion algorithm for prostate cancer radiotherapy.

PURPOSE: To evaluate the accuracy of measuring volumes using three-dimensional ultrasound (3D US), and to verify the feasibility of the replacement of CT-MR fusion images with CT-3D US in radiotherapy treatment planning. METHODS: Phantoms, consisting of water, contrast agent, and agarose, were manufactured. The volume was measured using 3D US, CT, and MR devices. A CT-3D US and MR-3D US image fusion software was developed using the Insight Toolkit library in order to acquire three-dimensional fusion images. The quality of the image fusion was evaluated using metric value and fusion images. RESULTS: Volume measurement, using 3D US, shows a 2.8 ± 1.5% error, 4.4 ± 3.0% error for CT, and 3.1 ± 2.0% error for MR. The results imply that volume measurement using the 3D US devices has a similar accuracy level to that of CT and MR. Three-dimensional image fusion of CT-3D US and MR-3D US was successfully performed using phantom images. Moreover, MR-3D US image fusion was performed using human bladder images. CONCLUSIONS: 3D US could be used in the volume measurement of human bladders and prostates. CT-3D US image fusion could be used in monitoring the target position in each fraction of external beam radiation therapy. Moreover, the feasibility of replacing the CT-MR image fusion to the CT-3D US in radiotherapy treatment planning was verified.

Long-term results of forward intensity-modulated radiation therapy for patients with early-stage breast cancer.

PURPOSE: To observe long-term clinical outcomes for patients with early-stage breast cancer treated with forward intensity-modulated radiation therapy (IMRT), including local control and clinical toxicities. MATERIALS AND METHODS: We retrospectively analyzed a total of 214 patients with stage I-II breast cancer who were treated with breast conserving surgery followed by adjuvant breast radiation therapy between 2001 and 2008. All patients were treated using forward IMRT. The whole breast was irradiated to a dose of 50 to 50.4 Gy followed by an 8 to 12 Gy electron boost to the surgical bed. RESULTS: The median age was 46 years (range, 21 to 82 years) and the medial follow-up time was 7.3 years (range, 2.4 to 11.7 years). Stage T1 was 139 (65%) and T2 was 75 (35%), respectively. Ipsilateral breast recurrence was observed in 3 patients. The 5- and 10-year local control rates were 99.1% and 97.8%, respectively. The cosmetic outcome was evaluated according to the Harvard scale and 89.4% of patients were scored as excellent or good. CONCLUSION: The whole breast radiation therapy as an adjuvant treatment using a forward IMRT technique showed excellent long-term local control as well as favorable outcomes of toxicity and cosmesis.

Serum levels of zinc, calcium, and iron are associated with the risk of preeclampsia in pregnant women.

Evidence indicates that various elements, including antioxidant minerals, might play an important role in preeclampsia (PE). This study was carried out to investigate the relationship between serum mineral levelsof zinc, calcium, iron, and selenium and the risk of preeclampsia in Korean women. Twenty-nine normal controls and 30 women with preeclampsia were recruited for the study. Preeclampsia was defined as having high blood pressure (≥ 140/90 mm Hg after 20 weeks gestation) and proteinuria (≥300 mg/24 hours). Serum mineral content was determined by instrumental neutron activation analysis. Serum zinc (P < .0001) and calcium (P = .0188) levels were lower in women with preeclampsia than those of normal women, while serum iron was significantly higher in women with preeclampsia (P = .0045). The odds ratio for preeclampsia was lower in women with higher serum zinc levels than those with lower levels after adjustment for age, height, and weight before delivery (P < .0001). The adjusted odds ratio for preeclampsia also decreased across tertiles of serum calcium concentration (P = .0452). However, there was an increased adjusted odds ratio for preeclampsia across tertiles of serum iron level (P = .0104). These results suggest that levels of serum minerals such as zinc, calcium, and iron may be associated with the risk of preeclampsia in Korean pregnant women.

Dose optimization with first-order total-variation minimization for dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT).

PURPOSE: A new treatment scheme coined as dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT) has recently been proposed to bridge the gap between IMRT and VMAT. By increasing the angular sampling of radiation beams while eliminating dispensable segments of the incident fields, DASSIM-RT is capable of providing improved conformity in dose distributions while maintaining high delivery efficiency. The fact that DASSIM-RT utilizes a large number of incident beams represents a major computational challenge for the clinical applications of this powerful treatment scheme. The purpose of this work is to provide a practical solution to the DASSIM-RT inverse planning problem. METHODS: The inverse planning problem is formulated as a fluence-map optimization problem with total-variation (TV) minimization. A newly released L1-solver, template for first-order conic solver (TFOCS), was adopted in this work. TFOCS achieves faster convergence with less memory usage as compared with conventional quadratic programming (QP) for the TV form through the effective use of conic forms, dual-variable updates, and optimal first-order approaches. As such, it is tailored to specifically address the computational challenges of large-scale optimization in DASSIM-RT inverse planning. Two clinical cases (a prostate and a head and neck case) are used to evaluate the effectiveness and efficiency of the proposed planning technique. DASSIM-RT plans with 15 and 30 beams are compared with conventional IMRT plans with 7 beams in terms of plan quality and delivery efficiency, which are quantified by conformation number (CN), the total number of segments and modulation index, respectively. For optimization efficiency, the QP-based approach was compared with the proposed algorithm for the DASSIM-RT plans with 15 beams for both cases. RESULTS: Plan quality improves with an increasing number of incident beams, while the total number of segments is maintained to be about the same in both cases. For the prostate patient, the conformation number to the target was 0.7509, 0.7565, and 0.7611 with 80 segments for IMRT with 7 beams, and DASSIM-RT with 15 and 30 beams, respectively. For the head and neck (HN) patient with a complicated target shape, conformation numbers of the three treatment plans were 0.7554, 0.7758, and 0.7819 with 75 segments for all beam configurations. With respect to the dose sparing to the critical structures, the organs such as the femoral heads in the prostate case and the brainstem and spinal cord in the HN case were better protected with DASSIM-RT. For both cases, the delivery efficiency has been greatly improved as the beam angular sampling increases with the similar or better conformal dose distribution. Compared with conventional quadratic programming approaches, first-order TFOCS-based optimization achieves far faster convergence and smaller memory requirements in DASSIM-RT. CONCLUSIONS: The new optimization algorithm TFOCS provides a practical and timely solution to the DASSIM-RT or other inverse planning problem requiring large memory space. The new treatment scheme is shown to outperform conventional IMRT in terms of dose conformity to both the targetand the critical structures, while maintaining high delivery efficiency.

Efficient IMRT inverse planning with a new L1-solver: template for first-order conic solver.

Intensity modulated radiation therapy (IMRT) inverse planning using total-variation (TV) regularization has been proposed to reduce the complexity of fluence maps and facilitate dose delivery. Conventionally, the optimization problem with L-1 norm is solved with quadratic programming (QP), which is time consuming and memory expensive due to the second-order Newton update. This study proposes to use a new algorithm, template for first-order conic solver (TFOCS), for fast and memory-efficient optimization in IMRT inverse planning. The TFOCS utilizes dual-variable updates and first-order approaches for TV minimization without the need to compute and store the enlarged Hessian matrix required for Newton update in the QP technique. To evaluate the effectiveness and efficiency of the proposed method, two clinical cases were used for IMRT inverse planning: a head and neck case and a prostate case. For comparison, the conventional QP-based method for the TV form was adopted to solve the fluence map optimization problem in the above two cases. The convergence criteria and algorithm parameters were selected to achieve similar dose conformity for a fair comparison between the two methods. Compared with conventional QP-based approach, the proposed TFOCS-based method shows a remarkable improvement in computational efficiency for fluence map optimization, while maintaining the conformal dose distribution. Compared with QP-based algorithms, the computational speed using TFOCS for fluence optimization is increased by a factor of 4 to 6, and at the same time the memory requirement is reduced by a factor of 3 to 4. Therefore, TFOCS provides an effective, fast and memory-efficient method for IMRT inverse planning. The unique features of the approach should be particularly important in inverse planning involving a large number of beams, such as in VMAT and dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT).

Scatter correction in cone-beam CT via a half beam blocker technique allowing simultaneous acquisition of scatter and image information.

PURPOSE: X-ray scatter incurred to detectors degrades the quality of cone-beam computed tomography (CBCT) and represents a problem in volumetric image guided and adaptive radiation therapy. Several methods using a beam blocker for the estimation and subtraction of scatter have been proposed. However, due to missing information resulting from the obstruction of the blocker, such methods require dual scanning or dynamically moving blocker to obtain a complete volumetric image. Here, we propose a half beam blocker-based approach, in conjunction with a total variation (TV) regularized Feldkamp-Davis-Kress (FDK) algorithm, to correct scatter-induced artifacts by simultaneously acquiring image and scatter information from a single-rotation CBCT scan. METHODS: A half beam blocker, comprising lead strips, is used to simultaneously acquire image data on one side of the projection data and scatter data on the other half side. One-dimensional cubic B-Spline interpolation/extrapolation is applied to derive patient specific scatter information by using the scatter distributions on strips. The estimated scatter is subtracted from the projection image acquired at the opposite view. With scatter-corrected projections where this subtraction is completed, the FDK algorithm based on a cosine weighting function is performed to reconstruct CBCT volume. To suppress the noise in the reconstructed CBCT images produced by geometric errors between two opposed projections and interpolated scatter information, total variation regularization is applied by a minimization using a steepest gradient descent optimization method. The experimental studies using Catphan504 and anthropomorphic phantoms were carried out to evaluate the performance of the proposed scheme. RESULTS: The scatter-induced shading artifacts were markedly suppressed in CBCT using the proposed scheme. Compared with CBCT without a blocker, the nonuniformity value was reduced from 39.3% to 3.1%. The root mean square error relative to values inside the regions of interest selected from a benchmark scatter free image was reduced from 50 to 11.3. The TV regularization also led to a better contrast-to-noise ratio. CONCLUSIONS: An asymmetric half beam blocker-based FDK acquisition and reconstruction technique has been established. The proposed scheme enables simultaneous detection of patient specific scatter and complete volumetric CBCT reconstruction without additional requirements such as prior images, dual scans, or moving strips.

Improved compressed sensing-based cone-beam CT reconstruction using adaptive prior image constraints.

Volumetric cone-beam CT (CBCT) images are acquired repeatedly during a course of radiation therapy and a natural question to ask is whether CBCT images obtained earlier in the process can be utilized as prior knowledge to reduce patient imaging dose in subsequent scans. The purpose of this work is to develop an adaptive prior image constrained compressed sensing (APICCS) method to solve this problem. Reconstructed images using full projections are taken on the first day of radiation therapy treatment and are used as prior images. The subsequent scans are acquired using a protocol of sparse projections. In the proposed APICCS algorithm, the prior images are utilized as an initial guess and are incorporated into the objective function in the compressed sensing (CS)-based iterative reconstruction process. Furthermore, the prior information is employed to detect any possible mismatched regions between the prior and current images for improved reconstruction. For this purpose, the prior images and the reconstructed images are classified into three anatomical regions: air, soft tissue and bone. Mismatched regions are identified by local differences of the corresponding groups in the two classified sets of images. A distance transformation is then introduced to convert the information into an adaptive voxel-dependent relaxation map. In constructing the relaxation map, the matched regions (unchanged anatomy) between the prior and current images are assigned with smaller weight values, which are translated into less influence on the CS iterative reconstruction process. On the other hand, the mismatched regions (changed anatomy) are associated with larger values and the regions are updated more by the new projection data, thus avoiding any possible adverse effects of prior images. The APICCS approach was systematically assessed by using patient data acquired under standard and low-dose protocols for qualitative and quantitative comparisons. The APICCS method provides an effective way for us to enhance the image quality at the matched regions between the prior and current images compared to the existing PICCS algorithm. Compared to the current CBCT imaging protocols, the APICCS algorithm allows an imaging dose reduction of 10-40 times due to the greatly reduced number of projections and lower x-ray tube current level coming from the low-dose protocol.