徐旭,王紅改,李謀思
(1.武漢大學測繪學院,湖北 武漢 430079; 2.四川電力設計咨詢有限責任公司,四川 成都 610041;3.成都市規(guī)劃信息技術中心,四川 成都 610041; 4.武漢市勘察設計有限公司,湖北 武漢 430000)
隨著越來越多電力鐵塔的建設與運營,鐵塔安全狀態(tài)受到了廣泛關注。提取鐵塔運營過程中的結構特征,計算鐵塔傾斜度,實現鐵塔安全監(jiān)測是十分重要的問題?,F階段主要利用高精度免棱鏡全站儀對鐵塔結構進行測量,但該方法測量特征點數據較少,難以全面反映鐵塔結構特征,可靠性較低。地面三維激光掃描技術的推廣,為鐵塔安全監(jiān)測提供了新的技術手段,它突破了傳統(tǒng)的單點測量方法,具有高精度、高效率等優(yōu)勢,能通過海量點云數據獲取被測物體表面幾乎全部特征信息[1,2]。
梅文勝、李必軍分別利用地面三維激光掃描儀進行了精細地形圖測繪[3]與建筑物特征提取[4],電力行業(yè)謝雄耀等人利用三維激光掃描儀實現了電力桿塔掃描檢測[5];劉求龍等人利用三維激光掃描儀點云數據建立了惠泉變電站三維模型[6]。由于電力鐵塔結構形式多變,且掃描過程中受外界環(huán)境影響粗差點較多,點云后處理困難;針對這些問題,本文利用成都某變電站監(jiān)測項目鐵塔點云數據建立塔身獨立坐標系,采用基于鐵塔結構幾何特征的RANSAC算法提取鐵塔結構特征點云,計算鐵塔傾斜度,實現了鐵塔運營階段安全監(jiān)測。
地面三維激光掃描儀采用儀器內部坐標系,坐標原點為儀器中心,X、Y軸位于橫向掃描面內,Z軸垂直于橫向掃描面[7]。掃描儀發(fā)射激光射向目標點P,激光經反射后由儀器內部接收器接收,通過激光發(fā)射和接收時間差T計算目標點到儀器中心距離r,并記錄橫向掃描角φ和縱向掃描角θ,獲取被測點球面坐標P(r,φ,θ),根據掃描儀球面坐標系與直角坐標系轉換關系(公式1),得到被測點P的儀器內部空間直角坐標P(X,Y,Z):
(1)
圖1 掃描儀內部坐標系
地面三維激光掃描儀獲取的點云數據采用儀器內部坐標系,本文針對獨立設站鐵塔掃描點云,為便于后期鐵塔結構點云提取、建模分析,利用鐵塔點云數據根據鐵塔結構幾何特征建立塔身獨立坐標系。圖2為采用Maptek I-Site 8820掃描儀在成都某變電站監(jiān)測項目中獨立設站采集的鐵塔點云數據。
圖2 鐵塔原始點云數據
圖3 建立塔身坐標系流程圖
圖4 塔身坐標系坐標原點與XY軸
原始點云粗差剔除、特征提取一般采用最小二乘法、特征值法,這類方法是通過給定的目標方程,計算出最優(yōu)模型參數進行特征提取,其對于初始參數與目標方程的設定要求很高,算法不具有穩(wěn)健性,因此本文在鐵塔結構點云提取中采用能從包含異常數據的樣本中通過迭代方式自動提取出最優(yōu)參數模型的隨機抽樣一致性算法。
隨機抽樣一致性(Random Sampling Consensus,RANSAC)算法由Fishier和Bolles在1981年提出,該算法原理:對于一個數據量為m的樣本S,利用隨機最小數據量n求出其參數模型M,例如直線、平面模型中n分別為2和3;選擇合適的指標(例如平面點云中為點到平面模型的距離)對樣本中其他數據逐個在參數模型M中進行判定,并通過設定的閾值t進行篩選,若小于閾值t,則為“內點”,反之則相反,記錄該模型內點個數N。通過在樣本S中迭代多次選擇最小數據量n,計算不同參數模型M中獲取含內點數N最多的模型作為該樣本最優(yōu)參數模型。RANSAC算法通過最少的限定條件,利用點云自身樣本數據進行最優(yōu)參數模型選擇,該算法提取平面點云流程如下:
(1)從原始點云集合S中選取任意3個不共線點{P1,P2,P3},計算其初始平面模型Ax+By+Cz+D=0。
(2)遍歷點云集合S剩下點到上述初始平面模型的距離:
(3)設定合適的閾值t,即平面厚度。當di (4)重復上述(1)~(3)步驟K次,即迭代次數;RANSAC算法中定義K滿足公式f=1-(1-εn)K,其中f為至少能獲取一個有效參數模型的概率,ε為原始點云集合S中有效點云所占的比例(預估值),n為確定參數模型最小數據量,在平面中值為3。 (5)選擇K次迭代后內點個數N值最大的參數模型作為最優(yōu)平面擬合模型,該參數模型對應的內點即為剔除粗差的平面點云數據。 (2)CP(a)/=1,考慮商群G=G/Φ(P)=P/Φ(P)(〈a〉Φ(P))/Φ(P),則由文獻[10]中定理其中l(wèi)=|P|.進而有 G=P Φ(P)QΦ(P)Qh2··· Φ(P)Qhl.又Φ(P)=CP(a)/=1,頂點集連通且連接P中與Φ(P)相連的頂點,故P?(G)的連通分支個數為k(P?(G))=s1(P)-s1(Φ(P))+1. 由上節(jié)2.1可知,利用RANSAC算法提取平面點云過程中主要涉及閾值t、置信度f、有效數據比例ε、迭代次數K這4個參數的選擇。本文基于鐵塔鋼梁結構幾何特征進行上述參數合理地選取與設定,以獲取最優(yōu)結構平面點云集合。 如圖5所示,白色線條區(qū)域內點云為鐵塔導線橫擔橫梁點云數據,可以看出橫梁上邊緣點云極不規(guī)則,粗差點較多,而靠近下邊緣點云排列整齊。鐵塔結構中該白色點云部分長約 2.9 m,寬約 0.12 m,根據儀器掃描分辨率與下邊緣結構點云可判斷每相鄰掃描兩點間距約為 0.03 m,因此可計算得出該結構橫梁有效點云約為500點,而統(tǒng)計白色部分點云個數約為 1 200點,因此有效數據比例ε=0.42,至少能獲取一個有效參數模型的置信度f,根據經驗公式設為0.99[8];根據f=1-(1-εn)K可計算得出迭代次數K為70次。由于鐵塔鋼梁結構表面光滑平整,為有效剔除粗差點,設置平面厚度閾值t為 0.01 m。利用上述設定參數通過RANSAC迭代算法提取獨立坐標系下導線橫擔鋼梁點云數據,即圖6中白色線條區(qū)域內點云,該組點云約450點,且排列規(guī)則整齊,兩點間距離約為 0.03 m,清晰準確地反映了該導線橫擔鋼梁結構特征,驗證了利用該改進算法能有效地剔除粗差,提取鐵塔點云結構。 圖5 導線橫擔結構點云 圖6 基于鐵塔結構的RANSAC算法點云提取 通過改進的RANSAC算法剔除鐵塔掃描點云粗差,現提取塔身獨立坐標系中XZ平面內鐵塔結構點云數據,獲取該鐵塔順電力線路前進方向各橫擔、塔身等模型如圖7所示: 圖7 XZ平面內鐵塔模型 對圖7模型中部分距離、角度等值進行標注,并從塔身與橫擔夾角、各橫擔長度、左右兩側橫擔差值、左右兩側塔身鋼材寬度等進行比較分析,結果如表1、表2所示。 塔身模型距離比對 表1 塔身模型傾角比對 表2 通過表1、表2結果可以得到以下信息:①通過建立塔身坐標系、RANSAC算法結構點云提取、點云建模等流程獲取的鐵塔模型具有高精度與可靠性,鐵塔左右兩側同一比對項目差值達到毫米級;②該型轉角塔左右兩側橫擔的差值約為 34 cm,地線橫擔長為 14.6 m,判斷為2號轉角塔;③該鐵塔塔身投影傾角約為87.8°,與該類型鐵塔設計傾角相符。通過與該型鐵塔初始設計參數比對,驗證了該鐵塔建模成果的準確性。 傾斜度檢測是鐵塔結構檢測的一項重要指標,通過塔身中心軸的鉛垂度判斷鐵塔傾斜值。目前常用有鉛錘法、經緯儀法、平面鏡法,同時一些重要線路上建立了自動化傾斜監(jiān)測裝置[9]。本文利用地面三維激光掃描獲取的鐵塔點云數據,計算該鐵塔傾斜度。圖8為該鐵塔地線橫擔塔架切片在塔身坐標系XY平面內的點云,利用塔架結構點云切片求出該投影面幾何中心(即圖8中標注幾何中心),作為該鐵塔中心軸線上一點。 圖8 塔身切片幾何中心 由于鐵塔橫擔結構特征清晰明顯,點云數據量大,因此從上至下選取圖9中白色部分九組截面數據在XY投影面內求取鐵塔塔架中心坐標作為鐵塔軸線,結果如表3所示。 圖9 白色截面點云為待分析點云 塔架中心軸坐標 表3 續(xù)表3 由于鐵塔傾斜實際是從下至上逐步增大,為便于分析,設定第9組坐標為基準點,所在的XY平面為基準面。繪制如圖10、圖11所示。 圖10 塔身坐標系中X軸方向偏移趨勢 圖11 塔身坐標系中Y軸(順電力線路方向)偏移趨勢 通過表3和圖10、圖11可以得出: (1)該鐵塔中心軸在塔身坐標系X軸上最大偏移約為 8 mm,變化趨勢不明顯; (2)在Y軸(順電力線路前進方向)上隨著距離基準面高度的增加,偏移值加大,偏移趨勢約呈一定的線性關系(圖11中的黑色趨勢線)。 根據電力桿塔傾斜度計算公式: G=E/H×100% (2) 其中G為傾斜度,E為偏移值,H為距離基礎高度。該鐵塔高 40 m,計算該9組截面傾斜度。得出地線橫擔處(即鐵塔頂端)的傾斜度最大,約為0.48%,與表4正常桿塔傾斜允許最大值1%差值較大。得出該鐵塔順電力線路方向傾斜度處于閾值范圍內,處于安全運營狀態(tài)。 正常桿塔傾斜允許最大值 表4 作為新一代測繪儀器,地面三維激光掃描儀已逐步應用于各個領域,該類儀器的非接觸式、海量點云、全景獲取等優(yōu)勢十分明顯。本文利用地面三維激光掃描儀對成都某變電站附近鐵塔掃描觀測,主要做了如下工作: (1)利用鐵塔點云結構明顯的幾何特征建立了塔身獨立坐標系,將任意設站掃描獲取的內部坐標系轉換為基于塔身結構的獨立坐標系,便于后期鐵塔結構提取分析。 (2)采用RANSAC算法處理原始點云,開創(chuàng)性地基于鐵塔結構特征選定RANSAC算法參數,進行鐵塔掃描點云粗差剔除、結構提取。 (3)通過鐵塔點云建模與精度分析肯定了三維激光掃描儀的掃描精度與可靠性,并獲取了鐵塔運營階段中的各結構狀態(tài)參數。 (4)利用掃描點云對鐵塔進行傾斜度檢測,定量確定了鐵塔處于安全運營狀態(tài)。在今后工作中還需要進一步研究鐵塔傾斜原理與傾斜規(guī)律,定量定性實現鐵塔安全監(jiān)測。3.2 基于鐵塔結構的RANSAC算法經驗參數選取與點云提取
4 點云建模與精度分析
5 鐵塔傾斜度檢測
6 結 語