Reconstructing the cross sections of coronary arteries from biplane angiograms

An algorithm that reconstructs the cross sections of the lumens of coronary arteries from two mutually orthogonal X-ray projections is described. The algorithm accommodates the possibility of elliptical, crescent, or star shapes. It represents each biplane projection of a transverse slice of the arterial lumen as a binary-valued image. The single-coordinate moments of these two projection images are equal to those of the slice. Since the cross-coordinate moments of the slice are not available from the projections, an algorithm to estimate these moments based on assumptions of smoothness and connectivity is developed. Once all the missing moments are estimated, the image of the slice can be estimated by inverting these moments, using the uniqueness theorem governing the relation between an image and its moments. Preliminary tests of the algorithm on synthetic data, on hardware phantoms and on a segment of a barium-enhanced in vitro coronary artery are reported.