The method can be applied to vector measurements with full variance-covariance matrix. The rotation matrix is estimated in two steps both of which are optimal in their own. The solution is analytical. Special case with only 2 non-collinear vector measurements is treated. Error analysis of the solution is provided.