Beam Orientation Optimization for Intensity Modulated Radiation Therapy using Adaptive l1 Minimization

Beam orientation optimization (BOO) is a key component in the process of IMRT treatment planning. It determines to what degree one can achieve a good treatment plan quality in the subsequent plan optimization process. In this paper, we have developed a BOO algorithm via adaptive l_1 minimization. Specifically, we introduce a sparsity energy function term into our model which contains weighting factors for each beam angle adaptively adjusted during the optimization process. Such an energy term favors small number of beam angles. By optimizing a total energy function containing a dosimetric term and the sparsity term, we are able to identify the unimportant beam angles and gradually remove them without largely sacrificing the dosimetric objective. In one typical prostate case, the convergence property of our algorithm, as well as the how the beam angles are selected during the optimization process, is demonstrated. Fluence map optimization (FMO) is then performed based on the optimized beam angles. The resulted plan quality is presented and found to be better than that obtained from unoptimized (equiangular) beam orientations. We have further systematically validated our algorithm in the contexts of 5-9 coplanar beams for 5 prostate cases and 1 head and neck case. For each case, the final FMO objective function value is used to compare the optimized beam orientations and the equiangular ones. It is found that, our BOO algorithm can lead to beam configurations which attain lower FMO objective function values than corresponding equiangular cases, indicating the effectiveness of our BOO algorithm.Comment: 19 pages, 2 tables, and 5 figure

Full modelling of high-intensity focused ultrasound and thermal heating in the kidney using realistic patient models

Objective: High-intensity focused ultrasound (HIFU) therapy can be used for non-invasive treatment of kidney (renal) cancer, but the clinical outcomes have been variable. In this study, the efficacy of renal HIFU therapy was studied using nonlinear acoustic and thermal simulations in three patients. Methods: The acoustic simulations were conducted with and without refraction in order to investigate its effect on the shape, size and pressure distribution at the focus. The values for the attenuation, sound speed, perfusion and thermal conductivity of the kidney were varied over the reported ranges to determine the effect of variability on heating. Furthermore, the phase aberration was studied in order to quantify the underlying phase shifts using a second order polynomial function. Results: The ultrasound field intensity was found to drop on average 11.1 dB with refraction and 6.4 dB without refraction. Reflection at tissue interfaces was found to result in a loss less than 0.1 dB. Focal point splitting due to refraction significantly reduced the heating efficacy. Perfusion did not have a large effect on heating during short sonication durations. Small changes in temperature were seen with varying attenuation and thermal conductivity, but no visible changes were present with sound speed variations. The aberration study revealed an underlying trend in the spatial distribution of the phase shifts. Conclusion: The results show that the efficacy of HIFU therapy in the kidney could be improved with aberration correction. Significance: A method is proposed by which patient specific pre-treatment calculations could be used to overcome the aberration and therefore make ultrasound treatment possible.Comment: Journal paper, IEEE Transactions on Biomedical Engineering (2018

Proton tracking in a high-granularity Digital Tracking Calorimeter for proton CT purposes

Radiation therapy with protons as of today utilizes information from x-ray CT in order to estimate the proton stopping power of the traversed tissue in a patient. The conversion from x-ray attenuation to proton stopping power in tissue introduces range uncertainties of the order of 2-3% of the range, uncertainties that are contributing to an increase of the necessary planning margins added to the target volume in a patient. Imaging methods and modalities, such as Dual Energy CT and proton CT, have come into consideration in the pursuit of obtaining an as good as possible estimate of the proton stopping power. In this study, a Digital Tracking Calorimeter is benchmarked for proof-of-concept for proton CT purposes. The Digital Tracking Calorimeteris applied for reconstruction of the tracks and energies of individual high energy protons. The presented prototype forms the basis for a proton CT system using a single technology for tracking and calorimetry. This advantage simplifies the setup and reduces the cost of a proton CT system assembly, and it is a unique feature of the Digital Tracking Calorimeter. Data from the AGORFIRM beamline at KVI-CART in Groningen in the Netherlands and Monte Carlo simulation results are used to in order to develop a tracking algorithm for the estimation of the residual ranges of a high number of concurrent proton tracks. The range of the individual protons can at present be estimated with a resolution of 4%. The readout system for this prototype is able to handle an effective proton frequency of 1 MHz by using 500 concurrent proton tracks in each readout frame, which is at the high end range of present similar prototypes. A future further optimized prototype will enable a high-speed and more accurate determination of the ranges of individual protons in a therapeutic beam.Comment: 21 pages, 8 figure

A Monte Carlo study of organ and effective doses of cone beam computed tomography (CBCT) scans in radiotherapy

Cone-beam CT (CBCT) scans utilized for image guided radiation therapy (IGRT) procedures have become an essential part of radiotherapy. The aim of this study was to assess organ and effective doses resulting from new CBCT scan protocols (head, thorax, and pelvis) released with a software upgrade of the kV on-board-imager (OBI) system. Influence of the scan parameters that were changed in the new protocols on the patient dose was also investigated. Organ and effective doses for protocols of the new software (V2.5) and a previous version (V1.6) were assessed using Monte Carlo (MC) simulations for the International Commission on Radiological Protection (ICRP) adult male and female reference computational phantoms. The number of projections and the mAs values were increased and the size of the scan field was extended in the new protocols. Influence of these changes on organ and effective doses of the scans was investigated. The OBI system was modelled in EGSnrc/BEAMnrc, and organ doses were estimated using EGSnrc/DOSXYZnrc. The MC model was benchmarked against experimental measurements. Organ doses resulting from the V2.5 protocols were higher than those of V1.6 for organs that were partially or fully inside the scans fields, and increased by (3 to 13)%, (10 to 77)%, and (13 to 21)% for the head, thorax, and pelvis protocols for both phantoms, respectively. As a result, effective doses rose by 14%, 17%, and 16% for the male phantom, and 13%, 18%, and 17% for the female phantom for the three scan protocols, respectively. The scan field extension for the V2.5 protocols contributed significantly in the dose increases, especially for organs that were partially irradiated such as the thyroid in head and thorax scans and colon in the pelvic scan. The contribution of the mAs values and projection numbers was minimal in the dose increases, up to 2.5%. The field size extension plays a major role in improving the treatment output by including more markers in the field of view to match between CBCT and CT images and hence setting up the patient precisely. Therefore, a trade-off between the risk and benefits of CBCT scans should be considered, and the dose increases should be monitored. Several recommendations have been made for optimization of the patient dose involved for IGRT procedures

Bounded perturbation resilience of projected scaled gradient methods

We investigate projected scaled gradient (PSG) methods for convex minimization problems. These methods perform a descent step along a diagonally scaled gradient direction followed by a feasibility regaining step via orthogonal projection onto the constraint set. This constitutes a generalized algorithmic structure that encompasses as special cases the gradient projection method, the projected Newton method, the projected Landweber-type methods and the generalized Expectation-Maximization (EM)-type methods. We prove the convergence of the PSG methods in the presence of bounded perturbations. This resilience to bounded perturbations is relevant to the ability to apply the recently developed superiorization methodology to PSG methods, in particular to the EM algorithm.Comment: Computational Optimization and Applications, accepted for publicatio