偏振成像入射角歧义解消除方法及应用与流程

本发明属于偏振成像技术领域，尤其涉及由反射光的偏振度反演入射光的入射角时，入射角歧义解消除方法及应用。

背景技术：

近年来，偏振成像技术得到迅速的发展，利用光电偏振成像不仅可获取目标场景的辐射强度信息，而且可解算得到目标场景的偏振度、偏振角和辐射率等信息，增加在复杂背景下对目标的探测和识别能力，在目标探测和分类、水面波纹测量和空间遥感探测等领域展现出广泛的应用前景。海洋波纹检测等是目前偏振成像技术研究广泛关注的方向，其以太阳光或天空光作为无偏的辐射光源，通过偏振成像系统对水面成像求得反射光的斯托克斯矢量，进而得到反射光的偏振度和偏振角信息，并通过菲涅耳公式反演出入射光的入射角以及水面法向量，复原出海洋的波面分布。然而，依据菲涅尔公式，由偏振度求解入射角时存在双值歧义解问题，难以唯一确定入射角。

对于如何去除偏振成像求解入射角歧义解已成为人们研究重要内容。2002年daisukemiyazaki提出分别采集旋转前后被测物体的偏振图像，分析被测曲面反射光偏振度的一阶微分来消除歧义解的方法。1995年mikepartridge和2002年daisukemiyazaki均研究了分别在可见波段和红外波段采集目标的偏振图像，由两条不同折射率的偏振度与入射角关系来获得唯一确定的入射角。2005年garyaatkinson和2012年alihmahmoud研究了利用漫反射光来消除入射角的歧义解。2012年christophestolz提出测量不同可见光波段的偏振度数据，分析偏振度对波长的微分来消除入射角的歧义值。这些方法及测量装置较复杂，对于基于海面波纹检测或水下目标探测等应用并不实用。

因此需要设计一种方法及测量装置较简单，实用于基于海面波纹检测或水下目标探测等应用的入射角歧义解消除方法。方法既要测量装置简单，又要避免对目标多次测量，以适应动态海面波纹测量需求。

技术实现要素：

为解决现有入射角歧义解消除方法中方法及测量装置较复杂，不适用于基于海面波纹检测或水下目标探测等应用的问题。本发明公开的偏振成像入射角歧义解消除方法及应用的目的是：实现偏振成像入射角歧义解消除，得到真实的偏振成像入射角信息，避免对同一目标的多次测量，提高测量效率，此外，由于本发明单次测量即能够消除偏振成像入射角歧义解，尤其通过单次测量得到动态海面的真实的偏振成像入射角信息，进而能适应动态海面波纹测量需求，实现对动态海面波纹测量。本发明具有测量装置简单的优点。

通过本发明公开的偏振成像入射角歧义解消除方法测量得到的动态海面波纹信息，通过进一步分析能够应用于水下运动目标的探测。

本发明目的是通过下述技术方案实现的。

本发明公开的偏振成像入射角歧义解消除方法，通过可见光偏振成像系统对目标的反射光成像，计算目标反射光偏振度dolp；根据菲涅尔公式，由偏振度dolp反演目标反射光的入射角θi，所述反演得到的入射角有两个解，其中一个入射角解为真实入射角信息，另一入射角解为歧义解；通过反射光光强的投影关系，求得反射光垂直方向强度ir⊥；测量反射界面处入射的自然光强，入射光垂直分量强度ii⊥为入射的自然光强iin的一半，垂直分量的反射率r⊥为反射光垂直方向强度ir⊥与入射光垂直分量强度ii⊥之比；根据垂直分量的反射率r⊥，通过垂直分量反射率r⊥和入射角θi的关系，反演出入射光的入射角，得到参考入射角θir。选取两个候选入射角θis和θim中更靠近参考入射角θir的入射角作为真实入射角θi信息。

本发明公开的偏振成像入射角歧义解消除方法，包括如下步骤：

步骤1、通过可见光偏振成像系统对目标的反射光成像，计算目标反射光偏振度dolp。

为了便于测量和降低噪声对可见光偏振成像系统测量精度的影响，作为优选，通过可见光偏振成像系统对反射光成像，选取检偏方向分别为0°、45°、90°、135°的四通道强度计算反射光偏振度，计算公式为：

其中，dolp为线偏振度，i0为检偏方向为0°的强度，i45为检偏方向为45°的强度，i90为检偏方向为90°的强度，i135为检偏方向为135°的强度。

步骤2、根据菲涅尔公式，由偏振度dolp反演目标反射光的入射角θi，所述反演得到的入射角有两个解，其中一个入射角解为真实入射角信息，另一入射角解为歧义解。

根据反射光的偏振度dolp，通过偏振度和入射角的关系，反演出入射光的入射角θi，得到两个候选入射角θis和θim。偏振度dolp和入射角θi的关系为：

步骤3、通过反射光光强的投影关系，求得反射光垂直方向强度ir⊥。

作为优选，相机灰度级与其接收辐射强度的事先标定，定义0°检偏方向为偏振成像系统的参考方向。利用0°和90°的光强i0、i90，以及偏振成像系统的参考方向与反射面的夹角αsys，由下式求得反射光垂直方向强度ir⊥：

ir⊥＝i0sin²αsys+i90cos²αsys

步骤4、测量反射界面处入射的自然光强iin，入射光垂直分量强度ii⊥为入射的自然光强iin的一半，垂直分量的反射率r⊥为反射光垂直方向强度ir⊥与入射光垂直分量强度ii⊥之比。

步骤4具体实现方法为：

测量反射界面处入射的自然光强iin，得到入射光垂直分量强度ii⊥＝iin/2，再由下式获得垂直分量的反射率r⊥：

r⊥＝ir⊥/ii⊥

步骤5、根据步骤4得到的垂直分量的反射率r⊥，通过垂直分量反射率r⊥和入射角θi的关系，反演出入射光的入射角θi，得到参考入射角θir。垂直分量反射率和入射角的关系为：

步骤6、选取两个候选入射角θis和θim中更靠近参考入射角θir的入射角作为真实入射角信息。所述两个候选入射角θis和θim已由步骤2得到。所述的真实入射角信息即为消除歧义解唯一的入射角。

步骤6具体实现方法为：

步骤2已得到两个候选入射角θis和θim，通过下式判断条件得到更靠近参考入射角θir的入射角作为真实入射角信息，所述的真实入射角信息即为消除歧义解唯一的入射角。

还包括步骤7：利用步骤6获得的消除歧义解唯一的入射角θi，计算动态海面每点微面元法向量n，重构动态海面波纹。

当所述偏振成像入射角歧义解消除方法应用于测量动态海面波纹时，通过步骤6得到动态海面的真实的偏振成像入射角信息θi，通过调整可见光偏振成像系统姿态，使反射光线方向平行于偏振成像系统光轴方向。在反射光沿光轴传播且反射光只存在s分量条件下，俯仰角θ等于入射角θi，通过如下公式关系得到动态海面波纹，实现对动态海面波纹测量。

n＝[cos(aop±90°)sinθi,sin(aop±90°)sinθi,cosθi]

其中，n为动态海面微面元法向量，aop为反射光偏振角。

还包括步骤8：由于水下运动目标与水体的作用，使得其在水面会形成特定形状的波纹，其中水下运动目标所在位置形成的水面凸起称为“伯努利水丘”，而目标后面所形成的尾迹称为“开尔文尾迹”。对“伯努利水丘”和“开尔文尾迹”的探测成为水下运动目标探测的信息源。通过分析步骤7测量得到的动态海面波纹信息，实现对水下运动目标的探测，进一步，利用步骤7测量得到的动态海面波纹信息，判断动态海面波纹信息中是否存在“伯努利水丘”和“开尔文尾迹”，进而判断是否在该区域存在水下运动目标。进一步，若动态海面波纹信息中存在“伯努利水丘”或“开尔文尾迹”，还可通过分析“伯努利水丘”或“开尔文尾迹”大小、形状、高度等特征，获取水下运动目标的大小、深度、运动速度等信息。

有益效果：

1、本发明公开的偏振成像入射角歧义解消除方法，通过可见光偏振成像系统对目标的反射光成像，计算目标反射光偏振度dolp；根据菲涅尔公式，由偏振度dolp反演目标反射光的入射角，所述反演得到的入射角θi有两个解，其中一个入射角解为真实入射角信息，另一入射角解为歧义解；通过反射光光强的投影关系，求得反射光垂直方向强度ir⊥；测量反射界面处入射的自然光强，入射光垂直分量强度ii⊥为入射的自然光强iin的一半，垂直分量的反射率r⊥为反射光垂直方向强度ir⊥与入射光垂直分量强度ii⊥之比；根据垂直分量的反射率r⊥，通过垂直分量反射率r⊥和入射角θi的关系，反演出入射光的入射角，得到参考入射角θir。选取两个候选入射角θis和θim中更靠近参考入射角θir的入射角作为真实入射角信息θi。在确保入射角测量精度及求解稳定性的同时，准确的去除偏振度计算入射角的歧义解。

2、本发明公开的偏振成像入射角歧义解消除方法，当所述偏振成像入射角歧义解消除方法用于测量动态海面波纹时，通过步骤6得到动态海面的真实的偏振成像入射角信息，通过所述入射角信息得到动态海面波纹，实现对动态海面波纹测量。

3、本发明公开的偏振成像入射角歧义解消除方法，能满足基于海面波纹检测或水下目标探测等应用的需求，对基于偏振成像的透明介质表面和海面波纹等成像检测具有理论意义和应用意义。

4、本发明公开的偏振成像入射角歧义解消除方法，通过分析测量得到的动态海面波纹信息，实现对水下运动目标的探测。

附图说明

图1是本发明公开的偏振成像入射角歧义解消除方法及应用的流程图。

图2是偏振度和入射角的关系；

图3是反射光振动方向、反射面与偏振成像参考方向的关系示意图；

图4是垂直分量反射率r⊥与入射角的关系；

图5是相机灰度级与其接收照度的关系；

图6是具体实施例1采用本发明后去除歧义解效果图，其中6a为平板玻璃的偏振特性检测实验原理，6b为候选入射角，6c为去除歧义解条件，6d为去除歧义解后的入射角；

图7是具体实施例2采用本发明后去除歧义解效果图，其中7a为水面的偏振特性检测实验原理，7b为候选入射角，7c为去除歧义解条件，7d为去除歧义解后的入射角。

具体实施方式

为使本发明的技术方案更加清楚明白，下面结合实例对本发明做进一步的详细说明。

实例1：以平板玻璃偏振特性检测为例:

选取平板玻璃(折射率n＝1.5163，布儒斯特角θb＝56.625°)作为示例，搭建平板玻璃偏振特性成像检测系统(图6a)，通过调整平板玻璃支架得到不同的入射角(入射角等于玻璃板倾角θ)，采集了10°，20°…80°八组入射倾角，且每个倾角采集0°、45°、90°和135°四个检偏方向的图像。

本实施例公开的偏振成像入射角歧义解消除方法，具体实现步骤如下：

步骤1、通过可见光偏振成像系统对目标的反射光成像，计算目标反射光偏振度。

通过可见光偏振成像系统对反射光成像，选取检偏方向分别为0°、45°、90°、135°的四通道强度计算反射光偏振度，计算公式为：

其中，dolp为线偏振度，i0为检偏方向为0°的强度，i45为检偏方向为45°的强度，i90为检偏方向为90°的强度，i135为检偏方向为135°的强度。

10°，20°…80°八组入射倾角测量的偏振度dolp分别为0.0263、0.1479、0.3526、0.5598、0.8718、0.9749、0.6711、0.2879。

步骤2、根据菲涅尔公式，由偏振度dolp反演目标反射光的入射角θi，所述反演得到的入射角有两个解，其中一个入射角解为真实入射角信息，另一入射角解为歧义解。

根据反射光的偏振度dolp，通过偏振度和入射角的关系，反演出入射光的入射角θi，得到两个候选入射角θis和θim(如图6b)。偏振度dolp和入射角θi的关系为：

步骤3、通过反射光光强的投影关系，求得反射光垂直方向强度ir⊥。

相机灰度级与其接收辐射强度的事先标定(如图5)，定义0°检偏方向为偏振成像系统的参考方向(如图3)。利用0°和90°的光强i0、i90，以及偏振成像系统的参考方向与反射面的夹角αsys，由下式求得反射光垂直方向强度ir⊥：

ir⊥＝i0sin²αsys+i90cos²αsys

步骤4、测量反射界面处入射的自然光强iin，入射光垂直分量强度ii⊥为入射的自然光强iin的一半，垂直分量的反射率r⊥为反射光垂直方向强度ir⊥与入射光垂直分量强度ii⊥之比。

测量反射界面处入射的自然光强iin，得到入射光垂直分量强度ii⊥＝iin/2，再由下式获得垂直分量的反射率r⊥：

r⊥＝ir⊥/ii⊥

10°，20°…80°八组入射倾角测量的垂直分量的反射率r⊥分别为0.0613、0.0659、0.0672、0.0971、0.0896、0.1820、0.3576、0.4728。

步骤5、根据步骤4得到的垂直分量的反射率r⊥，通过垂直分量反射率r⊥和入射角θi的关系，反演出入射光的入射角θi，得到参考入射角θir。垂直分量反射率和入射角的关系为：

10°，20°…80°八组入射倾角测量的参考入射角θir分别为30.497°、33.284°、34.009°、45.348°、43.117°、59.980°、72.765°、77.602°。

步骤6、选取两个候选入射角θis和θim中更靠近参考入射角θir的入射角作为真实入射角信息。所述两个候选入射角θis和θim已由步骤2得到。所述的真实入射角信息即为消除歧义解唯一的入射角。

步骤2已得到两个候选入射角θis和θim，通过下式判断条件得到更靠近参考入射角θir的入射角作为真实入射角信息(如图6c)，所述的真实入射角信息即为消除歧义解唯一的入射角(如图6d)。

还包括步骤7：利用步骤6获得的消除歧义解唯一的入射角θi，计算动态海面每点微面元法向量n，重构动态海面波纹。

当所述偏振成像入射角歧义解消除方法应用于测量动态海面波纹时，通过步骤6得到动态海面的真实的偏振成像入射角信息θi，通过调整可见光偏振成像系统姿态，使反射光线方向平行于偏振成像系统光轴方向。在反射光沿光轴传播且反射光只存在s分量条件下，俯仰角θ等于入射角θi，通过如下公式关系得到动态海面波纹，实现对动态海面波纹测量。

n＝[cos(aop±90°)sinθi,sin(aop±90°)sinθi,cosθi]

其中，n为动态海面微面元法向量，aop为反射光偏振角。

还包括步骤8：由于水下运动目标与水体的作用，使得其在水面会形成特定形状的波纹，其中水下运动目标所在位置形成的水面凸起称为“伯努利水丘”，而目标后面所形成的尾迹称为“开尔文尾迹”。对“伯努利水丘”和“开尔文尾迹”的探测成为水下运动目标探测的信息源。通过分析步骤7测量得到的动态海面波纹信息，实现对水下运动目标的探测，进一步，利用步骤7测量得到的动态海面波纹信息，判断动态海面波纹信息中是否存在“伯努利水丘”和“开尔文尾迹”，进而判断是否在该区域存在水下运动目标。进一步，若动态海面波纹信息中存在“伯努利水丘”或“开尔文尾迹”，还可通过分析“伯努利水丘”或“开尔文尾迹”大小、形状、高度等特征，获取水下运动目标的大小、深度、运动速度等信息。

实例2：以水面偏振特性检测为例:

本实施例选取水面(折射率n＝1.33，布儒斯特角θb＝53.089°)作为示例，搭建了水面偏振特性成像检测系统(图7a)，通过调整偏振成像系统姿态得到不同的入射角，采集10°，20°…80°八组入射倾角，且每个倾角采集0°、45°、90°和135°四个检偏方向的图像。操作步骤完全同实施例1，此处不再重复。

以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。