【技术实现步骤摘要】
【专利摘要】本专利技术公开了,属于大地测量
。该方法对于全球海洋扰动重力计算,通过将观测数据补充至全球区域可在每一纬度圈进行FFT计算,而后在经度方向上进行求和运算,这种处理方式保证了FFT计算后的结果与原解析算法是一致的。对于局部海洋扰动重力计算,通过对核函数和观测数据补充一定量数据从而使得核函数构成了循环矩阵,进而利用FFT计算后保证了结果与原解析算法是一致的。该专利技术为海洋测高卫星数据计算扰动重力提供了计算公式,同时给出了精确快速方法,本专利技术也可以应用于物理大地测量中大部分积分问题的求解如求解大地水准面、垂线偏差、地形改正等。【专利说明】
本专利技术涉及精确快速计算全球及局部海洋扰动重力的方法,属于大地测量
。
技术介绍
传统的海洋重力场确定一般是由卫星测高数据计算获得海洋区域的重力异常,随着海洋学、空间大地测量学的发展,扰动重力较重力异常表现出更多优势。首先海洋区域由于海面地形的量级较小,扰动重力的确定可以避免重力异常计算过程中的归算问题,因此扰动重力的确定要比重力异常精确,其次,从应用角度看,相同精度情况下,利用扰动重力计算重力场扰动场元的准确性要比重力异常高。基于上述原因,越来越多的学者将扰动重力作为地球重力场的基础数据。Molodensky (1960)给出了利用垂线偏差数据计算重力异常的公式,C.Hwang(1998)利用不同的推导方法也给出了利用垂线偏差数据计算重力异常的公式,然而到目前为止,没有公开文献给出利用垂线偏差数据计算扰动重力的公式。目前,随着卫星测高技术发展,观测数据的数量和精度大幅提升,计 ...
【技术保护点】
全球海洋扰动重力的精确快速计算方法,包括如下步骤:步骤一:采用卫星测高技术获得海洋区域的大地水准面高数据H海,采用重力场模型获得陆地区域的大地水准面高H陆。式中,GM表示引力常数,γ表示计算点的正常重力,a表示参考椭球长半轴,r表示地心向径,Nmax表示模型最高阶数,λ、θ分别表示计算点的经度和极距,表示重力场模型的位系数。将H海和H陆组合形成全球大地水准面高数据H,H的维数是(N+1)×(M+1)。步骤二:根据H按照下两式得到全球N×M格网的垂线偏差南北分量ξij和东西分量ηij,M=360°/Δλ,表示南北方向格网边长,Δλ表示东西方向格网边长。ηij=-Hi,j+1-Hi,jdsλ---(3)]]>式中,dsλ分别表示南北方向和东西方向的格网边长。步骤三:k=1;步骤四:按照下式计算第1至第N个纬度圈的核函数序列IVξ,ij、IVη,ij,i=1,2,…N,j=1,2,…M。IVξ,ij=K'·cosα (4)IVη,ij=K'·sinα (5)式中各符号含义如下:K′=-cosψ22sin2ψ2-cosψ2(1+2sin ...
【技术特征摘要】
【专利技术属性】
技术研发人员:翟振和,孙中苗,王兴涛,李迎春,
申请(专利权)人:西安测绘研究所,
类型:发明
国别省市:陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。