An algorithm has been developed to correct the image distortion introduced by the image intensifier - CCD camera chain in our volume tomographic angiography imaging system. First the distorted image is partitioned into small quadrilateral regions, then each partition is mapped individually to the corrected image by its corresponding bilinear functions. The parameters of these bilinear functions can be acquired by a special calibration procedure and stored as a lookup table. The maximum errors before and after the distortion correction were measured to be 5.96 mm and 0.03 mm. The method is further applied to the 3D image reconstruction of a simulated vascular phantom and live animal studies.