主权项 |
一种基于可见光偏振成像的实时非接触水面波纹测量方法,其特征在于,具体实施步骤如下:步骤一、利用偏振成像系统获取水面反射光的多偏振态图像,根据多偏振态图像计算反射光的偏振度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>·</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>θ</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>n</mi><mi>W</mi></msub><mo>·</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>±</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>·</mo><mo>|</mo><msub><mi>n</mi><mi>AO</mi></msub><mo>-</mo><msub><mi>n</mi><mi>i</mi></msub><mo>|</mo><mo>≥</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>求解水面波纹坡度,实现水面波纹非接触实时测量。 |