Method for calculating gravity anomaly high-order radial derivative by using band-limited thought

The invention relates to a method for calculating a gravity anomaly high-order radial derivative by using a band-limited thought, which is mainly technically characterized by comprising the following steps of representing an analytic kernel function of an earth external gravity anomaly Poisson integral formula as a spherical harmonic series expansion formula, and obtaining a gravity anomaly high-order radial derivative spherical harmonic series expression by using a derivation method, and truncating the spherical harmonic series expression of the gravity anomaly high-order radial derivative into a band-limited summation formula consistent with the spectral range of the gravity observation value by using a removal-recovery technology, compensating a far-region truncation error of the gravity anomaly high-order radial derivative by using a global gravity potential model, and calculating to obtain the gravity anomaly high-order radial derivative. The method is reasonable in design, the gravity anomaly high-order