In order to satisfy the requirement from the polarization navigation for the spatial position of feature points in the sky,a method is introduced to calculate the solar spatial position with the Rayleigh atmosphere polarization pattern.First,the distribution pattern of the polarized sky light is established on the basis of the Rayleigh scattering theory of the atmosphere optics.Then a method using nonlinear least square optimization is introduced to solve the matter of finding the solar spatial position with finite sampling information from the atmosphere polarization pattern.The Gauss-Newton method couple with global search for initial value guarantee less steady residual error.Experimental simulation results indicate that the error of the solar zenith angle and azimuth angle is respectively less than 10-5 degree and less than 10-6 degree with sampling the degree of polarization,and the errors of the two angles are all less than 10-6 with sampling the angle of the degree.And the method has shown good precision under different sampling types.As a result,the method is proved to be precise and adapt various sampling types,the polarization pattern can be efficiently processed into the solar spatial position.