CN115331488B - 通用航空器的数据处理方法 - Google Patents

Abstract

本发明提供了一种通用航空器的数据处理方法，包括：在检测到通用航空器在当前时刻启动时，获取通用航空器的初始信息；在通用航空器的航电系统的运行过程中，获取通用航空器的实时飞行信息；将获取的通用航空器的实时飞行信息进行存储，得到通用航空器的飞行信息统计表；对飞行信息统计表中的数据的异常值进行处理；基于处理后的数据获取航空器在本次飞行过程中的起飞时刻、结束时刻和架次；基于本次飞行过程中的起飞时刻、结束时刻和架次得到航空器的本次飞行时间和飞行架次。本发明能够基于航空器的飞行信息自动计算航空器的每次飞行时间和飞行架次。

Description

通用航空器的数据处理方法

技术领域

本发明涉及数据处理领域，特别是涉及一种通用航空器的数据处理方法。

背景技术

通用航空作为民用航空运输体系的“两翼”之一，在满足人民日益增长的美好生活方面发挥了重要作用。近年来，通用航空飞行时间不断增加，它作为一个关键指标，对通用航空公司运行管理中飞行质量管理、航空器状态监控等方面发挥了重要作用。通用航空飞行时间不仅是衡量飞行服务站服务能力的一个要素，也是衡量通用航空公司安全运行的重点。

但目前国内通用航空公司飞行时间统计数据多采取传统人工登记方法，这种人工登记方法存在以下缺点：（1）费时费力，效率低下；（2）易出现人为差错。基于此，会出现：（1）无法快速精准了解全部学员飞行学习情况，机务人员无法精准判断飞机时寿件到寿时间；（2）无法动态监管通用航空公司的安全运行。

发明内容

针对上述技术问题，本发明采用的技术方案为：

本发明实施例提供一种通用航空器的数据处理方法，包括以下步骤：

S100，获取任一通用航空器i在第m次飞行时的初始状态信息，包括初始航向角θ^im ₀、初始位置L^im ₀和初始速度V^im ₀；i的取值为1到n；m的取值为1到N（i），N（i）为通用航空器i的航电系统在设定时间段内的启动次数；n为通用航空器的数量；L^im ₀=（lon^im ₀,lat^im ₀,h^im ₀），lon^im ₀、lat^im ₀、h^im ₀分别为通用航空器i在启动时刻对应的初始经度、初始纬度和初始高度。

S200，获取通用航空器i在第m次飞行过程中通过通信方式r在任一采集时刻t^im-r _s1的采集的飞行状态信息集（t^im-r _s1，L^im-r _s1），其中，L^im-r _s1为通过通信方式r在采集时刻t^im-r _s1获取到的通用航空器i的飞行位置，L^im-r _s1=（lon^im-r _s1，lat^im-r _s1，h^im-r _s1），lon^im-r _s1，lat^im-r _s1，h^im-r _s1分别为通用航空器i在t^im-r _s1时对应的经度、纬度和高度；s1的取值为1到p（i，m，r），p（i，m，r）为第m次飞行过程中通过通信方式r采集信息对应的采集时刻的数量；r的取值为1到N，N为通信方式的数量。

S300，对飞行状态信息进行融合处理，得到融合后的飞行状态信息。

S400，基于数据融合后的飞行状态信息获取通用航空器i的飞行航向角和飞行速度，形成通用航空器i的飞行信息统计表，其中，飞行信息统计表的第j行包括（t^im _j，θ^im _j，L^im _j， V^im _j）；θ^im _j，L^im _j， V^im _j分别为通用航空器i在采集时刻t^im _j时获取的飞行航向角、飞行位置和飞行速度，j的取值为1到max{p（i，m，r）}。

S500，对飞行信息统计表进行修正，得到修正后的飞行信息统计表。

S600，基于修正后的飞行信息统计表和初始状态信息获取通用航空器i在第m次飞行过程中的起飞时刻T^im _s、结束时刻T^im _e和飞行架次。

S700，获取通用航空器i在第m次飞行过程中的总飞行时间T^im=T^im _e-T^im _s和总飞行架次。

本发明至少具有以下有益效果：

本发明实施例提供的通用航空器的数据处理方法，能够基于航空器的飞行信息自动计算航空器的在设定时间内的飞行时间和飞行架次，精准预测飞机时寿件到寿时间，及时提醒机务人员更换时寿件，缩短维修周期，从而能够降低通用航空器在维修工作中的成本。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的通用航空器的数据处理方法的流程图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例提供的通用航空器的数据处理方法的流程图。如图1所示，本发明实施例提供一种通用航空器的数据处理方法，包括以下步骤：

S100，获取任一通用航空器i在第m次飞行时的初始状态信息，包括初始航向角θ^im ₀、初始位置L^im ₀和初始速度V^im ₀；i的取值为1到n；m的取值为1到N（i），N（i）为通用航空器i的航电系统在设定时间段内的启动次数；n为通用航空器的数量；L^im ₀=（lon^im ₀,lat^im ₀,h^im ₀），lon^im ₀、lat^im ₀、h^im ₀分别为通用航空器i在启动时刻对应的初始经度、初始纬度和初始高度。

在本发明实施例中，通用航空器中安装有机载设备，该机载设备可支持N种通信方式，在一个优选实施例中，N可等于3，例如可为北斗、ads-b、5G/4G等通信方式。通用航空器上的机载设备开启后，一次包含完成的起飞和降落时间段，即一次启动和结束之间的时间段视为通用航空器的一次飞行时间。在本发明实施例中，设定时间段可基于实际需要设置，可以日历日、日历月、日历年为单位。

在本发明实施例中，在与多个通用航空器通信连接的远程服务器检测到任一通用航空器i的机载设备第m次启动时，会通过N种通信方式中的任一种或者所有通信方式获取通用航空器i的初始状态信息。在通过N种通信方式获取初始状态信息时，优选，以最先接收到的初始状态信息为准。具体地，通用航空器的初始时间和初始位置（包括初始高度和初始经纬度信息）会通过机载设备发送给远程服务器。

S200，获取通用航空器i在第m次飞行过程中通过通信方式r在任一采集时刻t^im-r _s1的采集的飞行状态信息集（t^im-r _s1，L^im-r _s1），其中，L^im-r _s1为通过通信方式r在采集时刻t^im-r _s1获取到的通用航空器i的飞行位置，L^im-r _s1=（lon^im-r _s1，lat^im-r _s1，h^im-r _s1），lon^im-r _s1，lat^im-r _s1，h^im-r _s1分别为通用航空器i在t^im-r _s1时对应的经度、纬度和高度；s1的取值为1到p（i，m，r），p（i，m，r）为第m次飞行过程中通过通信方式r采集信息对应的采集时刻的数量；r的取值为1到N，N为通信方式的数量。

通过S200，能够获取到通过N种通信方式在第m次飞行过程中采集到的所有飞行状态信息集。

在本发明实施例中，飞行过程是指通用航空器从起飞到落地之间的过程。

S300，对飞行状态信息进行融合处理，得到融合后的飞行状态信息。

进一步地，S300可具体包括：

S310，将通过N种通信方式获取的飞行状态信息按照时间顺序进行排序，得到排序后的飞行状态信息集组D1^im=（D^im ₁，D^im ₂，…，D^im _s，…，D_E ^im），D^im _s=（t^im _s，L^im _s），L^im _s为在采集时刻t^im _s采集到的飞行位置，s的取值为1到∑^N _r=1 p（i，m，r）。

S312，从D1^im中获取q个等时间信息集组（I^im ₁，I^im ₂，…，I_u ^im，…，I_q ^im），其中，第u个等时间信息集组I_u ^im=（D_u1 ^im，D_u2 ^im，…，D_uv ^im，…，D_ub（u） ^im），D_uv ^im=（t^im _uv，L^im _uv），t^im _uv，L^im _uv为飞行状态信息集D_uv ^im中的采集时刻和飞行位置；u的取值为1到q，v的取值为1到b（u），b（u）为I_u ^im中的飞行状态信息集的数量；其中，t^im _u1=t^im _u2=…=t^im _uv=…=t^im _{u b（u）}，即每个等时间信息组为远程服务器在同一采集时刻同时从多个通信方式渠道获取到的飞行状态信息集。

S314，对I_u ^im中的相同的飞行状态信息集进行合并，得到中间等时间信息集组I^* _u ^im。

在该步骤中，相同的飞行状态信息集是指状态集中的元素均相同，即采集时刻、飞机x轴、y轴和z轴上的坐标、飞行高度和飞行经纬度信息均相同，对于这些相同的飞行状态信息集，进行合并得到一个飞行状态信息集即可。

S316，从I^* _u ^im中获取目标飞行状态信息集（t^im* _u，L^im* _u），L^im* _u基于t^im* _u、t^im* _u-f1、G^im* _u-f1和L^im* _u-f2确定，其中，t^im* _u-f1为t^im* _u的前一个采集时刻，L^im* _u-f1和L^im* _u-f2分别为t^im* _u的前一个采集时刻和前两个采集时刻对应的飞行位置。

在一个示意性实施例中，L^im* _u通过如下步骤获取：

S3160，获取在t^im* _u时的预测位置G^im* _u-0。

本领域技术人员知晓，G^im* _u-0可基于前一采集时刻和前两个采集时刻的空间位置来获取。在一个示意性实施例中，G^im* _u-0=△G^im* _u-0+G^im* _u-f1。

其中，位移变化量

；G^im* _u-0= （x^im* _u-0，y^im* _u-0，z^im* _u-0），x^im* _u-0，y^im* _u-0，z^im* _u-0分别为在t^im* _u时预测的通用航空器i在x轴、y轴 和z轴上的坐标。△G^im* _u-0为预测位置差。G^im* _u-f1为L^im* _u-f1对应的空间位置，G^im* _u-f1=（x^im* _u-f1， y^im* _{u- f1}，z^im* _{u- f1}），x^im* _u-f1，y^im* _{u- f1}，z^im* _{u- f1}分别为基于L^im* _u-f1得到的通用航空器i在x轴、y 轴和z轴上的坐标；G^im* _u-f2为L^im* _u-f2对应的空间位置，G^im* _u-f2=（x^im* _u-f2，y^im* _{u- f2}，z^im* _{u- f2}）， x^im* _u-f2，y^im* _{u- f2}，z^im* _{u- f2}分别为基于L^im* _u-f2得到的通用航空器i在x轴、y轴和z轴上的坐标。 在本发明实施例中，“→”表示向量。

本领域技术人员知晓，基于经纬度和高度获取对应的空间坐标可为现有技术。在本发明实施例中，空间坐标可为世界坐标系中的坐标。

S3162，获取位置差集△G^im* _u=（△G^im* _u1，△G^im* _u2，…，△G^im* _ua，…，△G^im* _uc（u）），其中，△G^im* _ua为I^* _u ^im中第a个飞行状态信息集中的飞行位置与G^im* _u-0之间的距离，a的取值为1到c（u），c（u）为I_u ^im*中的飞行状态信息集的数量。在一个示意性实施例中，距离可为欧式距离、余弦距离等。优选为余弦距离。

S3164，如果min（∣△G^im* _ua∣）对应的数量等于1，则G^im* _u= min（△G^im* _ua）+G^im* _u-0，如果min（∣△G^im* _ua∣）∣对应的数量大于1，则说明存在多个满足与G^im* _u-0之间的距离最短的飞行位置，此时，设置G^im* _u=G^im* _u-0。min（∣△G^im* _ua∣）表示∣△G^im* _u1∣，∣△G^im* _u2∣，…，∣△G^im* _ua∣，…，∣△G^im* _uc（u）∣中的最小值。min（△G^im* _ua）表示min（∣△G^im* _ua∣）对应的距离。

S3165，基于G^im* _u获取L^im* _u。本领域技术人员知晓，基于空间坐标获取对应的经纬度和高度可为现有技术。

在另一个示意性实施例中，S316被替换为：

S317，从I^* _u ^im中获取目标飞行状态信息集（t^im* _u，G^im* _u），其中，G^im* _u基于G^im* _u-f1和G^im* _u-f2确定，其中，G^im* _u-f1和G^im* _u-f2分别为t^im* _u的前一个采集时刻和前两个采集时刻对应的飞行位置。

进一步地，S317可包括：

S3170，获取预测位置G^im* _u-0=2*G^im* _u-f1-G^im* _u-f2。

S3172，获取位置差集△G^im* _u=（△G^im* _u1，△G^im* _u2，…，△G^im* _ua，…，△G^im* _uc（u）），其中，△G^im* _ua为I^* _u ^im中第a个飞行状态信息集中的飞行位置与G^im* _u-0之间的距离，a的取值为1到c（u），c（u）为I_u ^im*中的飞行状态信息集的数量。在一个示意性实施例中，距离可为欧式距离、余弦距离等。优选为余弦距离。

S3174，如果min（∣△G^im* _ua∣）对应的数量等于1，则G^im* _u= min（△G^im* _ua）+G^im* _u-0，如果min ∣△G^im* _ua∣）对应的数量大于1，则G^im* _u=G^im* _u-0。

S318，将D1^im中I_u ^im替换为目标飞行状态信息集，得到目标飞行状态信息集组。

具体地，在D1^im中，将I_u ^im中除目标飞行状态信息集之外的飞行状态信息集删除，就可得到目标飞行状态信息组。

进一步，在本发明实施例中，在S318之后还可包括：

S319, 遍历目标飞行状态信息集组，对于目标飞行状态信息集组中相邻的两个飞行状态信息集中的采集时刻，如果两个采集时刻之间的差值大于设定阈值，则基于线性插值法在该相邻的两个飞行状态信息集之间添加飞行状态信息集。

在本发明实施例中，设定阈值可基于实际情况进行确定。例如，可为k*△t_max，k为大于1的自然数，△t_max为N中通信方式中数据发送间隔最大的通信方式对应的数据发送时间间隔。在一个示例中，优选，k可为3~5，优选，k可为3。在一个具体示例中，设定阈值可为5分钟。

本领域技术人员知晓，基于线性插值法在该相邻的两个飞行状态信息集之间添加飞行状态信息集可采用现有技术。

S319的技术效果在于，若在山区，5G/4G无通信基站，北斗更新速度慢以及ADS-B有遮蔽，影响通信时，可根据数据序列完好性，进行线性插值，填补无通信记录时刻的数据，从而能够使得获取的数据更加完整。

S300的技术效果在于，通过多种通信方式获取飞行信息并进行融合，能够获取尽可能多的飞行信息，能够确保数据完整性和可靠性。

S400，基于数据融合后的飞行状态信息获取通用航空器i的飞行航向角和飞行速度，形成通用航空器i的飞行信息统计表，其中，飞行信息统计表的第j行包括（t^im _j，θ^im _j，L^im _j， V^im _j）；θ^im _j，L^im _j， V^im _j分别为通用航空器i在采集时刻t^im _j时获取的飞行航向角、飞行位置和飞行速度，j的取值为1到max{p（i，m，r）}，max{p（i，m，r）}为p（i，m，1），p（i，m，2），…，p（i，m，r），…，p（i，m，N）中的最大值。

本领域技术人员知晓，基于通用航空器的飞行时间和飞行位置获取对应的飞行航向角和飞行速度可为现有技术。

S500，对飞行信息统计表进行修正，得到修正后的飞行信息统计表。

在本发明实施例中，S500可具体包括：

对于修正后的飞行信息统计表中第j行数据，如果V^im _j>Vⁱ _max或θ^im _j>360°或h^im _j<Hⁱ _A，则删除第j行数据，其中，Vⁱ _max为通用航空器i的最大飞行速度，Hⁱ _A为通用航空器i对应的机场的标高。

在本发明实施例中，Vⁱ _max可基于通用航空器i的任务性质、性能等获取。S500的技术效果在于，通过对实时航空信息中异常值进行修正，能够确保数据准确。

S600，基于修正后的飞行信息统计表和初始状态信息获取通用航空器i在第m次飞行过程中的起飞时刻T^im _s、飞行架次和结束时刻T^im _e。

具体地，S600可具体包括：

遍历修正后的飞行信息统计表，如果θ^im _k≠θ^im ₀，V^im _k+1>V^im _k>V^im _k-1>V^im ₀，h^im _k+1>h^im _k>h^im _k-1=h^im ₀，并且，lon^im _k和lat^im _k位于通用航空器i的飞行跑道对应的经纬度区域内，说明航空器i在起飞，则设置通用航空器i在第m次飞行过程中的起飞时刻T^im _s=t^im _k。如果θ^im _k+1=θ^im _k=θ^im _k-1=θ^im _S，V^im _k>V^im _k-1>V^im ₀，h^im _k-1>h^im _k =h^im ₀，并且，lon^im _k和lat^im _k位于通用航空器i的飞行跑道对应的经纬度区域内，说明通用航空器i已经完成一次飞行，则设置C_im=C_im+1，C_im为通用航空器i在第m次飞行过程中的架次计数器，初始值为0。θ^im _S为通用航空器i在t^im _k时的起飞航向角。

如果θ^im _k=θ^im _k-1=θ^im ₀，V^im ₀=V^im _k<V^im _k-1，h^im ₀=h^im _k=h^im _k-1，lon^im _k=lon^im ₀，lat^im _k=lat^im ₀，说明通用航空器i结束飞行，则设置通用航空器i在第m次飞行过程中的结束时刻T^im _e=t^im _k；k的取值为1到f（im），f（im）为修正后的飞行信息统计表中的行数。

在本发明实施例中，通用航空器i的飞行跑道对应的经纬度区域可根据通用航空器i的飞行跑道在地图中的经纬度信息网格化得到，具体实现方法可为现有方法。

S700，获取通用航空器i在第m次飞行过程中的总飞行时间T^im=T^im _e-T^im _s和总飞行架次。总飞行架次为S600中得到的C_im。

在本发明实施例中，飞行时间可采用60进制计算。

本发明实施例提供的通用航空器的数据处理方法，能够基于航空器的飞行信息自动计算航空器的每次飞行时间和飞行架次，能够提高飞行信息的获取速度和减少运行人员工作量。

进一步地，本发明实施例提供的数据处理方法还可包括：

如果T^im 大于设定预警时间和/或C_im大于设定预警架次，则进行预警。

设定预警时间和设定预警架次可基于实际需要进行设置。在一个具体实施例中，如果T^im和/或C_im达到50的倍数时，即，设定预警时间可为50*M小时，设定预警架次可为50*M架次，M为大于等于1的整数，例如，设定预警时间可为50小时,100小时，150小时，…。设定预警架次可为50架次，100架次，150架次，…。

预警信息可进行可视化显示，例如通过声光预警方式进行显示，以提醒运行人员对时寿件进行管理控制。

进一步地，本发明实施例提供的数据处理方法还可包括：

基于获取的每个通用航空器在每次飞行过程中的总飞行时间和总飞行架次，获取任一机型在指定时间段内的总飞行时间、总飞行架次和机型利用率。

由于数据库中事先存储有每个通用航空器的机型，因此，可基于获取的每个通用 航空器在每次飞行过程中的总飞行时间和总飞行架次，获取任一机型在指定时间段内的总 飞行时间、总飞行架次和机型利用率。例如，对于机型a的通用航空器一共有f（a）架，在指定 时间段TB内的总飞行时间

，T^a _d为f（a）架中的第d架通用航空器在指定时间段 内的总飞行时间，可基于指定时间段内的总飞行次数确定，每次飞行的总飞行时间可基于 S100~S700所示出的方法获取。同理，总飞行架次与总飞行时间类似。机型a在指定时间段内 的机型利用率等于

，T^a _d1为指定时间段TB内，f（a）架中通用航空器中可以使 用的f1（a）架通用航空器中的第d1架通用航空器的可利用时间，每架通用航空器的可利用 时间可基于实际情况进行确定。

在一个示意性实施例中，指定时间段可为一天。

进一步地，本发明实施例提供的数据处理方法还可包括：

基于获取的每个通用航空器在每次飞行过程中的总飞行时间和总飞行架次，获取任一通航公司的航空器的总飞行时间、总飞行架次、航空器利用率，以及指定时间段内的执行任务空域的利用率。

在本发明实施例中，数据库中会事先存储有每个通用航空器所属的通航公司和执行任务的机场，因此，可基于获取的每个通用航空器在每次飞行过程中的总飞行时间和总飞行架次，获取任一通航公司的航空器在指定时间段内的总飞行时间、总飞行架次、航空器利用率，以及指定时间段内的执行任务空域的利用率。

在一个示例中，例如，对于通航公司b，所拥有通用航空器一共有f（b）架，在指定时间段TB内的总飞行时间T^b _TB=∑^f(b) _d=1 T_d ^b, T^b _d为f（b）架中的第d架通用航空器在指定时间段内的总飞行时间，可基于指定时间段内的总飞行次数确定，每次飞行的总飞行时间可基于S100~S700所示出的方法获取。同理，总飞行架次与总飞行时间类似。通航公司b在指定时间段内的机型利用率等于T^b _TB/(∑^f1(b) _d1=1 T_d1 ^b),T^b _d1为指定时间段TB内，f（b）架中通用航空器中可以使用的f1（b）架通用航空器中的第d1架通用航空器的可利用时间，每架通用航空器的可利用时间可基于实际情况进行确定。

在一个示例中，例如，对于机场c，指定时间段内的执行任务空域的利用率等于∑^Tc _d2=1 T_d2 ^c和∑^Tc _d3=1 T_d3 ^c的比值，T^c _d2为机场c在第d2天的使用时间，等于在第d2天的机场飞行最晚结束时间和机场最早开飞时间之间的差值，Tc为指定时间段所包含的天数，T^c _d3为机场c在第d3天的可用于飞行的开放时间，即可利用时间。

本发明的实施例还提供了一种非瞬时性计算机可读存储介质，该存储介质可设置于电子设备之中以保存用于实现方法实施例中一种方法相关的至少一条指令或至少一段程序，该至少一条指令或该至少一段程序由该处理器加载并执行以实现上述实施例提供的方法。

本发明的实施例还提供了一种电子设备，包括处理器和前述的非瞬时性计算机可读存储介质。

本发明的实施例还提供一种计算机程序产品，其包括程序代码，当所述程序产品在电子设备上运行时，所述程序代码用于使该电子设备执行本说明书上述描述的根据本发明各种示例性实施方式的方法中的步骤。

虽然已经通过示例对本发明的一些特定实施例进行了详细说明，但是本领域的技术人员应该理解，以上示例仅是为了进行说明，而不是为了限制本发明的范围。本领域的技术人员还应理解，可以对实施例进行多种修改而不脱离本发明的范围和精神。本发明公开的范围由所附权利要求来限定。

Claims (7)

1.一种通用航空器的数据处理方法，其特征在于，包括以下步骤：

S100，获取任一通用航空器i在第m次飞行时的初始状态信息，包括初始航向角θ^im ₀、初始位置L^im ₀和初始速度V^im ₀；i的取值为1到n；m的取值为1到N（i），N（i）为通用航空器i的航电系统在设定时间段内的启动次数；n为通用航空器的数量；L^im ₀=（lon^im ₀,lat^im ₀,h^im ₀），lon^im ₀、lat^im ₀、h^im ₀分别为通用航空器i在启动时刻对应的初始经度、初始纬度和初始高度；

S200，获取通用航空器i在第m次飞行过程中通过通信方式r在任一采集时刻t^im-r _s1采集的飞行状态信息集（t^im-r _s1，L^im-r _s1），其中，L^im-r _s1为通过通信方式r在采集时刻t^im-r _s1获取到的通用航空器i的飞行位置，L^im-r _s1=（lon^im-r _s1，lat^im-r _s1，h^im-r _s1），lon^im-r _s1，lat^im-r _s1，h^im-r _s1分别为通用航空器i在t^im-r _s1时对应的经度、纬度和高度；s1的取值为1到p（i，m，r），p（i，m，r）为第m次飞行过程中通过通信方式r采集信息对应的采集时刻的数量；r的取值为1到N，N为通信方式的数量；

S300，对飞行状态信息集中的飞行状态信息进行融合处理，得到融合后的飞行状态信息；

S400，基于融合后的飞行状态信息获取通用航空器i的飞行航向角和飞行速度，形成通用航空器i的飞行信息统计表，其中，飞行信息统计表的第j行包括（t^im _j，θ^im _j，L^im _j，V^im _j）；θ^im _j，L^im _j，V^im _j分别为通用航空器i在采集时刻t^im _j时获取的飞行航向角、飞行位置和飞行速度，j的取值为1到max{p（i，m，r）}；

S500，对飞行信息统计表进行修正，得到修正后的飞行信息统计表；

S600，基于修正后的飞行信息统计表和初始状态信息获取通用航空器i在第m次飞行过程中的起飞时刻T^im _s、结束时刻T^im _e和飞行架次；

S700，获取通用航空器i在第m次飞行过程中的总飞行时间T^im=T^im _e-T^im _s和总飞行架次；

S300包括：

S310，将通过N种通信方式获取的飞行状态信息按照时间顺序进行排序，得到排序后的飞行状态信息集组D1^im=（D^im ₁，D^im ₂，…，D^im _s，…，D_E ^im），D^im _s=（t^im _s，L^im _s），L^im _s为在采集时刻t^im _s采集到的飞行位置，s的取值为1到∑^N _r=1 p（i，m，r）；

S312，从D1^im中获取q个等时间信息集组（I^im ₁，I^im ₂，…，I_u ^im，…，I_q ^im），其中，第u个等时间信息集组I_u ^im=（D_u1 ^im，D_u2 ^im，…，D_uv ^im，…，D_ub（u） ^im），D_uv ^im=（t^im _uv，L^im _uv），t^im _uv，L^im _uv分别为飞行状态信息集D_uv ^im中的采集时刻和飞行位置；u的取值为1到q，v的取值为1到b（u），b（u）为I_u ^im中的飞行状态信息集的数量；其中，t^im _u1=t^im _u2=…=t^im _uv=…=t^im _{u b（u）}；

S314，对I_u ^im中的相同的飞行状态信息集进行合并，得到中间等时间信息集组I^* _u ^im；

S316，从I^* _u ^im中获取目标飞行状态信息集（t^im* _u，L^im* _u），L^im* _u为在采集时刻t^im* _u采集到的飞行位置，L^im* _u基于t^im* _u、t^im* _u-f1、G^im* _u-f1和L^im* _u-f2确定，其中，t^im* _u-f1为t^im* _u的前一个采集时刻，L^im* _u-f1和L^im* _u-f2分别为t^im* _u的前一个采集时刻和前两个采集时刻对应的飞行位置，G^im* _u-f1为L^im* _u-f1对应的空间位置；

S318，将D1^im中I_u ^im替换为目标飞行状态信息集，得到目标飞行状态信息集组；

在S316中，L^im* _u通过如下步骤获取：

S3160，获取在t^im* _u时的预测位置G^im* _u-0；

S3162，获取位置差集△G^im* _u=（△G^im* _u1，△G^im* _u2，…，△G^im* _ua，…，△G^im* _uc（u）），其中，△G^im* _ua为I^* _u ^im中第a个飞行状态信息集中的飞行位置与G^im* _u-0之间的距离，a的取值为1到c（u），c（u）为I^* _u ^im中的飞行状态信息集的数量；

S3164，如果min（∣△G^im* _ua∣）对应的数量等于1，则G^im* _u= min（△G^im* _ua）+G^im* _u-0，如果min（∣△G^im* _ua∣）对应的数量大于1，则G^im* _u=G^im* _u-0；

S3165，基于G^im* _u获取L^im* _u。

2.根据权利要求1所述的方法，其特征在于，S600进一步包括：

遍历修正后的飞行信息统计表，如果θ^im _k≠θ^im ₀， V^im _k+1>V^im _k>V^im _k-1>V^im ₀，h^im _k+1>h^im _k>h^im _k-1=h^im ₀，并且，lon^im _k和lat^im _k位于通用航空器i的飞行跑道对应的经纬度区域内，则设置通用航空器i在第m次飞行过程中的起飞时刻T^im _s=t^im _k；如果θ^im _k+1=θ^im _k=θ^im _k-1=θ^im _S，V^im _k>V^im _k-1>V^im ₀，h^im _k-1>h^im _k =h^im ₀，并且，lon^im _k和lat^im _k位于通用航空器i的飞行跑道对应的经纬度区域内，则设置C_im=C_im+1，C_im为通用航空器i在第m次飞行过程中的架次计数器，初始值为0；θ^im _S为通用航空器i在t^im _k时的起飞航向角；如果θ^im _k=θ^im _k-1=θ^im ₀，V^im ₀=V^im _k<V^im _k-1，h^im ₀=h^im _k=h^im _k-1，lon^im _k=lon^im ₀，lat^im _k=lat^im ₀，则设置通用航空器i在第m次飞行过程中的结束时刻T^im _e=t^im _k；k的取值为1到f（im），f（im）为修正后的飞行信息统计表中的行数。

3.根据权利要求1所述的方法，其特征在于，还包括：

S319, 遍历目标飞行状态信息集组，对于目标飞行状态信息集组中相邻的两个飞行状态信息集中的采集时刻，如果两个采集时刻之间的差值大于设定阈值，则基于线性插值法在该相邻的两个飞行状态信息集之间添加飞行状态信息集。

4.根据权利要求1所述的方法，其特征在于，S500具体包括：

对于ST_im中第j行数据，如果V^im _j>Vⁱ _max或θ^im _j>360°或h^im _j<Hⁱ _A，则删除第j行数据，其中，Vⁱ _max为通用航空器i的最大飞行速度，Hⁱ _A为通用航空器i对应的机场的标高。

5.根据权利要求2所述的方法，其特征在于，还包括：

如果T^im 大于设定预警时间和/或C_im大于设定预警架次，则进行预警。

6.根据权利要求1所述的方法，其特征在于，还包括：

基于获取的每个通用航空器在每次飞行过程中的总飞行时间和总飞行架次，获取任一机型在指定时间段内的总飞行时间、总飞行架次和机型利用率。

7.根据权利要求1所述的方法，其特征在于，还包括：

基于获取的每个通用航空器在每次飞行过程中的总飞行时间和总飞行架次，获取任一通航公司的航空器的总飞行时间、总飞行架次、航空器利用率，以及指定时间段内的执行任务空域的利用率。

Images