CN103363953A - 一种提高卫星测高交叉点精度的方法 - Google Patents
一种提高卫星测高交叉点精度的方法 Download PDFInfo
- Publication number
- CN103363953A CN103363953A CN2013102823174A CN201310282317A CN103363953A CN 103363953 A CN103363953 A CN 103363953A CN 2013102823174 A CN2013102823174 A CN 2013102823174A CN 201310282317 A CN201310282317 A CN 201310282317A CN 103363953 A CN103363953 A CN 103363953A
- Authority
- CN
- China
- Prior art keywords
- point
- prime
- track
- latitude
- longitude
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种提高卫星测高交叉点精度的方法,包括如下步骤:计算各条轨迹上相邻数据点间的距离、数据密度;选取所有轨迹中数据密度最小的轨迹作为参考轨迹,其他轨迹为待平均轨迹;为每条待平均轨迹选取插值多项式的次数;利用次多项式拟合各条待平均轨迹;计算各条拟合曲线上距各参考点最近点的经度、纬度、海面高;计算平均轨迹数据点的经度、纬度、海面高;循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高:将两轨迹间的交叉点判断转换为任一轨迹中两相邻测高数据点连线段与其余轨迹中两相邻数据点连线段是否产生交点;若产生交点,则此交点为交叉点。较传统方法,本方法具有精度高、计算简单、适用范围广、运算效率高等特点。
Description
技术领域
本发明涉及卫星测高技术领域,特别是涉及一种提高卫星测高交叉点精度的方法。
背景技术
共线平均是根据周期卫星的重复性而设计的一种减小卫星轨道误差的方法。研究沿共线轨迹的海面高值,可以发现海面高的长波变化,这些变化的产生主要是径向轨道误差和海面的时变部分造成的。共线平均之后,可以计算交叉点不符值,其反映的也是轨道径向误差,它的确定对卫星测高中的交叉点平差至关重要,其精度主要由交叉点位置决定。
传统方法中,交叉点位置的确定由概略位置及精确位置共同计算,上升及下降弧段的曲线拟合是求解交叉点概略位置的关键,基本方法为分段拟合法和多项式拟合法,其计算工作量大,且存在的拟合误差将会直接影响交叉点位置的精度。
发明内容
发明目的:为了解决上述问题,本发明提出了一种提高卫星测高交叉点精度的方法。
技术方案:一种提高卫星测高交叉点精度的方法,包括如下步骤:
计算各条轨迹上相邻数据点间的距离、数据密度计算各条轨迹上相邻数据点间的距离、数据密度,其中数据密度ρ的计算公式为:
式中,di为第i段相邻数据点间距离,Ri为第i段相邻数据点间距离占总距离的比值。
选取所有轨迹中数据密度最小的轨迹作为参考轨迹,其他轨迹为待平均轨迹;
为每条待平均轨迹选取插值多项式的次数n;利用n次多项式拟合各条待平均轨迹,拟合方程为:
式中,xi为卫星测高数据点经度,yi为卫星测高数据点纬度;
计算各条拟合曲线上距各参考点最近点的经度、纬度、海面高;
计算平均轨迹数据点的经度、纬度、海面高;
循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高:将两轨迹间的交叉点判断转换为任一轨迹中两相邻测高数据点连线段与其余轨迹中两相邻数据点连线段是否产生交点;若产生交点,则此交点为交叉点。
其中,所述为每条待平均轨迹选取插值多项式的次数n的方法是:
从1开始,以1为增幅逐步增加n值,并计算拟合曲线的相关系数rn,其中,相关系数r的计算公式为:
当n=N时,若rN-rN-1的值小于等于10-8,停止计算且取n=N为插值多项式的次数。
所述的计算拟合曲线上距参考点的最近的点的经度X',纬度Y',海面高H'的方法为:
确定参考轨迹点至曲线的距离,计算公式为:
式中,X0,Y0分别为参考点的经度与纬度。
将曲线拟合方程y=Pn(x)代入上式,有:
令 即:S=G(x)。
分别对S=G(x)两边求导得到:S'=G'(x)。
计算X'使其满足:
可求得最近点为:(X',Y'),其中Y'=Pn(X')。
式中,HP、HQ分别为最近点的两个相邻点的海面高;YP、YQ分别为最近点的两个相邻点的纬度值。
所述的计算平均轨迹数据点海面高H的步骤为:
计算权Pi:
式中,X,Y为平均坐标的经度与纬度值,X′i,Y′i为参考点周围第i个最近点的经度、纬度,Si为平均坐标与周围第i个最近点间的距离。
计算平均轨迹数据点海面高:
式中,H′i为参考点周围第i个最近点的海面高。
所述循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高的方法是:对所有轨迹进行编号:1,2,…,M,编号方式任意,将所有轨迹分为1条判断轨迹与M-1条待判断轨迹,其中判断轨迹按编号顺序循序选择,判断轨迹需与每条待判断轨迹进行计算;
对所有轨迹中所有相邻数据点间线段由下至上进行编号,其中判断轨迹中相邻数据点间线段称为判断线段,待判断轨迹中相邻数据点间线段称为待判断线段;
按判断线段编号顺序将每条判断线段与待判断线段进行计算,将判断线段纬度值较小的点命名为A点,另一点则为B点,待判断线段两点命名为C点与D点,C与D的位置任意。计算的基本原理为对两线段所构成的图形作如下转换:顺时针旋转α度,并向X轴偏移长度为YA的距离,向Y轴偏移长度为XA的距离,直至A点与原点重合,线段AB与X轴重合,此时A点与B点的坐标由(XA,YA),(XB,YB)转换为(0,0),(X′B,Y′B)C点与D点的坐标由(XC,YC),(XD,YD)转换为(X′C,Y′C),(X′D,Y′D)。通过以上变换,将传统的求解方程组的问题转换为求解方程的问题,以减少计算量。
计算X′B,X′C,Y′C,X′D,Y′D,α的值,公式如下:
式中,XA,YA为A点经度与纬度,XB,YB为B点经度与纬度,XC,YC为C点经度与纬度,XD,YD为D点经度与纬度;若Y′C×Y′D>0,则必不存在交叉点,跳过这两条轨迹的剩余判断;若Y′C×Y′D≤0,且则存在交叉点,且交叉点的经度,纬度及两个海面高值为:
若产生交叉点则同样跳过这两条轨迹的剩余判断,循环计算直到所有判断进行完为止。
本发明采用上述技术方案,具有以下有益效果:本发明的精确求解卫星测高交叉点的方法,提出了新的参考轨迹的选取方法,该方法所选取的参考轨迹数据点密度大,且分布均匀;提出了新的平均轨迹确定方法,该方法所求得的平均轨迹更为准确,更能代表测高卫星的轨迹,可以有效消除径向轨道误差;提出了新的交叉点判断方法,该方法简单易操作,并且可以用于常规方法所不适用的大范围交叉点的确定;提出了新的交叉点精确位置的计算方法,该方法无需计算交叉点概略位置,且避免了轨迹拟合所带来的误差。相较传统交叉点求解方法,本方案具有精度高、计算简单、适用范围广、运算效率高等特点。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
如图1所示,本发明的提高卫星测高交叉点精度的方法,包括以下步骤:
(1)分别计算各条轨迹上相邻数据点间的距离di,计算公式如下:
式中,xi,xi+1为数据点经度,yi,yi+1为数据点纬度,N为各条轨迹上数据点个数。
(2)计算各条轨迹的数据密度ρ,计算公式如下:
式中,di为第i段相邻数据点间距离,Ri为第i段相邻数据点间距离占总距离的比值。
(3)选取所有轨迹中数据密度最小的轨迹作为参考轨迹,其他轨迹为待平均轨迹;
(4)为每条待平均轨迹选取插值多项式的次数n;
(5)利用n次多项式拟合各条待平均轨迹,拟合方程为:
式中,xi为卫星测高数据点经度,yi为卫星测高数据点纬度。
(6)计算各条拟合曲线上距各参考点最近的点的经度X',纬度Y'。
(7)计算最近点的海面高H':
式中,HP、HQ分别为最近点的两个相邻点的海面高;YP、YQ分别为最近点的两个相邻点的纬度值。
(8)计算平均轨迹数据点的经度X,纬度Y:
式中,X′i为最近点经度,Y′i为最近点纬度,m为参考点周围所相应的最近点个数。
按照上式计算所有平均轨迹数据点的经度,纬度。
(9)计算所有平均轨迹数据点海面高H;
(10)卫星测高轨迹交叉点的判断及计算,其基本方法为:将两轨迹间的交叉点判断转换为任一轨迹中两相邻测高数据点连线段与其余轨迹中两相邻数据点连线段是否产生交点。若产生交点,则此交点为交叉点。
步骤(4)中选取插值多项式的次数n的方法为:
(401)从1开始,以1为增幅逐步增加n值,并计算拟合曲线的相关系数rn,其中,相关系数r的计算公式为:
(402)当n=N时,若rN-rN-1的值小于等于10-8,停止计算且取n=N为插值多项式的次数。
步骤(6)计算拟合曲线上距参考点最近的点的经度X',纬度Y'的方法为:
(601)确定参考轨迹点至曲线的距离,计算公式为:
式中,X0,Y0分别为参考点的经度与纬度。
(602)将曲线拟合方程y=Pn(x)代入上式,有:
(603)令 即:S=G(x)。
(604)分别对S=G(x)两边求导得到:S'=G'(x)。
(605)计算X'使其满足:
可求得最近点为:(X',Y'),其中Y'=Pn(X')。
步骤(9)计算平均轨迹数据点海面高H的步骤为:
(901)计算权Pi:
式中,X,Y为平均坐标的经度与纬度值,X′i,Y′i为参考点周围第i个最近点的经度、纬度,Si为平均坐标与周围第i个最近点间的距离。
(902)计算平均轨迹数据点海面高:
式中,H′i为参考点周围第i个最近点的海面高。
步骤(10)卫星测高轨迹交叉点的判断及计算具体为:
(1001)对所有轨迹进行编号:1,2,…,M,编号方式任意。将所有轨迹分为1条判断轨迹与M-1条待判断轨迹,其中判断轨迹按编号顺序循序选择。判断轨迹需与每条待判断轨迹进行计算。
(1002)对所有轨迹中所有相邻数据点间线段由下至上进行编号。其中判断轨迹中相邻数据点间线段称为判断线段,待判断轨迹中相邻数据点间线段称为待判断线段。
(1003)按判断线段编号顺序将每条判断线段与待判断线段进行计算。将判断线段纬度值较小的点命名为A点,另一点则为B点,待判断线段两点命名为C点与D点,C与D的位置任意。计算的基本原理为对两线段所构成的图形作如下转换:顺时针旋转α度,并向X轴偏移长度为YA的距离,向Y轴偏移长度为XA的距离,直至A点与原点重合,线段AB与X轴重合,此时A点与B点的坐标由(XA,YA),(XB,YB)转换为(0,0),(X′B,Y′B)C点与D点的坐标由(XC,YC),(XD,YD)转换为(X′C,Y′C),(X′D,Y′D)。通过以上变换,将传统的求解方程组的问题转换为求解方程的问题,以减少计算量。
计算X′B,X′C,Y′C,X′D,Y′D,α的值,公式如下:
式中,XA,YA为A点经度与纬度,XB,YB为B点经度与纬度,XC,YC为C点经度与纬度,XD,YD为D点经度与纬度。若Y′C×Y′D>0,则必不存在交叉点,跳过这两条轨迹的剩余判断;若Y′C×Y′D≤0,且则存在交叉点,且交叉点的经度,纬度及两个海面高值为:
若产生交叉点则同样跳过这两条轨迹的剩余判断。循环计算直到所有判断进行完为止。
Claims (5)
1.一种提高卫星测高交叉点精度的方法,其特征在于,包括如下步骤:
计算各条轨迹上相邻数据点间的距离、数据密度,其中数据密度ρ的计算公式为:
式中,di为第i段相邻数据点间距离,Ri为第i段相邻数据点间距离占总距离的比值;
选取所有轨迹中数据密度最小的轨迹作为参考轨迹,其他轨迹为待平均轨迹;
为每条待平均轨迹选取插值多项式的次数n;利用n次多项式拟合各条待平均轨迹,拟合方程为:
式中,xi为卫星测高数据点经度,yi为卫星测高数据点纬度;
计算各条拟合曲线上距各参考点最近点的经度、纬度、海面高;
计算平均轨迹数据点的经度、纬度、海面高;
循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高:将两轨迹间的交叉点判断转换为任一轨迹中两相邻测高数据点连线段与其余轨迹中两相邻数据点连线段是否产生交点;若产生交点,则此交点为交叉点。
3.根据权利要求1所述的一种精确求解卫星测高交叉点的方法,其特征在于:计算拟合曲线上距参考点的最近的点的经度X',纬度Y',海面高H'的方法为:
确定参考轨迹点至曲线的距离,计算公式为:
式中,X0,Y0分别为参考点的经度与纬度;
将曲线拟合方程y=Pn(x)代入上式,有:
令 即:S=G(x);
分别对S=G(x)两边求导得到:S'=G'(x);
计算X'使其满足:
可求得最近点为:(X',Y'),其中Y'=Pn(X');
式中,HP、HQ分别为最近点的两个相邻点的海面高;YP、YQ分别为最近点的两个相邻点的纬度值。
4.根据权利要求1所述的一种精确求解卫星测高交叉点的方法,其特征在于计算平均轨迹数据点海面高H的步骤为:
计算权Pi:
式中,X,Y为平均坐标的经度与纬度值,X′i,Y′i为参考点周围第i个最近点的经度、纬度,Si为平均坐标与周围第i个最近点间的距离;
计算平均轨迹数据点海面高:
式中,H′i为参考点周围第i个最近点的海面高。
5.根据权利要求1所述的一种提高卫星测高交叉点精度的方法,其特征在于:所述循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高的方法是:
对所有轨迹进行编号:1,2,…,M,编号方式任意,将所有轨迹分为1条判断轨迹与M-1条待判断轨迹,其中判断轨迹按编号顺序循序选择,判断轨迹需与每条待判断轨迹进行计算;
对所有轨迹中所有相邻数据点间线段由下至上进行编号,其中判断轨迹中相邻数据点间线段称为判断线段,待判断轨迹中相邻数据点间线段称为待判断线段;
按判断线段编号顺序将每条判断线段与待判断线段进行计算,将判断线段纬度值较小的点命名为A点,另一点则为B点,待判断线段两点命名为C点与D点,C与D的位置任意;计算的方法为对两线段所构成的图形作如下转换:顺时针旋转α度,并向X轴偏移长度为YA的距离,向Y轴偏移长度为XA的距离,直至A点与原点重合,线段AB与X轴重合,此时A点与B点的坐标由(XA,YA),(XB,YB)转换为(0,0),(X′B,Y′B)C点与D点的坐标由(XC,YC),(XD,YD)转换为(X′C,Y′C),(X′D,Y′D)。通过以上变换,将传统的求解方程组的问题转换为求解方程的问题,以减少计算量;计算X′B,X′C,Y′C,X′D,Y′D,α的值,公式如下:
式中,XA,YA为A点经度与纬度,XB,YB为B点经度与纬度,XC,YC为C点经度与纬度,XD,YD为D点经度与纬度;若Y′C×Y′D>0,则必不存在交叉点,跳过这两条轨迹的剩余判断;若Y′C×Y′D≤0,且则存在交叉点,且交叉点的经度,纬度及两个海面高值为:
若产生交叉点则同样跳过这两条轨迹的剩余判断,循环计算直到所有判断进行完为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310282317.4A CN103363953B (zh) | 2013-07-05 | 2013-07-05 | 一种提高卫星测高交叉点精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310282317.4A CN103363953B (zh) | 2013-07-05 | 2013-07-05 | 一种提高卫星测高交叉点精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103363953A true CN103363953A (zh) | 2013-10-23 |
CN103363953B CN103363953B (zh) | 2015-06-17 |
Family
ID=49365842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310282317.4A Expired - Fee Related CN103363953B (zh) | 2013-07-05 | 2013-07-05 | 一种提高卫星测高交叉点精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103363953B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105180876A (zh) * | 2015-08-26 | 2015-12-23 | 国家电网公司 | 复杂情况下架空输电线路交跨距离精确测量方法 |
CN106970385A (zh) * | 2017-03-30 | 2017-07-21 | 福建农林大学 | 一种异源卫星高度计观测资料的数据融合方法 |
CN108279440A (zh) * | 2018-02-01 | 2018-07-13 | 中国国土资源航空物探遥感中心 | 一种航空重力测网交叉点非遍历搜索方法 |
CN110186423A (zh) * | 2019-06-11 | 2019-08-30 | 中国石油大学(华东) | 一种卫星测高数据共线处理方法 |
CN110470414A (zh) * | 2019-08-20 | 2019-11-19 | 宏人仁医医疗器械设备(东莞)有限公司 | 一种体温测量系统以及用于体温测量系统的校正方法 |
CN111680354A (zh) * | 2020-04-20 | 2020-09-18 | 北京航空航天大学 | 近地回归轨道卫星星下点和摄影点轨迹自交点的计算方法 |
CN111982062A (zh) * | 2020-07-13 | 2020-11-24 | 国家海洋信息中心 | 一种基于最小距离法确定测高卫星轨迹交叉点位置的方法 |
CN114882105A (zh) * | 2022-05-18 | 2022-08-09 | 武汉大学 | 一种改进的卫星地面轨迹交叉点提取方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070194980A1 (en) * | 2005-08-23 | 2007-08-23 | Mitchell Douglas A | System and Method for Removal of Sea-State Bias in Satellite Altimetry Data |
CN101825473A (zh) * | 2010-04-02 | 2010-09-08 | 北京四维图新科技股份有限公司 | 一种导航方法及导航系统 |
CN101887128A (zh) * | 2010-07-09 | 2010-11-17 | 中国科学院测量与地球物理研究所 | 确定全球卫星导航系统导航卫星频间偏差的方法 |
US7962253B1 (en) * | 2007-01-24 | 2011-06-14 | Rockwell Collins, Inc. | Integrated barometric altitude and satellite altitude-based vertical navigation system |
-
2013
- 2013-07-05 CN CN201310282317.4A patent/CN103363953B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070194980A1 (en) * | 2005-08-23 | 2007-08-23 | Mitchell Douglas A | System and Method for Removal of Sea-State Bias in Satellite Altimetry Data |
US7962253B1 (en) * | 2007-01-24 | 2011-06-14 | Rockwell Collins, Inc. | Integrated barometric altitude and satellite altitude-based vertical navigation system |
CN101825473A (zh) * | 2010-04-02 | 2010-09-08 | 北京四维图新科技股份有限公司 | 一种导航方法及导航系统 |
CN101887128A (zh) * | 2010-07-09 | 2010-11-17 | 中国科学院测量与地球物理研究所 | 确定全球卫星导航系统导航卫星频间偏差的方法 |
Non-Patent Citations (3)
Title |
---|
ANANDA PASCUAL ECT.: "Improved description of the ocean mesoscale variability by combing four satillte altimeters", 《GEOPHYSICAL RESEARCH LETTERS》 * |
岳建平等: "交会法测定水平位移的精度分析", 《测绘工程》 * |
金涛勇等: "基于多源卫星测高数据的新一代全球平均海面高模型", 《测绘学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105180876A (zh) * | 2015-08-26 | 2015-12-23 | 国家电网公司 | 复杂情况下架空输电线路交跨距离精确测量方法 |
CN105180876B (zh) * | 2015-08-26 | 2018-06-15 | 国家电网公司 | 复杂情况下架空输电线路交跨距离精确测量方法 |
CN106970385A (zh) * | 2017-03-30 | 2017-07-21 | 福建农林大学 | 一种异源卫星高度计观测资料的数据融合方法 |
CN108279440A (zh) * | 2018-02-01 | 2018-07-13 | 中国国土资源航空物探遥感中心 | 一种航空重力测网交叉点非遍历搜索方法 |
CN110186423B (zh) * | 2019-06-11 | 2020-04-28 | 中国石油大学(华东) | 一种卫星测高数据共线处理方法 |
CN110186423A (zh) * | 2019-06-11 | 2019-08-30 | 中国石油大学(华东) | 一种卫星测高数据共线处理方法 |
CN110470414A (zh) * | 2019-08-20 | 2019-11-19 | 宏人仁医医疗器械设备(东莞)有限公司 | 一种体温测量系统以及用于体温测量系统的校正方法 |
CN110470414B (zh) * | 2019-08-20 | 2021-01-29 | 宏人仁医医疗器械设备(东莞)有限公司 | 一种体温测量系统以及用于体温测量系统的校正方法 |
CN111680354A (zh) * | 2020-04-20 | 2020-09-18 | 北京航空航天大学 | 近地回归轨道卫星星下点和摄影点轨迹自交点的计算方法 |
CN111680354B (zh) * | 2020-04-20 | 2022-10-21 | 北京航空航天大学 | 近地回归轨道卫星星下点和摄影点轨迹自交点的计算方法 |
CN111982062A (zh) * | 2020-07-13 | 2020-11-24 | 国家海洋信息中心 | 一种基于最小距离法确定测高卫星轨迹交叉点位置的方法 |
CN111982062B (zh) * | 2020-07-13 | 2022-02-18 | 国家海洋信息中心 | 一种基于最小距离法确定测高卫星轨迹交叉点位置的方法 |
CN114882105A (zh) * | 2022-05-18 | 2022-08-09 | 武汉大学 | 一种改进的卫星地面轨迹交叉点提取方法 |
CN114882105B (zh) * | 2022-05-18 | 2024-09-20 | 武汉大学 | 一种改进的卫星地面轨迹交叉点提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103363953B (zh) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103363953B (zh) | 一种提高卫星测高交叉点精度的方法 | |
US9068849B2 (en) | Method and system for reducing shape points in a geographic data information system | |
Hidayat et al. | Discharge estimation in a backwater affected meandering river | |
CN102879820B (zh) | 基于三角网格的三维表层模型构建方法 | |
Le Heron et al. | Calculating ice volumes and ice flux to constrain the dimensions of a 440 Ma North African ice sheet | |
CN104567906A (zh) | 一种基于北斗的城市路网车辆路径规划方法及装置 | |
CN102561118A (zh) | 一种基于轨迹偏差的高铁无砟轨道快速精调方法 | |
CN104631392A (zh) | 一种基于河相关系的航道整治设计方法 | |
Wain et al. | Propagation and dissipation of the internal tide in upper Monterey Canyon | |
US8392103B2 (en) | Database for efficient storage of track geometry and feature locations | |
CN114882105B (zh) | 一种改进的卫星地面轨迹交叉点提取方法 | |
CN107504974A (zh) | 地形分块与地形测点加权的地形匹配定位方法 | |
CN103590363A (zh) | 下切性河段水位降幅的控制节点定位防护方法 | |
Suursaar et al. | Cyclicity in ridge patterns on the prograding coasts of Estonia | |
CN106123866A (zh) | 基于航测技术的铁路既有线里程放样方法 | |
CN111982062B (zh) | 一种基于最小距离法确定测高卫星轨迹交叉点位置的方法 | |
CN103870676A (zh) | 一种适用于中国近海区域的高度计海浪平均波周期反演方法 | |
CN106054252A (zh) | 一种叠前时间偏移的方法及装置 | |
CN107807386A (zh) | 一种预测点坝砂体位置的方法及装置 | |
CN110412312B (zh) | 涌潮流速垂向分布的一种计算方法 | |
CN102930176B (zh) | 大地水准面模型九网格五样条插值方法 | |
Hwang et al. | SGOTL: A Computer Program for Modeling High-Resolution, Height-Dependent Gravity Effect of Ocean Tide Loading. | |
CN106125136A (zh) | 一种起伏地表的共中心点大道集生成方法 | |
Fine et al. | A wavefront orientation method for precise numerical determination of tsunami travel time | |
CN103399362B (zh) | 一种基于三角高程网的大气折光系数反演方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150617 Termination date: 20180705 |
|
CF01 | Termination of patent right due to non-payment of annual fee |