2021, Vol. 44
2. 中国民航大学 飞行技术学院, 天津 300300;
3. 中国民航大学 中欧航空工程学院, 天津 300300
2. The Flight Technology College, Civil Aviation University of China, Tianjin 300300, P. R. China;
3. Sino-European institute of Aviation Engineering, Civil Aviation University of China, Tianjin 300300, P. R. China
民航历年发生的空中交通管理不安全事件中，管制员人因差错导致的不安全事件数至少占总数的80%[1]。可见，对管制过程开展人因研究对于提升民航运行安全水平具有重要意义。
人机交互分析是人因研究的重要方法。人的感觉器官(如眼睛和耳朵)接受信息，传到大脑进行处理、判断和决策，然后大脑发出的指令传递到口、手、脚等部分并执行[2]，管制指挥的关键在于准确及时地获取空域内运行航空器的信息。而管制员主要通过眼睛获取航空器的信息[3]，因此注视时间和扫视路径等眼动指标能在一定程度上反映被试选择信息的倾向性，目前眼动追踪技术已经能准确记录眼球运动数据。在航空领域，视觉搜索策略[4-6]、工作负荷[7-10]与管制人因差错[11-14]3个方面是国外国内相关研究的主要方面。此外，眼动追踪技术正成为新兴管制技术的重要安全评估手段，如远程塔台界面评估[15]等。综上可知，利用眼动追踪技术捕捉眼球运动信息已经成为国内外研究人员分析管制人机交互的重要手段，通过处理眼动数据，可以提取有意义的眼球运动指标，这些量化指标往往能准确反映管制员的信息搜索和信息获取过程。
雷达界面中航空器目标往往是管制员重点关注对象，这些注视点集中的区域也经常被划分为兴趣区来研究，因此兴趣区提取可以作为一种雷达目标标定的方法。目前兴趣区划分方法主要利用热图法手动划分，该方法虽然快捷简单，但是手动划分一方面缺乏切割兴趣区的理论依据与数据支撑，二是这种粗略的划分方法缺乏一定的精确性。而聚类算法因具有理论数据支撑和一定的准确度而进入人们的视野。但在使用传统的具有噪声的基于密度的聚类算法(DBSCAN，density-based spatial clustering of applications with noise)处理眼动数据时，输入参数不合理会导致聚类效果不理想[16]。因此有必要改进传统DBSCAN算法以避免在处理眼动数据时的问题。
为了提高管制人机交互的效率与安全水平，有必要提出一种能客观准确提取管制雷达界面航空器的目标方法，本研究从眼动数据处理入手进行雷达界面航空器目标标定。利用雷达管制模拟机、眼动仪等搭建模拟管制实验平台，招募8名民航管制学员为实验被试，记录被试在实验过程中的眼动数据，从原始眼动数据中提取注视点信息，对传统DBSCAN算法进行自适应获取参数改进，最后对相关改进进行验证。
1 模拟管制实验 1.1 实验被试总计8名航空管制学员参加本实验，年龄范围为20~25岁。所有被试在实验前均参与过相关管制培训课程，具备独立完成实验任务的能力。所有管制学员视力正常2.0±0.3。
1.2 实验平台实验平台硬件部分包括：
1) Tobii X2-60眼动仪：采样率60 Hz，精度0.34°；
2) 雷达管制模拟机；
3) 14英寸笔记本电脑：分辨率为1 366×768。
1.3 实验程序具体实验程序如下：
1) 连接好所有设备的电源，并将眼动仪与外界设备相连，待连接完成后开启眼动仪，让被试处于设备屏幕前方约0.6 m位置；
2) 待眼动仪开启后，先校准眼动仪，保证仪器能准确记录眼睛观察的位置；
3) 启动雷达管制模拟机，设定初始航空器数目为7；
4) 截取模拟实验过程中某一时刻的管制界面，让被试观察此时的雷达界面，并读取界面中的航空器信息，然后口头报告给主试，直到读取完所有航空器目标的数据，要求被试不能改变标牌的位置或形状，记录每名被试实验过程中的眼动数据，最后换下一被试看同样的雷达界面；
5) 关闭并整理所有实验器材。
1.4 实验数据预处理眼动仪无法记录被试眨眼或看向实验屏幕外时眼睛运动的数据，此时收集到的原始数据是空白数据，数据预处理时需将其删除。
一般认为当人的眼睛停留在某一位置的时间超过100 ms时会形成一次注视[17]。实验使用眼动仪的采样率为60 Hz，即1 s内记录60次数据。因此当眼动仪记录的数据点连续6次未改变位置时，表明眼睛至少在这一位置停留了100 ms，可以认为发生了一次注视过程，此过程眼睛对准的位置称为注视点。这是一种从眼动数据中提取注视点的方法。
2 基于DBSCAN算法的雷达目标标定 2.1 DBSCAN算法简介管制员在进行空中交通指挥时通过阅读飞机的标牌来获得相关的飞行参数，所以管制员的注视行为能直观反映出其关注的目标。雷达中的航空器目标往往受到相对频繁地关注，因此眼动数据能被用来对雷达界面中航空器目标进行标定。
图 1左图显示，雷达界面中有航班号分别为CCA1610、CCA1402和CCA937的3个航空器目标。图 1右图为所有被试在雷达界面上的眼动数据热点图，对照2图可以发现航空器附近的颜色较亮，说明注视点在航空器目标附近集中，根据热点图手动画出3个矩形框，得到3个航空器在雷达界面上的大概位置，显然这种划分的精确性较低，因此需要一种具有较高精确性的方法来标定航空器目标。
|
图 1 模拟雷达界面与注视点热图 Fig. 1 Radar simulation interface and gaze point heat map |
热点图表明在航空器周围分布的注视点密度大，而其他区域注视点密度小，因此可用基于密度的聚类算法进行处理，眼动数据中不可避免地会有一些由于头部运动等行为导致的噪声数据，因此在选择数据处理方法时需要考虑密度和噪声干扰，而DBSCAN算法能根据密度把数据划分为不同的簇，同时能在具有噪声的数据集中进行聚类[18]，因此引入DBSCAN算法处理眼动数据。
DBSCAN算法有2个重要的参数，ε邻域值和密度阈值(Minpts)。对于给定的数据点Q，到Q的距离不大于ε值的区域称为该点的ε邻域，单位是像素。密度阈值一般是一个固定值。DBSCAN算法中，首先统计所有数据点在给定ε值范围内其他数据点的个数，如果其他数据点个数大于给定的密度阈值，即满足建立数据簇的条件，然后遍历所有的数据簇，将密度可达数据簇合并，当搜索不到新的可达数据簇时，结束迭代并建立类别[19]。
DBSCAN算法取不同参数的计算结果对比如图 2，图中以不同形状区分不同类别的数据点，其中的黑空心圆点为噪声点。图 2表明，设定的ε邻域数值或者密度阈值都不能合理划分航空器目标和飞机标牌。原因可能是上述2个参数对DBSCAN算法的聚类结果影响较大，如果输入过大的ε邻域值，无法有效过滤噪声点，更容易建立数据簇，从而降低算法的精确性；输入过小的ε邻域值，很多有效数据点被错误地分类为噪声点，降低算法准确性，导致目标区域标定过小；而密度阈值的选取也为全局参数，但不同注视点数据簇的分布密度差异较大，设定单一的密度阈值处理眼动数据并没有取得理想的效果。
|
图 2 用取不同参数的DBSCAN算法处理CCA1610目标附近眼动数据的结果 Fig. 2 The results of using DBSCAN algorithm with different parameters to process the eye movement data near the target of CCA1610 |
因此，为了避免不合理设定参数导致错误的结果，提高DBSCAN算法的准确性与精确性，将从自适应获取ε邻域值和变密度阈值设定等方面对传统DBSCAN算法进行改进。
2.2 DBSCAN算法改进 2.2.1 自适应获取ε邻域值传统DBSCAN算法需要用户合理设定参数，图 2聚类效果对比图表明参数值的变化对聚类效果的影响很大，如果用户随机输入参数很可能使不符合预期结果。因此，研究将K-最临近距离分布图与统计分析相结合来自适应获取ε邻域值。
在选取DBSCAN算法推荐的邻域值时，Ester M[20]曲线拐点通过设定密度阈值为4来对比选取ε值。因此为了得到推荐邻域值，设定K-最临近距离分布法中的K=4，即对眼动数据集的一个数据点Q，找出距离Q第4近的数据点P，计算Q点与P点之间的距离，访问所有数据点计算出距离，将距离值放入集合DistMinpts中，然后按距离值大小将数据点升序排序。最后画出DistMinpts排序图，如图 3，数据点序号作为横轴，序号对应的距离值为纵轴。一般将排序图中曲线拐点的纵坐标设定为推荐的ε值[21]。图 3箭头对准的位置是曲线的拐点，其纵坐标值约为12.5 px，所以推荐ε值为12.5 px。
|
图 3 DistMinpts排序图 Fig. 3 Sorting chart of DistMinpts |
DistMinpts排序图推荐值是为达到给定密度阈值，大多数的数据点需要满足的最小ε值。应用概率论中分布概率密度函数的特点，如果能拟合出DistMinpts距离分布的概率密度函数，计算出概率密度函数的极值点，并令ε等于极值点能够达到同样的效果[22]。
利用不同随机分布的概率密度函数对DistMinpts分布进行拟合，结果如图 4。图中虚线曲线、实线曲线分别为正态分布的概率密度曲线以及含有尺度参数和位置参数的t分布(下文简称为含参t分布)概率密度曲线。图 4表明，含参t分布的曲线与DistMinpts的分布情况最接近。同时含参t分布的概率密度函数存在唯一的极值点，结合眼动数据，采用最大似然法对参数求解，算得极值点为7.8，因此对DistMinpts分布进行拟合确定的邻域值为7.8。
|
图 4 分布的概率密度函数与拟合曲线 Fig. 4 Probability density function and fitting curves of DistMinpts distribution |
此数值小于DistMinpts排序图法的推荐值，如果直接利用拟合法获得的邻域值会导致大量的数据点被误删，为保证后续聚类过程的效果，需要将拟合法求出的极值进行放大，放大系数m应取1.6，为作图法推荐值与拟合函数极值点的比值。
2.2.2 设定变密度阈值人在观察图片时对不同对象的观察深度是不同的，其根源是不同对象所拥有的信息量和信息重要度迥异，因此不同对象周围的注视点密度可能相差很大。由于DBSCAN算法使用了固定的全局参数[16]，此算法对密度差别较大的数据集(如眼动数据集)聚类效果较差，图 2聚类结果表明，设定全局的密度阈值来处理眼动数据可能会产生不理想的分类效果。因此为了提高该算法处理眼动数据的可靠性与精确性，改进算法还应能自适应地改变密度阈值。
首先，根据前文描述的方法获取邻域值后，统计眼动数据集D中所有数据点的ε邻域内存在的其他数据点个数N。将数据点按照对应的N值大小进行降序排序。显然，这时该数据集中序号越靠前的数据点所在数据簇中注视点密度越大。
然后计算各数据点对应N值的均值N，将对应的N值不小于N的数据点都放入数据集D1，以N作为密度阈值进行聚类。剩余的数据点放入数据集D2，重新计算D2中数据点N值的平均值N1，重复上述过程。由于确定ε邻域值时令K=4，聚类过程结束条件设定为：如果第i次重新计算后的数据点的均值Ni小于4，结束聚类过程。最终输出聚类结果。
2.2.3 改进算法整体流程为了提高DBSCAN算法处理眼动数据的可靠性和准确性，对DBSCAN算法进行了改进，即根据数据集自适应获取ε邻域值和变密度阈值设计。对于眼动数据集D，使用改进算法进行处理，最终输出是带有类别标签的聚类分析结果，改进算法整体流程如图 5。
|
图 5 改进DBSCAN算法流程图 Fig. 5 Flow chart for the improved DBSCAN algorithm |
首先是对自适应获取ε邻域方法的适用性验证。在改进的DBSCAN算法中，由于拟合分布概率函数极值ε0与作图法推荐的ε邻域值存在较大的差异，于是设定了一个放大系数，当对其他航空器目标采用一样的放大处理时，需要对这种放大处理和放大系数的适用性进行检验。
利用8名被试的全部眼动数据做出热点图，在热图中注视点密度较大区域对未标定的另外2架航空器进行粗标。截取出粗标后的眼动数据，处理得到注视点，先利用作图法求出推荐邻域值，再进行分布拟合求出分布函数极值，计算结果对比图如图 6。
|
图 6 改进方法的ε邻域估计值 Fig. 6 Estimated ε value of the improved method |
图 6中的修正邻域值是将分布函数极值放大1.6倍后的结果，即是另外2个航空器目标采用自适应法得到的ε邻域值，与作图法推荐值相差分别为8.6%和10%，差异较小，表明对于其他位置的航空器目标，放大系数值是合理的。
最后是对改进聚类算法的标定结果进行检验，改进算法应能够准确划分航空器目标和航空器标牌。在上面检验中确定的ε值的基础上，基于MATLAB工具使用改进算法对目标区域的眼动数据进行聚类，聚类结果如图 7(a)所示。
|
图 7 眼动数据聚类结果和目标区域提取结果 Fig. 7 Eye movement data clustering results and target region extraction results |
图 7(a)表明改进后的算法能准确划分航空器目标和航空器标牌，并且对雷达界面上其他区域的航空器目标依然能很好地完成聚类。
在3架航空器目标聚类结果的基础上，删除掉噪声点，并使用最小包络矩形法提取出航空器的目标区域，实现对不同航空器目标的标定。提取效果如图 7(b)，图中左边方框中是航空器区域，右边方框是标牌区域。显然，提出的目标区域标定方法较热图法更为精确，且对于雷达界面中不同的航空器目标都具有良好的适用性。
3 结论眼动数据的热点图直观反映出航空器目标附近的注视点密度较大，表明航空器目标是雷达界面中管制员重点关注的目标，在热点图上能粗略标定航空器目标，但仍需进一步处理眼动数据以提高目标提取精度。
研究引入DBSCAN算法并进行改进，提出了一种基于眼动数据的雷达目标标定方法。改进算法能自动处理眼动数据。通过对模拟雷达管制界面中其他目标使用改进算法处理，发现对于不同的航空器，自适应算法获取的ε值与推荐的ε值相差仅为8.6%和10%，证明采用拟合分布函数极值点结合放大系数确定ε值的改进措施是合理的。聚类和提取结果表明，研究提出的雷达目标标定方法，对雷达管制界面中不同航空器目标具有一定的准确性和普适性。
| [1] |
王永刚, 董保健. 管制员人为差错影响因素及指标权重分析[J]. 中国安全生产科学技术, 2011, 7(2): 28-33. Wang Y G, Dong B J. Analysis of human error influence factors and evaluation indexes weights determination for air traffic controller[J]. Journal of Safety Science and Technology, 2011, 7(2): 28-33. (in Chinese) |
| [2] |
刘亚威. 管制疲劳的眼动指标研究[D]. 天津: 中国民航大学, 2018. Liu Y W. The research of eye movement index for detecting air traffic controllers' fatigue[D]. Tianjin: Civil Aviation University of China, 2018. (in Chinese) |
| [3] |
李斌. 随机场景人眼检测实时追踪及其应用技术研究[D]. 北京: 中国科学院大学(中国科学院西安光学精密机械研究所), 2018. Li B. Researches on real time eye tracking in random scenes and application technology[D]. Beijing: University of Chinese Academy of Sciences, 2018. (in Chinese) |
| [4] |
McClung S N, Kang Z. Characterization of visual scanning patterns in air traffic control[J]. Comput Intell Neurosci, 2016, 2016: 8343842. |
| [5] |
Gal Z. Gaze behavior and visual attention: a review of eye tracking studies in aviation[J]. The International Journal of Aviation Psychology, 2016, 26(3/4): 75-104. |
| [6] |
Manske P G, Schier S L. Visual scanning in an air traffic control tower-A simulation study[J]. Procedia Manufacturing, 2015, 3: 3274-3279. DOI:10.1016/j.promfg.2015.07.397 |
| [7] |
Peißl S, Wickens C D, Baruah R. Eye-tracking measures in aviation: a selective literature review[J]. The International Journal of Aerospace Psychology, 2018, 28(3/4): 98-112. |
| [8] |
靳慧斌, 洪远, 蔡亚敏. 基于交互指标的空中交通管制员工作负荷实时测量方法研究[J]. 安全与环境工程, 2015, 22(3): 147-150, 157. Jin H B, Hong Y, Cai Y M. Research on real-time measuring method of ATC workload based on interaction indexes[J]. Safety and Environmental Engineering, 2015, 22(3): 147-150, 157. (in Chinese) |
| [9] |
Rodríguez S, Sánchez L, López P, et al. Pupillometry to assess air traffic controller workload through the mental workload model[C]//Proceedings of the 5th International Conference on Application and Theory of Automation in Command and Control Systems -ATACCS'15. September 30-October 2, 2015. Toulouse, France. New York: ACM Press, 2015: 95-104.
|
| [10] |
杨新湦, 王茜. 眼动指标、管制负荷及航空器动态特征之间的相关性分析[J]. 科学技术与工程, 2018, 18(15): 333-340. Yang X S, Wang Q. The correlation analysis of eye movement indicators, controler's workload and aircraft dynamics[J]. Science Technology and Engineering, 2018, 18(15): 333-340. (in Chinese) DOI:10.3969/j.issn.1671-1815.2018.15.052 |
| [11] |
Edwards T, Homola J, Mercer J, et al. Multifactor interactions and the air traffic controller: The interaction of situation awareness and workload in association with automation[J]. IFAC-Papers OnLine, 2016, 49(19): 597-602. DOI:10.1016/j.ifacol.2016.10.623 |
| [12] |
Giraudet L, Imbert J P, Tremblay S, et al. High rate of inattentional deafness in simulated air traffic control tasks[J]. Procedia Manufacturing, 2015, 3: 5169-5175. DOI:10.1016/j.promfg.2015.07.555 |
| [13] |
di Stasi L L, McCamy M B, Catena A, et al. Microsaccade and drift dynamics reflect mental fatigue[J]. European Journal of Neuroscience, 2013, 38(3): 2389-2398. DOI:10.1111/ejn.12248 |
| [14] |
Friedrich M, Biermann M, Gontar P, et al. The influence of task load on situation awareness and control strategy in the ATC tower environment[J]. Cognition, Technology & Work, 2018, 20(2): 205-217. |
| [15] |
Li W C, Kearney P, Braithwaite G, et al. How much is too much on monitoring tasks? Visual scan patterns of single air traffic controller performing multiple remote tower operations[J]. International Journal of Industrial Ergonomics, 2018, 67: 135-144. DOI:10.1016/j.ergon.2018.05.005 |
| [16] |
冯振华, 钱雪忠, 赵娜娜. Greedy DBSCAN: 一种针对多密度聚类的DBSCAN改进算法[J]. 计算机应用研究, 2016, 33(9): 2693-2696, 2700. Feng Z H, Qian X Z, Zhao N N. Greedy DBSCAN: an improved DBSCAN algorithm on multi-density clustering[J]. Application Research of Computers, 2016, 33(9): 2693-2696, 2700. (in Chinese) |
| [17] |
马勇. 基于眼动分析的汽车驾驶员视觉搜索模式研究[D]. 西安: 长安大学, 2006. Ma Y. Study on driver' visual search pattern based on analysis of eye movements[D]. Xi'an: Changan University, 2006. (in Chinese) |
| [18] |
Wang L M, Li M Y, Han X M, et al. An improved density-based spatial clustering of application with noise[J]. International Journal of Computers and Applications, 2018, 40(3): 1-7. |
| [19] |
Mahesh Kumar K, Rama Mohan Reddy A. A fast DBSCAN clustering algorithm by accelerating neighbor searching using Groups method[J]. Pattern Recognition, 2016, 58: 39-48. DOI:10.1016/j.patcog.2016.03.008 |
| [20] |
Ester M, Kriegel H P, Sander J, et al. A density-based algorithm for discovering clusters in large spatial databases with noise[C]//Proceedings of the Second International Conference on Knowledge Discovery and Data Mining. USA: AAAI Press, 1996: 226-231.
|
| [21] |
夏鲁宁, 荆继武. SA-DBSCAN: 一种自适应基于密度聚类算法[J]. 中国科学院研究生院学报, 2009, 26(4): 530-538. Xia L N, Jing J W. SA-DBSCAN: a self-adaptive density-based clustering algorithm[J]. Journal of the Graduate School of the Chinese Academy of Sciences, 2009, 26(4): 530-538. (in Chinese) |
| [22] |
霍百明. 基于眼动追踪的雷达动态目标标定方法研究[D]. 天津: 中国民航大学, 2019. Huo B M. Research of radar dynamic target calibration method based on eye tracking[D]. Tianjin: Civil Aviation University of China, 2019. (in Chinese) |