The method for celestial bodies’ center of mass position relative to their figures determination on the basis of harmonic analysis of the expansion in spherical functions in order to refine the physical libration parameters