Automatic correction of biplane projection imaging geometry

A novel method is presented for correcting errors in measurements of biplane projection imaging geometry without prior identification of corresponding points in the two images. For imaged objects that project onto both images, a constraint equation is obtained that relates weighted integrals along corresponding epipolar lines. The integrals are computed to first order in the angular beamwidth, which is assumed to be small. Starting from measured or estimated values, geometrical parameters are computed iteratively in order to maximize the correlation between epipolar line integrals in the two images. Improvement in the computation of corresponding epipolar lines is demonstrated on images of a wire phantom. The root mean square distance of the epipolar lines from the corresponding reference points is improved from 15 pixel widths to less than 4 pixel widths (1.3 mm). Convergence is demonstrated on phantom images for individual parameter variations up to 70% in relative magnification, a relative shift of the imaging planes by 50 pixels, or a relative rotation of at least 35 degrees around either of two axes. Applicability to clinical images is demonstrated by using a biplane angiogram of a pig to align corresponding points determined from images of a Perspex cube acquired with the same geometry. (C) 1996 American Association of Physicists in Medicine.