Shape Rectangularization Problems in Intensity-Modulated Radiation Therapy

  • Authors:
  • Nikhil Bansal;Danny Z. Chen;Don Coppersmith;Xiaobo S. Hu;Shuang Luan;Ewa Misiołek;Baruch Schieber;Chao Wang

  • Affiliations:
  • IBM T.J. Watson Research Center, P.O. Box 218, 10598, Yorktown Heights, NY, USA;University of Notre Dame, Department of Computer Science and Engineering, 46556, Notre Dame, IN, USA;IDA Center for Communications Research, 805 Bunn Drive, 08540, Princeton, NJ, USA;University of Notre Dame, Department of Computer Science and Engineering, 46556, Notre Dame, IN, USA;University of New Mexico, Department of Computer Science, 87131-0001, Albuquerque, NM, USA;Saint Mary’s College, Mathematics Department, 46556, Notre Dame, IN, USA;IBM T.J. Watson Research Center, P.O. Box 218, 10598, Yorktown Heights, NY, USA;University of Notre Dame, Department of Computer Science and Engineering, 46556, Notre Dame, IN, USA

  • Venue:
  • Algorithmica
  • Year:
  • 2011

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper, we present a theoretical study of several shape approximation problems, called shape rectangularization (SR), which arise in intensity-modulated radiation therapy (IMRT). Given a piecewise linear function f such that f(x)≥0 for any x∈ℝ, the SR problems seek an optimal set of constant window functions to approximate f under a certain error criterion, such that the sum of the resulting constant window functions equals (or well approximates) f. A constant window function W(⋅) is defined on an interval I such that W(x) is a fixed value h0 for any x∈I and is 0 otherwise. A constant window function can be viewed as a rectangle (or a block) geometrically, or as a vector with the consecutive a’s property combinatorially. The SR problems find applications in setup time and beam-on time minimization and dose simplification of the IMRT treatment planning process. We show that the SR problems are APX-Hard, and thus we aim to develop theoretically efficient and provably good quality approximation SR algorithms. Our main contribution is to present algorithms for a key SR problem that achieve approximation ratios better than 2. For the general case, we give a $\frac{24}{13}$ -approximation algorithm. For unimodal input curves, we give a $\frac{9}{7}$ -approximation algorithm. We also consider other variants for which better approximation ratios are possible. We show that an important SR case that has been studied in medical literature can be formulated as a k-MST(k-minimum-spanning-tree) problem on a certain geometric graph G; based on a set of geometric observations and a non-trivial dynamic programming scheme, we are able to compute an optimal k-MST in G efficiently.