A three-dimensional image registration method also uses an algorithm and applies it to the solution of such problems in 3D CT DSA. The method can deal with incomplete or partial volumetric data in registration and correct patient global motion prior to subtraction even when it is coupled with local unconscious/nonrigid movements. Experimental results demonstrate the effectiveness of this algorithm on several clinical spiral CT data.