本发明属于座舱照明人机功效范畴,涉及一种对飞行器座舱眩光进行评估的方法。
背景技术:
飞行器座舱中由于照明灯或外部光线在座舱风挡玻璃或显示器玻璃屏幕上发生反射进入人眼将产生眩光,影响飞行员对显示器上信息和舱外情况的判断。因此,需要对眩光进行评估,计算眩光位置、判断眩光源。目前通常采用软件仿真或者人为评估的方法对眩光位置以及眩光源进行评估,采用软件仿真的方法建模及分析过程耗时长、工作量大,与实际可能存在偏差,采用人为评估的方法费时费力,存在主观性,并对眩光位置、眩光源缺乏量化评估。
为克服上述困难进而实现座舱眩光位置的快速计算,以及眩光源的判断,需要一种能够实现对飞行器座舱眩光环境进行量化评估的方法。
技术实现要素:
本发明的目的是:提供一种驾驶舱眩光源位置计算方法,针对现有的采用软件仿真以及人为评价对座舱眩光进行评估的方法费时费力、工作强度大等问题进行改善,实现对眩光进行高效、科学评估的目标。通过建立光源信息库、采集人眼视角的图像信息,并借助图像处理技术,解算眩光坐标、判断眩光源位置,提高评估的效率,提升评估结果的准确性。
本发明的技术方案是:
一种驾驶舱眩光源位置计算方法,包括:
根据眩光点在图像坐标系中的坐标、相机参数以及水平像差解算眩光源的空间坐标,根据眩光点空间坐标,光源坐标、反光面信息判断眩光源位置,最终得到眩光评估结果。
包括:
步骤一:采集飞行员左眼、右眼眼位处光环境图像信息。
步骤二:建立座舱照明光源、显示器信息库。
步骤三:将光环境图像进行灰度处理。
步骤四:提取眩光点在光环境图像上的平面坐标。
步骤五:解算眩光点空间坐标集glarecor={(x1,y1,z1),(x2,y2,z2)…(xn,yn,zn)}。
步骤六:判断眩光源。
步骤一具体为:
在左、右眼位点处设置像机,相机焦距f,保持垂直相差为零,光轴间距等于标准瞳距dise,采集座舱光环境图像。
步骤二具体为:
建立光源及显示器反光面数据库。已知驾驶舱中发光灯具数目为nl,光源的空间坐标集表示为lampcor={(xlamp1,ylamp1,zlamp1),(xlamp2,ylamp2,zlamp2)…(xlampnl,ylampnl,zlampnl)},已知驾驶舱中显示器发光平面数目为np,反光平面集表示为displayplane={(a1,b1,c1),(a2,b2,c2)…(anp,bnp,cnp)};
其中,(xlamp1,ylamp1,zlamp1)表示第一个光源的空间坐标,依此类推,(xlampnl,ylampnl,zlampnl)表示第nl个光源的空间坐标,(a1,b1,c1)表示第一个显示器反光平面的平面坐标参数,对应第一个显示器反光平面的平面方程a1x+b1y+c1z=1,依此类推,(anp,bnp,cnp)表示第np个显示器反光平面的平面坐标参数,对应第np个显示器反光平面的平面方程anpx+bnpy+cnpz=1,
步骤三具体为:
将左眼、右眼图像转换成灰度图像,图像矩阵分别标记为picl、picr。
步骤四具体为:
设左眼、右眼图像中眩光点数目分别为ml、mr,通过对图像矩阵元素顺序历遍以及阈值比较的方法,求解眩光点平面坐标;
计算过程包括:设置最佳截取阈值threshold,对图像矩阵picl、picr进行二值化处理,得到眩光矩阵piclg、picrg,对piclg、picrg进行顺序历遍,按历遍顺序得到眩光点中心坐标点集piclcor={(xleft1,yleft1),(xleft2,yleft2)…(xleftml,yleftml)},picrcor={(xright1,yright1),(xright2,yright2)…(xrightmr,yrightmr)};
其中,(xleft1,yleft1)表示左眼图像中第1个眩光点的平面坐标,依此类推,(xleftml,yleftml)表示左眼图像中第ml个眩光点的平面坐标,类似地,(xright1,yright1)表示左眼图像中第1个眩光点的平面坐标,依此类推,(xrightmr,yrightmr)表示左眼图像中第mr个眩光点的平面坐标。
步骤五具体为:
利用左眼、右眼图像中某眩光点平面坐标,计算该眩光点在显示器反光面上的空间坐标;
以第i个眩光点空间坐标的计算为例,设第i个眩光点空间坐标为(xi,yi,zi),坐标可用如下公式计算:
xi=dise·xlefti/d
yi=dise·ylefti/d
zi=dise·f/d
d=xlefti-xrighti
其中,dise表示双目标准瞳距,f表示相机焦距,xlefti、ylefti分别表示第i个眩光点在左眼图像中的平面坐标,xrighti、yrighti分别表示第i个眩光点在右眼图像中的平面坐标,由于垂直相差为零,因此ylefti与yrighti相等,d表示眩光点在左、右眼图像中的平面横坐标的差值。
步骤六具体为:
通过顺序历遍光源及显示器反光面数据库,寻找满足判据的光源坐标,即为该显示器上某一眩光点的眩光源;
依然以第i个眩光点的眩光源坐标求解为例,计算过程包括:
首先,计算第i个眩光点对应的像平面坐标,设第i个眩光点空间坐标glarecor(i)=(xi,yi,zi),在左像平面上的平面坐标为piclcor(i)=(xlefti,ylefti),对应的灰度等级为grayl,在右像平面上的平面坐标为picrcor(i)=(xrighti,yrighti),对应的灰度等级为grayr,通过如下公式比较grayl与grayr得到piccor(i):
piccor(i)=(grayl>grayr)?piclcor(i):picrcor(i)
记piccor(i)=(xpi,ypi);
然后,计算第i个眩光点所在的显示器平面坐标。设piccor(i)所在显示器平面参数为displayplane(j)=(aj,bj,cj),表示该眩光点所在反光平面坐标为ajx+bjy+cjz=1,将眩光点空间坐标代入反光平面坐标,满足等式,则该眩光点在此反光平面上,判据可表示为glarecor(i)与displayplane(j)内积为1:
最后,计算第i个眩光点对应的眩光源坐标。眩光源、显示器上的眩光点以及像平面上的眩光点构成反射路径平面,设光源j空间坐标为lampcor(j)=(xlampj,ylampj,zlampj),设反射路径平面参数为reflineplane=(p1,p2,p3),其平面方程p1x+p2y+p3z=1;
其中,
p2=(a-xpic)/(ypia–xpib)
p1=(1-p2ypi)/xpi
p3=((xi-xlampj)-p2(xiylampj-zixlampj))/(xizlampj-zixlampj)
a=xlampjzi–xizlampj
b=ylampjzi–yizlampj
c=zi–zlampj
对光源信息库进行历遍,通过如下判据,满足则光源j为眩光点i的眩光源:
reflineplane·displayplane(i)=0。
本发明的优点是:通过本专利中的评估方法,可实现座舱中眩光的快速评估。该评估方法模拟人眼视角,并借助图像处理技术,解算眩光坐标、判断眩光源位置,提高评估的效率,提升评估结果的准确性,评估效率高,可给出客观、量化的评估结果,可提升座舱光环境评估的效率和精度。
附图说明
图1是本发明对眩光源空间坐标进行识别的流程。
具体实施方式
如图1所示,下面分别以某型飞机座舱中眩光、眩光源的位置识别为例,作进一步详细说明。
一种驾驶舱眩光源位置计算方法,包括:
根据眩光点在图像坐标系中的坐标、相机参数以及水平像差解算眩光源的空间坐标,根据眩光点空间坐标,光源坐标、反光面信息判断眩光源位置,最终得到眩光评估结果。
步骤包括:
步骤一:采集飞行员左眼、右眼眼位处光环境图像信息。
在左、右眼位点处设置像机,相机焦距f,保持垂直相差为零,光轴间距等于标准瞳距dise,采集座舱光环境图像。
步骤二:建立座舱照明光源、显示器信息库。
建立光源及显示器反光面数据库。已知驾驶舱中发光灯具数目为nl,光源的空间坐标集表示为lampcor={(xlamp1,ylamp1,zlamp1),(xlamp2,ylamp2,zlamp2)…(xlampnl,ylampnl,zlampnl)},已知驾驶舱中显示器发光平面数目为np,反光平面集表示为displayplane={(a1,b1,c1),(a2,b2,c2)…(anp,bnp,cnp)};
其中,(xlamp1,ylamp1,zlamp1)表示第一个光源的空间坐标,依此类推,(xlampnl,ylampnl,zlampnl)表示第nl个光源的空间坐标,(a1,b1,c1)表示第一个显示器反光平面的平面坐标参数,对应第一个显示器反光平面的平面方程a1x+b1y+c1z=1,依此类推,(anp,bnp,cnp)表示第np个显示器反光平面的平面坐标参数,对应第np个显示器反光平面的平面方程anpx+bnpy+cnpz=1,
步骤三:将光环境图像进行灰度处理。
将左眼、右眼图像转换成灰度图像,图像矩阵分别标记为picl、picr。
步骤四:提取眩光点在光环境图像上的平面坐标。
设左眼、右眼图像中眩光点数目分别为ml、mr,通过对图像矩阵元素顺序历遍以及阈值比较的方法,求解眩光点平面坐标。
计算过程包括:设置最佳截取阈值threshold,对图像矩阵picl、picr进行二值化处理,得到眩光矩阵piclg、picrg,对piclg、picrg进行顺序历遍,按历遍顺序得到眩光点中心坐标点集piclcor={(xleft1,yleft1),(xleft2,yleft2)…(xleftml,yleftml)},picrcor={(xright1,yright1),(xright2,yright2)…(xrightmr,yrightmr)}。
其中,(xleft1,yleft1)表示左眼图像中第1个眩光点的平面坐标,依此类推,(xleftml,yleftml)表示左眼图像中第ml个眩光点的平面坐标,类似地,(xright1,yright1)表示左眼图像中第1个眩光点的平面坐标,依此类推,(xrightmr,yrightmr)表示左眼图像中第mr个眩光点的平面坐标。
步骤五:解算眩光点空间坐标集glarecor={(x1,y1,z1),(x2,y2,z2)…(xn,yn,zn)}。
利用左眼、右眼图像中某眩光点平面坐标,计算该眩光点在显示器反光面上的空间坐标。
以第i个眩光点空间坐标的计算为例,设第i个眩光点空间坐标为(xi,yi,zi),坐标可用如下公式计算:
xi=dise·xlefti/d
yi=dise·ylefti/d
zi=dise·f/d
d=xlefti-xrighti
其中,dise表示双目标准瞳距,f表示相机焦距,xlefti、ylefti分别表示第i个眩光点在左眼图像中的平面坐标,xrighti、yrighti分别表示第i个眩光点在右眼图像中的平面坐标,由于垂直相差为零,因此ylefti与yrighti相等,d表示眩光点在左、右眼图像中的平面横坐标的差值。
步骤六:判断眩光源。
通过顺序历遍光源及显示器反光面数据库,寻找满足判据的光源坐标,即为该显示器上某一眩光点的眩光源。
依然以第i个眩光点的眩光源坐标求解为例,计算过程包括:
首先,计算第i个眩光点对应的像平面坐标,设第i个眩光点空间坐标glarecor(i)=(xi,yi,zi),在左像平面上的平面坐标为piclcor(i)=(xlefti,ylefti),对应的灰度等级为grayl,在右像平面上的平面坐标为picrcor(i)=(xrighti,yrighti),对应的灰度等级为grayr,通过如下公式比较grayl与grayr得到piccor(i):
piccor(i)=(grayl>grayr)?piclcor(i):picrcor(i)
记piccor(i)=(xpi,ypi)。
然后,计算第i个眩光点所在的显示器平面坐标。设piccor(i)所在显示器平面参数为displayplane(j)=(aj,bj,cj),表示该眩光点所在反光平面坐标为ajx+bjy+cjz=1,将眩光点空间坐标代入反光平面坐标,满足等式,则该眩光点在此反光平面上,判据可表示为glarecor(i)与displayplane(j)内积为1:
最后,计算第i个眩光点对应的眩光源坐标。眩光源、显示器上的眩光点以及像平面上的眩光点构成反射路径平面,设光源j空间坐标为lampcor(j)=(xlampj,ylampj,zlampj),设反射路径平面参数为reflineplane=(p1,p2,p3),其平面方程p1x+p2y+p3z=1。
其中,
p2=(a-xpic)/(ypia–xpib)
p1=(1-p2ypi)/xpi
p3=((xi-xlampj)-p2(xiylampj-zixlampj))/(xizlampj-zixlampj)
a=xlampjzi–xizlampj
b=ylampjzi–yizlampj
c=zi–zlampj
对光源信息库进行历遍,通过如下判据,光源j是否为眩光点i的眩光源:reflineplane·displayplane(i)=0。