Estimating the 3D skeletons and transverse areas of coronary arteries from biplane angiograms

A method for estimating the three-dimensional (3D) skeletons and transverse areas of the lumens of coronary arteries from digital X-ray angiograms is described. The method is based on the use of a 3D generalized cylinder (GC) consisting of a series of 3D elliptical disks transverse to and centered on a 3D skeleton (medial axis) of the coronary arteries. The estimates of the transverse areas are based on a nonlinear least-squares-error estimation technique described by D.W. Marquardt (1963). This method exploits densitometric profiles, boundary estimates, and the orientation of the arterial skeleton in 3-space and includes an automatic artery tracking procedure. It applies an adaptive window to the densitometric profile data that are used in the parameter estimation. Preliminary experimental tests of the procedure on angiograms of in vivo human coronaries and on synthetic images yield encouraging results.