一种针对飞机飞行轨迹数据跳变震荡的处理方法与流程

技术特征：

1.一种针对飞机飞行轨迹数据跳变震荡的处理方法，其特征是，包括对飞行轨迹数据中高度跳点的处理，且包括以下步骤：

S1、输入含飞机经纬度、高度、水平速度、竖直速度、时间信息的离散点集合；

S2、对离散点集合中的每一个点进行处理：利用前一个有效点与当前点比较，用速度和时间计算两点间距离，与读入的原始数据对比，若相差过大，则更正或删除；

S3、输出无高度跳点的离散点集合。

2.如权利要求1所述的一种针对飞机飞行轨迹数据跳变震荡的处理方法，其特征是，包括以下步骤：

S1、输入离散点；

其中：(x,y)：经纬度，height：高度，sv：水平速度，hv：竖直速度，time：时间；

S2、读入一个离散点Pi(∈Original−Cdata)如果经纬度数据缺失或者高度为负，移除；否则：

若为第一个数据，Cdata^∗.pushback(Pi)；

若非第一个数据，与前一个离散点exP(∈Cdata^∗)进行比较：

如果两点时间间隔超过δt，Cdata^∗.pushback(Pi)；

否则：

∗realDistance：用三维坐标计算两点之间平行于地球表面的距离；

其中R=6371004米为地球半径，P1(x1,y1,h1)，P2(x2,y2,h2)为给定两点三维坐标；

∗computedDistance：用水平速度和时间计算两点之间的距离；

∗realHeight=Pi.height：原始数据Pi的高度；∗computedHeight：用前一个数据的竖直速度和时间计算的高度；

如果：d<ds且h<dv，Cdata^∗.pushback(Pi)；

如果：d<ds且h≥dv，用速度和时间计算出的高度computedHeight更新Cdata的高度数据，Cdata^∗.pushback(Pi)；

否则移除；

∗这里取δt=10，ds=0.25，ds=0.1；即，间隔超过10秒不做比较直接存入，距离误差与高度误差分别小于25%和10%；

S3、输出新的离散点；

使得每个Pi^∗数据尽可能完整，且无高度跳点。

3.如权利要求1或2所述的一种针对飞机飞行轨迹数据跳变震荡的处理方法，其特征是，还包括对飞行轨迹数据中经纬度震荡的处理，且包括以下步骤：

S′1、输入含飞机经纬度信息的离散点集合；

S′2、对离散点集合中的每一个点进行处理：给定顺序相邻三个点：...,Pi,Pi+1,Pi+2,...；

若以中间点Pi+1为顶点，Pi+1Pi,Pi+1Pi+2为边的角度过小，则交换Pi+1和Pi+2，...,Pi,Pi+2,Pi+1,...；再检查这三个点之前的一个点和这三个点交换以后的前两个点，即Pi−1,Pi,Pi+2；反之，三点顺序不变，检查这三个点的后两个点和接下来的一个点，即Pi+1,Pi+2,Pi+3；

S′3、输出无经纬度震荡的离散点集合。

4.如权利要求3所述的一种针对飞机飞行轨迹数据跳变震荡的处理方法，其特征是，包括以下步骤：

S′1、输入离散点；

其中，(x,y)：经纬度；

S′2、先对数据进行分类，将经纬度相同的点按时间顺序放在一起，，

Qi={Pi|所有Pi经纬度相同}；

按顺序对相邻三个集合Qi,Qi+1,Qi+2进行分析，从三个集合中各取一个点，Pl∈Qi，

Pj∈Qi+1，Pk∈Qi+2：

计算每两个点之间的距离dlj=|Pj−Pl|,dlk=|Pk−Pl|,djk=|Pk−Pj|给定两点Pj(xj,yj)，Pk(xk,yk)，距离djk计算公式：

φ=arccos[cos(90−yk)cos(90−yj)+sin(90−yk)sin(90−yj)cos(xj−xk)]

其中R=6371004米为地球半径；

计算PjPl与PjPk的夹角ϕl：

•若ϕi过小，Qi+1与Qi+2交换位置，对Qi−1,Qi,Qi+2进行分析，转到（a）（b）；

•若Qi+1与Qi+2不需要交换位置，对Qi+1,Qi+2,Qi+3进行分析，转到（a）（b）；

S′3、输出离散点为对下标i重新排序，使得这组数据Cdata^∗在经纬度坐标下不会来回跳。