張紅偉,朱運東,林雪竹,李麗娟,劉濤
(長春理工大學 光電工程學院,長春 130022)
大尺寸測量主要指幾米至幾百米范圍內物體的空間位置、尺寸、形狀、運動軌跡等測量。當前,大尺寸測量技術在航空、航海和工業建設等眾多領域發揮著重要作用,視覺三維測量技術以其高效率、高精度等優勢在測量領域得到迅速發展[1-2]。
在測量領域,如何進行網絡規劃是測量過程中的難點問題,規劃布局的好壞直接影響到整體的測量精度和效率[3]。Olague等人[4]提出了全局性網絡優化設計思想,利用離散化組合對參數解空間進行求解;林雪竹等人[5]提出了基于iGPS測量系統的大部件對接位姿測量優化設計方法,通過建立iGPS系統測量模型和不確定度模型對接測量網絡,優化設計了對接測量網絡iGPS多發射器的布站方式;趙子越等人[6]詳細闡述了基于高精度控制場的全局組網定向方法的數學原理,包括約束方程建立以及迭代初值的獲取并通過實驗進行了驗證;劉虹[7]提出了三維掃描測頭精確跟蹤的攝影測量方法,并通過雙目傳感器驗證了該方法的可靠性;王文祥[8]通過對三維掃描測頭跟蹤系統的研究,利用多傳感器跟蹤解決了跟蹤視場盲區問題;劉曉利等人[9]闡述了基于條紋相位結構光的多節點三維傳感器測量網絡的網絡構造與規劃。以上研究均未涉及視覺測量網絡構建的研究,缺乏對傳感器站位設置方面的討論,實際工業測量中網絡規劃憑經驗布局,造成數據出現分層現象,導致測量失敗。同時開展視覺測量網絡規劃研究對我國大尺寸視覺掃描技術的發展和應用具有重要的意義[10]。
綜上,為了組建高精度視覺測量網絡,減小站位測量誤差,本文提出了基于雙目視覺測量網絡的站位布局規劃方法。對系統結構及測量原理進行闡述;根據視覺測量原理,構建雙目視覺測量模型,通過理論分析建立誤差模型;將影響參數根據作用機理將其劃分為與傳感器相對位置相關的內部結構參數和與測量特征相關的外部測量參數,仿真分析各參數變化對測量結果的影響,為網絡布局提供初值參數;經過區域劃分和站位優化,成功實現測量網絡構建;最后,通過實驗進行了驗證。
大尺寸視覺測量系統組成如圖1所示,由視覺跟蹤標準架(手持式三維形貌掃描儀)、視覺傳感器、計算機、控制盒及其他附件組成。視覺跟蹤標準架是一種其上按一定規則分布著若干個控制點,且具有線掃描功能的手持式設備。視覺跟蹤標準架測量被測工件時,按下相應按鈕,雙目視覺傳感器開始采集視覺跟蹤標準架上控制點圖像,通過相應算法可得到控制點成像中心坐標[11]。視覺跟蹤標準架作為一個剛體,通過控制點成像中心坐標實現視覺跟蹤標準架在雙目視覺傳感器站位下的實時跟蹤定位。最后,將視覺跟蹤標準架測得的點云數據進行基于世界坐標系的點云拼接,完成整體測量。為了實現大尺寸工件整體測量,需要組建多站位測量網絡,同時為了保證測量精度,避免點云拼接出現分層問題,需要合理規劃測量網絡中各個站位位置,如圖2所示,為由四個傳感器構成的雙站位雙目視覺測量網絡,虛線表示其中一組雙目視覺的有效測量區域。

圖1 系統組成框圖

圖2 視覺測量網絡
大尺寸視覺測量掃描系統中坐標系有視覺跟蹤標準架坐標系Os-XsYsZs、世界坐標系Ow-XwYwZw、視覺傳感器坐標系Ov-XvYvZv,如圖2所示。假設測量范圍內任一點P,在視覺跟蹤標準架坐標系下坐標矢量為Ps,在世界坐標系下坐標矢量為Pw,則二者的關系表示為:

測量點P在視覺跟蹤標準架坐標系Os-XsYsZs下的矢量Ps由標定好的視覺跟蹤標準架測量得到,測量時需要標定視覺跟蹤標準架坐標系與視覺傳感器坐標之間的Tvs轉換關系和視覺傳感器坐標系與世界坐標系之間的Twv轉換關系。至此,可進行整體的測量。
以對稱雙目視覺測量系統為出發點,通過理論分析,建立誤差傳遞模型,分析站位參數對測量結果的影響,為參數取值提供理論依據。
雙目視覺測量模型如圖3所示,視覺傳感器焦距分別為f1,f2,Ov1、Ov2為心,C1-x1y1,C2-x2y2為兩傳感器像平面坐標系,C1、C2為兩傳感器成像平面中心,傳感器坐標系Ov1-Xv1Yv1Zv1,Ov2-Xv2Yv2Zv2的Xv1,Yv2軸分別與其對應的像平面坐標系x1,x2軸平行,Zv1,Zv2垂直鏡頭指向傳感器背面,根據右手準則Yv1,Yv2垂直向上。設Ov1-Xv1Yv1Zv1坐標系為世界坐標系Ow-XwYwZw,C1、Ov1連線與C2、Ov2連線夾角稱為兩傳感器光軸夾角2β,兩傳感器光心之間的距離稱為基線B。

圖3 雙目視覺測量原理
設傳感器焦距f1=f2=f,空間任意一點P(X,Y,Z),在兩傳感器成像平面上的投影點坐標分別為p1(x1,y1),p2(x2,y2),對應其傳感器光心的水平視場角和垂直視場角γ1,γ2,θ1,θ2表示為:

根據雙目視覺成像原理,解得空間點P世界坐標(X,Y,Z)表示為:

為方便描述,將式(3)寫為:

根據誤差分配與合成理論[12]得空間點在各個坐標軸的測量誤差為:

其中,i分別表示各個變量;Δi表示各變量誤差,則空間點綜合測量誤差表示為:


通過以上分析,影響測量結果的參數主要有水平視場角、垂直視場角、基線距離、光軸夾角,根據各個參數對系統作用機制的不同,將參數劃分為兩部分,即由基線、光軸夾角構成的雙傳感器內部結構參數,和由水平視場角、垂直視場角構成的外部測量參數。利用綜合誤差模型分別對各個參數進行分離變量仿真,分析參數對系統測量精度影響趨勢,確定有效測量范圍。
2.3.1 內部結構參數對測量結果的影響分析
(1)基線對測量結果的因素分析
固定光軸夾角,分析基線距離對測量結果的影響。通過MATLAB仿真,結果如圖4所示,隨著基線距離的增大,測量誤差成線性增長。

圖4 基線對測量結果的影響
(2)光軸夾角對測量結果的因素分析
光軸夾角表示雙傳感器之間的相對姿態,由圖5可知,當基線距離一定時,隨著光軸夾角的增大,測量誤差呈先減小后增大的趨勢。

圖5 光軸夾角對測量結果的影響
基線和光軸夾角共同決定雙傳感器之間的相對位姿,綜合分析基線和光軸夾角雙因素對測量結果的影響,如圖6所示。

圖6 基線與光軸夾角對測量結果的影響
由圖6可知,隨著基線距離的增大,測量誤差增長的越快,隨著光軸夾角的增大,測量誤差減小的趨勢趨于平緩。然而并不意味著基線距離設置越小越好,基線設置越小,相應的測量范圍就越小,測量同一工件,需要的站位數目就多。
因此,需要結合實際測量需求,設置合理的基線距離和光軸夾角。在基線距離和光軸夾角確定的情況下,一個穩定的雙目視覺測量系統就構建完成了。
2.3.2 外部測量參數對測量結果的影響分析
(1)水平視場角對測量結果的因素分析
水平視場角的正負反映了測量點處于兩傳感器光軸夾角內部還是外部區域,保持垂直視場角一定,分析水平視場角對站位精度影響。
圖7為分析γ1,γ2對測量結果的影響,當水平視場角γ1,γ2同為正角度時,測量誤差較小;γ1,γ2同為負角度時,且隨著角度增大,測量誤差成遞增趨勢。因此,測量時,應盡量保證測量特征處于視場中心。

圖7 水平視場角對測量結果的影響
(2)垂直視場角對測量結果的因素分析
在水平視場角確定的情況下,垂直視場角可以確定空間點的具體位置。垂直視場角決定了成像點在成像芯片上的縱向距離,垂直視場角對測量結果的影響如圖8所示。

圖8 垂直視場角對測量結果的影響
由圖8可以看出,隨著垂直視場角的增大,測量誤差逐漸增大,即像點在成像平面的縱向距離和水平中線距離越遠,誤差越大。為了控制綜合測量誤差在測量范圍內均勻分布,保證像點縱向距離沿水平中心對稱分布,即測量特征沿水平中線對稱分布。
相同視場角的情況下,不同的工作距離測量結果不同。為了便于分析工作距離對測量結果的影響,設圖3中特征點相對傳感器基線距離為工作距離Z,在傳感器光軸方向的投影點記為X,設:

由圖9可知k=0.8-1.7內誤差變化相對平穩,k=1.4系統產生的誤差最小。實際測量中基線距的選取受到測量空間、成本等限制,需要結合實際測量環境進行設置。

圖9 基線距與物距關系對測量結果的影響
通過外部測量參數影響分析,對雙目視覺系統的相應測量區域設置提供了參考依據。
大尺寸物體需要設置多站位才能實現對所有測量特征的有效測量,且同一測量區域,會因為站位位置的不同超出約束導致站位失效或站位不合理的情況出現。因此,在滿足測量約束的條件下規劃站位位置,獲取較優的測量站位尤為重要。
針對大尺寸物體外形測量,結合視覺傳感器測量范圍提出區域劃分思想,每個子區域對應一個測量站位點。
區域劃分實質是對測量特征進行聚類劃分,單元子區域內測量特征在相應規劃的站位下可完成測量。通過對視覺測量系統工作原理分析,構建視覺測量多約束模型,在視覺測量多約束模型約束的作用下,實現測量特征的區域劃分。視覺測量多約束模型如下:
同時為了保證區域劃分結果均勻性與合理性,定義區域劃分調整參數為:

式中,yt代表第t個子區域測量的特征值。
其中,

Wi為排除規劃點后重新排列數據中的第i個點,VSN為當前站位的位置,Y為所有測量點數目,k為當前未規劃點數目。
結合區域優化參數,面向大尺寸物體進行區域劃分步驟如下:
(1)導入待規劃測量特征;
(2)生成初始站位位置信息;
(3)判斷測量特征是否符合視覺測量多約束模型;
(4)符合,將其劃分到子區域組,反之,將其規劃到不可測區域組;
(5)統計子區域特征值,即區域劃分調整參數值,根據其函數值調整站位位置,直到其值不再增加時停止;
(6)記錄結果,包括站位位置、子區域組特征信息和不可測區域組特征信息;
(7)判斷不可測區域組是否為空,如果是,執行下一步;反之,重新生成初始站位信息,執行步驟(2);
(8)輸出結果;
(9)結束。
流程圖如圖10所示。

圖10 區域劃分流程
站位規劃約束條件越多,復雜度越高,求解過程越復雜[13]。經過計算,每個子區域特征都滿足視覺傳感器特定約束。
大尺寸物體同一測量區域,不同站位與各測量特征距離、角度的不同,導致有效測量區域的差異。
因此,區域劃分完成之后,需要對各子區域規劃建立最優測量站位。
在第二部分中對雙目視覺測量系統進行了測量模型構建,并對其進行了測量結果精度實驗分析,視覺測量系統精度主要與投影角、光軸夾角、基線有關。根據雙目視覺系統測量精度實驗數據分析結果,運用前文構建的視覺測量系統精度數學模型對測量站位位置進行測量精度判斷,在站位信息迭代過程中以綜合誤差函數值為計算依據,確定各分區域下的較優測量站位位置。固定光軸夾角、基線距離,以子區域內測量特征相對于傳感器光心的水平視場角、垂直視場角為站位模型優化判斷參數,對子區域內所有特征進行測量約束參數遍歷計算,判斷是否超出測量約束值。綜上,大尺寸物體組網測量網絡優化驟如下:
(1)導入待規劃測量特征,將全局測量特征劃分為j組子區域;
(2)生成初始站位位置,遍歷各子區域特征是否符合當前站位多約束模型;
(3)符合,繼續下一步;反之,細化測量區域,重新進行區域劃分;
(4)計算綜合誤差值;
(5)重新迭代計算站位位置;
(6)輸出綜合誤差最小值站位位置;
(7)結束。
優化模型目標是為了尋找較優站位,使綜合誤差最小。考慮到優化模型的多參數和多維度等問題,提出利用遺傳算法[14-15]求解。遺傳算法(GA)有高度并行性、自適應和自學習等特征,能夠較好地解決具有較高組合性能的優化問題,其利用編碼基因表示求解變量,通過遺傳算子作用多個基因構成的染色體模擬自然進化過程,經過多次迭代獲得優化解。算法執行主要涉及編碼種類選擇、適應度函數fitness構建、遺傳算子參數值設定等步驟。編碼方法主要有:二進制、十進制、格雷碼等,編碼種類直接影響到算法效率,二進制和格雷碼均需要解碼,遺傳算子作用時涉及到復雜變換,為了節約效率,采用十進制編碼方式;適應度函數值體現了解的優良程度,適應度值越高,意味著有較大的概率進入下一步迭代,而綜合誤差值越小代表結果越好,因此進行適應度函數設置時需要對綜合誤差函數表達式作適當改變,令fitness=Fmax-Δ,Fmax為較大常數,遠遠大于Δ;遺傳算子包括選擇、交叉和變異算子,采用輪盤賭進行選擇操作,交叉和變異概率依次設為0.8和0.2,模擬自然界的進化過程;當全局適應度值連續迭代60次而未發生變化時,認為此時已達到最優解。
故基于遺傳算法的視覺傳感器站位優化模型求解步驟如下:
(1)導入待規劃測量特征,將全局測量特征劃分為j組子區域;
(2)種群初始化,生成初始站位位置,遍歷各子區域特征是否符合當前站位多約束模型;
(3)符合,繼續下一步;反之,增加區域劃分個數,重新進行劃分;
(4)遺傳算子作用,進行選擇、交叉和變異操作,產生新的站位信息;
(5)計算適應度函數值;
(6)是否滿足迭代截止條件;
(7)滿足,繼續下一步;反之,執行步驟(4),遺傳算子作用繼續進化;
(8)輸出適應度函數最大值即綜合誤差最小值對應的站位模型信息;
(9)結束。
根據上述算法建立的站位規劃流程圖如圖11所示。

圖11 站位優化流程圖
計算機型號HP-ZBOOK15,兩個工業相機,每秒可實時記錄30幅3D圖像,分辨率為7912×5432 pixel,實驗中通過云臺實現對相機水平角和俯仰角的控制,三角架實現對相機高度位置的調節。控制盒一端與相機和視覺跟蹤標準架相連,一端與計算機相連,計算機與視覺跟蹤標準架通過網線連接。
系統主要針對大尺寸物體進行多站位組網測量,實驗過程較為復雜,且需要進行坐標系轉換,綜上情況需要對視覺測量場進行精度驗證。如圖12所示,固定雙傳感器相對姿態,在雙傳感器優化測量范圍內隨機分布13個磁力靶座,以其中一個磁力靶座作為測量中心,測量另外12個磁力靶座到測量中心的距離。為減小隨機誤差影響,每組距離均重復測量5次,取5次測量結果的平均值作為測量結果。雙傳感器測量結果視為測量值,激光跟蹤儀測量的結果視為真值,計算相對誤差,并作為判斷標準,數據如表1所示。

圖12 測量環境

表1 實驗數據
由表1結果可知,該視覺測量系統的長度測量相對誤差均在0.03 mm以內,符合精度要求。
某車型如圖13所示,長3 600 mm,寬1 600 mm,高1 400 mm。設車幾何中心為坐標系原點,沿車身長方向為X方向,且正方向指向車頭,平行水平面為XOY,根據右手準則,Z軸垂直向上。根據以上分析結合車體的實際待測量區域,設置雙傳感器基線為1 500 mm,光軸夾角為40°,將車身劃分為五個子區域,站位信息如表2所示,俯視示意圖如圖14所示。

圖13 測量車

表2 站位數據

圖14 俯視圖
將在各站位得到的數據進行去噪、濾波,并將站位之間的共同測量數據進行偏差比較,結果如表3所示,1和2站位之間公共測量數據最小最大偏差分別為0.017 mm、0.032 mm,2和3站位之間公共測量數據最小最大偏差分別為0.017 mm、0.033 mm,3和4站位之間公共測量數據最小最大偏差分別為0.019 mm、0.032 mm,4和5站位之間公共測量數據最小最大偏差分別為0.017 mm、0.034 mm,1和5站位之間公共測量數據最小最大偏差分別為0.018 mm、0.035 mm,2和4站位之間公共測量數據最小最大偏差分別為0.018 mm、0.036 mm。由結果可知各站位數據偏差均在0.04 mm以內。

表3 站位共同區域測量數據偏差
最后將各個站位數據進行拼接得到整體測量數據,如圖15所示,各站位數據之間均沒有分層,實現了高精度拼接,滿足了工業測量精度要求。

圖15 整體數據
研究了基于雙目視覺的大尺寸測量網絡站位規劃方法,提高了站位測量精度,保證了整體測量精度。通過雙目視覺測量模型,基于理論分析建立誤差傳遞函數;根據參數影響機制將其劃分為內部結構參數和外部參數兩類,通過仿真得出各參數對測量結果的影響趨勢,為網絡布局優化過程中雙傳感器相對位置設定和區域劃分提供理論基礎;實驗結果表明,經過規劃后的測量網絡各站位測量數據相對誤差均在0.04 mm以內,實現了多站位數據之間的高精度拼接,滿足了工業測量要求。