Matching and reconstruction of brachytherapy seeds using the Hungarian algorithm (MARSHAL)
- Department of Computer Science, Johns Hopkins University, Baltimore, Maryland 21218 (United States) and Department of Mechanical Engineering, State University of New York at Stony Brook, Stony Brook, New York 11794 (United States) and Department of Mechanical Engineering, Johns Hopkins University, Baltimore, Maryland 21218 (United States) and Acoustic MedSystems, Inc., 2110 Clearlake Boulevard, Champaign, Illinois 61822 (United States)
Intraoperative dosimetric quality assurance in prostate brachytherapy critically depends on discerning the three-dimensional (3D) locations of implanted seeds. The ability to reconstruct the implanted seeds intraoperatively will allow us to make immediate provisions for dosimetric deviations from the optimal implant plan. A method for seed reconstruction from segmented C-arm fluoroscopy images is proposed. The 3D coordinates of the implanted seeds can be calculated upon resolving the correspondence of seeds in multiple x-ray images. We formalize seed-matching as a combinatorial optimization problem, which has salient features: (a) extensively studied solutions by the computer science community; (b) proof for the nonexistence of any polynomial time exact algorithm; and (c) a practical pseudo-polynomial algorithm that mostly runs in O(N{sup 3}) time using any number of images. We prove that two images are insufficient to correctly match the seeds, while a third image renders the matching problem to be of nonpolynomial complexity. We utilize the special structure of the problem and propose a pseudopolynomial time algorithm. Using three presegmented images, matching and reconstruction of brachytherapy seeds using the Hungarian algorithm achieved complete matching in simulation experiments; and 98.5% in phantom experiments. 3D reconstruction error for correctly matched seeds has a mean of 0.63 mm, and 0.9 mm for incorrectly matched seeds. The maximum seed reconstruction error in each implant was typically around 1.32 mm. Both on synthetic data and in phantom experiments, matching rate and reconstruction error achieved using presegmented images was found to be sufficient for prostate brachytherapy. The algorithm is extendable to deal with arbitrary number of images without any loss in speed or accuracy. The algorithm is sufficiently generic to provide a practical solution to any correspondence problem, across different imaging modalities and features.
- OSTI ID:
- 20726912
- Journal Information:
- Medical Physics, Vol. 32, Issue 11; Other Information: DOI: 10.1118/1.2104087; (c) 2005 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA); ISSN 0094-2405
- Country of Publication:
- United States
- Language:
- English
Similar Records
BrachyView: Proof-of-principle of a novel in-body gamma camera for low dose-rate prostate brachytherapy
SU-D-BRF-07: Ultrasound and Fluoroscopy Based Intraoperative Image-Guidance System for Dynamic Dosimetry in Prostate Brachytherapy