分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法技术

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

【技术实现步骤摘要】
分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法
本专利技术涉及数据融合
,尤其涉及多源海洋重力数据融合的应用领域,具体说是一种分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法。
技术介绍
高精度高分辨率的海洋重力场有助于研究地球内部地质构造,同时也是研究海洋动力环境变化的重要基础,此外,高精度高分辨率的重力异常对于勘探开发海洋矿产资源、保障航天器和远程武器的发射也具有重要意义。传统意义上的海洋重力场探测主要依靠海底布设或海面舰船搭载的重力测量仪器,其精度高但成本高昂且作业周期长,重力异常数据的分辨率较低,且由于各种因素的制约,无法获得敏感区域的重力异常数据。随着航天技术的发展,利用卫星测高反演海洋重力异常技术已经成熟,同时其也成为获取大范围高分辨率海洋重力异常的唯一手段,但其不足之处在于反演重力异常精度较低。在此情形下,如何对船测高精度数据和卫星测高反演的高分辨率数据进行融合成为一大焦点性问题。目前对于两种数据的融合主要有统计法和解析法两种,统计法主要是基于统计学原理,应用最小二乘配置法进行数据融合,但这种方式需要大量观测数据构建协方差函数,且本文档来自技高网...

【技术保护点】
1.一种分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法,其特征在于,包括以下步骤:(1)确定n个船测重力数据gk(1≤k≤n)所在的卫星测高重力网格,利用网函数将网格角点的卫星测高重力数据Gi(1≤i≤4)插值到船测重力数据gk(1≤k≤n)处,得到船测重力点的卫星测高重力插值结果Gk'(1≤k≤n);(2)计算卫星测高重力插值结果Gk'(1≤k≤n)与船测重力数据gk(1≤k≤n)的差值ΔGk(1≤k≤n);(3)对ΔGk(1≤k≤n)进行聚类分析,将船测重力点区域划分成t个子区域,每个子区域有s条纵向网格边、r条横向网格边和s·r个船测重力点;(4)在第p个子区域分别沿s...

【技术特征摘要】
1.一种分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法,其特征在于,包括以下步骤:(1)确定n个船测重力数据gk(1≤k≤n)所在的卫星测高重力网格,利用网函数将网格角点的卫星测高重力数据Gi(1≤i≤4)插值到船测重力数据gk(1≤k≤n)处,得到船测重力点的卫星测高重力插值结果Gk'(1≤k≤n);(2)计算卫星测高重力插值结果Gk'(1≤k≤n)与船测重力数据gk(1≤k≤n)的差值ΔGk(1≤k≤n);(3)对ΔGk(1≤k≤n)进行聚类分析,将船测重力点区域划分成t个子区域,每个子区域有s条纵向网格边、r条横向网格边和s·r个船测重力点;(4)在第p个子区域分别沿s条纵向网格边和r条横向网格边对ΔGk(1≤k≤s·r)进行分形插值,得到第p个子区域中每条纵向和横向网格边上的重力数据改正值ΔGAp'(1≤p≤t);(5)根据步骤(4)中得到的第p个子区域每条纵向和横向网格边上的重力数据改正值ΔGAp'(1≤p≤t),利用1-网函数插值计算第p个子区域内的重力数据改正值ΔGBp'(1≤p≤t),得到第p个子区域所有重力数据改正值ΔGp'={ΔGAp',ΔGBp'};(6)以第p个子区域重力数据改正值ΔGp'(1≤p≤t)为中心,卫星测高重力网格分辨率的2~3倍为搜索半径,确定第p个子区域在搜索范围内的u个卫星测高重力数据Gpi(1≤p≤t,1≤i≤u);(7)根据第p个子区域重力数据改正值ΔGp'和步骤(6)中确定的卫星测高重力数据Gpi(1≤p≤t,1≤i≤u)得到第p个子区域数据融合结果(8)根据每个子区域中的数据融合结果得到t个子区域中所有的数据融合产品2.根据权利要求1所述的一种分形插值和网函数插值结合的卫星测高重力数据与船测重力数据融合方法,其特征在于:所述步骤(1)中网函数插值类型为1-网函数插值,公式为其中G1(L0,B0),G2(L1,B0),G3(L1,B1),G4(L0,B1)是船测重力点gk(1≤k≤n)所在卫星测高重力网格的四个数据点,船测重力点gk(1≤k≤n)坐标为(L,B),Qi(i=1,2,3,4)是过船测重力点gk(1≤k≤n)与卫星测高重力网格平行的两条直线在网格边界上截得的四个点,Ai(i=1,2,3,4)是两条直线将网格分成四个小矩形的面积,A是网格总面积,A=A1+A2+A3+A4,其中A1=(L-L0)(B-B0)...

【专利技术属性】
技术研发人员:刘善伟赵吉祥万剑华赵华刚于晓东
申请(专利权)人:中国石油大学华东山东天元信息技术股份有限公司
类型:发明
国别省市:山东,37

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1