主权项 |
一种基于倾斜视角下透镜成像模型的河流水面流场定标方法,其特征在于,包括以下五个主要步骤:(1)设置断面标志杆,在测量断面1上设置两根和水面相交的标志杆2、3,用于指示断面方向,并将标志杆2和水面的交点作为水位参考点4;(2)架设成像测量装置,首先将一台工业相机和一台激光测距仪固连构成成像测量装置5,并架设在河岸一侧;其次调节三脚架的位置及云台的水平,使得单独改变俯仰角后激光测距仪的靶点能够依次照准标志杆2、3,以确保相机光轴6位于测量断面1上;然后用激光测距仪测量其到水位参考点4的斜距D和俯仰角δ,并由内置的三角测量程序计算测距仪到水面的垂直距离:H<sub>L</sub>=D·sinδ接下来选择焦距f合适的光学镜头并调节拍摄俯仰角,使得相机视场覆盖完整的测量断面并具有尽可能高的空间分辨率,进而根据激光测距仪当前的俯仰角α′求出工业相机光轴的俯仰角:α=α′+Δα其中,Δα表示预先标定的二者间的角度差值;最后计算工业相机到水面的垂直距离:H=H<sub>L</sub>‑d·cosα其中,d表示激光测距仪测量中心到工业相机成像中心的距离;(3)采集图像序列,成像测量装置稳定后,以Δt为时间间隔连续拍摄总时长为T的N幅图像用于流速测量;(4)标定测速线起点距,首先读取图像序列中的第一幅图像,并根据测速需求在图像中沿断面方向设置I条与顺流方向平行且大小为M×1的测速线L<sub>i</sub>(1≤i≤I),用(x<sub>i</sub>,y<sub>i</sub>)表示测速线中点的图像坐标;然后在图像中提取水位参考点4的图像坐标(x<sub>0</sub>,y<sub>0</sub>)作为计算测速线起点距的参考零点;最后根据倾斜视角下的透镜成像模型,利用图像中纵坐标j从y<sub>0</sub>到y<sub>i</sub>的各像素点在断面Y方向的物像尺度因子ΔY(j)标定测速线的起点距:<maths num="0001"><math><![CDATA[<mrow><msub><mi>Dist</mi><mi>i</mi></msub><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>j</mi><mo>=</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><msub><mi>y</mi><mi>i</mi></msub></munderover><mi>Δ</mi><mi>Y</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>j</mi><mo>=</mo><msub><mi>y</mi><mn>0</mn></msub></mrow><msub><mi>y</mi><mi>i</mi></msub></munderover><mi>H</mi><mo>·</mo><mo>{</mo><mrow><mn>1</mn><mo>/</mo><mi>tan</mi><mo>[</mo><mi>α</mi><mo>+</mo><mi>arctan</mi><mfrac><mrow><mo>(</mo><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo><mi>s</mi></mrow><mi>f</mi></mfrac><mo>]</mo><mo>-</mo><mn>1</mn><mo>/</mo><mi>tan</mi><mo>[</mo><mi>α</mi><mo>+</mo><mi>arctan</mi><mfrac><mrow><mo>(</mo><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><mi>j</mi><mo>)</mo><mi>s</mi></mrow><mi>f</mi></mfrac><mo>]</mo></mrow><mo>}</mo></mrow>]]></math><img file="FDA0001005389410000021.GIF" wi="1926" he="159" /></maths>其中,m×n表示相机的图像分辨率,s表示图像传感器的像元尺寸;(5)标定测速线流速值,对于每条测速线L<sub>i</sub>,首先以测速线上像素的空间为横坐标、每帧图像的时间为纵坐标,建立大小为M×N的时空图像;然后采用时空图像测速法估计测速线上运动矢量的大小:<maths num="0002"><math><![CDATA[<mrow><msub><mi>v</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>s</mi><mi>i</mi></msub><mi>T</mi></mfrac><mo>=</mo><mfrac><msub><mi>s</mi><mi>i</mi></msub><mrow><mi>N</mi><mo>·</mo><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>t</mi><mi>a</mi><mi>n</mi><mi>θ</mi></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow>]]></math><img file="FDA0001005389410000022.GIF" wi="464" he="119" /></maths>其中,s<sub>i</sub>表示时间T内目标沿顺流方向运动的像素距离,θ表示时空图像的纹理主方向;最后根据倾斜视角下的透镜成像模型,利用测速线上像素点在顺流X方向的物像尺度因子ΔX(y<sub>i</sub>)标定测速线上物理流速的大小:<maths num="0003"><math><![CDATA[<mrow><msub><mi>V</mi><mi>i</mi></msub><mo>=</mo><msub><mi>v</mi><mi>i</mi></msub><mo>·</mo><mi>Δ</mi><mi>X</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>v</mi><mi>i</mi></msub><mo>·</mo><mfrac><mrow><mi>H</mi><mi>s</mi></mrow><msqrt><mrow><msup><mrow><mo>[</mo><mrow><mo>(</mo><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>s</mi><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>f</mi><mn>2</mn></msup></mrow></msqrt></mfrac><mo>/</mo><mi>t</mi><mi>a</mi><mi>n</mi><mo>[</mo><mi>α</mi><mo>+</mo><mi>arctan</mi><mfrac><mrow><mo>(</mo><mi>n</mi><mo>/</mo><mn>2</mn><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo><mi>s</mi></mrow><mi>f</mi></mfrac><mo>]</mo></mrow>]]></math><img file="FDA0001005389410000023.GIF" wi="1478" he="179" /></maths>其中,V<sub>i</sub>的正负反映了测速线上流速的方向。 |