发明名称 基于可见光偏振成像的实时非接触水面波纹测量方法
摘要 一种基于可见光偏振成像的实时非接触水面波纹测量方法,首先利用偏振成像系统获取水面反射光的多偏振态图像,根据多偏振态图像计算反射光的偏振度和偏振方位角;然后利用偏振度,结合菲涅尔定律计算出入射至水面光线的入射角;最后基于偏振方位角和入射角,计算水面微面元的法向量。本发明方法为非接触成像测量,避免测量方法影响水面波纹,测量时不需要事先在水下放置靶标或模板等,适用于机载水下成像图像校正的实际需求;同时,能够实时给出水面波纹的二维空间坡度分布,便于机载水下成像的图像校正。
申请公布号 CN103389075B 申请公布日期 2015.07.22
申请号 CN201310309050.3 申请日期 2013.07.22
申请人 北京理工大学 发明人 金伟其;刘敬;王霞;裘溯;陈振跃;管文
分类号 G01C13/00(2006.01)I 主分类号 G01C13/00(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 仇蕾安;高燕燕
主权项 一种基于可见光偏振成像的实时非接触水面波纹测量方法,其特征在于,具体实施步骤如下:步骤一、利用偏振成像系统获取水面反射光的多偏振态图像,根据多偏振态图像计算反射光的偏振度DOLP和偏振方位角Ψ;步骤二、利用偏振度DOLP,结合菲涅尔定律计算出入射至水面光线的入射角θ<sub>i</sub>;步骤三、基于偏振方位角Ψ和入射角θ<sub>i</sub>,计算水面微面元的法向量n<sub>W</sub>;该步骤的具体过程为:(1)以偏振成像系统的光学主点为坐标原点O,以偏振成像系统中CCD的行方向和列方向分别为X轴和Y轴,并以光轴方向为Z轴,建立偏振成像系统坐标系O‑XYZ;设定入射光线由水面微面元进行反射;定义入射光和反射光所在平面为反射面,且入射光线与反射光线在水面微面元上的交点为A;定义反射光线与偏振成像系统上像面的交点为D,定义反射面与像面的交线与像面下边界线的交点为B;(2)计算直线AO的方向向量n<sub>AO</sub>,基于偏振方位角Ψ计算直线BD的方向向量n<sub>BD</sub>,然后根据n<sub>AO</sub>和n<sub>BD</sub>计算反射面的法向量n<sub>R</sub>;该步骤的具体过程为:根据点D(u,v,‑f)和O(0,0,0)的坐标,计算直线AO的方向向量n<sub>AO</sub>如式(6):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>n</mi><mi>AO</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>x</mi><mi>AO</mi></msub><mo>,</mo><msub><mi>y</mi><mi>AO</mi></msub><mo>,</mo><msub><mi>z</mi><mi>AO</mi></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>=</mo><msup><mrow><mo>[</mo><mfrac><mi>u</mi><msqrt><msup><mi>u</mi><mn>2</mn></msup><mo>+</mo><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mo>-</mo><mi>f</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>,</mo><mfrac><mi>v</mi><msqrt><msup><mi>u</mi><mn>2</mn></msup><mo>+</mo><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mo>-</mo><mi>f</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>,</mo><mfrac><mrow><mo>-</mo><mi>f</mi></mrow><msqrt><msup><mi>u</mi><mn>2</mn></msup><mo>+</mo><msup><mi>v</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mo>-</mo><mi>f</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>]</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000721448000000011.GIF" wi="1840" he="219" /></maths>根据步骤一中求得的偏振方位角Ψ求得直线BD的方向向量n<sub>BD</sub>如式(7):n<sub>BD</sub>=[x<sub>BD</sub>,y<sub>BD</sub>,z<sub>BD</sub>]<sup>T</sup>=[cosΦ,sinΦ,0]<sup>T</sup>=[‑sinΨ,cosΨ,0]<sup>T</sup>   (7)根据n<sub>AO</sub>和n<sub>BD</sub>计算反射面的法向量n<sub>R</sub>,如式(8):n<sub>R</sub>=n<sub>AO</sub>×n<sub>BD</sub>   (8)(3)根据入射角θ<sub>i</sub>、直线AO的方向向量n<sub>AO</sub>、反射面的法向量n<sub>R</sub>和入射光的方向向量n<sub>i</sub>,根据式(14)计算出水面微面元的法向量n<sub>W</sub>;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>n</mi><mi>W</mi></msub><mo>&CenterDot;</mo><msub><mi>n</mi><mi>AO</mi></msub><mo>=</mo><mo>|</mo><msub><mi>n</mi><mi>W</mi></msub><mo>|</mo><mo>|</mo><msub><mi>n</mi><mi>AO</mi></msub><mo>|</mo><mi>cos</mi><msub><mi>&theta;</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>n</mi><mi>W</mi></msub><mo>&CenterDot;</mo><msub><mi>n</mi><mi>R</mi></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mo>|</mo><msub><mi>n</mi><mi>W</mi></msub><mo>|</mo><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><msub><mi>n</mi><mi>AO</mi></msub><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>=</mo><mo>&PlusMinus;</mo><msub><mi>n</mi><mi>W</mi></msub><mo>|</mo><msub><mi>n</mi><mi>AO</mi></msub><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>|</mo></mtd></mtr><mtr><mtd><msub><mi>z</mi><mi>i</mi></msub><mo>=</mo><msub><mi>z</mi><mi>AO</mi></msub><mover><mo>+</mo><mo>-</mo></mover><msub><mi>z</mi><mi>W</mi></msub><mo>&CenterDot;</mo><mo>|</mo><msub><mi>n</mi><mi>AO</mi></msub><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>|</mo><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000721448000000021.GIF" wi="1430" he="424" /></maths>其中,z<sub>i</sub>为入射光的方向向量沿光轴方向的分量,z<sub>AO</sub>为直线AO的方向向量沿光轴方向的分量,z<sub>W</sub>为水面微面元法向量沿光轴方向的分量;步骤四、根据步骤三计算的水面微面元的法向量n<sub>W</sub>求解水面波纹坡度,实现水面波纹非接触实时测量。
地址 100081 北京市海淀区中关村南大街5号