This invention discloses a method for calculating the sub-
satellite point based on the SGP4 orbital model, relating to the field of
satellite navigation and
orbit prediction technology. The polar spherical
projector of this invention eliminates
longitude jumps in
latitude >75°, ensuring continuous trajectory output for scenarios such as polar
scientific expeditions and
glacier monitoring. By directly mapping three-dimensional radial information to generate
latitude, it avoids the iterative
convergence problem of the ellipsoidal model, ensuring natural symmetry of the North and South Poles coordinates and continuous
longitude domain. Second-order differential compensation of the mean anomaly angle, combined with dynamic correction of the Earth's rotational
angular velocity, strictly aligns the UTC
timestamp with the orbital position. It suppresses the cumulative error generated by traditional linear interpolation at perigee. It enhances adaptability to complex perturbation environments;
atmospheric drag vectorization compensation captures real-time aerodynamic direction through
relative velocity cross terms, combined with continuous
processing of the density index model, significantly improving the orbital attenuation prediction accuracy of low-Earth
orbit satellites during geomagnetic storms.