发明名称 一种提高卫星测高交叉点精度的方法
摘要 本发明公开了一种提高卫星测高交叉点精度的方法,包括如下步骤:计算各条轨迹上相邻数据点间的距离、数据密度;选取所有轨迹中数据密度最小的轨迹作为参考轨迹,其他轨迹为待平均轨迹;为每条待平均轨迹选取插值多项式的次数;利用次多项式拟合各条待平均轨迹;计算各条拟合曲线上距各参考点最近点的经度、纬度、海面高;计算平均轨迹数据点的经度、纬度、海面高;循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高:将两轨迹间的交叉点判断转换为任一轨迹中两相邻测高数据点连线段与其余轨迹中两相邻数据点连线段是否产生交点;若产生交点,则此交点为交叉点。较传统方法,本方法具有精度高、计算简单、适用范围广、运算效率高等特点。
申请公布号 CN103363953B 申请公布日期 2015.06.17
申请号 CN201310282317.4 申请日期 2013.07.05
申请人 河海大学 发明人 岳建平;梁子亮;吕志才;顾景强;郑应新;董杰
分类号 G01C5/00(2006.01)I 主分类号 G01C5/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种提高卫星测高交叉点精度的方法,其特征在于,包括如下步骤:计算各条轨迹上相邻数据点间的距离、数据密度,其中数据密度ρ的计算公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>R</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>d</mi><mi>i</mi></msub><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>d</mi><mi>i</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mi>&rho;</mi><mo>=</mo><msqrt><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><msub><mi>R</mi><mi>i</mi></msub><mn>2</mn></msup></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></mfrac></msqrt></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow>]]></math><img file="FDA0000658592430000011.GIF" wi="627" he="445" /></maths>式中,d<sub>i</sub>为第i段相邻数据点间距离,R<sub>i</sub>为第i段相邻数据点间距离占总距离的比值;选取所有轨迹中数据密度最小的轨迹作为参考轨迹,其他轨迹为待平均轨迹;为每条待平均轨迹选取插值多项式的次数n;利用n次多项式拟合各条待平均轨迹,拟合方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mn>0</mn></msub><mo>+</mo><msub><mi>a</mi><mn>1</mn></msub><mi>x</mi><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>a</mi><mi>n</mi></msub><msup><mi>x</mi><mi>n</mi></msup></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mi>P</mi><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>y</mi><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000658592430000012.GIF" wi="588" he="164" /></maths>式中,x<sub>i</sub>为卫星测高数据点经度,y<sub>i</sub>为卫星测高数据点纬度;计算各条拟合曲线上距各参考点最近点的经度、纬度、海面高;计算平均轨迹数据点的经度、纬度、海面高;循环判断是否产生交叉点并计算交叉点的经度、纬度、海面高:将两轨迹间的交叉点判断转换为任一轨迹中两相邻测高数据点连线段与其余轨迹中两相邻数据点连线段是否产生交点;若产生交点,则此交点为交叉点。
地址 211100 江苏省南京市江宁开发区佛城西路8号