缔合Legendre函数的快速插值计算及其在六边形网格模型重力异常计算中的应用

扫码查看

原文链接

NETL
NSTL
万方数据

中文摘要：鉴于局部区域非等纬度分布点的超高阶次球谐综合计算效率较低,本文深入研究了缔合Legendre函数的插值算法,结合模型重力异常求解的展开式,提出了超高阶模型空间重力异常插值快速计算方法,为了验证本文方法在东西狭长分布的点集更具有应用优势,采用球谐旋转变换技术进一步提升了计算效率.试验结果表明,利用2160阶次EGM2008模型计算日本地区同一高度的30 303个非等纬度分布的六边形网格点处模型重力异常,插值快速计算方法相比逐点计算,在误差不超±0.005 mGal水平下,计算耗时从3 669.41 s缩减到98.05 s,同时经过球谐旋转变换,将南北狭长的日本区域点集旋转为东西狭长分布,使得上述相应计算内容的耗时进一步由98.05 s缩减到19.06 s,计算效率相比最初方法提速比达到近200倍,有效解决了非等纬度分布点的超高阶模型重力异常快速解算的效率难题,且验证了该方法在东西狭长分布的情况下具有更高的提速比.

外文标题：A fast method for interpolation of the associated Legendre functions and its application to the calculation of local hexagonal grid point gravity anomalies using an ultra-high-degree gravity field model

外文摘要：In view of the low efficiency of the calculation of ultra-high-degree spherical harmonic synthesis for non-equal latitude distributed points of regional area,this article has conducted a deep research on the interpolation algorithm of associated Leg-endre functions.Combined with the harmonic expansion of the model gravity anomaly,a fast calculation method using interpo-lation was proposed for ultra-high-degree model free-air gravity anomaly.In order to verify the advantages of this new method over point sets with narrow distribution from east to west,the efficiency of the new method was further improved by adopting the spherical harmonic rotation(SHR)technology.The experimental results showed that the proposed method using interpola-tion reduced the calculation time from 3 669.41 s to 98.05 s compared with point-by-point approaches,with errors not excee-ding±0.005 mGal,when we used EGM2008 model up to 2160 degree and order to achieve the model free-air gravity anomalies of 30 303 hexagonal grid points at the same height level in Japan with non-equal latitude distribution.Meanwhile,by applying SHR and rotating these points distributed in the north-south direction to the east-west direction,the time consumption for the corresponding solution was further reduced from 98.05 s to 19.06 s,with a speed-up ratio of nearly 200 times to the original method.The method proposed in this article effectively solves the problem of low efficiency in solving ultra-high-degree model free-air gravity anomalies for regional non-equal latitude distributed points,and has a higher computing speed-up effect in the case of east-west elongated distribution.

外文关键词：

Earth gravitational modelspherical harmonic synthesisLegendre functionhexagonal gridspherical harmonic transformation under rotationinterpolation method

作者：

李新星、范昊鹏、万宏发、范雕、冯进凯

展开 >

作者单位：

信息工程大学地理空间信息学院,河南郑州 450001

地理信息工程国家重点实验室,陕西西安 710054

智慧地球重点实验室,北京 100020

关键词：

地球重力场模型球谐综合勒让德函数六边形网格球谐旋转变换插值方法

基金：

国家自然科学基金国家自然科学基金国家自然科学基金国家自然科学基金地理信息工程国家重点实验室基金中国地质大学(武汉)中央高校基本科研业务费

项目编号：

U23A2028414040204217400742074013SKLGIE2023-Z-1-12024XLB13

出版年：

2024

DOI：

10.11947/j.AGCS.2024.20230110

测绘学报

中国测绘学会

测绘学报

CSTPCD北大核心

影响因子：1.602

ISSN：1001-1595

年,卷(期)：2024.53(9)

参考文献量10