本发明属于民航飞行安全领域，更具体地说，涉及一种飞机在飞行过程中所遇环境的检测方法，可用于评价飞行环境的飞行操作难度。
背景技术：
近三十年来，各国航空公司的飞机已逐渐加装了qar(quickaccessrecorder，快速存储记录器)系统[1,2]，其记录数据包括飞机着陆时的磁航向角、风向、风速、气温等，为研究foqa(flightoperationsqualityassurance，飞行操作品质监控)提供了数据基础，具有实时性、普适性、易用性的特点。
多年的研究发现，飞行事故是由飞行员、环境、飞机三者共同作用[3]所引发的。飞机的飞行过程是人(飞行员操作)，机(飞机运行状态)，环(航线环境)综合作用的结果。因此，如何定量评价航线环境具有非常重要的实用价值。更重要的是，对飞行环境的评价将有助于实时评价空域以及机场的着陆条件，据此能够在同等飞行条件下对飞行员的驾驶水平给出更客观的评价，指导航空公司改进飞行员培养方案。
目前的飞行环境检测方法主要基于气象卫星、无人机[4]、机场雷达[5,6]采集环境信息，监测成本高，缺乏针对性。目前绝大多数民航飞机上都装有qar和众多传感器，其自身一直在实时采集传感器数据。
本发明借用飞机已有的传感器和qar采集环境数据，使用成本比无人机更低。本发明监测飞机飞行途中所遇环境，可应用于起飞着陆跑道环境，航路环境，为后来的飞机提供参考。监测范围缩小到跑道，比机场气象雷达更有针对性；对机场设备不做要求，更具有普适性。不仅如此，本发明率先引入信息熵[7]来衡量环境的混乱度，以环境熵作为一段航路中环境混乱度的定量度量，这是其他环境检测方法所忽视的测量指标。
[1]kimjh,lyouj.enhancedqarflightdataencodinganddecodingalgorithmforcivilaircraft[c]sice-icase,2006.internationaljointconference.ieee,2007:5169-5173.
[2]walkerg.redefiningtheincidentstolearnfrom:safetyscienceinsightsacquiredonthejourneyfromblackboxestoflightdatamonitoring[j].safetyscience,2017.
[3]mcfaddenkl,toweller.aviationhumanfactors:aframeworkforthenewmillennium[j].journalofairtransportmanagement,1999,5(4):177-184.
[4]甘肃中茂环保科技有限公司.小型飞行环境监测仪:中国,cn204264448u[p].2015-04-15.
[5]呼曦,魏航科.飞行环境监视系统综合监控方法及声光报警装置:中国,cn103019216a[p].2013-04-03.
[6]t·豪尔德森.用于在飞行器上测量空气湍流且还用于机场和风力发电厂的方法和激光雷达系统:德国,cn101268384[p].2008-09-17.
[7]shannonce.amathematicaltheoryofcommunication[j].bellsystemtechnicaljournal,1948,27(3):379-423.
[8]sanderj,esterm,kriegelhp,etal.density-basedclusteringinspatialdatabases:thealgorithmgdbscananditsapplications[j].dataminingandknowledgediscovery,1998,2(2):169-194.
技术实现要素：
本发明的技术解决问题：克服现有技术的不足，解决了目前卫星、雷达气象报告中缺乏航班针对性的问题，提供一种飞机飞行过程环境的检测方法，能够实时评价每个航班飞行环境，并且能够准确实时地反应飞行航线上的环境情况。
本发明的技术解决方案：一种飞机飞行过程环境的检测方法，包括环境熵计算方法、环境编码方法及其使用方法。(1)环境熵为飞机在飞行过程中所遇环境的无序性的一种度量指标，以相关qar(quickaccessrecorder)数据为输入，对飞机飞行过程中所遇环境进行编码，最后计算环境编码的信息熵。(2)环境编码是对飞机在飞行过程中所遇环境的编码，不同的环境被赋予不同的编码。(3)环境熵表明了飞机飞行过程中所遇环境的混乱程度，同时也表明了该环境的飞行操作难度，可用来评价飞机途径该环境时的飞行操作难度。本发明填补了业内定量检测飞行过程环境及驾驶操作方法的空白，为执行foqa(flightoperationsqualityassurance，飞行操作品质监控)提供了检测指标和方法。
具体步骤如下：
步骤一：从快速存储记录器qar设备中采集飞机在飞行过程中每一时间段的环境数据，所述环境数据包括以下时序参数：航向、风向、风速、气温，将所有环境数据标准化，得到环境数据序列；
步骤二：按时间段遍历各个环境数据序列，在当任意环境数据序列的累积变化幅度大于环境状态划分阈值时，将环境数据序列划分入不同的环境状态，并以该时间段内各环境参数的统计特征表示该环境状态；
步骤三：对飞机在某飞行过程中所有时间段内的环境状态聚类，得到各环境状态的类别标记；
步骤四：以飞机在飞行过程中所遇的环境状态的类别标记的顺序序列作为飞机在该飞行过程中的环境编码；
步骤五：根据环境编码出现的频率，计算环境编码的熵值，即环境熵值，环境熵是指飞机在飞行过程中所遇环境的无序性的一种度量指标值，是飞机在飞行过程中所遇环境的混乱程度。
所述的环境数据可根据实际情况简化采集，当该飞行过程中某项环境数据没有明显变化时，则可以忽略该项数据，从而减少环境数据的采集难度，提升计算环境熵值的运算速度。
所述步骤二中所述环境划分阈值设为2。
所述步骤二中的时间段内各环境参数的统计特征取均值、中位数、众数，优选取均值。
所述步骤三中的聚类选用基于密度的聚类方法，所述聚类方法选择为dbscan，同时选择环境划分阈值的1/2作为dbscan算法[8]所需的聚类半径。
所述步骤五具体实现如下：统计各环境状态类别标记出现的频率，如标记i出现的频率为pi，计算得到环境编码的熵值s：
s＝-∑ipilogpi。
环境熵为0时，说明环境稳定；当环境极度不稳定时，环境熵可达到4甚至5，说明该飞行过程中的环境状态非常混乱，飞行员难以保持稳定的飞行。
本发明与现有技术相比的优点在于：本发明首先提出了环境熵这一新概念，用于描述某架飞机飞行过程的大气环境混乱程度；并给出了环境熵的计算方法，计算数据主要来自于qar，具有普适性、标准性和可计算性，由于使用的是每架飞机自身的传感器测量参数，因此能够定量描述每个航班的航线环境情况。环境稳定时环境熵为0；当环境极度不稳定时，环境熵可达到4甚至5，说明该飞行过程中的环境状态非常混乱，飞行员难以保持稳定的飞行。
附图说明
图1为本发明一种飞机飞行过程环境的检测方法流程图。
具体实施方式
下面结合附图及实施例对本发明进行详细说明。
如图1所示，本发明具体实现如下：
(1)初始化聚类半径r，环境状态划分阈值a，初始化上一状态区间结束时刻last为0,状态数组m为空；
(2)从qar设备读取所需检测时段内的环境数据序列，包括航向、风向、风速、气温，并将以上输入的序列标准化为0均值、单位标准差的序列，每种参数的标准化处理公式为：
其中，input为输入序列，mean为求均值，std为求标准差；
(3)遍历各参数标准化后的序列,若该序列包含n种环境参数，且每种参数的记录时长为t秒，将该序列存储在t行n列的二维数组x中，当
(4)采用dbscan算法，以聚类半径r对状态数组m聚类,相近的环境状态被聚为一类，同一类环境持有同样的标记,即第0类环境状态标记0，第1类环境状态标记1，以此类推，第i类环境状态标记为i，以这些标记作为飞机在飞行过程中所遇到环境的环境编码。
(5)统计每类环境编码的出现频率pi，计算得到环境编码的熵值s：
s＝-∑ipilogpi。
当计算结果，环境熵为0时，说明环境稳定；环境熵可达到4甚至5时，环境极度不稳定时，说明该飞行过程中的环境状态非常混乱，飞行员难以保持稳定的飞行。
提供以上实施例仅仅是为了描述本发明的目的，而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改，均应涵盖在本发明的范围之内。