[发明专利]分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法

摘要

本发明公开了一种分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法，基本步骤为：利用网函数插值对卫星测高重力数据插值；计算插值结果与船测重力数据的差值；对差值进行聚类分析，将整体区域划分成多个子区域；在每个子区域各纵向和横向网格边上对差值进行分形插值，在子区域内部进行网函数插值，得到重力数据改正值；以重力数据改正值为中心，根据卫星测高重力值融合计算得到每个子区域数据融合产品，最终得到整个区域数据融合产品。本发明提供的方法科学合理，实现过程简单易懂，在对船测重力数据和卫星测高重力数据进行融合后能够得到高精度高分辨率的重力数据，同时充分保持了多源重力数据局部特征与整体特征相似性。

主权项

1.一种分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法，其特征在于，包括以下步骤：(1)确定n个船测重力数据g_k(1≤k≤n)所在的卫星测高重力网格，利用网函数将网格角点的卫星测高重力数据G_i(1≤i≤4)插值到船测重力数据g_k(1≤k≤n)处，得到船测重力点的卫星测高重力插值结果G_k'(1≤k≤n)；(2)计算卫星测高重力插值结果G_k'(1≤k≤n)与船测重力数据g_k(1≤k≤n)的差值ΔG_k(1≤k≤n)；(3)对ΔG_k(1≤k≤n)进行聚类分析，将船测重力点区域划分成t个子区域，每个子区域有s条纵向网格边、r条横向网格边和s·r个船测重力点；(4)在第p个子区域分别沿s条纵向网格边和r条横向网格边对ΔG_k(1≤k≤s·r)进行分形插值，得到第p个子区域中每条纵向和横向网格边上的重力数据改正值ΔGA_p'(1≤p≤t)；(5)根据步骤(4)中得到的第p个子区域每条纵向和横向网格边上的重力数据改正值ΔGA_p'(1≤p≤t)，利用1‑网函数插值计算第p个子区域内的重力数据改正值ΔGB_p'(1≤p≤t)，得到第p个子区域所有重力数据改正值ΔG_p'＝{ΔGA_p'，ΔGB_p'}；(6)以第p个子区域重力数据改正值ΔG_p'(1≤p≤t)为中心，卫星测高重力网格分辨率的2～3倍为搜索半径，确定第p个子区域在搜索范围内的u个卫星测高重力数据G_pi(1≤p≤t,1≤i≤u)；(7)根据第p个子区域重力数据改正值ΔG_p'和步骤(6)中确定的卫星测高重力数据G_pi(1≤p≤t,1≤i≤u)得到第p个子区域数据融合结果(8)根据每个子区域中的数据融合结果得到t个子区域中所有的数据融合产品