1
|
Schröder J, Censor Y, Süss P, Küfer KH. Immunity to Increasing Condition Numbers of Linear Superiorization versus Linear Programming. ARXIV 2024:arXiv:2407.18709v1. [PMID: 39108290 PMCID: PMC11302679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Figures] [Subscribe] [Scholar Register] [Indexed: 08/13/2024]
Abstract
Given a family of linear constraints and a linear objective function one can consider whether to apply a Linear Programming (LP) algorithm or use a Linear Superiorization (LinSup) algorithm on this data. In the LP methodology one aims at finding a point that fulfills the constraints and has the minimal value of the objective function over these constraints. The Linear Superiorization approach considers the same data as linear programming problems but instead of attempting to solve those with linear optimization methods it employs perturbation resilient feasibility-seeking algorithms and steers them toward feasible points with reduced (not necessarily minimal) objective function values. Previous studies compared LP and LinSup in terms of their respective outputs and the resources they use. In this paper we investigate these two approaches in terms of their sensitivity to condition numbers of the system of linear constraints. Condition numbers are a measure for the impact of deviations in the input data on the output of a problem and, in particular, they describe the factor of error propagation when given wrong or erroneous data. Therefore, the ability of LP and LinSup to cope with increased condition numbers, thus with ill-posed problems, is an important matter to consider which was not studied until now. We investigate experimentally the advantages and disadvantages of both LP and LinSup on examplary problems of linear programming with multiple condition numbers and different problem dimensions.
Collapse
Affiliation(s)
- Jan Schröder
- Optimization Department, Fraunhofer - ITWM, Kaiserslautern 67663, Germany
| | - Yair Censor
- Department of Mathematics, University of Haifa, Mt. Carmel, Haifa 3498838, Israel
| | - Philipp Süss
- Optimization Department, Fraunhofer - ITWM, Kaiserslautern 67663, Germany
| | - Karl-Heinz Küfer
- Optimization Department, Fraunhofer - ITWM, Kaiserslautern 67663, Germany
| |
Collapse
|
2
|
Barkmann F, Censor Y, Wahl N. Superiorization of projection algorithms for linearly constrained inverse radiotherapy treatment planning. Front Oncol 2023; 13:1238824. [PMID: 38033492 PMCID: PMC10685292 DOI: 10.3389/fonc.2023.1238824] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Accepted: 09/18/2023] [Indexed: 12/02/2023] Open
Abstract
Objective We apply the superiorization methodology to the constrained intensity-modulated radiation therapy (IMRT) treatment planning problem. Superiorization combines a feasibility-seeking projection algorithm with objective function reduction: The underlying projection algorithm is perturbed with gradient descent steps to steer the algorithm towards a solution with a lower objective function value compared to one obtained solely through feasibility-seeking. Approach Within the open-source inverse planning toolkit matRad, we implement a prototypical algorithmic framework for superiorization using the well-established Agmon, Motzkin, and Schoenberg (AMS) feasibility-seeking projection algorithm and common nonlinear dose optimization objective functions. Based on this prototype, we apply superiorization to intensity-modulated radiation therapy treatment planning and compare it with (i) bare feasibility-seeking (i.e., without any objective function) and (ii) nonlinear constrained optimization using first-order derivatives. For these comparisons, we use the TG119 water phantom, the head-and-neck and the prostate patient of the CORT dataset. Main results Bare feasibility-seeking with AMS confirms previous studies, showing it can find solutions that are nearly equivalent to those found by the established piece-wise least-squares optimization approach. The superiorization prototype solved the linearly constrained planning problem with similar dosimetric performance to that of a general-purpose nonlinear constrained optimizer while showing smooth convergence in both constraint proximity and objective function reduction. Significance Superiorization is a useful alternative to constrained optimization in radiotherapy inverse treatment planning. Future extensions with other approaches to feasibility-seeking, e.g., with dose-volume constraints and more sophisticated perturbations, may unlock its full potential for high performant inverse treatment planning.
Collapse
Affiliation(s)
- Florian Barkmann
- Institute for Machine Learning, Department of Computer Science, ETH Zürich, Zurich, Switzerland
- Department of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Baden-Württemberg, Germany
- Heidelberg Institute for Radiation Oncology (HIRO) and National Center for Radiation Research in Oncology (NCRO), Heidelberg, Baden-Württemberg, Germany
| | - Yair Censor
- Department of Mathematics, Faculty of Natural Sciences, University of Haifa, Haifa, Israel
| | - Niklas Wahl
- Department of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg, Baden-Württemberg, Germany
- Heidelberg Institute for Radiation Oncology (HIRO) and National Center for Radiation Research in Oncology (NCRO), Heidelberg, Baden-Württemberg, Germany
| |
Collapse
|
3
|
Yu S, Censor Y, Jiang M, Luo G. Per-RMAP: Feasibility-Seeking and Superiorization Methods for Floorplanning with I/O Assignment. INTERNATIONAL SYMPOSIUM OF ELECTRONICS DESIGN AUTOMATION 2023; 2023:286-291. [PMID: 38988948 PMCID: PMC11236085 DOI: 10.1109/iseda59274.2023.10218694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/12/2024]
Abstract
The feasibility-seeking approach provides a systematic scheme to manage and solve complex constraints for continuous problems, and we explore it for the floorplanning problems with increasingly heterogeneous constraints. The classic legality constraints can be formulated as the union of convex sets. However, the convergence of conventional projection-based algorithms is not guaranteed when the constraints sets are non-convex, which is the case with unions of convex sets. In this work, we propose a resetting strategy to greatly eliminate the divergence issue of the projection-based algorithm for the feasibility-seeking formulation. Furthermore, the superiorization methodology (SM), which lies between feasibility-seeking and constrained optimization, is firstly applied to floorplanning. The SM uses perturbations to steer the iterates of a feasibility-seeking algorithm to a feasible solution with shorter total wirelength. The proposed algorithmic flow is extendable to tackle various constraints and variants of floorplanning problems, e.g., floorplanning with I/O assignment problems. We have evaluated the proposed algorithm on the MCNC benchmarks. We can obtain legal floorplans only two times slower than the branch-and-bound method in its current prototype using MATLAB, with only 3% wirelength inferior to the optimal results. We evaluate the effectiveness of the algorithmic flow by considering the constraints of I/O assignment, and our algorithm achieves 8% improvement on wirelength.
Collapse
Affiliation(s)
- Shan Yu
- Department of Information and Computational Sciences, School of Mathematical Sciences, Peking University, Beijing, China
| | - Yair Censor
- Department of Mathematics, University of Haifa, Haifa, Israel
| | - Ming Jiang
- Department of Information and Computational Sciences, School of Mathematical Sciences, Peking University, Beijing, China
| | - Guojie Luo
- National Key Laboratory for Multimedia Information Processing, School of Computer Science, Peking University, Beijing, China
- Center for Energy-efficient Computing and Applications, Peking University, Beijing, China
| |
Collapse
|
4
|
DeJongh DF, DeJongh EA. An Iterative Least Squares Method for Proton CT Image Reconstruction. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2022; 6:304-312. [PMID: 36061217 PMCID: PMC9432481 DOI: 10.1109/trpms.2021.3079140] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Clinically useful proton Computed Tomography images will rely on algorithms to find the three-dimensional proton stopping power distribution that optimally fits the measured proton data. We present a least squares iterative method with many features to put proton imaging into a more quantitative framework. These include the definition of a unique solution that optimally fits the protons, the definition of an iteration vector that takes into account proton measurement uncertainties, the definition of an optimal step size for each iteration individually, the ability to simultaneously optimize the step sizes of many iterations, the ability to divide the proton data into arbitrary numbers of blocks for parallel processing and use of graphical processing units, and the definition of stopping criteria to determine when to stop iterating. We find that it is possible, for any object being imaged, to provide assurance that the image is quantifiably close to an optimal solution, and the optimization of step sizes reduces the total number of iterations required for convergence. We demonstrate the use of these algorithms on real data.
Collapse
Affiliation(s)
- Don F. DeJongh
- ProtonVDA LLC, 1700 Park St Ste 208, Naperville, IL 60563 USA
| | | |
Collapse
|
5
|
Censor Y, Schubert KE, Schulte RW. Developments in Mathematical Algorithms and Computational Tools for Proton CT and Particle Therapy Treatment Planning. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2022. [DOI: 10.1109/trpms.2021.3107322] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
6
|
Meyer S, Pinto M, Parodi K, Gianoli C. The impact of path estimates in iterative ion CT reconstructions for clinical-like cases. Phys Med Biol 2021; 66. [PMID: 33765672 DOI: 10.1088/1361-6560/abf1ff] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Accepted: 03/25/2021] [Indexed: 11/11/2022]
Abstract
Ion computed tomography (CT) promises to mitigate range uncertainties inherent in the conversion of x-ray Hounsfield units into ion relative stopping power (RSP) for ion beam therapy treatment planning. To improve accuracy and spatial resolution of ion CT by accounting for statistical multiple Coulomb scattering deflection of the ion trajectories from a straight line path (SLP), the most likely path (MLP) and the cubic spline path (CSP) have been proposed. In this work, we use FLUKA Monte Carlo simulations to investigate the impact of these path estimates in iterative tomographic reconstruction algorithms for proton, helium and carbon ions. To this end the ordered subset simultaneous algebraic reconstruction technique was used and coupled with a total variation superiorization (TVS). We evaluate the image quality and dose calculation accuracy in proton therapy treatment planning of cranial patient anatomies. CSP and MLP generally yielded nearly equal image quality with an average RSP relative error improvement over the SLP of 0.6%, 0.3% and 0.3% for proton, helium and carbon ion CT, respectively. Bone and low density materials have been identified as regions of largest enhancement in RSP accuracy. Nevertheless, only minor differences in dose calculation results were observed between the different models and relative range errors of better than 0.5% were obtained in all cases. Largest improvements were found for proton CT in complex scenarios with strong heterogeneities along the beam path. The additional TVS provided substantially reduced image noise, resulting in improved image quality in particular for soft tissue regions. Employing the CSP and MLP for iterative ion CT reconstructions enabled improved image quality over the SLP even in realistic and heterogeneous patient anatomy. However, only limited benefit in dose calculation accuracy was obtained even though an ideal detector system was simulated.
Collapse
Affiliation(s)
- Sebastian Meyer
- Department of Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States of America.,Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Garching b. München, Germany
| | - Marco Pinto
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Garching b. München, Germany
| | - Katia Parodi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Garching b. München, Germany.,Shared senior authorship
| | - Chiara Gianoli
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Garching b. München, Germany.,Shared senior authorship
| |
Collapse
|
7
|
Kaser S, Bergauer T, Birkfellner W, Burker A, Georg D, Hatamikia S, Hirtl A, Irmler C, Pitters F, Ulrich-Pur F. First application of the GPU-based software framework TIGRE for proton CT image reconstruction. Phys Med 2021; 84:56-64. [PMID: 33848784 DOI: 10.1016/j.ejmp.2021.03.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Revised: 02/26/2021] [Accepted: 03/03/2021] [Indexed: 11/28/2022] Open
Abstract
In proton therapy, the knowledge of the proton stopping power, i.e. the energy deposition per unit length within human tissue, is essential for accurate treatment planning. One suitable method to directly measure the stopping power is proton computed tomography (pCT). Due to the proton interaction mechanisms in matter, pCT image reconstruction faces some challenges: the unique path of each proton has to be considered separately in the reconstruction process adding complexity to the reconstruction problem. This study shows that the GPU-based open-source software toolkit TIGRE, which was initially intended for X-ray CT reconstruction, can be applied to the pCT image reconstruction problem using a straight line approach for the proton path. This simplified approach allows for reconstructions within seconds. To validate the applicability of TIGRE to pCT, several Monte Carlo simulations modeling a pCT setup with two Catphan® modules as phantoms were performed. Ordered-Subset Simultaneous Algebraic Reconstruction Technique (OS-SART) and Adaptive-Steepest-Descent Projection Onto Convex Sets (ASD-POCS) were used for image reconstruction. Since the accuracy of the approach is limited by the straight line approximation of the proton path, requirements for further improvement of TIGRE for pCT are addressed.
Collapse
Affiliation(s)
- Stefanie Kaser
- Institute of High Energy Physics, Austrian Academy of Sciences, Vienna, Austria.
| | - Thomas Bergauer
- Institute of High Energy Physics, Austrian Academy of Sciences, Vienna, Austria
| | - Wolfgang Birkfellner
- Center for Medical Physics and Biomedical Engineering, Medical University of Vienna, Vienna, Austria
| | | | - Dietmar Georg
- Department of Radiation Oncology, Medical University of Vienna, Vienna, Austria; MedAustron Ion Therapy Center, Wiener Neustadt, Austria
| | - Sepideh Hatamikia
- Center for Medical Physics and Biomedical Engineering, Medical University of Vienna, Vienna, Austria; Austrian Center for Medical Innovation and Technology, Wiener Neustadt, Austria
| | | | - Christian Irmler
- Institute of High Energy Physics, Austrian Academy of Sciences, Vienna, Austria
| | - Florian Pitters
- Institute of High Energy Physics, Austrian Academy of Sciences, Vienna, Austria
| | - Felix Ulrich-Pur
- Institute of High Energy Physics, Austrian Academy of Sciences, Vienna, Austria
| |
Collapse
|
8
|
Abstract
Objective Particle radiobiology has contributed new understanding of radiation safety and underlying mechanisms of action to radiation oncology for the treatment of cancer, and to planning of radiation protection for space travel. This manuscript will highlight the significance of precise physical and biologically effective dosimetry to this translational research for the benefit of human health. This review provides a brief snapshot of the evolving scientific basis for, and the complex current global status, and remaining challenges of hadron therapy for the treatment of cancer. The need for particle radiobiology for risk planning in return missions to the Moon, and exploratory deep-space missions to Mars and beyond are also discussed. Methods Key lessons learned are summarized from an impressive collective literature published by an international cadre of multidisciplinary experts in particle physics, radiation chemistry, medical physics of imaging and treatment planning, molecular, cellular, tissue radiobiology, biology of microgravity and other stressors, theoretical modeling of biophysical data, and clinical results with accelerator-produced particle beams. Results Research pioneers, many of whom were Nobel laureates, led the world in the discovery of ionizing radiations originating from the Earth and the Cosmos. Six radiation pioneers led the way to hadron therapy and the study of charged particles encountered in outer space travel. Worldwide about 250,000 patients have been treated for cancer, or other lesions such as arteriovenous malformations in the brain between 1954 and 2019 with charged particle radiotherapy, also known as hadron therapy. The majority of these patients (213,000) were treated with proton beams, but approximately 32,000 were treated with carbon ion radiotherapy. There are 3500 patients who have been treated with helium, pions, neon or other ions. There are currently 82 facilities operating to provide ion beam clinical treatments. Of these, only 13 facilities located in Asia and Europe are providing carbon ion beams for preclinical, clinical, and space research. There are also numerous particle physics accelerators worldwide capable of producing ion beams for research, but not currently focused on treating patients with ion beam therapy but are potentially available for preclinical and space research. Approximately, more than 550 individuals have traveled into Lower Earth Orbit (LEO) and beyond and returned to Earth. Conclusion Charged particle therapy with controlled beams of protons and carbon ions have significantly impacted targeted cancer therapy, eradicated tumors while sparing normal tissue toxicities, and reduced human suffering. These modalities still require further optimization and technical refinements to reduce cost but should be made available to everyone in need worldwide. The exploration of our Universe in space travel poses the potential risk of exposure to uncontrolled charged particles. However, approaches to shield and provide countermeasures to these potential radiation hazards in LEO have allowed an amazing number of discoveries currently without significant life-threatening medical consequences. More basic research with components of the Galactic Cosmic Radiation field are still required to assure safety involving space radiations and combined stressors with microgravity for exploratory deep space travel. Advances in knowledge The collective knowledge garnered from the wealth of available published evidence obtained prior to particle radiation therapy, or to space flight, and the additional data gleaned from implementing both endeavors has provided many opportunities for heavy ions to promote human health.
Collapse
|
9
|
Dickmann J, Rit S, Pankuch M, Johnson RP, Schulte RW, Parodi K, Dedes G, Landry G. An optimization algorithm for dose reduction with fluence‐modulated proton CT. Med Phys 2020; 47:1895-1906. [DOI: 10.1002/mp.14084] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Revised: 01/30/2020] [Accepted: 02/05/2020] [Indexed: 01/12/2023] Open
Affiliation(s)
- J. Dickmann
- Department of Medical Physics Faculty of Physics Ludwig‐Maximilians‐Universität München Am Coulombwall 1 85748 Garching b. München Germany
| | - S. Rit
- Univ Lyon INSA‐Lyon Université Claude Bernard Lyon 1 UJM‐Saint Étienne CNRS, Inserm CREATIS UMR 5220 U1206 F‐69373 Lyon France
| | - M. Pankuch
- Northwestern Medicine Chicago Proton Center Warrenville IL 60555 USA
| | - R. P. Johnson
- Department of Physics University of California Santa Cruz Santa Cruz CA 95064 USA
| | - R. W. Schulte
- Division of Biomedical Engineering Sciences Loma Linda University Loma Linda CA 92354 USA
| | - K. Parodi
- Department of Medical Physics Faculty of Physics Ludwig‐Maximilians‐Universität München Am Coulombwall 1 85748 Garching b. München Germany
| | - G. Dedes
- Department of Medical Physics Faculty of Physics Ludwig‐Maximilians‐Universität München Am Coulombwall 1 85748 Garching b. München Germany
| | - G. Landry
- Department of Medical Physics Faculty of Physics Ludwig‐Maximilians‐Universität München Am Coulombwall 1 85748 Garching b. München Germany
- Department of Radiation Oncology University Hospital, LMU Munich 81377 Munich Germany
- German Cancer Consortium (DKTK) 81377 Munich Germany
| |
Collapse
|