The method proposed is based on the principal axes and the central moments theories, which enable the geometrical registration of two volumes which are differently oriented in the space, with respect to the fixed Cartesian global coordinate system  whose versors are (x, y, z).

The proposed approach is an iterative method which avoid rotations in the 3D space, but performs three sequential rotations in the coronal plane, in the axial plane, and finally in the sagittal plane. At each step the central inertia matrix of the rigid body is computed from where the rotation angle is estimated. After each rotation, the eigenvectors of the volume are changed then they have to be recomputed by means of a new inertia matrix.  The results are more accurate since the interpolation errors are reduced by restricting the rotation side effects on a single plane at a time.

This method is automatic and it does not require any input of parameters.