Randomized algorithms for mixed matching and covering in hypergraphs in 3D seed reconstruction in brachytherapy

Brachytherapy is a radiotherapy method for cancer. In its low dose radiation (LDR) variant a number of radioactive implants, so-called seeds, are inserted into the affected organ through an operation. After the implantation, it is essential to determine the locations of the seeds in the organ. A common method is to take three X-ray photographs from different angles; the seeds show up on the X-ray photos as small white lines. In order to reconstruct the three-dimensional configuration from these X-ray photos, one has to determine which of these white lines belong to the same seed. We model the problem as a mixed packing and covering hypergraph optimization problem and present a randomized approximation algorithm based on linear programming. We analyse the worst-case performance of the algorithm by discrete probabilistic methods and present results for data of patients with prostate cancer from the university clinic of Schleswig-Holstein, Campus Kiel. These examples show an almost optimal performance of the algorithm which presently cannot be matched by the theoretical analysis.

[1]  M. Okamoto.  Some inequalities relating to the partial sum of binomial probabilities , 1959 .

[2]  P J Biggs,et al.  Geometric reconstruction of seed implants using a three-film technique. , 1983, Medical physics.

[3]  M D Altschuler,et al.  Rapid, accurate, three-dimensional location of multiple sees in implant radiotherapy treatment planning. , 1983, Physics in medicine and biology.

[4]  W. T. Gowers,et al.  RANDOM GRAPHS (Wiley Interscience Series in Discrete Mathematics and Optimization) , 2001 .

[5]  D Tubic,et al.  Automated seed detection and three-dimensional reconstruction. II. Reconstruction of permanent prostate implants using simulated annealing. , 2001, Medical physics.

[6]  Egon Balas,et al.  An Algorithm for the Three-Index Assignment Problem , 1991, Oper. Res..

[7]  Prabhakar Raghavan,et al.  Randomized rounding: A technique for provably good algorithms and algorithmic proofs , 1985, Comb..

[8]  W. Hoeffding.  Probability Inequalities for sums of Bounded Random Variables , 1963 .

[9]  Anand Srivastav,et al.  Derandomization in Combinatorial Optimization , 2001 .

[10]  R. Lougee-Heimer,et al.  The Common Optimization INterface for Operations Research: Promoting open-source software in the operations research community , 2003 .

[11]  M S Rosenthal,et al.  An automatic seed identification technique for interstitial implants using three isocentric radiographs. , 1983, Medical physics.

[13]  R L Siddon,et al.  Two-film brachytherapy reconstruction algorithm. , 1985, Medical physics.

[14]  Leslie G. Valiant,et al.  Fast probabilistic algorithms for hamiltonian circuits and matchings , 1977, STOC '77.

[15]  H. Chernoff.  A Measure of Asymptotic Efficiency for Tests of a Hypothesis Based on the sum of Observations , 1952 .

[16]  R. Marks,et al.  Three-dimensional seed reconstruction from an incomplete data set for prostate brachytherapy. , 2004, Physics in medicine and biology.

[17]  H I Amols,et al.  A three-film technique for reconstruction of radioactive seed implants. , 1981, Medical physics.

[18]  W. Brogan.  Algorithm for Ranked Assignments with Applications to Multiobject Tracking , 1989 .

[19]  M Zaider,et al.  Operator-free, film-based 3D seed reconstruction in brachytherapy. , 2002, Physics in medicine and biology.

[20]  Anand Srivastav,et al.  Three-dimensional reconstruction of seed implants by randomized rounding and visual evaluation. , 2007, Medical physics.

[21]  Noga Alon,et al.  The Probabilistic Method , 2015, Fundamentals of Ramsey Theory.

[22]  D. Ash,et al.  ESTRO/EAU/EORTC recommendations on permanent seed implantation for localized prostate cancer. , 2000, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[23]  Frank-André Siebert,et al.  The design and testing of a solid phantom for the verification of a commercial 3D seed reconstruction algorithm. , 2005, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[24]  R. J. Gallagher,et al.  Treatment planning for brachytherapy: an integer programming model, two computational approaches and experiments with permanent prostate implant planning. , 1999, Physics in medicine and biology.