CN111178706A - 一种基于眼动数据的管制员注意力分配评估方法 - Google Patents
一种基于眼动数据的管制员注意力分配评估方法 Download PDFInfo
- Publication number
- CN111178706A CN111178706A CN201911307937.2A CN201911307937A CN111178706A CN 111178706 A CN111178706 A CN 111178706A CN 201911307937 A CN201911307937 A CN 201911307937A CN 111178706 A CN111178706 A CN 111178706A
- Authority
- CN
- China
- Prior art keywords
- eye movement
- duration
- coordinates
- point
- controller
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000004424 eye movement Effects 0.000 title claims abstract description 62
- 238000011156 evaluation Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000008859 change Effects 0.000 claims abstract description 24
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000001035 drying Methods 0.000 claims abstract description 11
- 238000007619 statistical method Methods 0.000 claims abstract description 3
- 239000006185 dispersion Substances 0.000 claims description 14
- 230000004434 saccadic eye movement Effects 0.000 claims description 14
- 230000001186 cumulative effect Effects 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 210000001747 pupil Anatomy 0.000 claims description 3
- 230000004397 blinking Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000006399 behavior Effects 0.000 description 6
- 230000007423 decrease Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000004962 physiological condition Effects 0.000 description 1
- 230000001711 saccadic effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
- G06F3/011—Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
- G06F3/013—Eye tracking input arrangements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Human Computer Interaction (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明公开了一种基于眼动数据的管制员注意力分配评估方法,该方法包括:采集管制员在执行管制过程中不同航空器数量场景中的眼动数据并进行去燥处理;利用眼动算法区分去燥后的眼动数据;选择反映管制员注意力分配的眼动指标;利用统计分析和多重分形去趋势波算法分析选取的眼动指标的变化特征对管制员注意力分配进行评估。本发明通过管制员在面对管制场景时眼动数据的变化和相关算法进行结合,得出了较优的管制员注意力相关指标,可以有效评估管制员注意力分配特征。
Description
技术领域
本发明涉及一种基于眼动数据的管制员注意力分配评估方法,属于注意力评估技术领域。
背景技术
民用航空运输凭借其安全、快捷、舒适、不受地形限制的优点,在交通运输业中占有独特的地位。随着航空制造业的快速发展,民用航空运输需求的逐年增长,民用航空运输的环境日益复杂,由此造成的空中交通管制任务难度逐渐增加。管制员作为管制扇区的管理者,其行为直接影响着整个空中交通系统的运行。空中交通流量的增加和管制场景的复杂度会增加管制员的工作难度,增强管制员犯错的概率。近年来,由于管制员工作差错造成的空管严重事故征候事件略有增加,例如,2016年在上海虹桥国际机场发生的跑道入侵事件,塔台管制员遗忘飞机动态造成人为原因严重事故征候。通常,为保证管制工作的准确性和管制的安全性,雷达管制员在监视雷达屏幕时,要合理的分配注意力,避免错、忘、漏现象,快速、准确的发现航空器之间潜在的冲突。因此,评估管制员的注意力分配行为对于空中交通运行安全具有重要的理论意义,是预测管制员获取信息的方式和管制意图的前提,也是实现管制自动化的研究基础。
视觉搜索是管制员快速找到目标、获取航空信息最直接的方式。管制员利用扫视快速掌握管制场景的动态信息、凭借精准的注视获取目标信息。研究表明,眼动行为是心理活动的外在的表现,眼动指标的时空特性能够表征管制员信息搜索过程中的行为,反映管制员的管制情景意识及注意力分配状况。
因此,眼动数据可以反映管制员心理和生理状况。通过分析管制员在管制场景中的各项眼动指标,挖掘各项指标特征,可以深入理解管制员的注意力分配行为,正确评估注意力分配。
发明内容
本发明所要解决的技术问题是:提供一种基于眼动数据的管制员注意力分配评估方法,基于管制场景变化对管制员产生心理-生理变化会表现为眼动数据变化的原理,分析眼动数据波形变化,从而判定不同场景下注意力分配指标的特征变化,为评估注意力分配提供一种新的判定方法。
本发明为解决上述技术问题采用以下技术方案:
一种基于眼动数据的管制员注意力分配评估方法,包括如下步骤:
步骤1,采集管制员在执行管制过程中不同航空器数量场景下的眼动数据,包括瞳孔直径、眨眼频率、视线与屏幕的交点坐标;
步骤2,对步骤1采集的眼动数据进行去燥处理,删除存在缺失或异常的数据;
步骤3,采用眼动算法将经过去燥处理后的数据区分为注视组和扫视组,并计算得到每个注视组的注视中心坐标、注视开始时间和注视持续时间,计算得到每个扫视组的扫视持续时间、扫视开始时间和扫视角速度;
步骤4,从注视组和扫视组中选择反映管制员注意力分配的眼动指标,即注视持续时间、扫视持续时间、扫视角速度,并多重分形去趋势波算法分析选取的眼动指标的变化特征,根据变化特征对管制员注意力分配进行评估。
作为本发明的一种优选方案,步骤1所述采集管制员在执行管制过程中不同航空器数量场景下的眼动数据,是通过faceLAB眼动仪进行采集。
作为本发明的一种优选方案,所述步骤3的具体过程为:
步骤31,选取经过去燥处理后的视线与屏幕的交点坐标作为眼动算法的数据,设置两个空间参数,包括定界距离d1和聚簇距离d2,设置一时间参数,为最小注视持续时间mt,其中,d1=0.01m,sxa、sya分别表示经过去燥处理后的视线与屏幕的交点坐标中所有横、纵坐标的标准差,mt=50ms;
步骤32,将经过去燥处理后的视线与屏幕的交点坐标记为(x1,y1),(x2,y2),…,(xM,yM),设置第一个坐标(x1,y1)为第一个聚类簇的起点和均值点;
步骤33,对于后续每个坐标,计算坐标与当前均值点之间的欧几里得距离,当欧几里得距离<定界距离d1时,则将该坐标归到当前聚类簇中,并对均值点坐标进行更新,新的均值点坐标的横、纵轴分别对应为当前聚类簇中所有坐标的横、纵轴的均值;当欧几里得距离≥定界距离d1时,当前聚类簇聚类完成,转入步骤34;
步骤34,将欧几里得距离≥定界距离d1对应的坐标设置为下一个聚类簇的起点和均值点,对后续每个坐标重复步骤33-34,直至遍历完所有坐标,则完成初步区分;
步骤35,根据初步区分结果,将每个聚类簇视为一个注视组,将簇中所有坐标的横、纵轴均值作为该簇均值点坐标的横、纵轴,计算簇中各坐标与均值点坐标之间的欧几里得距离,将欧几里得距离>聚簇距离d2对应的坐标去除后,计算该簇中心点坐标,中心点坐标的横、纵轴即为去除欧几里得距离>聚簇距离d2对应的坐标后剩余坐标的横、纵轴的均值,用(xf,yf,st_timef,dur_timef)表示每个聚类簇,其中,xf、yf分别为注视组中心点坐标的横、纵轴,st_timef表示注视开始时间,dur_timef表示注视持续时间;
步骤36,若某注视组的注视持续时间<最小注视持续时间mt,则将该注视组转为扫视组,用(sac_times,st_times,ang_velocitys)表示每个扫视组,其中,sac_times表示扫视持续时间,st_times表示扫视开始时间,ang_velocitys表示扫视角速度,区分结束。
作为本发明的一种优选方案,所述步骤4的具体过程为:
步骤41,选择注视持续时间、扫视持续时间、扫视角速度,作为反映管制员注意力分配的眼动指标进行统计分析,对每个指标都进行步骤42-49的操作;
步骤42,以注视持续时间指标为例,计算注视持续时间指标序列的均值,公式为:
步骤43,构造去均值的和序列,确定累积离差,公式为:
其中,Y(i)表示累积离差,i=1,…,N;
步骤44,将Y(i)序列划分成Ns个互不相同的区间,每个区间包含s个数据,则Ns=int(N/s);
步骤45,对每个区间内的每个数据进行多项式拟合:
yv(j)=a1jk+a2jk-1+…+akj+ak+1,j=1,…,s
其中,v=1,…,Ns,yv(j)表示划分的第v个区间中第j个数据的多项式,k=1,…,6,a1、a2、…、ak、ak+1为多项式系数;
步骤46,计算去趋势后的均方误差,公式为:
其中,F2(s,v)表示均方误差,Y[(v-1)s+j]表示累积离差;
步骤47,计算波动函数,公式为:
当q=0时,波动函数为
当q≠0时,波动函数为
其中,Fq(s)表示波动函数,h(q)表示广义Hurst指数;
步骤48,根据波动函数Fq(s)取对数,求得其变化斜率即h(q),从而得到qmax和qmin,计算Δh(q)=h(qmin)-h(qmax),根据h(q)计算奇异指数α和多重分形谱宽f(α);
步骤49,当Δh(q)>0.5时记1分,否则记0分;当α<1时记1分,否则记0分;当f(α)>1.3时记1分,否则记0分;
步骤410,统计注视持续时间指标序列、扫视持续时间指标序列、扫视角速度指标序列对应的上述三个参数的分值,得到总分值,总分值在0-9分之间。
作为本发明的一种优选方案,步骤43所述累积离差的修正公式为:
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1、本发明首次提出以管制员眼动数据的注意力分配评估方法,基于管制场景变化对管制员产生心理-生理变化会表现为眼动数据变化的原理,分析眼动数据波形变化,从而判定不同场景下注意力分配指标的特征变化,为评估注意力分配提供一种新的判定方法。
2、本发明通过管制员在面对管制场景时眼动数据的变化和相关算法进行结合,得出了较优的管制员注意力相关指标,可以有效评估管制员注意力分配特征。
附图说明
图1是本发明一种基于眼动数据的管制员注意力分配评估方法的流程图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
如图1所示,一种基于眼动数据的管制员注意力分配评估方法,包括如下步骤:
S1,采集管制员在执行管制过程中不同航空器数量场景中的眼动数据;
使用faceLAB眼动仪采集管制员眼动数据,包括瞳孔直径、眨眼频率、视线与屏幕的交点坐标。
S2,对眼动数据进行去燥处理;
对采集到的眼动数据进行去噪处理,删除缺少相关指标的数据和异常数据。
S3,利用眼动算法区分去燥后的眼动数据;
设置空间约束和时间约束,即计算过程设置三个基本参数:两个空间参数(定界距离d1、聚簇距离d2)和一个时间参数(最小注视持续时间mt)。
根据去燥后的眼动数据,从中按时间提取视线与屏幕的交点坐标序列,记为(x1,y1),(x2,y2),…,(xM,yM),其中,(x1,y1)、(x2,y2)、…、(xM,yM)是视线与屏幕的交点坐标,将视线与屏幕的交点坐标序列作为眼动算法的原始数据。
S31:将去燥后的眼动数据第一个记录点(x1,y1)设置为第一个聚类簇的起点及该簇的均值点坐标(mean_x1,mean_y1);
S32:计算后续记录点与均值点之间的欧几里德距离dist:
若dist<d1,则将该点归为当前的聚类簇中,并重新计算均值点的坐标:
式中(x1t,y1t)为当前聚类簇中第t个点的坐标,T为当前的聚类簇中点的个数;
S33:重复S32的计算过程,直到dist≥d1,则聚类簇聚类完成,并将该簇定义为一个注视组,将当前的均值点坐标定义为注视组的均值点坐标;将当前的记录点设置为下一个聚类簇的起点,并作为该簇新的均值点坐标;
S34:重复S32、S33的计算过程,直到遍历完成该段时间序列数据,则初步的区分结果已经完成;
S35:根据S34的初步区分结果,计算每一个注视组中各点与该组均值点之间的欧几里德距离,如果某点与均值点的距离大于d2,则该点不用于计算所在注视组的中心点坐标,中心点坐标为该注视组所有点坐标的均值。根据以上的筛选结果,用(xf,yf,st_timef,dur_timef)四个参数表示每个注视组,其中xf,yf表示注视中心的坐标,st_timef表示注视开始时间,dur_timef表示注视持续时间;
S36:若dur_timef<mt,则该数据簇不作为注视组;
S37:两注视组之间的数据参数视为扫视组,根据注视组的区分结果,用(sac_times,st_times,ang_velocitys)三个参数表示扫视组,其中sac_times表示扫视持续时间,st_times表示扫视开始时间,ang_velocitys表示扫视角速度。
S4,选择反映管制员注意力分配的眼动指标,利用多重分形去趋势波算法分析选取的眼动指标的变化特征对管制员注意力分配进行评估。
选择反映管制员注意力分配的眼动数据包括注视持续时间、扫视持续时间和扫视角速度。
计算选取数据的均值,即
其中:N为时间序列的长度;xk为时间序列;
构造去均值的和序列,确定其累积离差,即
分割数据,即
将Y(i)序列划分成Ns个互不相同的区间,每个区间含有s个数据,其中Ns=int(N/s);
多项式拟合,即
yv(j)=a1jk+a2jk-1+…+akj+ak+1j=1,…,s,k=1,…,6
其中:v(v=1,…,Ns)为划分的每个小区间;k为阶拟合多项式的阶数;a1、a2、…、ak、ak+1为各项系数;
计算去趋势后的均方误差,即
计算去趋势后的波动函数,即
当q=0时,波动函数为
当q≠0时,波动函数为
其中:Fq(s)是序列数据长度s和阶数q的函数,表示不同的q描述信号不同波动程度对Fq(s)的影响,Fq(s)随着s的增大呈幂律增长;根据上述两个公式可得Fq(s)、s与h(q)的计算其关系为Fq(s)∝sh(q)
其中:h(q)称为广义Hurst指数,即为logFq(s)~log函数式中的变化斜率;
修正MF-DFA算法,即
当h(q)趋近于0时需要修正MF-DFA算法,需要在MF-DFA算法过程的第一步求两次累积离差:
重复步骤分割数据、多项式拟合、计算均方误差和计算去趋势后的均方误差值和波动函数,计算广义波动函数,即
对于大标度有
计算广义Hurst指数和多重分形理论中的标度指数,即
τ(q)=qh(q)-1
计算α和多重分形谱,即
f(α)=q[α-h(q)]+1
其中:上述两个公式根据勒让德变换α=τ′(q)和f(α)=qα-τ(q)得到。
在本实施例中,分析各个指标的多重分形特征时选取航空器数量分别为4、8、12、16时的各个指标信号数据,MF-DFA方法中每个场景中波动函数取值范围为[-20,20]。
表1不同数量下注视持续时间MF-DFA参数
当q从-20增加到+20时,注视持续时间的广义Hurst指数h(q)曲线均依赖q随之下降,不为常数,说明注视持续时间存在明显的多重分形特征,且其波动函数符合幂律分布。表1所示,随着航空器数量的增加,Δh(q)值越来越大,说明注视持续时间的波形随着数量的增加多重分形越来越强,当航空器数量增加至16时,且整体呈现长程正相关,该变化并不是一个随机过程,过去的变化特性影响现在和未来序列的特性变化。多重分形谱宽f(α)值逐渐增加,而分形维数差Δf逐渐减小,且Δf<0,说明注视持续时间局部变化越强烈,波动幅度越大,振动信号波谷数目多于波峰数目,且随着数量的增加分布复杂性增加,多重分形越显著。
表2不同数量下扫视角速度MF-DFA参数
随着阶数q的增加,h(q)逐渐减小,并不是常数,说明扫视角速度序列的波动函数符合幂律分布。结合表2所示,随着航空器数量的增加,Δh(q)逐渐减小,说明扫视角速度的多重分形越来越弱;表明扫视角速度序列总体呈现长程正相关。因此,可以表明扫视角速度在一定时间段内会继续保持目前的变化状态,随着数量的增加多重分形强度减弱。根据表4.5显示的多重分形谱参数,随着航空器数量的增加,f(α)基本呈现递减趋势,而Δf都小于零,说明扫视角速度信号多处于波谷,随数量的增加,奇异性减弱,信号逐渐平滑,复杂性减小。
表3不同数量下扫视持续时间MF-DFA参数
随着q的增加,广义Hurst指数h(q)逐渐减小,并不是常数,Δh(q)逐渐减小,说明扫视持续时间序列的波动函数符合幂律分布,该序列多重分形特性逐渐减弱。由表3可知,表示扫视持续时间呈整体正相关变化;多重分形谱宽f(α)随航空器数量的增加而减小,说明扫视持续时间信号的概率分布均匀,奇异性减小,波动程度较平缓,多重分形减弱,复杂性减小,且未来一段时间内变化趋势与现在一致。Δf值随着航空器数量的增加而减小,说明扫视持续时间的信号波形波峰多于波谷的数目逐渐减少,当航空器数量增加至16时,波谷数目远大于波峰数目。
(1)Hurst指数
自相似性是指某时间序列从不同尺度来看,其波动具有相似性;长程相关性表明了时间序列信号在某时间间隔内两数据具有统计相关性,反映了信号内部的波动特性。对某眼动数据时间序列:
Δh(q)=h(qmin)-h(qmax),其中h(qmin)表示阶数q为最小值时所对应的Hurst指数,h(qmax)表示阶数q为最大值时所对应的Hurst指数,Δh(q)可用来度量时间序列信号多重分形特征的强弱程度,Δh(q)值越大,表明其多重分形特征性越强。
(2)奇异指数α和多重分形谱宽f(α)
奇异指数α反映了具有分形序列在某局部概率测试分布上的不均匀程度,多重分形谱宽f(α)是奇异指数α的分维函数,形状一般为单峰拱形。由α与f(α)得出描述多重分析谱的参数,这些参数能够精确地表达多重分析时间序列动力学行为:
奇异指数α,α越小说明信号局部奇异性越强,即信号局部变化越强烈,相反,α越大说明该信号局部奇异性越弱,信号越平滑,因而α在一定程度上也反映了信号的局部变化剧烈程度;
多重分形谱宽f(α)=αmax-αmin,f(α)表示了信号多重分形性的强弱,f(α)越大则时间序列的概率测度分布得越不均匀,振动的越剧烈,可以反映眼动及生理信号波动的剧烈程度和分布的复杂性。
根据MF-DFA算法计算的三个眼动指标,每个眼动指标都有以上三个参数Δh(q)、奇异指数α、多重分形谱宽f(α)。对于三个参数都有Δh(q)大于0.5记1分,小于等于0.5记0分;奇异指数α小于1记1分,大于等于1记0分;多重分形谱宽f(α)大于1.3记1分,小于等于1.3记0分。管制员注意力分配评估得分最低0分,最高9分,根据结果最后用得分量化管制员注意力分配情况,得分越高说明管制员注意力分配情况越好。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (5)
1.一种基于眼动数据的管制员注意力分配评估方法,其特征在于,包括如下步骤:
步骤1,采集管制员在执行管制过程中不同航空器数量场景下的眼动数据,包括瞳孔直径、眨眼频率、视线与屏幕的交点坐标;
步骤2,对步骤1采集的眼动数据进行去燥处理,删除存在缺失或异常的数据;
步骤3,采用眼动算法将经过去燥处理后的数据区分为注视组和扫视组,并计算得到每个注视组的注视中心坐标、注视开始时间和注视持续时间,计算得到每个扫视组的扫视持续时间、扫视开始时间和扫视角速度;
步骤4,从注视组和扫视组中选择反映管制员注意力分配的眼动指标,即注视持续时间、扫视持续时间、扫视角速度,并多重分形去趋势波算法分析选取的眼动指标的变化特征,根据变化特征对管制员注意力分配进行评估。
2.根据权利要求1所述基于眼动数据的管制员注意力分配评估方法,其特征在于,步骤1所述采集管制员在执行管制过程中不同航空器数量场景下的眼动数据,是通过faceLAB眼动仪进行采集。
3.根据权利要求1所述基于眼动数据的管制员注意力分配评估方法,其特征在于,所述步骤3的具体过程为:
步骤31,选取经过去燥处理后的视线与屏幕的交点坐标作为眼动算法的数据,设置两个空间参数,包括定界距离d1和聚簇距离d2,设置一时间参数,为最小注视持续时间mt,其中,d1=0.01m,sxa、sya分别表示经过去燥处理后的视线与屏幕的交点坐标中所有横、纵坐标的标准差,mt=50ms;
步骤32,将经过去燥处理后的视线与屏幕的交点坐标记为(x1,y1),(x2,y2),…,(xM,yM),设置第一个坐标(x1,y1)为第一个聚类簇的起点和均值点;
步骤33,对于后续每个坐标,计算坐标与当前均值点之间的欧几里得距离,当欧几里得距离<定界距离d1时,则将该坐标归到当前聚类簇中,并对均值点坐标进行更新,新的均值点坐标的横、纵轴分别对应为当前聚类簇中所有坐标的横、纵轴的均值;当欧几里得距离≥定界距离d1时,当前聚类簇聚类完成,转入步骤34;
步骤34,将欧几里得距离≥定界距离d1对应的坐标设置为下一个聚类簇的起点和均值点,对后续每个坐标重复步骤33-34,直至遍历完所有坐标,则完成初步区分;
步骤35,根据初步区分结果,将每个聚类簇视为一个注视组,将簇中所有坐标的横、纵轴均值作为该簇均值点坐标的横、纵轴,计算簇中各坐标与均值点坐标之间的欧几里得距离,将欧几里得距离>聚簇距离d2对应的坐标去除后,计算该簇中心点坐标,中心点坐标的横、纵轴即为去除欧几里得距离>聚簇距离d2对应的坐标后剩余坐标的横、纵轴的均值,用(xf,yf,st_timef,dur_timef)表示每个聚类簇,其中,xf、yf分别为注视组中心点坐标的横、纵轴,st_timef表示注视开始时间,dur_timef表示注视持续时间;
步骤36,若某注视组的注视持续时间<最小注视持续时间mt,则将该注视组转为扫视组,用(sac_times,st_times,ang_velocitys)表示每个扫视组,其中,sac_times表示扫视持续时间,st_times表示扫视开始时间,ang_velocitys表示扫视角速度,区分结束。
4.根据权利要求1所述基于眼动数据的管制员注意力分配评估方法,其特征在于,所述步骤4的具体过程为:
步骤41,选择注视持续时间、扫视持续时间、扫视角速度,作为反映管制员注意力分配的眼动指标进行统计分析,对每个指标都进行步骤42-49的操作;
步骤42,以注视持续时间指标为例,计算注视持续时间指标序列的均值,公式为:
步骤43,构造去均值的和序列,确定累积离差,公式为:
其中,Y(i)表示累积离差,i=1,…,N;
步骤44,将Y(i)序列划分成Ns个互不相同的区间,每个区间包含s个数据,则Ns=int(N/s);
步骤45,对每个区间内的每个数据进行多项式拟合:
yv(j)=a1jk+a2jk-1+…+akj+ak+1,j=1,…,s
其中,v=1,…,Ns,yv(j)表示划分的第v个区间中第j个数据的多项式,k=1,…,6,a1、a2、…、ak、ak+1为多项式系数;
步骤46,计算去趋势后的均方误差,公式为:
其中,F2(s,v)表示均方误差,Y[(v-1)s+j]表示累积离差;
步骤47,计算波动函数,公式为:
当q=0时,波动函数为
当q≠0时,波动函数为
其中,Fq(s)表示波动函数,h(q)表示广义Hurst指数;
步骤48,根据波动函数Fq(s)取对数,求得其变化斜率即h(q),从而得到qmax和qmin,计算Δh(q)=h(qmin)-h(qmax),根据h(q)计算奇异指数α和多重分形谱宽f(α);
步骤49,当Δh(q)>0.5时记1分,否则记0分;当α<1时记1分,否则记0分;当f(α)>1.3时记1分,否则记0分;
步骤410,统计注视持续时间指标序列、扫视持续时间指标序列、扫视角速度指标序列对应的上述三个参数的分值,得到总分值,总分值在0-9分之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911307937.2A CN111178706A (zh) | 2019-12-18 | 2019-12-18 | 一种基于眼动数据的管制员注意力分配评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911307937.2A CN111178706A (zh) | 2019-12-18 | 2019-12-18 | 一种基于眼动数据的管制员注意力分配评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111178706A true CN111178706A (zh) | 2020-05-19 |
Family
ID=70652050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911307937.2A Pending CN111178706A (zh) | 2019-12-18 | 2019-12-18 | 一种基于眼动数据的管制员注意力分配评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111178706A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112785890A (zh) * | 2020-12-24 | 2021-05-11 | 中国人民解放军空军工程大学 | 一种塔台飞行指挥模拟训练管制员注意力评估系统及其方法 |
CN112949404A (zh) * | 2021-02-01 | 2021-06-11 | 吉林大学 | 一种夜间驾驶人视觉搜索能力量化方法 |
CN113778070A (zh) * | 2020-07-17 | 2021-12-10 | 北京京东振世信息技术有限公司 | 机器人的控制方法和装置 |
CN114648102A (zh) * | 2022-05-24 | 2022-06-21 | 华东交通大学 | 火灾报警方法、系统、可读存储介质及计算机设备 |
CN115844404A (zh) * | 2023-03-01 | 2023-03-28 | 中国民航大学 | 一种基于眼动数据的管制员注意力特征评价方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106251065A (zh) * | 2016-07-28 | 2016-12-21 | 南京航空航天大学 | 一种基于眼动行为指标体系的管制效率评估方法 |
US20180053103A1 (en) * | 2016-08-19 | 2018-02-22 | Massachusetts Institute Of Technology | Systems and methods for providing visual allocation management |
-
2019
- 2019-12-18 CN CN201911307937.2A patent/CN111178706A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106251065A (zh) * | 2016-07-28 | 2016-12-21 | 南京航空航天大学 | 一种基于眼动行为指标体系的管制效率评估方法 |
US20180053103A1 (en) * | 2016-08-19 | 2018-02-22 | Massachusetts Institute Of Technology | Systems and methods for providing visual allocation management |
Non-Patent Citations (3)
Title |
---|
刘银鑫: "基于不同场景的管制员眼动行为分析研究" * |
张美玉,孙轶,王艳军: "航班离场时隙跳变时空特征分析" * |
王艳军,刘银鑫: "空中交通管制专家与新手眼动行为的统计分析研究" * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113778070A (zh) * | 2020-07-17 | 2021-12-10 | 北京京东振世信息技术有限公司 | 机器人的控制方法和装置 |
CN113778070B (zh) * | 2020-07-17 | 2024-04-19 | 北京京东振世信息技术有限公司 | 机器人的控制方法和装置 |
CN112785890A (zh) * | 2020-12-24 | 2021-05-11 | 中国人民解放军空军工程大学 | 一种塔台飞行指挥模拟训练管制员注意力评估系统及其方法 |
CN112949404A (zh) * | 2021-02-01 | 2021-06-11 | 吉林大学 | 一种夜间驾驶人视觉搜索能力量化方法 |
CN112949404B (zh) * | 2021-02-01 | 2022-04-19 | 吉林大学 | 一种夜间驾驶人视觉搜索能力量化方法 |
CN114648102A (zh) * | 2022-05-24 | 2022-06-21 | 华东交通大学 | 火灾报警方法、系统、可读存储介质及计算机设备 |
CN115844404A (zh) * | 2023-03-01 | 2023-03-28 | 中国民航大学 | 一种基于眼动数据的管制员注意力特征评价方法及装置 |
CN115844404B (zh) * | 2023-03-01 | 2023-05-12 | 中国民航大学 | 一种基于眼动数据的管制员注意力特征评价方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111178706A (zh) | 一种基于眼动数据的管制员注意力分配评估方法 | |
CN111951637B (zh) | 一种任务情景相关联的无人机飞行员视觉注意力分配模式提取方法 | |
CN113743471B (zh) | 一种驾驶评估方法及其系统 | |
CN108053110B (zh) | 一种基于pmu数据的变压器状态在线诊断方法 | |
CN104503916A (zh) | 一种系统界面可用性定量评估方法 | |
CN115828140A (zh) | 邻域互信息与随机森林相融合故障检测方法、系统及应用 | |
CN106650250B (zh) | 基于等效评价因子的红外多目标威胁排序方法 | |
CN110619482A (zh) | 基于逻辑回归和单层次分析赋权法的驾驶行为评分方法 | |
CN112215120B (zh) | 视觉搜索区域的确定方法、装置及驾驶模拟器 | |
CN106155876A (zh) | 一种金丝雀分析方法、应用和计算设备 | |
Kuravsky et al. | Assessing the aircraft crew actions with the aid of a human factor risk model | |
CN109343062B (zh) | 一种径向干扰回波和降水回波的识别方法及系统 | |
CN117629122B (zh) | 一种水坝位移的监测与预警方法及系统 | |
CN117406027A (zh) | 一种配电网故障测距方法及其系统 | |
CN113311858A (zh) | 一种基于无人机的飞控系统及方法 | |
CN114469138B (zh) | 基于时频域的脑电爆发抑制模式的检测方法、系统及介质 | |
Xiao et al. | Mental workload prediction based on attentional resource allocation and information processing | |
CN104391599B (zh) | 一种平衡精确度与灵敏度的触屏点过滤方法 | |
CN110928752A (zh) | 空管台站健康度评估方法、装置及设备 | |
US12012219B2 (en) | Aircraft buffet detection | |
CN108986554B (zh) | 一种基于模糊综合评判的空域扇区拥挤程度动态识别方法 | |
CN111091910B (zh) | 基于画钟测试的智能评估系统 | |
CN111966220A (zh) | 一种基于眼动数据的管制员情景意识水平分析方法 | |
CN113125439B (zh) | 一种弹簧端面检测系统及其检测方法 | |
CN115546339B (zh) | 利用景观生态指标模型的城市森林制图方法、装置及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |