Treatment Planning Optimization for VMAT, Tomotherapy and Cyberknife
|
|
- Amberlynn Cornelia Chandler
- 5 years ago
- Views:
Transcription
1 Treatment Planning Optimization for VMAT, Tomotherapy and Cyberknife Kerem Akartunalı Department of Management Science Strathclyde Business School Joint work with: Vicky Mak-Hau and Thu Tran 14 July 2015
2 Outline 1 Introduction 2 Solution Methods Improving the Formulation Methodologies 3 Computational Results
3 Outline 1 Introduction 2 Solution Methods Improving the Formulation Methodologies 3 Computational Results
4 VMAT: Gantry
5 VMAT: Multi-Leaf Collimator
6 3-D Imaging (Voxels) and Dose
7 IMRT vs. VMAT/Tomotherapy/Cyberknife: Overview Intensity-Modulated Radiation Therapy (IMRT) Technology widely used and been around for a while. Lots of literature, lots of solution methods. Beaming limited to a handful of selected angles only. Volumetric-Modulated Arc Therapy (VMAT), Tomotherapy and Cyberknife More recent technologies and limited literature. VMAT: Gantry rotates 360 degrees-single arc or dual arcs in co-planar manner. Beam on all the time. MLC leaf movement linear to gantry rotation. Tomotherapy: Rotation in a helical manner; simple binary MLC. Cyberknife: Source of radiation is mounted on a robotic arm; various collimators (such as circular).
8 Overview of Decision Problem A number of physical limitations. Shapes between consecutive leaves. Changes in shapes and dosage between consecutive angles. The aim is to optimize treatment plans Minimize treatment time, maximize radiation on tumor, etc. Not too much radiation on sensitive tissue/organs. But also need a minimum level of radiation to tumor. Motivation: Software currently in use do not even generate feasible treatment plans. Focus on minimizing error. Aim: To optimize the fluence and MLC apertures simultaneously. Incorporate all volume limits for dose upper bounds for organs at risk (OAR) and dose lower bounds for planning target volumes (PTV) as hard constraints, but allow to relax either if necessary.
9 Decisions to Make... There are a number of decisions to consider. When do we change shapes of MLC? How do the shapes of MLC look like? How does units of radiation applied change? Discretization of angles/snapshots. This eliminates the questions of when, given sufficient number of snapshots. Matrix-discretization of MLC. Rows already natural, columns also added.
10 Notation: Sets and Parameters - I = {1,..., m} be the index set of the MLC rows. - J = {1,..., n} be the index set of the MLC columns. Each cell (i, j) is called a beamlet or a bixel. Define J = {0, n + 1} J. - V be the index set of all voxels V t targeted tumor. V o organs/sensitive tissue. - L = {(l, r) l, r J, l < r}. - K be the index set of beam-angles. - Dijv k already defined. - L v prescription dose for voxel for v V t. - U v maximum dose allowed for voxel v V o.
11 Notation: Variables - yi(l,r) k {0, 1} 1 if the leaf position in the k-th snapshot and i-th row is (l, r). 0 otherwise. - z k : Number of radiation units used for snapshot k. Physical limitation of beamer: M z k 0 - D v : Total dose of radiation received by voxel v. More to come later about how we calculate this. - x v {0, 1} 1 if tumor voxel v receives at least target amount d of radiation. Note L v is a bare minimum, d is preferable. 0 otherwise.
12 Problem Constraints (l,r) L y k i,(l,r) = 1 i I, k K For each row on each snapshot, we can only choose one pair of right-left leaf positions. r δ 1 r=0 r 1 l=0 y k+1 + n+1 i( l, r) r 1 r=r+δ+1 l=0 i I, (l, r) L, k K y k+1 i( l, r) 1 y k i(l,r) Right-leaf position can move at most δ units between consecutive snapshots. Same is also true for left-leaf position...
13 Problem Constraints (cont d) n+1 r=l+1 r 1 yi(l, r) k + l=0 r=1 l y k (i+1)( l, r) 1 i I, l J, k K
14 Problem Constraints (cont d) r 1 n+1 y k i( l,r) + l=0 r= l+1 l=r n y k (i+1)( l, r) 1 i I, r J, k K
15 Problem Constraints (cont d) z k z k+1 k K Number of unit radiation applied can be changed at most between consecutive snapshots. d v L v d v U v d v dx v v V t v V o v V t Lower- and upper-bounds on dose received by tumor and organs, resp.; is d achieved for tumor voxels?
16 Problem Objective max v V t x v Simply maximize number of tumor cells that receive at least d units of radiation. A more successful treatment plan Lesser number of treatments. There are other objectives used in practice. Minimize total treatment time (to comfort patients + use limited resources as much as possible). Minimize total number of shapes (to minimize technical effort). More to read on these in Ernst, Mak and Mason [2009,2010].
17 Defining Dose Received d v = k K i I j J 1 2 zk Dijv k (l,r) L l<j<r yi(l,r) k + zk+1 D k+1 ijv (l,r) L l<j<r y k+1 i(l,r)
18 Defining Dose Received d v = k K i I j J 1 2 zk Dijv k (l,r) L l<j<r yi(l,r) k + zk+1 D k+1 ijv (l,r) L l<j<r y k+1 i(l,r)
19 Outline 1 Introduction 2 Solution Methods Improving the Formulation Methodologies 3 Computational Results
20 Improving the Formulation: Linearizing D v d v = k K i I j J 1 2 zk Dijv k (l,r) L l<j<r yi(l,r) k + zk+1 D k+1 ijv (l,r) L l<j<r y k+1 i(l,r) All nonlinear terms are bilinear. Define a new variable z ij k to indicate radiation amount for beam angle k and (i, j) th location of MLC. d v = 1 ( ) Dijv k z k ij + D k+1 ijv z k+1 ij 2 k K i I j J
21 Improving the Formulation: Linearizing D v (cont d) To finalize the linearization, we add the following constraints: z ij k M k K, i I, j J (l,r) L l<j<r z ij k z k y k i(l,r) z ij k M 1 + z k ij 0 (l,r) L l<j<r yi(l,r) k + zk k K, i I, j J k K, i I, j J k K, i I, j J
22 Improving the Formulation: Valid Inequalities Proposition d v L v ( d L v )x v (1) is valid and dominates the constraint d v dx v.
23 Improving the Formulation: Valid Inequalities (cont d) Proposition Let Di,l+1,r 1,v k = r 1 j=l+1 Dk ijv. Then, the inequality z k z k ij (l,r) L l j r j min{ M, U v D k i,l+1,r 1,v is valid for the VMAT problem and dominates z ij k M 1 + (l,r) L l<j<r yi(l,r) k + zk }y k i(l,r) (2)
24 Improving the Formulation: Few Notes Linearization simplifies the problem significantly. Added constraints are nice. The proposed inequalities have no additional cost. They are tighter constraints and hence provide more compact formulation. They are facet-defining for some important subproblems (e.g., single-voxel, single-row subproblem). However, even the small problems might be still challenging to solve. So we need practical tools in addition to these theoretical tools.
25 Variable Neighborhood Search Heuristic (VNHS) Generate feasible shapes (y) for snapshots. Randomly assign radiation dosages (z) on each snapshot. Define neighborhoods around y and z. Penalize dose UB/LB violations. Apply guided-search: If # overdose > # underdose; then reduce z. If # underdose > # overdose; then increase z. Also can vary neighborhood size, e.g. proportional to % violated voxels. Random restarts after 50 iterations to avoid solution being trapped. Advantage: Very fast!
26 Lagrangian Relaxation A complex structure, a number of potential candidates, different expectations. LR 1 Relax constraints on max distance, dose lower/upper bounds, and max weight change K subproblems for each snapshot. LR 2 Relax constraints on interleaf, dose lower/upper bounds, linearization provides equal bound to the LP relaxation. LR 3 Relax linearization constraints two subproblems, one for y, and one for z and x. Subgradient optimization on its own is known to be computationally not very efficient. But maybe we can use these for some heuristics?
27 A Heuristic Related to LR 1 (Heur1) Output: K candidate solutions for k K do Solve the problem P1: min{ v V t (L v d v ) + + v V (dv Uv )+ (y, z, z) X rel } ; where X rel is the LR 1 subproblem, for only k. ; Fix shapes for k ; for k = k + 1 to k = K and k = k 1 to k = 1 do Solve the problem P1 with inter-snapshot constraints for only k. ; Fix shapes for k ; end Now all shapes fixed, solved the original problem. ; end
28 A Heuristic Related to LR 3 (Heur2) Idea: Combine subgradient optimization and heuristics. Initialize Lagrangian multipliers α k ij and β k ij using LPR duals. We obtain two separate problems. Pr.1 over x, z, z variables. Pr.2 over y variables. Apply subgradient optimization to update the multipliers. Also take the valid solution of Pr.2, fix all these shapes and run the original problem. Already know this heuristic runs fast.
29 A Centering-Based Heuristic (Heur3) Number of y variables and inter-leaf constraints prohibitive for real-size problems Idea 1: If we knew the location of the center of the opening in a row, we could have simply defined n binary variables (one for each location, either left or right). Idea 2: If we extend this to the whole snapshot (i.e., one center for the whole snapshot), then we can eliminate inter-leaf constraints. Advantage: A significant reduction in the number of variables and constraints (crucial for real-size problems). Disadvantage: The candidate opening patterns for a snapshot limited to only centered patterns (might cause infeasibilities).
30 Outline 1 Introduction 2 Solution Methods Improving the Formulation Methodologies 3 Computational Results
31 Generating Random Problems It is important to test different methodologies on as many different problems as possible. Real test cases are important, but they are scarce. Most papers published have only one or two cases. Fully exact methods suffer on any real-size problem. Optimization is practically impossible. Hence no insight on the performance of a method. Hence we created a random problem generator. We can solve some of these problems to optimality. We can run a lot of quick tests and comparisons. We can also see how difficult problems can get. Available at:
32 Random Problems: An Overview of Complexity p-mlc Size - # Voxels - # Snapshots d = 0.5 max{uv } d = max{uv } p LPR 0.02s LPR 0.02s Opt 0.05s Opt 2.11s p LPR 0.11s LPR 0.14s Opt 5.97s Opt s p LPR 0.52s LPR 12.3s Opt Opt > 8h p LPR 2.2s LPR 3.39s Opt s Opt NA In reality, MLCs discretization 5mm, voxels 4mm, and 5-8 degrees between snapshots. For a 10cm treatment field: MLC: 20x20; Treatment field: 25x25x25; Number of Snapshots:
33 Overall Comparison of Heuristic Methods Problem Set # best overall solutions (max MLC size) Heur1 Heur2 Heur3 GVNS Small (8 8) Medium (12 12) Large (15 15) Very Large Problem Set # instances with no solutions (max MLC size) Heur1 Heur2 Heur3 GVNS Small (8 8) Medium (12 12) Large (15 15) Very Large
34 GVNS: Huge Data Instances 1st sol Overall m n w h d K # Di,j,v k > 0 Time # ITER Time # ITER , 400, , 500, , , 500, , , 687, 500 7, , , 000, , , 000, , , 000, , , , 500, 000 5, , , 000, , , 400, , Best solution with original objective recorded 3 times before 1,000 iterations reached. In comparison to Peng et al. (2012), where highest # Di,j,v k > 0 is 114, 315, 187.
35 Conclusions A unified MIP formulation for VMAT, TomoTherapy, and CyberKnife. Heur 2 is very efficient for small problems (but fails otherwise). Heur3 and GVNS often fail for small problems. Heur1 and Heur3 have similar performance for medium and large problems. GVNS is powerful for larger problems (hence usable in real-life applications). A unified mixed-integer programming model for simultaneous fluence weight and aperture optimization in VMAT, Tomotherapy, and Cyberknife, Computers & Operations Research, Vol. 56, April 2015.
A Unified Mixed-Integer Programming Model for Simultaneous Fluence Weight and Aperture Optimization in VMAT, Tomotherapy, and CyberKnife
A Unified Mixed-Integer Programming Model for Simultaneous Fluence Weight and Aperture Optimization in VMAT, Tomotherapy, and CyberKnife Kerem Akartunalı Department of Management Science, University of
More informationRadiation therapy treatment plan optimization
H. Department of Industrial and Operations Engineering The University of Michigan, Ann Arbor, Michigan MOPTA Lehigh University August 18 20, 2010 Outline 1 Introduction Radiation therapy delivery 2 Treatment
More informationInteractive Treatment Planning in Cancer Radiotherapy
Interactive Treatment Planning in Cancer Radiotherapy Mohammad Shakourifar Giulio Trigila Pooyan Shirvani Ghomi Abraham Abebe Sarah Couzens Laura Noreña Wenling Shang June 29, 212 1 Introduction Intensity
More informationIterative regularization in intensity-modulated radiation therapy optimization. Carlsson, F. and Forsgren, A. Med. Phys. 33 (1), January 2006.
Iterative regularization in intensity-modulated radiation therapy optimization Carlsson, F. and Forsgren, A. Med. Phys. 33 (1), January 2006. 2 / 15 Plan 1 2 3 4 3 / 15 to paper The purpose of the paper
More information7/29/2017. Making Better IMRT Plans Using a New Direct Aperture Optimization Approach. Aim of Radiotherapy Research. Aim of Radiotherapy Research
Making Better IMRT Plans Using a New Direct Aperture Optimization Approach Dan Nguyen, Ph.D. Division of Medical Physics and Engineering Department of Radiation Oncology UT Southwestern AAPM Annual Meeting
More informationAvailable online at ScienceDirect. Procedia Computer Science 100 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 00 (206 ) 644 65 Conference on ENTERprise Information Systems / International Conference on Project MANagement / Conference
More informationAn optimization framework for conformal radiation treatment planning
An optimization framework for conformal radiation treatment planning Jinho Lim Michael C. Ferris Stephen J. Wright David M. Shepard Matthew A. Earl December 2002 Abstract An optimization framework for
More informationA hybrid framework for optimizing beam angles in radiation therapy planning
A hybrid framework for optimizing beam angles in radiation therapy planning Gino J. Lim and Laleh Kardar and Wenhua Cao December 29, 2013 Abstract The purpose of this paper is twofold: (1) to examine strengths
More informationTwo Effective Heuristics for Beam Angle Optimization in Radiation Therapy
Two Effective Heuristics for Beam Angle Optimization in Radiation Therapy Hamed Yarmand, David Craft Department of Radiation Oncology Massachusetts General Hospital and Harvard Medical School, Boston,
More informationOperations Research and Optimization: A Primer
Operations Research and Optimization: A Primer Ron Rardin, PhD NSF Program Director, Operations Research and Service Enterprise Engineering also Professor of Industrial Engineering, Purdue University Introduction
More informationMonaco VMAT. The Next Generation in IMRT/VMAT Planning. Paulo Mathias Customer Support TPS Application
Monaco VMAT The Next Generation in IMRT/VMAT Planning Paulo Mathias Customer Support TPS Application 11.05.2011 Background What is Monaco? Advanced IMRT/VMAT treatment planning system from Elekta Software
More informationLagrangean relaxation - exercises
Lagrangean relaxation - exercises Giovanni Righini Set covering We start from the following Set Covering Problem instance: min z = x + 2x 2 + x + 2x 4 + x 5 x + x 2 + x 4 x 2 + x x 2 + x 4 + x 5 x + x
More informationGPU applications in Cancer Radiation Therapy at UCSD. Steve Jiang, UCSD Radiation Oncology Amit Majumdar, SDSC Dongju (DJ) Choi, SDSC
GPU applications in Cancer Radiation Therapy at UCSD Steve Jiang, UCSD Radiation Oncology Amit Majumdar, SDSC Dongju (DJ) Choi, SDSC Conventional Radiotherapy SIMULATION: Construciton, Dij Days PLANNING:
More informationRadiotherapy Plan Competition TomoTherapy Planning System. Dmytro Synchuk. Ukrainian Center of TomoTherapy, Kirovograd, Ukraine
Radiotherapy Plan Competition 2016 TomoTherapy Planning System Dmytro Synchuk Ukrainian Center of TomoTherapy, Kirovograd, Ukraine Beam Geometry 6MV fan beam 3 jaw options 1.0, 2.5 and 5 cm 64 leaves binary
More informationMinimizing Setup and Beam-On Times in Radiation Therapy
Minimizing Setup and Beam-On Times in Radiation Therapy Nikhil Bansal 1, Don Coppersmith 2, and Baruch Schieber 1 1 IBM T.J. Watson Research Center, P.O. Box 218, Yorktown Heights, NY 10598, {nikhil,sbar}@us.ibm.com
More informationarxiv: v1 [physics.med-ph] 26 Oct 2009
arxiv:0910.4934v1 [physics.med-ph] 26 Oct 2009 On the tradeoff between treatment time and plan quality in rotational arc radiation delivery 1. Introduction David Craft and Thomas Bortfeld Department of
More informationAn experimental investigation on the effect of beam angle optimization on the reduction of beam numbers in IMRT of head and neck tumors
JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS, VOLUME 13, NUMBER 4, 2012 An experimental investigation on the effect of beam angle optimization on the reduction of beam numbers in IMRT of head and neck tumors
More information3 INTEGER LINEAR PROGRAMMING
3 INTEGER LINEAR PROGRAMMING PROBLEM DEFINITION Integer linear programming problem (ILP) of the decision variables x 1,..,x n : (ILP) subject to minimize c x j j n j= 1 a ij x j x j 0 x j integer n j=
More informationOPTIMIZATION METHODS IN INTENSITY MODULATED RADIATION THERAPY TREATMENT PLANNING
OPTIMIZATION METHODS IN INTENSITY MODULATED RADIATION THERAPY TREATMENT PLANNING By DIONNE M. ALEMAN A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT
More informationOPTIMIZATION MODELS FOR RADIATION THERAPY: TREATMENT PLANNING AND PATIENT SCHEDULING
OPTIMIZATION MODELS FOR RADIATION THERAPY: TREATMENT PLANNING AND PATIENT SCHEDULING By CHUNHUA MEN A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF
More informationONLY AVAILABLE IN ELECTRONIC FORM
MANAGEMENT SCIENCE doi 10.1287/mnsc.1070.0812ec pp. ec1 ec7 e-companion ONLY AVAILABLE IN ELECTRONIC FORM informs 2008 INFORMS Electronic Companion Customized Bundle Pricing for Information Goods: A Nonlinear
More informationA Fully-Automated Intensity-Modulated Radiation Therapy Planning System
A Fully-Automated Intensity-Modulated Radiation Therapy Planning System Shabbir Ahmed, Ozan Gozbasi, Martin Savelsbergh Georgia Institute of Technology Ian Crocker, Tim Fox, Eduard Schreibmann Emory University
More informationAn Application of the Extended Cutting Angle Method in Radiation Therapy
An Application of the Extended Cutting Angle Method in Radiation Therapy by Valentin Koch A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF Bachelor of Science Honours in The
More informationUsing a research real-time control interface to go beyond dynamic MLC tracking
in partnership with Using a research real-time control interface to go beyond dynamic MLC tracking Dr. Simeon Nill Joint Department of Physics at The Institute of Cancer Research and the Royal Marsden
More information15.082J and 6.855J. Lagrangian Relaxation 2 Algorithms Application to LPs
15.082J and 6.855J Lagrangian Relaxation 2 Algorithms Application to LPs 1 The Constrained Shortest Path Problem (1,10) 2 (1,1) 4 (2,3) (1,7) 1 (10,3) (1,2) (10,1) (5,7) 3 (12,3) 5 (2,2) 6 Find the shortest
More informationVolumetric Modulated Arc Therapy - Clinical Implementation. Outline. Acknowledgement. History of VMAT. IMAT Basics of IMAT
Volumetric Modulated Arc Therapy - Clinical Implementation Daliang Cao, PhD, DABR Swedish Cancer Institute, Seattle, WA Acknowledgement David M. Shepard, Ph.D. Muhammad K. N. Afghan, Ph.D. Fan Chen, Ph.D.
More informationURGENT IMPORTANT FIELD SAFETY NOTIFICATION
Subject: Product: Incorrect Rotation of the Collimator or Couch ERGO Scope: Sites affected will be those: 1. Running ERGO version 1.7.3 and higher and, 2. Using a Multileaf Collimator (MLC) device for
More informationOperations Research Methods for Optimization in Radiation Oncology
Journal of Radiation Oncology Informatics Operations Research Methods for Optimization in Radiation Oncology Review Article Matthias Ehrgott 1, Allen Holder 2 1 Department of Management Science, Lancaster
More informationHugues Mailleux Medical Physics Department Institut Paoli-Calmettes Marseille France. Sunday 17 July 2016
Hugues Mailleux Medical Physics Department Institut Paoli-Calmettes Marseille France Sunday 17 July 2016 AGENDA 1. Introduction 2. Material 3. Optimization process 4. Results 5. Comments 6. Conclusion
More informationMixed-Integer Programming Techniques for Decomposing IMRT Fluence Maps Using Rectangular Apertures
Mixed-Integer Programming Techniques for Decomposing IMRT Fluence Maps Using Rectangular Apertures Z. Caner Taşkın J. Cole Smith H. Edwin Romeijn March 31, 2008 Abstract We consider a matrix decomposition
More informationOutline. Column Generation: Cutting Stock A very applied method. Introduction to Column Generation. Given an LP problem
Column Generation: Cutting Stock A very applied method thst@man.dtu.dk Outline History The Simplex algorithm (re-visited) Column Generation as an extension of the Simplex algorithm A simple example! DTU-Management
More informationColumn Generation: Cutting Stock
Column Generation: Cutting Stock A very applied method thst@man.dtu.dk DTU-Management Technical University of Denmark 1 Outline History The Simplex algorithm (re-visited) Column Generation as an extension
More informationFundamentals of Integer Programming
Fundamentals of Integer Programming Di Yuan Department of Information Technology, Uppsala University January 2018 Outline Definition of integer programming Formulating some classical problems with integer
More informationHigh Throughput Computing and Sampling Issues for Optimization in Radiotherapy
High Throughput Computing and Sampling Issues for Optimization in Radiotherapy Michael C. Ferris, University of Wisconsin Alexander Meeraus, GAMS Development Corporation Optimization Days, Montreal, May
More informationAnalysis of RapidArc optimization strategies using objective function values and dose-volume histograms
JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS, VOLUME 11, NUMBER 1, WINTER 2010 Analysis of RapidArc optimization strategies using objective function values and dose-volume histograms Mike Oliver, a Isabelle
More informationA Fully-Automated Intensity-Modulated Radiation Therapy Planning System
A Fully-Automated Intensity-Modulated Radiation Therapy Planning System Shabbir Ahmed, Ozan Gozbasi, Martin Savelsbergh Georgia Institute of Technology Tim Fox, Ian Crocker, Eduard Schreibmann Emory University
More informationIMRT and VMAT Patient Specific QA Using 2D and 3D Detector Arrays
IMRT and VMAT Patient Specific QA Using 2D and 3D Detector Arrays Sotiri Stathakis Outline Why IMRT/VMAT QA AAPM TG218 UPDATE Tolerance Limits and Methodologies for IMRT Verification QA Common sources
More informationInstituto de Engenharia de Sistemas e Computadores de Coimbra Institute of Systems Engineering and Computers INESC - Coimbra
Instituto de Engenharia de Sistemas e Computadores de Coimbra Institute of Systems Engineering and Computers INESC - Coimbra Humberto Rocha Joana Matos Dias On the Optimization of Radiation Therapy Planning
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More informationGENERAL ASSIGNMENT PROBLEM via Branch and Price JOHN AND LEI
GENERAL ASSIGNMENT PROBLEM via Branch and Price JOHN AND LEI Outline Review the column generation in Generalized Assignment Problem (GAP) GAP Examples in Branch and Price 2 Assignment Problem The assignment
More informationIMSURE QA SOFTWARE FAST, PRECISE QA SOFTWARE
QA SOFTWARE FAST, PRECISE Software for accurate and independent verification of monitor units, dose, and overall validity of standard, IMRT, VMAT, SRS and brachytherapy plans no film, no phantoms, no linac
More informationTHE WIRELESS PHANTOM PERFORM ACCURATE PATIENT QA IN LESS TIME THAN EVER!
THE WIRELESS PHANTOM PERFORM ACCURATE PATIENT QA IN LESS TIME THAN EVER! Confidence in complex treatments Modern radiation therapy uses complex plans with techniques such as IMRT, VMAT and Tomotherapy.
More informationAlgorithms for Integer Programming
Algorithms for Integer Programming Laura Galli November 9, 2016 Unlike linear programming problems, integer programming problems are very difficult to solve. In fact, no efficient general algorithm is
More informationTomotherapy Physics. Machine Twinning and Quality Assurance. Emilie Soisson, MS
Tomotherapy Physics Machine Twinning and Quality Assurance Emilie Soisson, MS Tomotherapy at UW- Madison Treating for nearly 5 years Up to ~45 patients a day on 2 tomo units Units twinned to facilitate
More informationA column generation approach for evaluating delivery efficiencies of collimator technologies in IMRT treatment planning
A column generation approach for evaluating delivery efficiencies of collimator technologies in IMRT treatment planning M Gören and Z C Taşkın Department of Industrial Engineering, Boğaziçi University,
More informationGate Sizing by Lagrangian Relaxation Revisited
Gate Sizing by Lagrangian Relaxation Revisited Jia Wang, Debasish Das, and Hai Zhou Electrical Engineering and Computer Science Northwestern University Evanston, Illinois, United States October 17, 2007
More information8/4/2016. Emerging Linac based SRS/SBRT Technologies with Modulated Arc Delivery. Disclosure. Introduction: Treatment delivery techniques
Emerging Linac based SRS/SBRT Technologies with Modulated Arc Delivery Lei Ren, Ph.D. Duke University Medical Center 2016 AAPM 58 th annual meeting, Educational Course, Therapy Track Disclosure I have
More informationSolving lexicographic multiobjective MIPs with Branch-Cut-Price
Solving lexicographic multiobjective MIPs with Branch-Cut-Price Marta Eso (The Hotchkiss School) Laszlo Ladanyi (IBM T.J. Watson Research Center) David Jensen (IBM T.J. Watson Research Center) McMaster
More informationarxiv: v3 [physics.med-ph] 6 Dec 2011
Multicriteria VMAT optimization David Craft, Dualta McQuaid, Jeremiah Wala, Wei Chen, Ehsan Salari, Thomas Bortfeld December 5, 2011 arxiv:1105.4109v3 [physics.med-ph] 6 Dec 2011 Abstract Purpose: To make
More informationAn Introduction to Dual Ascent Heuristics
An Introduction to Dual Ascent Heuristics Introduction A substantial proportion of Combinatorial Optimisation Problems (COPs) are essentially pure or mixed integer linear programming. COPs are in general
More informationLINEAR PROGRAMMING FORMULATIONS AND ALGORITHMS FOR RADIOTHERAPY TREATMENT PLANNING
LINEAR PROGRAMMING FORMULATIONS AND ALGORITHMS FOR RADIOTHERAPY TREATMENT PLANNING ARINBJÖRN ÓLAFSSON AND STEPHEN J. WRIGHT Abstract. Optimization has become an important tool in treatment planning for
More informationFAST, precise. qa software
qa software FAST, precise Software for accurate and independent verification of monitor units, dose, and overall validity of standard, IMRT, rotational or brachytherapy plans no film, no phantoms, no linac
More informationOptimization with Multiple Objectives
Optimization with Multiple Objectives Eva K. Lee, Ph.D. eva.lee@isye.gatech.edu Industrial & Systems Engineering, Georgia Institute of Technology Computational Research & Informatics, Radiation Oncology,
More informationAn Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy
An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy Chenyang Xu 1, Siemens Corporate Research, Inc., Princeton, NJ, USA Xiaolei Huang,
More informationNetwork Flows. 7. Multicommodity Flows Problems. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 7. Multicommodity Flows Problems 7.2 Lagrangian Relaxation Approach Fall 2010 Instructor: Dr. Masoud Yaghini The multicommodity flow problem formulation: We associate nonnegative
More informationTEPZZ Z754_7A_T EP A1 (19) (11) EP A1 (12) EUROPEAN PATENT APPLICATION. (51) Int Cl.: A61N 5/10 ( )
(19) TEPZZ Z74_7A_T (11) EP 3 07 417 A1 (12) EUROPEAN PATENT APPLICATION (43) Date of publication: 0..16 Bulletin 16/ (1) Int Cl.: A61N / (06.01) (21) Application number: 16163147.8 (22) Date of filing:
More informationMonaco Concepts and IMRT / VMAT Planning LTAMON0003 / 3.0
and IMRT / VMAT Planning LTAMON0003 / 3.0 and Planning Objectives By the end of this presentation you can: Describe the cost functions in Monaco and recognize their application in building a successful
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
More informationLagrangean Methods bounding through penalty adjustment
Lagrangean Methods bounding through penalty adjustment thst@man.dtu.dk DTU-Management Technical University of Denmark 1 Outline Brief introduction How to perform Lagrangean relaxation Subgradient techniques
More informationCrew Scheduling Problem: A Column Generation Approach Improved by a Genetic Algorithm. Santos and Mateus (2007)
In the name of God Crew Scheduling Problem: A Column Generation Approach Improved by a Genetic Algorithm Spring 2009 Instructor: Dr. Masoud Yaghini Outlines Problem Definition Modeling As A Set Partitioning
More informationSUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING
Bulletin of Mathematics Vol. 06, No. 0 (20), pp.. SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING Eddy Roflin, Sisca Octarina, Putra B. J Bangun,
More informationMVE165/MMG631 Linear and integer optimization with applications Lecture 9 Discrete optimization: theory and algorithms
MVE165/MMG631 Linear and integer optimization with applications Lecture 9 Discrete optimization: theory and algorithms Ann-Brith Strömberg 2018 04 24 Lecture 9 Linear and integer optimization with applications
More information2. Modeling AEA 2018/2019. Based on Algorithm Engineering: Bridging the Gap Between Algorithm Theory and Practice - ch. 2
2. Modeling AEA 2018/2019 Based on Algorithm Engineering: Bridging the Gap Between Algorithm Theory and Practice - ch. 2 Content Introduction Modeling phases Modeling Frameworks Graph Based Models Mixed
More informationDose Calculation and Optimization Algorithms: A Clinical Perspective
Dose Calculation and Optimization Algorithms: A Clinical Perspective Daryl P. Nazareth, PhD Roswell Park Cancer Institute, Buffalo, NY T. Rock Mackie, PhD University of Wisconsin-Madison David Shepard,
More information56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997
56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997 Answer #1 and any five of the remaining six problems! possible score 1. Multiple Choice 25 2. Traveling Salesman Problem 15 3.
More informationDose-volume-based IMRT fluence optimization: A fast least-squares approach with differentiability
Dose-volume-based IMRT fluence optimization: A fast least-squares approach with differentiability Yin Zhang and Michael Merritt Technical Report TR06-11 Department of Computational and Applied Mathematics
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationMathematical Optimization in Radiotherapy Treatment Planning
1 / 35 Mathematical Optimization in Radiotherapy Treatment Planning Ehsan Salari Department of Radiation Oncology Massachusetts General Hospital and Harvard Medical School HST S14 May 13, 2013 2 / 35 Outline
More informationAdvanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras
Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras Lecture - 35 Quadratic Programming In this lecture, we continue our discussion on
More informationSimulation. Lecture O1 Optimization: Linear Programming. Saeed Bastani April 2016
Simulation Lecture O Optimization: Linear Programming Saeed Bastani April 06 Outline of the course Linear Programming ( lecture) Integer Programming ( lecture) Heuristics and Metaheursitics (3 lectures)
More informationWireless frequency auctions: Mixed Integer Programs and Dantzig-Wolfe decomposition
Wireless frequency auctions: Mixed Integer Programs and Dantzig-Wolfe decomposition Laszlo Ladanyi (IBM T.J. Watson Research Center) joint work with Marta Eso (The Hotchkiss School) David Jensen (IBM T.J.
More informationInstituto de Engenharia de Sistemas e Computadores de Coimbra Institute of Systems Engineering and Computers INESC - Coimbra
Instituto de Engenharia de Sistemas e Computadores de Coimbra Institute of Systems Engineering and Computers INESC - Coimbra Humberto Rocha Brígida da Costa Ferreira Joana Matos Dias Maria do Carmo Lopes
More informationPreviously Local sensitivity analysis: having found an optimal basis to a problem in standard form,
Recap, and outline of Lecture 20 Previously Local sensitivity analysis: having found an optimal basis to a problem in standard form, if the cost vectors is changed, or if the right-hand side vector is
More information15.083J Integer Programming and Combinatorial Optimization Fall Enumerative Methods
5.8J Integer Programming and Combinatorial Optimization Fall 9 A knapsack problem Enumerative Methods Let s focus on maximization integer linear programs with only binary variables For example: a knapsack
More informationAn Optimisation Model for Intensity Modulated Radiation Therapy
An Optimisation Model for Intensity Modulated Radiation Therapy Matthias Ehrgott Department of Engineering Science University of Auckland New Zealand m.ehrgott@auckland.ac.nz Abstract Intensity modulated
More informationarxiv: v1 [physics.med-ph] 11 Jan 2019
Dynamic fluence map sequencing using piecewise linear leaf position functions Matthew Kelly 1, Jacobus H.M. van Amerongen 2, Marleen Balvert 2,3,4, David Craft 5 1 Department of Mechanical Engineering,
More informationCurrent state of multi-criteria treatment planning
Current state of multi-criteria treatment planning June 11, 2010 Fall River NE AAPM meeting David Craft Massachusetts General Hospital Talk outline - Overview of optimization - Multi-criteria optimization
More informationMVE165/MMG630, Applied Optimization Lecture 8 Integer linear programming algorithms. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming algorithms Ann-Brith Strömberg 2009 04 15 Methods for ILP: Overview (Ch. 14.1) Enumeration Implicit enumeration: Branch and bound Relaxations Decomposition methods:
More informationSolutions for Operations Research Final Exam
Solutions for Operations Research Final Exam. (a) The buffer stock is B = i a i = a + a + a + a + a + a 6 + a 7 = + + + + + + =. And the transportation tableau corresponding to the transshipment problem
More informationThe MIP-Solving-Framework SCIP
The MIP-Solving-Framework SCIP Timo Berthold Zuse Institut Berlin DFG Research Center MATHEON Mathematics for key technologies Berlin, 23.05.2007 What Is A MIP? Definition MIP The optimization problem
More informationSurrogate Gradient Algorithm for Lagrangian Relaxation 1,2
Surrogate Gradient Algorithm for Lagrangian Relaxation 1,2 X. Zhao 3, P. B. Luh 4, and J. Wang 5 Communicated by W.B. Gong and D. D. Yao 1 This paper is dedicated to Professor Yu-Chi Ho for his 65th birthday.
More informationCS 473: Algorithms. Ruta Mehta. Spring University of Illinois, Urbana-Champaign. Ruta (UIUC) CS473 1 Spring / 36
CS 473: Algorithms Ruta Mehta University of Illinois, Urbana-Champaign Spring 2018 Ruta (UIUC) CS473 1 Spring 2018 1 / 36 CS 473: Algorithms, Spring 2018 LP Duality Lecture 20 April 3, 2018 Some of the
More informationDiscrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity
Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity Marc Uetz University of Twente m.uetz@utwente.nl Lecture 5: sheet 1 / 26 Marc Uetz Discrete Optimization Outline 1 Min-Cost Flows
More informationOptimization Techniques for Design Space Exploration
0-0-7 Optimization Techniques for Design Space Exploration Zebo Peng Embedded Systems Laboratory (ESLAB) Linköping University Outline Optimization problems in ERT system design Heuristic techniques Simulated
More informationADVANCING CANCER TREATMENT
The RayPlan treatment planning system makes proven, innovative RayStation technology accessible to clinics that need a cost-effective and streamlined solution. Fast, efficient and straightforward to use,
More informationLast topic: Summary; Heuristics and Approximation Algorithms Topics we studied so far:
Last topic: Summary; Heuristics and Approximation Algorithms Topics we studied so far: I Strength of formulations; improving formulations by adding valid inequalities I Relaxations and dual problems; obtaining
More informationIntegrating column generation in a method to compute a discrete representation of the non-dominated set of multi-objective linear programmes
4OR-Q J Oper Res (2017) 15:331 357 DOI 10.1007/s10288-016-0336-9 RESEARCH PAPER Integrating column generation in a method to compute a discrete representation of the non-dominated set of multi-objective
More informationA List Heuristic for Vertex Cover
A List Heuristic for Vertex Cover Happy Birthday Vasek! David Avis McGill University Tomokazu Imamura Kyoto University Operations Research Letters (to appear) Online: http://cgm.cs.mcgill.ca/ avis revised:
More informationDose Calculations: Where and How to Calculate Dose. Allen Holder Trinity University.
Dose Calculations: Where and How to Calculate Dose Trinity University www.trinity.edu/aholder R. Acosta, W. Brick, A. Hanna, D. Lara, G. McQuilen, D. Nevin, P. Uhlig and B. Slater Dose Calculations - Why
More informationAdvanced Operations Research Techniques IE316. Quiz 2 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 2 Review Dr. Ted Ralphs IE316 Quiz 2 Review 1 Reading for The Quiz Material covered in detail in lecture Bertsimas 4.1-4.5, 4.8, 5.1-5.5, 6.1-6.3 Material
More informationCOLUMN GENERATION IN LINEAR PROGRAMMING
COLUMN GENERATION IN LINEAR PROGRAMMING EXAMPLE: THE CUTTING STOCK PROBLEM A certain material (e.g. lumber) is stocked in lengths of 9, 4, and 6 feet, with respective costs of $5, $9, and $. An order for
More informationA Generic Separation Algorithm and Its Application to the Vehicle Routing Problem
A Generic Separation Algorithm and Its Application to the Vehicle Routing Problem Presented by: Ted Ralphs Joint work with: Leo Kopman Les Trotter Bill Pulleyblank 1 Outline of Talk Introduction Description
More informationThe MSKCC Approach to IMRT. Outline
The MSKCC Approach to IMRT Spiridon V. Spirou, PhD Department of Medical Physics Memorial Sloan-Kettering Cancer Center New York, NY Outline Optimization Field splitting Delivery Independent verification
More informationUnit.9 Integer Programming
Unit.9 Integer Programming Xiaoxi Li EMS & IAS, Wuhan University Dec. 22-29, 2016 (revised) Operations Research (Li, X.) Unit.9 Integer Programming Dec. 22-29, 2016 (revised) 1 / 58 Organization of this
More informationSUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING
ASIAN JOURNAL OF MATHEMATICS AND APPLICATIONS Volume 2014, Article ID ama0156, 11 pages ISSN 2307-7743 http://scienceasia.asia SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING
More informationParallel Branch & Bound
Parallel Branch & Bound Bernard Gendron Université de Montréal gendron@iro.umontreal.ca Outline Mixed integer programming (MIP) and branch & bound (B&B) Linear programming (LP) based B&B Relaxation and
More informationInteger Programming ISE 418. Lecture 7. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 7 Dr. Ted Ralphs ISE 418 Lecture 7 1 Reading for This Lecture Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Wolsey Chapter 7 CCZ Chapter 1 Constraint
More informationLP-Modelling. dr.ir. C.A.J. Hurkens Technische Universiteit Eindhoven. January 30, 2008
LP-Modelling dr.ir. C.A.J. Hurkens Technische Universiteit Eindhoven January 30, 2008 1 Linear and Integer Programming After a brief check with the backgrounds of the participants it seems that the following
More informationInteger Programming Chapter 9
1 Integer Programming Chapter 9 University of Chicago Booth School of Business Kipp Martin October 30, 2017 2 Outline Branch and Bound Theory Branch and Bound Linear Programming Node Selection Strategies
More informationDepartment of Mathematics Oleg Burdakov of 30 October Consider the following linear programming problem (LP):
Linköping University Optimization TAOP3(0) Department of Mathematics Examination Oleg Burdakov of 30 October 03 Assignment Consider the following linear programming problem (LP): max z = x + x s.t. x x
More information