Three-dimensional seed reconstruction from an incomplete data set for prostate brachytherapy

Intra-operative dosimetry in prostate brachytherapy requires 3D coordinates of the implanted, radioactive seeds. Since CT is not readily available during the implant operation, projection x-rays are commonly used for intra-operative seed localization. Three x-ray projections are usually used. The requirement of the current seed reconstruction algorithms is that the seeds must be identified on all three projections. However, in practice this is often difficult to accomplish due to the problem of heavily clustered and overlapping seeds. We have developed an algorithm that permits seed reconstruction from an incomplete data set. Instead of all three projections, the new algorithm requires only one of the three projections to be complete. Furthermore, even if all three projections are incomplete, it can reconstruct 100% of the implanted seeds depending on how the undetected seeds are distributed among the projections. The method utilizes the principles of epipolar imaging geometry and pseudo-matching of the undetected seeds. The algorithm was successfully applied to a large number of clinical cases where seeds imperceptibly overlap in some projections.