In this paper, the method to register multiple sets of skull CT images to absolute coordinate system is proposed. Contrary to correspondence paired mapping of previous techniques, four anatomical landmark points, three coplanar points and one non-coplanar point, compose three principal axes simple and unique for efficient registration by means of rigid body transformation. Throughout the numerical simulation with added random noises, the error performances in terms of different rotation and rounding-off of landmark points, and incorrect localization of anatomical landmark and target points are quantitatively analyzed to generalize the proposed technique. Experiments using real skull CT images demonstrate the feasibility for an efficient use in clinical practice.