A practical global distortion correction method for an image intensifier based x-ray fluoroscopy system

X-ray images acquired on systems with image intensifiers (II) exhibit characteristic distortion which is due to both external and internal factors. The distortion is dependent on the orientation of the II, a fact particularly relevant to II's mounted on C arms which have several degrees of freedom of motion. Previous descriptions of distortion correction strategies have relied on a dense sampling of the C-arm orientation space, and as such have been limited mostly to a single arc of the primary angle, alpha. We present a new method which smooths the trajectories of the segmented vertices of the grid phantom as a function of a prior to solving the two-dimensional warping problem. It also shows that the same residual errors of distortion correction could be achieved without fitting the trajectories of the grid vertices, but instead applying the previously described global method of distortion correction, followed by directly smoothing the values of the polynomial coefficients as functions of the C-arm orientation parameters. When this technique was applied to a series of test images at arbitrary alpha, the root-mean-square (RMS) residual error was 0.22 pixels. The new method was extended to three degrees of freedom of the C-arm motion: the primary angle, alpha; the secondary angle, beta; and the source-to-intensifier distance, lambda. Only 75 images were used to characterize the distortion for the following ranges: alpha, +/- 45 degrees (Deltaalpha = 22.5 degrees); beta, +/- 36 degrees (Deltabeta = 18 degrees); lambda, 98-118 cm (Deltalambda = 10 cm). When evaluated on a series of test images acquired at arbitrary (alpha, beta, lambda), the RMS residual error was 0.33 pixels. This method is targeted at applications such as guidance of catheter-based interventions and treatment planning for brachytherapy, which require distortion-corrected images over a large range of C-arm orientations.