隧道三維點云孔洞修復方法

錢伯至，藍秋萍

(河海大學 地球科學與工程學院，江蘇 南京 211100)

?

隧道三維點云孔洞修復方法

錢伯至，藍秋萍

(河海大學 地球科學與工程學院，江蘇 南京 211100)

針對隧道三維激光掃描點云中的孔洞，設計一種基于斷面輪廓的孔洞修復方法。隧道具有深長性與截面恒定性兩大特征，文中基于這兩大特征，利用間鄰穩定輪廓擬合完整隧道修復孔洞。該方法能夠對三維激光掃描點云進行圓柱投影，利用Delaunay三角構網，提出基于隧道表面三角網的孔洞識別方法，最后根據相鄰隧道輪廓相似的特點，利用相鄰的輪廓對孔洞進行修復。該方法能很好的修復孔洞并保留隧道結構細節，并以真實的隧道掃描數據驗證算法的有效性和可靠性。

三維點云；Delaunay 三角構網；孔洞識別；孔洞修補

近年來，隨著交通建設發展，隧道得到人們關注，而隧道的變形監測是隧道建設中必不可少的一項環節。由于隧道存在光線不足、狹長等問題，常規的測量難以適用于隧道，而三維激光掃描技術又稱為“實景復制技術”，它具有非接觸、掃描速度快、獲取信息量大、精度高且對環境要求較低等優點[1]。激光掃描通過自主發射光束，可以解決隧道光線不足的問題，適合隧道參數的測量。

然而因為隧道內交通不可中斷，且掃描儀的掃描垂直角存在一定的限度，致使隧道掃描點云中存在大量孔洞。從現有的文獻可見，目前已有許多孔洞修復的方法：如多視圖立體三維重建孔洞修復算法[2]，該方法基于視圖重建得到的三維點云與可見外殼點云融合算法[3-4]，在點云的稀疏度約束和曲率約束下修復點云，然而該方法對部分凹型區域也進行一定修復，對原有結構造成一定破壞；晏海平等提出基于徑向基函數的散亂點云孔洞修復算法[5]，通過內外邊界邊的識別進行孔洞識別，然后基于徑向函數的隱式曲面插值修復孔洞[6-7]，但是對于邊緣這類突變棱角，利用連續光滑曲面的約束條件重建缺失數據顯然有違于原始對象的幾何特征。

本文基于隧道輪廓的特征，利用相鄰隧道斷面輪廓修復孔洞。基于隧道斷面輪廓線檢測的孔洞識別最優截面模型，以中軸線為基礎復刻完整隧道，實現孔洞擬合修復。

1 截取斷面輪廓

1.1 中軸線計算及三角構網

圖1 角度與軸線距離匹配

表1平均邊長變化規律表

圖2 三角網網型

1.2 斷面輪廓截取及孔洞識別

基于隧道表面三角網，采用垂直于中軸線的截斷面與三角網進行空間求交獲得交點。Delaunay三角構網后所形成的隧道三角網網型較密且孔洞范圍內也存在三角網，在截取斷面輪廓線時難以區別孔洞。根據隧道孔洞附近三角網特點打斷孔洞位置的三角網邊，顯現孔洞。

圖3 點云孔洞

圖4 斷面截取

2 輪廓線修復

從隧道的特征出發，利用相鄰完整的斷面輪廓對孔洞位置進行修復。隧道內壁上存在吊燈、通風扇、管線等附屬物，使得隧道斷面輪廓有細微的差別，本文采用識別最優輪廓的方式，甄別出與需要填補的斷面輪廓最相似的斷面，并以此斷面為基準進行平移修補。最后利用補充中點的方式對隧道檢修道進行孔洞修復。

(1)

圖6中黑色點表示待修復斷面P，紅色點為某

圖5 孔洞修復

圖6 最優輪廓面

由于三維激光掃描的特點，導致隧道內檢修道面等曲率變化很小的結構在每個截斷面上的點都比較少，經過上述方法修復孔洞后該結構依舊會存在小孔洞，在已經保留隧道結構細節的前提下，對該結構利用插入中點的方式修復孔洞，使得隧道所有孔洞均被修復，且保證隧道內部的結構細節不被破壞。

3 實驗分析

實驗數據來自寧杭高速梯子山隧道的三維激光掃描點云，本次實驗采用該隧道其中某一測站的點云數據，點云個數為187 568個，該段隧道長度為108.785 m。隧道三維激光掃描過程中交通不可阻斷，所以隧道豎直墻面存在由汽車等物體遮擋所導致的孔洞，且由于掃描儀本身掃描垂直角存在一定限度，使得測站中心所在位置的拱頂與地面存在較大孔洞。實驗利用上述方法對隧道孔洞進行修復，并保留吊燈、通風扇、管線等附屬物。

對點云進行三角構網并截取斷面輪廓。依據斷面輪廓線的長度可識別孔洞的存在。隨后利用相鄰斷面輪廓進行孔洞修復(根據修復效果，可進行多次修復，以保證孔洞修復點云均勻)。最后通過插入中點的方式進行修復，完成孔洞的修復。

圖7(a)為隧道原始點云，(b)表示隧道拱頂存在的吊燈，(c)表示隧道拱頂和路面的孔洞，(d)表示隧道豎直墻面由于遮擋物遮擋導致的小孔洞。進行輪廓線修復后，效果見圖8。

圖7 隧道點云孔洞

圖8 修復后的點云

圖8是修復后隧道三維激光掃描點云，其中紅色為原始點云，綠色為修復后的點云，(a)中可以看出道路和拱頂的兩個孔洞被修復，(b)表示隧道兩側的小孔洞已被修復，(c)表示隧道吊燈在修復孔洞時被保留，(d)為修復過后隧道三維激光掃描點云沿中軸線方向觀察的形狀結構，圖中可以看出修復后隧道結構細節依舊保持完整。

4 結束語

針對隧道三維激光掃描點云孔洞提出一種可行的修復算法。基于點云平面擬合所得的法向量，實現高精度中軸線的獲取，同時討論“在利用圓柱投影構建隧道三角網過程中，投影后兩個坐標尺度差異對三維構網結果的影響”，并給出統一坐標尺度的方法，通過計算一個合適的半徑，將圓柱投影所得的弧度值轉換成弧度所對應的弧長，確保坐標尺度一致。利用Delaunay平面構網方式對點云進行三角構網，最后利用殘缺斷面的非殘缺點到參考斷面輪廓線的距離，優選出用于修復的參考斷面，并利用該斷面的截斷點修復孔洞，再加入插值孔洞修復法，以保證所有孔洞都已被修復。實驗結果表明，本文的算法不僅能完美地修復隧道點云的孔洞，且能確保隧道細節不被破壞。然而該方法目前只適用于直線型隧道，對于彎曲型隧道，由于中軸線隨隧道方向的變化而改變[15]，從而導致截斷面之間點的互相填補不能靠簡單的直線平移進行修復，有待進一步的研究。

[1] 崔磊，張鳳錄.面三維激光掃描系統內外業一體化研究[J].北京測繪，2013(1):47-49.

[2] 吳曉軍，文飛，溫佩芝.視圖立體三維重建中的孔洞修復[J].計算機輔助設計與圖形學報，2012，24(12)：1606-1613.

[3] SONG P, WU X J, WANG Y. A robust and accurate method for visual hull computation[C]//Proceedings of IEEE International conference on Information and Automation. Los Alamitos: IEEE Computer Society Press,2009:784-789.

[4] SONG P, WU X J, WANG M Y .Volumetric stereo and silhouette fusion for image-based modeling[J].The Visual Computer,2010,26(12):1435-1450.

[5] 晏海平，吳祿慎，陳華偉.基于徑向基函數的散亂點云孔洞修復算法[J].計算機工程與設計，2014，35(4)：1253-1257.

[6] 劉越華，廖文和，劉浩.逆向工程中散亂點云的K鄰域搜索算法研究[J].機械設計與制造，2012(2):256-258.

[7] 孫殿柱，朱昌志，李延瑞.散亂點云邊界特征快速提取算法[J].山東大學學報(工學版)，2009，39(1)：84-86.

[8] 藍秋萍，洪超，林歡，等.三維點云中自動提取隧道特征線[J].測繪工程，2015，24(10)：1-4.

[9] 曹詩卉，亓迎川.基于最小二乘法平面擬合的點云法矢算法[J].計算機技術，2016，30(1)：41-43.

[10] GREEN P J.SIBSON R.Computing Dirichlet tessellations in the plane[J].Computer Journal,1978,21(2):168-173.

[11] BRASSEL K E,REIF D.A procedure to generate Thiessen polygons[J].Geographical Analysis,1979,11(3):289-303.

[12] MCCAULLAGH M T,ROSSC G.Delaunay triangulation of a random data set isarithmic mapping[J].The Cartographic Journal,1980,17(2):93-99.

[13] MAUS A.Delaunay triangulation and the convex hull of n points in expected linear time[J].Bit Numberical Mathematics,1984,24(24):151-163.

[14] 余杰，呂品，鄭昌文.Delaunay三角網構建方法比較研究[J].中國圖像圖形學報，2010,15(8)：1158-1167.

[15] 汪子豪. 從隧道三維點云中自動截取斷面輪廓的方法[J].水利與建筑工程學報，2015(2)：47-52.

[責任編輯：張德福]

Repairing tunnel 3D point cloud hole

QIAN Bozhi, LAN Qiuping

(School of Earth Science and Engineering, Hohai University,Nanjing 211100，China)

Due to the hole in tunnel 3D laser scanning point cloud, this paper designs a hole-filling method based on cross-sections. Tunnel is narrow and the cross-sections are similar, so the nearby intact cross-sections can be used to repair the hole. First of all, this paper changes 3D laser scanning point cloud into 2D by cylindrical projections and then comes up with a method to recognize the hole after constructing the Delaunay triangles. Finally, because the nearby cross-sections are similar, the hole is repaired with the nearby cross-sections. This method can not only repair the hole effectively, but also keep tunnel detail structure without destroyed. The effectiveness and reliability of this method are verified with the scanning data form a real tunnel.

3D point cloud; Delaunay triangle construction; hole identification; hole repair

10.19349/j.cnki.issn1006-7949.2017.03.010

2016-05-18

國家自然科學基金資助項目(41301406)；江蘇省自然科學基金資助項目(BK20130829)

錢伯至(1994-)，男，本科生.

P231

A

1006-7949(2017)03-0046-05

引用著錄：錢伯至，藍秋萍.隧道三維點云孔洞修復方法[J].測繪工程，2017，26(3)：46-50，55.