蘇 芬,余 銳
(廣州市城市規(guī)劃勘測設(shè)計研究院,廣州 510060)
?
基于切片點云中心的形變監(jiān)測
蘇芬,余銳
(廣州市城市規(guī)劃勘測設(shè)計研究院,廣州 510060)
摘要:針對建筑物特征表面的特性,提出基于橫斷面的切片點云中心的變形監(jiān)測方法,利用建筑物的某個側(cè)面點云作為約束面,建立約束面和用于提取切片基準(zhǔn)面之間的關(guān)系,根據(jù)基準(zhǔn)面與特征表面點云法向量之間的關(guān)系模型,確定最佳基準(zhǔn)面搜索方法。利用得到的基準(zhǔn)面逐層地確定切片點云,對每個切片點云進行分區(qū)段重心的求取,通過計算整個區(qū)段重心的中心獲取切片的重心。利用無變形的建筑物對點云精度進行評價,并將其用于指導(dǎo)隧道的變形監(jiān)測。通過對比不同期切片重心的變化實現(xiàn)隧道的變形分析。
關(guān)鍵詞:變形監(jiān)測;三維激光;整體最小二乘;切片點云;基準(zhǔn)面
隨著三維激光掃描技術(shù)的不斷應(yīng)用,其在變形監(jiān)測領(lǐng)域的應(yīng)用得到了不斷的重視,主要是由于該技術(shù)不像傳統(tǒng)的測量技術(shù)那樣,只能對某個點進行掃描,該技術(shù)可以獲取高精度海量的點云數(shù)據(jù),并且在掃描過程中,不受天氣影響,因此,基于該技術(shù)的變形監(jiān)測也得到了很多學(xué)者的研究,如文獻[1]利用三維激光掃描儀對巖石壁進行掃描,通過提取特征點的變形實現(xiàn)變形監(jiān)測;文獻[2]將三維激光掃描技術(shù)應(yīng)用在了地鐵隧道的變形監(jiān)測中,并對變形監(jiān)測結(jié)果的精度進行分析;文獻[3]在考慮隧道內(nèi)壁特點的情況下,利用三維激光掃描儀對其內(nèi)壁進行掃描,并根據(jù)點云的變化狀況確定隧道頂部變形的區(qū)域,相似的方法在文獻[4]也有應(yīng)用。針對隧道內(nèi)壁特殊的結(jié)構(gòu),文獻[5]采用矢量等距離方法對隧道內(nèi)壁的點云進行切片提取,利用提取的切片可以實現(xiàn)隧道不同部位的變形分析[6]??傊梢园l(fā)現(xiàn),目前有較多的學(xué)者根據(jù)切片點云的變化來實現(xiàn)形變監(jiān)測。針對點云切片的提取,文獻[7-8]首先實現(xiàn)了點云的隨機采樣處理,計算采樣后點云的密度,根據(jù)不同區(qū)域的點云密度進行切片點位置的確定,該方法的缺點是點云特征的反映不是特別明顯,況且采樣過程中采用的Christofides算法時間復(fù)雜度達到了O(n4),其用在切片點云獲取上效率較低。針對此缺點,文獻[9-10]利用移動最小二乘的方法來跟蹤切片的提取,該方法實現(xiàn)了精度較高的切片點云提取,但是對于有重復(fù)投影狀況的點云無法進行使用。獲取切片后,文獻[11]采用初始點和初始方向來判斷切片細(xì)化的精度,以此來提高切片細(xì)化的精度,該算法對切片細(xì)化程度較好,但是不能應(yīng)用在多重輪廓散亂點云中,同時也無法快速獲取建筑物的切片。
為了能夠充分利用建筑物表面的點云,本文根據(jù)建筑物表面的特性,將點云法向量和基準(zhǔn)面的關(guān)系引入到切片點云的提取中。通過分析建筑物某一個面的擬合參數(shù)與用于切片提取的基準(zhǔn)面之間的關(guān)系來確定基準(zhǔn)面的初始參數(shù);再結(jié)合基準(zhǔn)面和點云法向量之間的關(guān)系,通過基準(zhǔn)面不斷搜索的方法得到了最佳的基準(zhǔn)面;利用逐層切片和基準(zhǔn)面互換的方法不斷提取切片點云,對提取的切片點云進行分區(qū)段的重心獲??;利用各個區(qū)段的重心點來獲取整個切片的重心,通過對比不同期切片重心的變化實現(xiàn)建筑物變形分析。
1基準(zhǔn)面與約束面函數(shù)關(guān)系的建立
(1)
假設(shè)掃描的建筑物點云序列為(XiYiZi),i=1,2,…,n,將其進行歸一化處理,從而將點云序列的坐標(biāo)原點移到建筑物區(qū)域,則點云序列變?yōu)?/p>
(2)
假設(shè)通過坐標(biāo)原點的基準(zhǔn)面為a1x+b1y+
c1z=0,根據(jù)任意約束面垂直于通過坐標(biāo)原點的基準(zhǔn)面,即可得到關(guān)系式
(3)
從而建立了基準(zhǔn)面和任意約束面之間的函數(shù)關(guān)系。
2利用法向量確定最佳基準(zhǔn)面
(4)
由于基準(zhǔn)面與建筑物整體方向垂直,則基準(zhǔn)面方程與法向量之間的關(guān)系為
(5)
將式(3)與式(5)組合得
(6)
根據(jù)式(6)可解得a1,b1,c1的關(guān)系為
(7)
將式(7)代入基準(zhǔn)面a1x+b1y+c1z=0,可以得到n個基準(zhǔn)面的表達式
(8)
(9)
(10)
以此類推,可以得到n個基準(zhǔn)面與建筑物表面點云法向量之間的夾角{cosθj1cosθj2…cosθjii=1,…,n,j=1,…,n},其中,i表示的是基準(zhǔn)面數(shù),j表示的是法向量數(shù)。
計算每個基準(zhǔn)面與法向量夾角之和,如式(11)所示。
(11)
搜索最小的sum(cosθji)所對應(yīng)的基準(zhǔn)面
3切片點云提取及中心的計算
由第2節(jié)根據(jù)掃描點云不同面的法向量,確定最佳基準(zhǔn)面,再由最佳基準(zhǔn)面計算點云到基準(zhǔn)面的距離,設(shè)定初始切片到基準(zhǔn)面距離,便可實現(xiàn)初始切片的提取,再將初始切片視為初始基準(zhǔn)面,再根據(jù)第2節(jié)進行第2層最佳基準(zhǔn)面的搜索,以此類推,便可準(zhǔn)確地獲取整個建筑物的切片及切片的重心,如圖1所示。
圖1 切片重心提取所需變量
計算每個切片的重心,傳統(tǒng)方法都是直接采用點云擬合的方法獲得,該方法雖然簡單,但對于噪聲抑制不夠明顯。本文將首先對切片點云進行分段,分別擬合各段的重心,再利用各段的重心擬合得到切片的重心,具體步驟如下:
1)將切片點云投影到XOY平面上,計算得到假重心點O′,以該點為原點建立極坐標(biāo)系。
3)由于視角的原因,使得點云掃描的過程中存在不均勻,則不同區(qū)段內(nèi)點云數(shù)會有差別,為了體現(xiàn)公平區(qū)段點云擬合原則,設(shè)置區(qū)段點云數(shù)的最小閾值為q,若某區(qū)段內(nèi)點云數(shù)量Pm(m∈N,1≤m≤n)大于閾值q,則對該區(qū)段點云進行重心擬合得到區(qū)段中心為Om。
4)若區(qū)段內(nèi)點云數(shù)Pm小于閾值q,則由相鄰區(qū)段Om-1,Om+1中心擬合得到該區(qū)段中心。
5)根據(jù)n個區(qū)段中心擬合得到該層切片點云的中心O。
4實例分析
4.1精度指標(biāo)分析
利用三維激光掃描儀對3S雕塑建筑物進行兩次掃描,由于該建筑物沒有變形,因此,通過兩次掃描的變化量可以確定掃描精度,利用本文方法提取切片,結(jié)果如圖2所示。
采用本文方法提取兩次掃描切片的重心,并根據(jù)重心差值確定由掃描誤差引起的變化量,結(jié)果如圖3所示。
圖2 目標(biāo)點云處理前后灰度圖
圖3 不同方向變化量
由圖3可知,兩期切片重心x,y,z方向變化量都沒超過2 mm。點位中誤差最大值為2.714 mm,最小值為0.87 mm,小于2 mm占了75%,因此,點位中誤差絕大部分都小于2 mm,且大于2 mm的全部在2 mm左右。而該掃描儀的點位及點云模型精度也是在2 mm左右,因此,在已知激光掃描儀精度的前提下,通過本文方法可進行相關(guān)的變形監(jiān)測。
4.2隧道變形分析
某隧道掌子面由于剛剛開挖,各種防護都沒有跟上,所以急需采用各種監(jiān)測手段來觀測其變形情況,以便確保身處第一線的工作人員的安全。為了對其進行區(qū)域變形監(jiān)測,使用Leica ScanStation C10(見圖4)對該隧道掌子面進行掃描(見圖5),采集不同期的點云數(shù)據(jù),由于多期之間的坐標(biāo)系統(tǒng)不統(tǒng)一,為了實現(xiàn)多期之間點云數(shù)據(jù)的統(tǒng)一,需要進行基準(zhǔn)的確定,為了實現(xiàn)不同期點云基準(zhǔn)的統(tǒng)一,需要進行控制點的布設(shè),在不同掃描的過程中,將用于點云基準(zhǔn)統(tǒng)一的標(biāo)靶布設(shè)在控制點上,從而可實現(xiàn)多期掃描點云數(shù)據(jù)統(tǒng)一在第1期掃描中。并將第1期作為其基準(zhǔn)數(shù)據(jù),利用本文提出的方法分析隧道內(nèi)壁的變形狀況。由于該掃描距離為10 m,因此,其掃描的特征點位精度為2 mm, 掃描的模型表面精度為2 mm,該掃描精度只是適用于掃描距離較近的情況(一般在50 m以內(nèi))。根據(jù)隧道工程變形監(jiān)測控制標(biāo)準(zhǔn),區(qū)間跨度<8 m時,隧道區(qū)間的地表沉降允許位移控制值為30 mm,拱頂沉降允許位移控制值為30 mm,因此,該方法完全滿足對隧道區(qū)間跨度變形的監(jiān)測。掃描的隧道點云影像圖,如圖6(a)所示。由于掃描的過程中,受外界很多因素的影響,使得隧道內(nèi)壁會產(chǎn)生大量的噪聲,為了能夠準(zhǔn)確地對隧道內(nèi)壁的變形狀況進行分析,需要對噪聲進行剔除,噪聲剔除后的點云,如圖6(b)所示。
圖4 Leica ScanStation C10及其使用的標(biāo)靶
圖5 隧道掃描現(xiàn)場
圖6 隧道點云
采用本文方法提取隧道內(nèi)壁的切片點云,如圖7所示。由圖7可知,利用本文方法提取的隧道內(nèi)壁切片相對于原始隧道沒有發(fā)生扭曲或變化,從而驗證了本文方法提取切片的可行性。計算不同期切片點云的重心,結(jié)果如表1所示。
根據(jù)表1不同期的切片重心,計算兩次重心差值,結(jié)果如圖8所示。
由圖8可發(fā)現(xiàn),兩期切片重心,x方向變形量最大值接近10 mm,最小值為1 mm左右;y方向變形量最大值為5 mm左右,最小值為2 mm左右,z方向變形量最大值為3 mm左右,最小值為1 mm左右,同時點位中誤差最大值為11 mm,最小值為2 mm左右,該儀器在50 m以內(nèi)掃描的特征點位精度在2 mm以內(nèi),模型精度在2 mm以內(nèi),由4.1節(jié)分析可知,兩期差值在2 mm左右時,可認(rèn)為無變形。由圖8(d)可知,第1個切片和第2個切片變化量在2 mm左右,該變化量主要由激光掃描儀誤差引起的,可認(rèn)為第1個和第2個切片沒有變形,而從第3個切片開始,變化量逐漸增大,且變化量超過了2 mm,說明從第3個切片開始隧道內(nèi)壁出現(xiàn)了變形,且變形量逐步增大。而第1個切片到第7個切片屬于開挖時間不同的區(qū)域,第1個切片屬于掌子面開挖比較久的區(qū)域,其隧道內(nèi)壁比較穩(wěn)定,而第7個切片屬于掌子面剛開挖的區(qū)域,隧道內(nèi)壁不穩(wěn)定,會出現(xiàn)脫落等現(xiàn)象。因此,越是開挖比較久的區(qū)域變化量越小,隧道內(nèi)壁越穩(wěn)定,而越是剛開挖的區(qū)域變化量越大,隧道內(nèi)壁越不穩(wěn)定。
圖7 隧道內(nèi)壁切片點云
切片數(shù)123456-1.01750-0.96650-1.04890-1.08710-1.13660-1.41730一期切片重心-4.81878-5.53008-6.01948-6.41468-6.91148-7.466980.956010.974910.957110.945410.970011.06211-1.01830-0.96750-1.04590-1.09050-1.13130-1.42230二期切片重心-4.82058-5.52808-6.01558-6.41768-6.91648-7.472280.956710.976410.958710.946410.969011.06431
圖8 不同方向變形分析結(jié)果
5結(jié)論
本文根據(jù)隧道內(nèi)壁的特點,提出了提取隧道內(nèi)壁橫斷面的切片點云算法,應(yīng)用該算法可以在不用標(biāo)靶張貼的情況下自動提取隧道的基準(zhǔn)面,并利用該基準(zhǔn)面,通過逐層切片和基準(zhǔn)面的互換提取較為準(zhǔn)確的切片點云。通過實例,驗證該算法可以快速、準(zhǔn)確地自動獲取隧道內(nèi)壁的切片點云,對于提取的切片點云,本文利用分段區(qū)域重心擬合的方法計算每個切片的重心,并基于對無變形建筑物的切片變化狀況,在已知掃描精度的前提下,驗證該算法用于變形監(jiān)測的可行性,并得到了變形監(jiān)測的精度。同時,經(jīng)過不同期掃描某個正在開挖的隧道內(nèi)壁點云的對比分析,通過對不同部位隧道內(nèi)壁點云切片的分析,驗證了不同開挖狀況區(qū)域的隧道內(nèi)壁變形狀況,利用不同期相同切片重心的變化狀況及掃描點云精度指標(biāo),實現(xiàn)了隧道的變形監(jiān)測。
參考文獻:
[1]Monserrat O,CROSETTO M.Deformation measurement using terrestrial laser scanning data and least squares 3D surface matching[J].ISPRS Journal of Photogrammetry & Remote Sensing 63,2008:142-154.
[2]康志忠,托雷,王保前,等.基于三維激光掃描的地鐵隧道連續(xù)形變監(jiān)測數(shù)據(jù)處理軟件系統(tǒng)[J].測繪工程,2013,22(5):82-88.
[3]簡驍,童鵬.基于地面激光雷達技術(shù)的隧道變形監(jiān)測方法研究[J].鐵道勘察,2011(6):19-21.
[4]李健,萬幼川,江夢華,等.基于地面激光技術(shù)的隧道變形監(jiān)測技術(shù)[J].地理空間信息,2012,10(1):14-17.
[5]呂震.反求工程CAD建模中的特征技術(shù)研究[D].杭州:浙江大學(xué),2002:70-102.
[6]孟娜.基于激光掃描點云的數(shù)據(jù)處理技術(shù)研究[D].濟南:山東大學(xué),2009.
[7]劉云峰,柯映林.反求工程中的混合切片技術(shù)[J].計算機輔助設(shè)計與圖形學(xué)學(xué)報,2003,15(3):741-745.
[8]柯映林,王青.反求工程中點云切片算法研究[J].計算機輔助設(shè)計與圖形學(xué)學(xué)報,2005,17(8):1798-1802.
[9]龔友平,陳國金,陳立平.基于切片方法截面數(shù)據(jù)處理[J].計算機輔助設(shè)計與圖形學(xué)學(xué)報,2008,20(3):321-325.
[10] 蓋紹彥,達飛鵬,雷明濤,等.三維重構(gòu)中的雜亂點云排序問題研究[J].計算機與現(xiàn)代化,2003,98(10):33-35.
[11] 龔友平,金濤,童水光.截面切片數(shù)據(jù)的自動細(xì)化算法[J].浙江大學(xué)學(xué)報(工學(xué)版),2008,42(2):337-340.
[12] 浮丹丹,周紹光,徐洋,等.基于主成分分析的點云平面擬合技術(shù)研究[J].測繪工程,2014,23(4):20-24.
[13] 張賢達.矩陣分析與應(yīng)用[M].北京:清華大學(xué)出版社,2004.
[14] HOPPE H,DEROSE T,DUCHAMP T,et al.Surface reconstruction from unorganized points[C]//ACM Siggraph 1992,1992:71-78.
[15] Huang H,Li D,Zhang H,et al.Consolidation of unorganized point clouds for surface reconstruction[J].ACM Trans Graph,2009,28(5):176.
[責(zé)任編輯:劉文霞]
DOI:10.19349/j.cnki.issn1006-7949.2016.08.013
收稿日期:2015-05-10
作者簡介:蘇芬(1978-),男,工程師.
中圖分類號:TU196
文獻標(biāo)識碼:A
文章編號:1006-7949(2016)08-0057-08
The deformation monitoring based on the center of point cloud of slice
SU Fen,YU Rui
(Guangzhou Urban Planning Survey Design and Research Institute,Guangzhou 510060,China)
Abstract:This paper proposes a deformation monitoring method,based on the center of point cloud of slice.First,the relationship between the constraint plane and the datum is constructed on one side of the building.The search method of best datum is determined according to the relationship between the datum and normal vectors of point cloud.Second,the point cloud of slice is determined by using the datum,and the local center of each sect1 of the slice is computed.The center of slice is determined according to the center of the whole local center.The building of non-deformation is used to evaluate the point cloud accuracy.Finally,the deformation monitoring of building is done according to the comparison between different centers of point cloud of slice.
Key words:deformation monitoring;3D laser;total least squares;point cloud of slice;datum plane