卫星测高数据共线处理方法.pdf

《卫星测高数据共线处理方法.pdf》由会员分享，可在线阅读，更多相关《卫星测高数据共线处理方法.pdf（9页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910501238.5 (22)申请日 2019.06.11 (71)申请人 中国石油大学 （华东） 地址 266580 山东省青岛市黄岛区长江西 路66号 (72)发明人 孙钦婷万剑华刘善伟 (51)Int.Cl. G01C 5/00(2006.01) (54)发明名称 一种卫星测高数据共线处理方法 (57)摘要 本发明公开了一种卫星测高数据共线处理 方法， 涉及卫星测高数据处理技术领域， 本发明 以卫星测高数据为基础， 利用二次多项式构建参 考轨道， 能够更加逼近离散数。

2、据； 然后利用卫星 沿轨采样间隔计算出参考轨道上正常点的经纬 度， 减少了数据的缺失； 最后利用距离加权法确 定各正常点的海面高， 从而将同一编号不严密重 复的多条轨道拟合成一条轨道。 权利要求书1页 说明书3页 附图4页 CN 110186423 A 2019.08.30 CN 110186423 A 1.一种卫星测高数据共线处理方法， 其特征在于， 包括以下步骤： (1)对卫星高度计的GDR数据进行预处理得到卫星测高数据点； (2)利用步骤(1)中在同一轨道不同周期的卫星测高数据点拟合出参考轨道； (3)计算参考轨道上第一个正常点的经纬度； (4)利用步骤(3)得到的参考轨道上第一个正常点。

3、的经纬度和卫星沿轨采样间隔S计算 得到参考轨道上所有正常点的经纬度； (5)以步骤(4)得到的各正常点为中心， 每个正常点获取搜索半径S范围内的个卫星 测高数据点， 并分别利用距离加权平均法确定各正常点的海面高， 实现卫星测高数据的共 线处理。 其中， S为卫星沿轨采样间隔。 2.根据权利要求1所述的一种卫星测高数据共线处理方法， 其特征在于： 所述步骤(2) 中拟合方法为二次多项式拟合， 具体如下所示： 式中， xi为各卫星测高数据点的经度， yi为各卫星测高数据点的纬度， a0、 a1、 a2为待求系 数， f(x)是二次多项式拟合的公式。 3.根据权利要求1所述的一种卫星测高数据共线处理。

4、方法， 其特征在于： 所述步骤(3) 中参考轨道上第一个正常点的纬度Y1为卫星测高数据点中纬度最大值； 经度X1根据公式(1) 计算得到。 4.根据权利要求1所述的一种卫星测高数据共线处理方法， 其特征在于： 所述步骤(4) 中正常点经纬度的计算公式为： 式中， R为地球平均半径6371km， S为卫星沿轨采样间隔(单位为千米)， Xi为正常点的经 度， Yi为正常点的纬度。 5.根据权利要求1所述的一种卫星测高数据共线处理方法， 其特征在于： 所述步骤(5) 中若0则删除该正常点， 若0则正常点的海面高计算公式为： 式中， xj为卫星测高数据点的经度， yj为卫星测高数据点的纬度， X为正常。

5、点的经度， Y为 正常点的纬度， dj为搜索半径内各卫星测高数据点到正常点的距离， sshj为搜索半径内各卫 星测高数据点的海面高， ssh为正常点的海面高。 权利要求书 1/1 页 2 CN 110186423 A 2 一种卫星测高数据共线处理方法 技术领域 0001 本发明涉及卫星测高数据处理技术领域， 具体涉及一种卫星测高数据共线处理方 法。 背景技术 0002 共线处理是卫星测高数据处理的首要步骤。 在理想状态下， 每个周期的卫星轨道 都是严密重合的。 但是由于轨道偏移、 仪器、 引力等因素的影响， 卫星的重复轨道并不能精 确共线， 其轨迹之间的最大距离可达1公里。 共线平差是根据具有。

6、重复周期卫星测高任务特 点而设计的一种消除卫星轨道误差并确定平均海平面及其变化的方法。 0003 使用传统共线处理方法时， 首先在相同轨道编号中确定一条观测数据质量最好、 数据量最多的轨道作为参考轨道。 并把参考轨道上的观测点作为正常点， 最后利用相同编 号临近的观测值通过线性内插计算出正常点海面高。 利用上述方法得到的参考轨道位置或 许位于所有轨道的边缘， 再进行插值计算误差较大， 并且会忽略许多位置上的观测点。 鉴于 此， 本发明利用二次多项式构建参考轨道， 能够更加逼近离散数据。 然后利用卫星沿轨采样 间隔计算出参考轨道上正常点的经纬度， 减少了数据的缺失。 最后利用距离加权法确定各 正。

7、常点的海面高， 从而将同一编号不严密重复的多条轨道拟合成一条轨道， 有利于后续卫 星测高数据交叉点平差的进行。 发明内容 0004 (一)要解决的技术问题 0005 为解决上述问题， 本发明提供了一种卫星测高数据共线处理方法。 0006 (二)技术方案 0007 本发明包含以下步骤： 0008 (1)对卫星高度计的GDR数据进行预处理得到卫星测高数据点； 0009 (2)利用步骤(1)中在同一轨道不同周期的卫星测高数据点拟合出参考轨道； 0010 (3)计算参考轨道上第一个正常点的经纬度； 0011 (4)利用步骤(3)得到的参考轨道上第一个正常点的经纬度和卫星沿轨采样间隔S 计算得到参考轨道。

8、上所有正常点的经纬度； 0012 (5)以步骤(4)得到的各正常点为中心， 每个正常点获取搜索半径S范围内的个 卫星测高数据点， 并分别利用距离加权平均法确定各正常点的海面高， 实现卫星测高数据 的共线处理。 0013 其中， S为卫星沿轨采样间隔。 0014 进一步， 所述步骤(1)中从卫星高度计的GDR读取卫星测高数据点的经度、 纬度和 海面高， 对经度和纬度预处理方法为： 0015 以0度经线为基准， 东经为正值， 西经为负值， 北纬取90减该位置纬度值， 南纬取90 加该位置纬度值。 说明书 1/3 页 3 CN 110186423 A 3 0016 进一步， 所述步骤(2)中拟合方法。

9、为二次多项式拟合， 具体如下所示： 0017 0018 式中， xi为各卫星测高数据点的经度， yi为各卫星测高数据点的纬度， a0、 a1、 a2为待 求系数， f(x)是二次多项式拟合的公式。 0019 进一步， 所述步骤(3)中参考轨道上第一个正常点的纬度Y_1为卫星测高数据点中 纬度最大值； 经度X_1根据公式(1)计算得到。 0020 进一步， 所述步骤(4)中正常点经纬度的计算公式为： 0021 0022 式中， R为地球平均半径6371km， S为卫星沿轨采样间隔(单位为千米)， Xi为正常点 的经度， Yi为正常点的纬度。 0023 进一步， 所述步骤(5)中若0则删除该正常点。

10、， 若0则正常点的海面高计算 公式为： 0024 0025 式中， xj为卫星测高数据点的经度， yj为卫星测高数据点的纬度， X为正常点的经 度， Y为正常点的纬度， dj为搜索半径内各卫星测高数据点到正常点的距离， sshj为搜索半径 内各卫星测高数据点的海面高， ssh为正常点的海面高。 0026 (三)有益效果 0027 本发明的优点体现在： 0028 本发明以卫星测高数据为基础， 能够利用二次多项式构建参考轨道， 能够更加逼 近离散数据。 然后利用卫星沿轨采样间隔计算出参考轨道上的正常点的经纬度， 减少了数 据的缺失。 最后利用距离加权法确定各正常点的海面高， 从而将同一编号不严密重。

11、复的多 条轨道拟合成一条轨道。 附图说明 0029 图1为本发明实施的步骤流程图， 0030 图2为拟合参考轨道示意图， 0031 图3为正常点经纬度计算示意图， 0032 图4为正常点海面高计算示意图。 具体实施方式 0033 为使本发明的目的、 内容、 和优点更加清楚， 下面结合附图和实施例， 对本发明的 说明书 2/3 页 4 CN 110186423 A 4 具体实施方式作进一步详细描述： 0034 参照图1， 本发明的具体实施步骤为： 0035 (1)对卫星高度计的GDR数据进行预处理得到卫星测高数据点； 0036 其中， 从卫星高度计的GDR读取卫星测高数据点的经度、 纬度和海面高。

12、， 对经度和 纬度预处理方法为： 0037 以0度经线为基准， 东经为正值， 西经为负值， 北纬取90减该位置纬度值， 南纬取90 加该位置纬度值。 0038 (2)参照图2， 利用步骤(1)中在同一轨道不同周期的卫星测高数据点拟合出参考 轨道； 0039 其中拟合方法为二次多项式拟合， 具体如下所示： 0040 0041 式中， xi为各卫星测高数据点的经度， yi为各卫星测高数据点的纬度， a0、 a1、 2为待 求系数， f(x)是二次多项式拟合的公式。 0042 (3)计算参考轨道上第一个正常点的经纬度； 0043 其中参考轨道上第一个正常点的纬度Y1为卫星测高数据点中纬度最大值； 经。

13、度X1 根据公式(1)计算得到。 0044 (4)参照图3， 利用步骤(3)得到的参考轨道上第一个正常点的经纬度和卫星沿轨 采样间隔S计算得到参考轨道上所有正常点的经纬度； 0045 其中正常点经纬度的计算公式为： 0046 0047 式中， R为地球平均半径6371km， S为卫星沿轨采样间隔(单位为千米)， Xi为正常点 的经度， Yi为正常点的纬度。 0048 (5)参照图4， 以步骤(4)得到的各正常点为中心， 每个正常点获取搜索半径S范围 内的个卫星测高数据点， 并分别利用距离加权平均法确定各正常点的海面高， 实现卫星 测高数据的共线处理。 0049 其中， S为卫星沿轨采样间隔。 。

14、若0则删除该正常点， 若0则正常点的海面 高计算公式为： 0050 0051 式中， xj为卫星测高数据点的经度， yj为卫星测高数据点的纬度， X为正常点的经 度， Y为正常点的纬度， dj为搜索半径内各卫星测高数据点到正常点的距离， sshj为搜索半径 内各卫星测高数据点的海面高， ssh为正常点的海面高。 说明书 3/3 页 5 CN 110186423 A 5 图1 说明书附图 1/4 页 6 CN 110186423 A 6 图2 说明书附图 2/4 页 7 CN 110186423 A 7 图3 说明书附图 3/4 页 8 CN 110186423 A 8 图4 说明书附图 4/4 页 9 CN 110186423 A 9 。