蔣海里 陳柳花 程效軍 朱明芳 李金濤
(1. 上海公路橋梁(集團)有限公司, 上海 200433; 2. 同濟大學測繪與地理信息學院, 上海 200092;3. 自然資源部現(xiàn)代工程測量重點實驗室, 上海 200092)
作為一種新興的建筑結構形式,鋼結構相比于混凝土結構具有強度高、施工周期短、低碳環(huán)保等特點,因此鋼結構越來越多地被應用于超高層建筑、機場、體育館、大型橋梁等項目中[1-2]。大型鋼結構一般采用分段制造現(xiàn)場拼接的方式。為控制鋼結構質量、保證各部分構件在現(xiàn)場能順利實現(xiàn)拼裝,分段制造的鋼結構在出廠前需要進行模擬預拼裝。
三維激光掃描技術可以快速獲取鋼結構構件的三維點云模型,具有速度快、精度高、非接觸、全天候等特點,在大型鋼結構模擬預拼裝領域得到了許多工程界學者的廣泛關注。文獻[3]和[4]通過布設公共標靶,從掃描點云數(shù)據(jù)中提取標靶坐標并人工指定對應匹配點對,實現(xiàn)兩構件點云的模擬預拼裝,該方法耗時耗力,且當兩構件相距較遠時公共靶標難以滿足兩測站共視條件,靶標難以布設;文獻[5]擬合點云中的螺栓、構件圓心等特征點,通過點云配準實現(xiàn)構件點云的模擬預拼裝,該方法自動化程度低,精度較差;文獻[6]在點云模型中人工提取關鍵點,通過點云配準實現(xiàn)鋼結構的模擬預拼裝,該方法提取的特征點精度低,拼接誤差大;文獻[7]采用點云與BIM模型對齊的方式將多段構件分別與整體設計模型對齊后,分析相鄰構件間的拼裝誤差,從而實現(xiàn)模擬預拼裝?,F(xiàn)有的利用三維激光掃描技術進行鋼結構模擬預拼裝的方法主要為基于特征點匹配的方式,存在自動化程度低、工作量大且精度差的問題。
為此本文提出了一種基于線特征匹配的鋼結構模擬預拼裝方法,該方法首先提取待拼接處兩構件的特征線,然后通過線特征自動匹配的方式實現(xiàn)兩構件的模擬預拼裝。
本文提出的基于線特征匹配的鋼結構模擬預拼裝方法主要包括特征線提取和線特征匹配兩個過程,算法流程如圖1所示。通過將提取的線特征進行自動匹配實現(xiàn)兩構建點云的模擬預拼裝。
圖1 算法流程圖
為實現(xiàn)鋼結構構件的模擬預拼裝,針對配準后的兩待拼裝構件點云數(shù)據(jù),首先需要人工指定要進行拼接的兩構件位置。該過程通過CloudCompare軟件中點云的移動操作實現(xiàn),在軟件中通過簡單的平移操作將某構件待拼接處大致靠近另一構件的待拼接處,指定用兩構件的該處進行拼接。
本文采取沿構件拼接處特征線方向進行切片的方式實現(xiàn)點云的局部加密,進而從切片中提取該方向上的線特征點。為有效提取線特征點,首先需要確定切片方向。大型鋼結構構件拼接處特征線主要沿三個相互垂直的方向,因此需要確定三個相互垂直的方向并將其作為切片方向。在進行點云數(shù)據(jù)采集時對三維激光掃描儀進行的整平操作可使點云數(shù)據(jù)坐標系的Z坐標軸垂直于地面,大型鋼結構構件放置時拼接面通常與地面垂直,因此Z坐標軸為切片方向,將其作為第一切片方向(圖2中S1),與構件主軸線一致垂直于拼接面的方向為另一切片方向。為確定構件主軸線方向,將構件點云投影到X、Y平面,構件的主軸線方向在X、Y平面內可采用主成分分析法(Principal Component Analysis, PCA)[8]計算,對所有的投影點構造一協(xié)方差矩陣D
圖2 切片方向
(1)
矩陣D的最大特征值所對應的特征向量即為構件的主軸線方向,將其作為第二切片方向(圖2中S2)。第三切片方向與第一、第二切片方向垂直正交,構成右手直角坐標系(圖2中S3)。
對構件點云分別沿第一、第二、第三3個切片方向進行切片,即采用一組與切片方向相垂直的平行平面對構件點云進行分割,并將構建的三維點云投影到對應的二維分割平面上。設原始構建點云的坐標范圍為(Xmin,Ymin,Zmin)~(Xmax,Ymax,Zmax),以沿第一切片方向S1進行切片為例,此時平面點集由一組坐標序列Zm構成
Zm=Zmin+m·Zthick
(2)
其中,m=0,1,…,|h/Zthick|;m表示切片層數(shù);h=Zmax-Zmin;Zthick表示切片厚度;|·|表示向上取整函數(shù)。根據(jù)鋼結構構件拼接面的復雜程度和掃描點的疏密程度,Zthick取為3~5倍的平均點間距。
橋梁構件經(jīng)1.2中的方式沿線特征的三個切片方向分別進行切片后,橋梁構件的線特征點在切片上均表現(xiàn)為端點或拐角點。由于原始點云在掃描時具有一定的掃描密度和掃描誤差,此處的端點和拐角點均指以該點為中心,以r為半徑的近鄰域內點的集合,根據(jù)原始點云掃描密度的不同,r可取為1~5 mm。圖3為某個沿第一切片方向S1的切片,切片上的端點和拐角點均為鋼結構構件中垂直于地面的線特征點。因此構件中特征線的提取轉換為切片平面上端點和拐角點的提取,提取方式如下:
(1)循環(huán)某個切片上的所有點Pi(i=1,2,…,I,I為對應切片上的點數(shù)),對每個點Pi,搜索Pi的R近鄰域點,用R近鄰域點做局部特征分析,R取為鋼結構構件中方向一致的相鄰直線間的最短距離。
(2)計算R鄰域內每個點到Pi點的距離,若不存在距離大于r的點,則將Pi看作孤立的噪聲點,不進行保留;若R鄰域內存在距離大于r的點,則執(zhí)行步驟3~5。
(3)尋找R鄰域內距離Pi點最遠的點,將其作為為參考點Pr,Pr如圖3中的虛線點所示。
(4)將PiPr方向作為參考方向,計算Pi到R鄰域內其他點Pj的方向與參考方向之間的角度θj。
(5)若除r鄰域點(圖3中淺色實點)外,不存在角度在135°~180°之間的點,則將Pi作為特征點(端點或拐角點),對其進行保留;若除r鄰域點外,存在角度在135°~180°之間的點,則將Pi作為非特征點(內部點),不對其進行保留。
圖3 切片上特征點提取
切片上提取的平面特征點對應該切片垂直方向上的線特征點,因此將切片平面中的端點和拐角點映射到三維空間即得該切片方向上的特征線。圖4為沿第一切片方向的某相鄰幾個切片及映射到三維空間的線特征點。
圖4 平面切片及映射回三維空間的線特征點
鋼結構構件間的模擬預拼裝實際上為剛體坐標轉換過程,通過空間坐標變換使兩個鋼結構構件在拼接處匹配到一起。如圖5所示,圖5(a)為兩構件在拼接處的特征線,以圖5(a)中左側特征線為基準,將右側特征線通過坐標變換統(tǒng)一到左側特征線所在的坐標系下,拼接后兩構件的特征線應最優(yōu)地重合在一起,如圖5(b)所示。為實現(xiàn)該坐標變換過程,以2.3中提取的線特征點為匹配對象,采用最鄰近點迭代法(Iterative Closest Point, ICP)[9-11]進行坐標變換參數(shù)的計算。
圖5 模擬預拼裝圖示
采用線特征匹配的方式計算出拼裝構件的坐標變換參數(shù)后,對拼裝構件按照計算出的參數(shù)進行坐標變換,將兩構件統(tǒng)一到同一坐標系下即可實現(xiàn)鋼結構構件的模擬預拼裝過程。
以某鋼結構加工廠生產(chǎn)的某大型橋梁構件點云作為實驗數(shù)據(jù)對本文提出的模擬預拼裝方法進行實驗驗證。相鄰兩個橋梁構件的掃描點如圖6所示,構件A和構件B的長度分別為17.72和26.86 m,點數(shù)量分別為15,561,955和40,375,521,待拼接處點云的平均點間距為3.9 mm。
圖6 橋梁鋼結構構件點云數(shù)據(jù)
以26.86 m的長構件為基準,將17.72 m的短構件拼接到長構件上。首先在CloudCompare軟件中通過平移短構件點云來指定兩構件的拼接位置,如圖7所示。然后提取待拼接處的特征線,提取結果如圖8所示。通過特征線匹配實現(xiàn)橋梁構建的模擬預拼裝,其整體拼接結果如圖9所示,圖10為拼接處四個拼接面的放大拼接效果。
圖7 指定待拼接位置
圖8 提取待拼接處特征線
圖9 整體拼接結果
圖10 局部拼接結果
從圖8中可以看出,本文提出的特征線提取方法能較完整地提取出鋼結構構件待拼接處的特征線。從圖9和圖10中可以看出,本文方法能夠實現(xiàn)較好的預拼接效果,對點云密度較小的構件底部平面,也能實現(xiàn)較完整地拼接。圖10中左右兩邊分別為對應的兩個鋼結構構件。左側對應拼接的基準構件(長構件),右側對應待拼接構件(短構件)。
為定量分析本文方法的拼接精度,將拼接后兩構件在拼接處特征線的重合程度作為拼接精度的衡量指標,該重合程度用拼接構件拼接處特征線上所有的點到基準構件拼接處特征線上最近點的距離的平均值表示。采用該方法計算得實驗結果的拼接精度為7.6 mm,約為兩倍的平均點間距,具有較高的拼接精度。
與大多數(shù)基于特征點匹配的方式利用點云數(shù)據(jù)進行大型鋼結構構件模擬預拼裝的方法不同,本文提出了一種基于線特征匹配的模擬預拼裝方法。該方法首先基于切片的方式提取構件待拼接處的特征線,然后通過特征線匹配采用ICP算法計算最優(yōu)坐標變換參數(shù)進而實現(xiàn)構件的模擬預拼裝。該方法大大提高了大型鋼結構構件模擬拼裝的自動化程度,減少了人工干預。且實驗表明,該方法具有較高的模擬預拼裝精度,可用于具體鋼結構模擬預拼裝實踐。該方法對利用三維激光掃描技術實現(xiàn)大型鋼結構自動模擬預拼裝具有一定的指導和借鑒意義。