林科 王紅 張英

摘 要： 目前，我國社會主義市場經(jīng)濟正處于飛速發(fā)展的階段，檢測行業(yè)在發(fā)展中也取得了顯著的成效，在對建筑項目的變形情況進行測量時，都是以傳統(tǒng)的測量方法為主，主要包括以下幾種方式：常規(guī)測量法、攝影測量和GPS測量等。雖然這些方式都得到了全面普及和推廣，但是在對數(shù)據(jù)進行采集時依舊存在一定的片面性，都是以單點式和影像式的方式為主。在使用前者對數(shù)據(jù)進行采集時，收集的資料比較局限，無法充分體現(xiàn)出檢測對象的整體情況。而采取后者的形式，雖然可以覆蓋比較大的面積，但是數(shù)據(jù)的準確度還有待提高。相比而言，在對建筑物表面的數(shù)據(jù)進行采集時，使用三維激光掃描術(shù)，可以確保數(shù)據(jù)的有效性和準確性，同時覆蓋的范圍也會更廣泛，對工作人員的要求比較低，并且可以充分滿足對數(shù)據(jù)的精確度要求，在此技術(shù)上得到建筑在形狀變化方面的數(shù)據(jù)。所以，對三維激光掃描技術(shù)進行研究具備長遠的意義，本文在研究的過程中，提出了切實可行的分析方法，在點云數(shù)據(jù)的基礎(chǔ)之上對隧道工程的形狀變化進行探究。

關(guān)鍵詞： 變形監(jiān)測;三維激光掃描技術(shù);點云壓縮

【中圖分類號】U452 ? ? 【文獻標識碼】A ? ? 【文章編號】1674-3733（2020）21-0032-04

1 三維激光掃描點云數(shù)據(jù)預(yù)處理

這些年來，伴隨著人類社會的不斷進步，全球工業(yè)發(fā)展水平正處于全面發(fā)展的階段，隨之誕生了大量的大型施工項目，并且建筑的整體框架和造型也日益繁瑣。建筑物因為受到外在重力和承載力的影響，形狀會發(fā)生一定的變化，嚴重的話甚至?xí)l(fā)生安全事故。所以，目前社會大眾都更加重視這些項目的安全程度。大型建筑物之所以形狀會發(fā)生變化，會受到很多因素的影響，因為自然環(huán)境在不斷變化，所以建筑物自身的承重力、整體框架和地理環(huán)境等都會破壞大型項目的穩(wěn)定性和安全性。因此不管是在建項目，還是已經(jīng)竣工的項目，都必須長時間進行監(jiān)督和測量，這樣就可以在第一時間對建筑物的變化進行掌握，并找到其中潛在的規(guī)律，一旦工程出現(xiàn)任何變化，就可以及時進行改善和優(yōu)化。以此確保大型工程項目的足夠安全，所以我們有必要對大型項目展開長期監(jiān)測和觀察。

使用三維激光掃描技術(shù)，是因為激光會對距離進行檢測，利用該原理，對檢測物體表面的立體坐標和紋理等信息進行觀察與記錄，同時通過先進的科學(xué)技術(shù)，對數(shù)據(jù)進行準確有效地處理，可以對多種不同的測量技術(shù)進行綜合。通過三維激光掃描技術(shù)可以在極短的時間內(nèi)對坐標信息進行采集，采集數(shù)量數(shù)以萬計，對于數(shù)量如此之多的坐標點來說，學(xué)術(shù)界稱之為“點云”。通過三維激光掃描技術(shù)，可以得到點云，這樣一來，就可以從不同的角度對檢測物體的形狀進行準確地描述。點云包含大量的數(shù)據(jù)和信息，并且排序是毫無規(guī)律的，其中包含很多不同的特征，比如坐標和系統(tǒng)不匹配，坐標易變形等特征。因此在實施檢測之前，要對相關(guān)數(shù)據(jù)進行配對，消除噪音，并對坐標系統(tǒng)進行歸類，同時對點云進行壓縮等，從而形成對應(yīng)的文件，為后期的分析奠定堅實的基礎(chǔ)，換句話說，就是對點云數(shù)據(jù)實施有效地處理。

1.1 三維激光掃描系統(tǒng)

1.1.1 系統(tǒng)分類

（1）地面型激光掃描系統(tǒng)

通過激光脈沖，可以準確有效地獲取地面物體的三維坐標，并對相關(guān)坐標進行測量，這就是地面型激光掃描系統(tǒng)。該系統(tǒng)具備很多不同的特征，比如覆蓋面廣、運行速度快、準確度高等等。

當(dāng)前，市場上有很多這種類型的掃描儀，用戶在使用的時候可以結(jié)合自身的需求進行選擇。

（2）機載型激光掃描系統(tǒng)

這種掃描系統(tǒng)是建立在無人機或直升機的基礎(chǔ)之上，該系統(tǒng)包括很多不同的部分，比如GPS系統(tǒng)、采集器和飛行引導(dǎo)系統(tǒng)等。通過該系統(tǒng)可以在最短的時間內(nèi)獲得大量的三維坐標數(shù)據(jù)和信息。

（3）便攜式激光掃描系統(tǒng)

這種掃描系統(tǒng)非常輕便、自身可以進行定位，并且速度非常靈敏，在對室內(nèi)坐標進行掃描時，可以使用該系統(tǒng)，一般情況下在對一些體積比較小的物體進行掃描時，都會使用這種系統(tǒng)。

1.1.2 系統(tǒng)工作原理

使用三維激光掃描儀進行掃描的時候，主要是按照以下標準來進行的，首先經(jīng)過相應(yīng)的發(fā)射裝置，將激光傳達到物體表面之后再按照原來的路線返回原地，運用時間計數(shù)器，從光束發(fā)射開始計算，到最后的接收，計算其中耗費的時間，緊接著將儀器到測量點的距離計算出來，在對距離進行測量時，該掃描儀不僅可以記錄水平方向上的角度，還可以記錄垂直方向上的角度，通過測量的距離和角度進行計算，得到相應(yīng)的三維坐標。

1.1.3 系統(tǒng)工作流程

系統(tǒng)在運行的時候，主要包含兩個階段，一個是外業(yè)工作階段，一個是內(nèi)業(yè)工作階段。

1.1.3.1 外業(yè)工作流程

（1）首先選好適當(dāng)?shù)臏y量點和標靶點，此時需要注意的是，在選取測量點的時候需要在掃描儀可以覆蓋的范圍之內(nèi)，要和前一個測量點的各個標靶保持通視，另外要保證測量站點的穩(wěn)定性。在選擇標靶點和測量站點的時候，要注意兩者不能位于同一條水平線上，并且要和測量站點保持通視，在選取標靶設(shè)置點的時候，最少要超過3個。

（2）安裝測量腳架和三腳架，將掃描儀有效放置之后，再對準中心保持水平，將儀器的高度進行測量。

（3）安裝標靶，隨后設(shè)置對應(yīng)的三腳架，再找到中心點和水平線的位置，將標靶的高度進行測量。

（4）將掃描儀和電源進行連接，如果掃描儀運行時間較長，則可以使用移動電源。

（5）將相關(guān)數(shù)據(jù)輸入進去，并在掃描儀中輸入相應(yīng)的儀器高度、標靶高度和測量站點的坐標等，根據(jù)掃描儀內(nèi)部現(xiàn)有的軟件對儀器進行控制，正式開始掃描。

1.1.3.2 內(nèi)業(yè)工作流程

在使用三維激光掃描儀進行掃描的時候，業(yè)內(nèi)工作流程占有極大的比例，且有著至關(guān)重要的作用。通過業(yè)內(nèi)工作流程，可以對大量的數(shù)據(jù)和信息進行處理和調(diào)整，其中包括點云輸出、消除噪音和壓縮等步驟。

1.2 點云導(dǎo)出、配準和去噪

三維激光掃描儀可以采集到大量數(shù)據(jù)，不過其中部分數(shù)據(jù)對于多余數(shù)據(jù)，不會給監(jiān)測分析帶來幫助，有時候甚至?xí)绊懛治龅捻樌M行。倘若不加區(qū)分地分析全部數(shù)據(jù)，便會增加數(shù)據(jù)誤差。比方說點云中存在超出掃描預(yù)計氛圍的點，或者存在監(jiān)測對象之外的點。此外，掃描儀的激光束會或多或少產(chǎn)生離散效應(yīng)，此時掃描儀接收器便會捕捉到大量雜亂反射光，這些反射光將形成大量噪聲點，并聚集在監(jiān)測對象邊緣。如果噪聲點過多的話，點云數(shù)據(jù)中便會包含龐大的冗余數(shù)據(jù)，導(dǎo)致監(jiān)測分析效率變低、準確度下降。所以在分析點云數(shù)據(jù)之前，必須有效去噪，減少冗余數(shù)據(jù)量。

1.3 點云坐標系歸化

在以各階段點云數(shù)據(jù)為基礎(chǔ)開展監(jiān)測分析的時候，倘若各階段拼接方案存在差異，將使得各階段點云數(shù)據(jù)坐標系出現(xiàn)明顯區(qū)別，也就是說出現(xiàn)相同監(jiān)測對象多個階段點云數(shù)據(jù)存在不止一個坐標系的狀況。此時應(yīng)當(dāng)進行點云坐標系歸化，也就是旋轉(zhuǎn)各階段點云數(shù)據(jù)的坐標系，讓各階段點云數(shù)據(jù)處在相同坐標系之下。本文主要通過以公共控制點為基礎(chǔ)的空間坐標系旋轉(zhuǎn)方法來進行點云坐標系歸化。

1.3.1 空間坐標系旋轉(zhuǎn)算法

假設(shè)原本的空間坐標系是OXYZ，在維持原點不變的情況下，將X Y Z軸旋轉(zhuǎn)至另一個空間坐標系OXYZ中之后，角度Q就是X Y Z軸和OX.OY.OZ軸之間的角度，α2β2γ2就是Ob軸和OX.OY.OZ軸之間的角度。α3β3γ3就是QZ軸和OX，OY，OZ軸之間的角度。坐標系OXYZ里面的P點坐標是（X，Y，Z），坐標系OX'Y'Z'里面的P點坐標是（X1，Y1，Z1）。

1.3.2 基于3個公共控制點的坐標系歸化

從各階段點云數(shù)據(jù)內(nèi)挑選出3個共有控制點，其中隧道兩端分別有2個控制點，剩下的1個控制點可自由挑選。以某階段點云數(shù)據(jù)坐標系歸化的計算過程為例，首先挑選出D1、D2與D3這3個共有控制點，隨后按照尺度因子1對隧道控制點坐標系OXYZ進行旋轉(zhuǎn)，產(chǎn)生全新的坐標系OX1Y1Z。

其中D1充當(dāng)全新坐標系的原點，D1、D2之間的連線則充當(dāng)全新坐標系的X軸走向，方向指向D2，設(shè)D1、D2兩點間的距離是S1。以D3為基礎(chǔ)明確XOY平面方位，接著形成全新坐標系的Y軸。設(shè)D1、D3兩點距離在X軸上的分量是S2。設(shè)D1、D3兩點距離在Y軸上的分量是S3。最后按照右手坐標系原則以及X軸、Y軸的方向，就能夠明確Z軸走向。

以D1、D2、D3三個共有控制點在兩個空間坐標系的坐標值為基礎(chǔ)，借助坐標系旋轉(zhuǎn)算法來計算出坐標系的轉(zhuǎn)換參數(shù)，然后旋轉(zhuǎn)全部點云數(shù)據(jù)的坐標系，從而實現(xiàn)坐標系歸化，為之后的監(jiān)測分析提供便利。

1.3.3 基于多個公共控制點的坐標系歸化

倘若各階段點云數(shù)據(jù)擁有3個以上的共有控制點，那么可以在原坐標系各點關(guān)系的基礎(chǔ)上進行空間坐標轉(zhuǎn)換，比方說點之間的距離、角度等等。然后基于全部限制條件完成附有限制條件的最小二乘平差，從而得出坐標系旋轉(zhuǎn)參數(shù)。最后計算出各個共有控制點旋轉(zhuǎn)前后的詳細坐標，這樣就能夠顯著降低坐標系轉(zhuǎn)換帶來的精度損失。

第一步：旋轉(zhuǎn)各個共有控制點的坐標系，采用的旋轉(zhuǎn)方法是以3個控制點坐標系規(guī)劃算法為基礎(chǔ)制定的，確保旋轉(zhuǎn)前后坐標精度不出現(xiàn)明顯變化。

第二步：以各個共有控制點的旋轉(zhuǎn)參數(shù)解算為基礎(chǔ)，并按照各個共有控制點在全新坐標系中的坐標以及標輪點坐標計算出旋轉(zhuǎn)參數(shù)。

1.4 點云數(shù)據(jù)壓縮

通過點云數(shù)據(jù)，可以對物體表面的各個點的信息進行詳細地記錄，其中包括每個點的坐標等，如果掃描的密度不同，那么物體表面經(jīng)過掃描之后，每個點的信息量也是截然不同的，在KB級和G級之間變動。使用G級技術(shù)對數(shù)據(jù)進行計算時，會導(dǎo)致計算機系統(tǒng)面臨極大的問題，并且會對工作效率造成嚴重影響，甚至?xí)档蛿?shù)據(jù)的有效性和準確性。所以我們需要結(jié)合工程的實際情況，在對點云數(shù)據(jù)進行處理的時候，如果出現(xiàn)信息量過大的現(xiàn)象，就要對數(shù)據(jù)進行壓縮處理。當(dāng)前，在對點云數(shù)據(jù)進行壓縮處理的時候，出現(xiàn)了很多不同類型的計算方法，但是從檢測工程的角度來說，有一些方法還存在一定的片面性，無法滿足檢測的實際需求。本文在研究的過程中，對過去常用的壓縮方法進行了詳細的說明，并對八叉樹分割方法和加權(quán)分層的方法進行比較，研究證實，因為加權(quán)分層的壓縮方式可以對壓縮的方向和比例進行把控，所以在對隧道等大型工程進行檢測時，可以使用這種方法。

1.4.1 基于八叉樹分割的點云壓縮算法

在對數(shù)據(jù)進行處理時，為了充分提高數(shù)據(jù)處理的速度和效率，可以將這些數(shù)據(jù)進行劃分，成為很多不同的立方體盒子，在對數(shù)據(jù)進行搜集時，可以將各個盒子的序號作為搜索的關(guān)鍵詞，將其中對應(yīng)的數(shù)據(jù)進行存儲，這樣一來就可以充分降低工作的難度和工作量。我們可以將點云的盒子看成是一個包圍盒，其中包含點云的所有信息和數(shù)據(jù)。

1.4.2 基于加權(quán)分層的點云壓縮算法

本文在研究的過程中，使用加權(quán)分層的方法，對點云進行壓縮，并將其分成不同的層次，在每一層搜索有效數(shù)據(jù)，這些數(shù)據(jù)之間的距離都滿足一定的比例和范圍，從而提高數(shù)據(jù)的準確性和有效性，同時可以將物體表面的形狀和特點進行有效保存，這種計算方法的有效性比較高，可以對壓縮的比例進行控制。

這種計算方法是將相關(guān)數(shù)據(jù)進行劃分，成為不同的層次，緊接著對每一層的點云進行連接。將測量物體進行劃分，從而形成一定規(guī)格的點云數(shù)據(jù)，呈片狀的形式，這就是點云切片。隨后將這些切片進行逐步累加，從而形成相對立體的物體形象。結(jié)合物體自身的特征，比如垂直方向和水平方向等，分成不同的層次，以此獲得充足的數(shù)據(jù)和信息，緊接著將數(shù)據(jù)進行投射，在垂直面上獲得物體的輪廓線，再在輪廓線的周圍提取相應(yīng)的數(shù)據(jù)。

第一次分層：特征方向是點云數(shù)據(jù)排序的過程，之后設(shè)計成包絡(luò)盒，以固定間隔h在這個方向上進行分層，（Xmax-XminYmax-Yminh）是每一層的容量，（Z，Z+h）是各層Z坐標區(qū)間，詳細記下不同層的各個下屬點。

第二次分層：在第一步完成之后，首先需要找出距等高線臨界點范圍內(nèi)的點，從各層每個包圍盒的點云數(shù)據(jù)之中，再與截面的方式，設(shè)△為臨界值，（Xmax-XminYmax-YminZ±1）是該截面區(qū)域點的空間范圍，點云成果也是各層提取的壓縮值。

切片厚度設(shè)置：切片厚度在提取時不能越過臨界值，點云數(shù)據(jù)的密度也與切片厚度有關(guān)。切片厚度在正常情況下是和點云密度成反比?？蓪設(shè)置為平均算法誤差。各個丹麥你的厚度臨界值=（eM）/（nLi）M用加權(quán)的方式來確定總的點個數(shù)，第一次分層之后第i層的點個數(shù)為Li，之后第一次分的層數(shù)是n。

1.4.3 兩種壓縮算法的對比分析

以上述研究成果能夠了解到，點云壓縮計算的方式較為復(fù)雜，而且如果數(shù)據(jù)量較為龐大的話，各領(lǐng)域點個數(shù)、標準偏差限差難以明確，另外，其中的壓縮比例也比較大，壓縮方向很難對整體壓縮形成操控，壓縮之后，在進行斷面切割的過程中，很容易便會使得斷面點數(shù)過低。而擬合誤差就會因為此種情況的出現(xiàn)而加大。當(dāng)進行經(jīng)驗測試壓縮后，各斷面點數(shù)如果較多，那么所得結(jié)果才可以幫助到之后的研究工作更好開展，才可以達成擬合精度標準。經(jīng)過對比，針對以加權(quán)分層作為基礎(chǔ)的點云壓縮方式而言，整體計算的方式較為簡便，對內(nèi)存率占用較少，而且在進行軟件編程的過程中，也會提供很多便利，壓縮過程的方向、比例均可以實際需求進行設(shè)定，在規(guī)模較大的項目中可以發(fā)揮巨大作用，比如說精簡隧道數(shù)據(jù)等等。至此，本篇文章為了使得擬合精度得到保障，探測變形水平有所強化，所以選擇加權(quán)分層的壓縮方式進行使用。

1.5 點云數(shù)據(jù)斷面切割

1.5.1 基于斜率收斂算法的斷面切割

對于斜率收斂算法而言，其計算原理便是以隧道斷面切割作為基礎(chǔ)，在數(shù)據(jù)后處理的過程中，假設(shè)難以直接得到隧道中心線，那么就要對中心線進行擬合，為了能夠保證中心線整體擬合的準確度，在截取斷面的過程中，一定要對點云厚度值以及斜率線限差值。如果算法可以實現(xiàn)迭代收斂，切面厚度值以及斜率限差值和斷面擬合誤差之間便會成正比，其精度則成反比，反之亦然。在測試后得知，如果點云厚度不超過3毫米，斜率限差等于0.1時，各斷面擬合半徑誤差小于8毫米。至此，若算法可以實現(xiàn)迭代收斂，那么就需要盡可能低的設(shè)置切面厚度以及斜率限差。

1.5.2 斷面點云數(shù)據(jù)擬合

對于隧道斷面而言，其形狀類型較多，以空間圓擬合形舉例，作為球、平面之間的交集，首先要對平面進行擬合，其次擬合球，最后，平面和球兩方程的交集則組成了斷面圓形方程。

以點云數(shù)據(jù)作為基礎(chǔ)，對平面進行擬合的方式眾多，像最小二乘法、特征數(shù)值法等等。在這之中，平面法向量作為單位向量的話，適用特征數(shù)值法，利用點至面距離平方以及最小條件標準得出ax+by+cz=d這一平面方程，此種方式將方向軸中的誤差全部考慮在內(nèi)?？墒?，對于點云數(shù)據(jù)中的噪聲點并沒有涉及，本篇文章所借助的最小二乘平面擬合法可以將三個坐標方向全部考慮在內(nèi)，此種方式是特征值法的衍生算法來。

1.5.3 斷面分析對象的提取

對于整體的分析對象提取而言，在進行傳統(tǒng)測量模式的時候，正常情況下，借助監(jiān)測點針對回歸平面進行擬合，從而得到各期平面重心位置偏移，這種方式是比較慣用于整體分析的，針對分析三維激光掃描技術(shù)，本篇文章選擇的研究辦法是以各期點云數(shù)據(jù)作為基礎(chǔ)，以隧道擬合中心作為標準，對斷面進行截取，之后擬合所得的斷面點云。得到斷面擬臺方程之后，需要對比檔次有差異，但斷面無差異情況下的半徑，之后得出差值，在這個基礎(chǔ)上，開展研究斷面收斂變形。

對于提取局部分析點而言，在以往的測量過程中，局部分析實際上就是比較各期監(jiān)測點的坐標，本篇文章進行斷面截取的辦法是以一個隧道擬合中心線位置作為基礎(chǔ)，從而對不同期數(shù)據(jù)所進行的。以零度作為斷面水平方向，以逆時針開展360°的截取，以a度作為截取點位的角度。得到斷面之后，對比不同期但是角度值一致的斷面點坐標，從而得出其中間距，最終得出此角度形變量。

1.6 精度影響定性與定量分析

1.6.1 定性分析精度影響因素

在整個數(shù)據(jù)處理過程中，影響后續(xù)數(shù)據(jù)處理精度的因素有：

（1）后續(xù)數(shù)據(jù)的分析的準確度與點云數(shù)據(jù)的密度有著必然聯(lián)系，斷面切割時的斷面的點個數(shù)與點云分布密度是相關(guān)聯(lián)的，擬合的準確度也和點數(shù)的數(shù)量存在聯(lián)系，也可以說密度與擬合的準確度成正比。

（2）之后數(shù)據(jù)分析的準確度與點云拼接存在聯(lián)系，兩站之間的標靶需要進行配對，無縫拼接在配對時很難做到，因此拼接也會發(fā)生偏差，這在測站的點云數(shù)據(jù)拼接時是常見現(xiàn)象。

（3）之后數(shù)據(jù)分析的準確度也會受到點云去噪的干擾，部分有效的點云數(shù)據(jù)會被點云去噪給去除掉。擬合質(zhì)量會隨著去除的有效數(shù)據(jù)的數(shù)量而衰減，特別是數(shù)據(jù)分析準確度這方面就是會過度下滑。

（4）之后數(shù)據(jù)分析的準確度與坐標系歸化有關(guān)，掃描時的標靶點架設(shè)，其中會有一些偏差，之后這些偏差還會借助標靶點和控制點一同參與坐標系歸化。

（5）之后數(shù)據(jù)分析的準確度會受到點云壓縮的干擾，許多的有效點云數(shù)據(jù)會在點云壓縮的過程中舍去掉，之后數(shù)據(jù)分析的準確度會下降，斷面的擬合效果會受到干擾，而且所受影響還沒辦法準確化。

（6）之后數(shù)據(jù)分析的準確度會受到斷面切割厚度的干擾，斷面擬合準確度會被斷面切割厚度所干擾，在斷面切割擬合時還會有著擬合偏差。

1.6.2 定量分析斷面點云厚度

本文的斷面厚度測試研究使用的是中等密度的點云數(shù)據(jù)?？晒烙嬎淼罃嗝孀髠?cè)是0度，到了頂部為90度，最右側(cè)是180度，之后每隔4度就是使用分析點和設(shè)計斷面來對比分析。對擬合斷面和設(shè)計斷面進行對比分析，針對1mm到9mm厚度。

分析1mm至9mm的厚度切片對比研究：將設(shè)計斷面和擬合的隧道斷面進行對比，斷面各個方向上的最大差值與厚度成正比。最大誤差增值在3mm內(nèi)，這要在厚度小于8mm時增加厚度才能實現(xiàn)，這也沒有超過隧道管徑測量的誤差范圍;隧道管徑測量的誤差范圍3mm會在厚度超過8mm時被超出。直線型且斷面呈圓或圓拱形的隧道收斂變形分析的方面，主要是針對其將點云切片厚度控制在8mm內(nèi)，使得之后斷面點云擬合誤差不會超出隧道結(jié)構(gòu)變形的誤差范圍，這樣也符合隧道收斂變形分析的準度要求。

1.7 本章小結(jié)

點云數(shù)據(jù)的預(yù)處理系統(tǒng)是本章節(jié)的主要研究內(nèi)容。優(yōu)先介紹了坐標系的歸化算法，與公共控制點有關(guān)，其中穿插了許多控制點進行參數(shù)轉(zhuǎn)換。本章節(jié)使用了各期的對比分析和兩種歸化算法，這也證明了歸化方式想要得到到更準確的結(jié)果就需要對多個公共控制點進行轉(zhuǎn)換參數(shù)，坐標系的全體轉(zhuǎn)換更直觀的完成。之后還用加權(quán)分層的點云壓縮進行介紹，還采取了兩種壓縮方式進行對比，即加權(quán)分層點云壓縮和八叉樹分割點云壓縮，從對比的數(shù)據(jù)可知加權(quán)分層點云壓縮方法更加適合隧道工程的點云數(shù)據(jù)壓縮。還將隧道斷面點云截取方法以斜率收斂的方式帶入，當(dāng)前方法使用斜率限差的迭代計算方式，將隧道的中心線完整擬合出來，對隧道斷面進行切分，使得斷面點云數(shù)據(jù)的價值得到保證。

2 三維激光掃描隧道超欠挖計算與檢測

本章節(jié)詳細陳述了各階段擬合斷面和設(shè)計斷面的對比，具體來說就是經(jīng)由點云數(shù)據(jù)預(yù)處理之后，能夠得到的斷面點云數(shù)據(jù)進行隧道的超欠挖計算與檢測。在正式計算前，應(yīng)當(dāng)首先詳細分解設(shè)計中心線，從而產(chǎn)生可以被計算機正確識別的中心線文件。

2.1 隧道橫斷面提取

提取出隧道橫斷面的方法多種多樣，本文主要介紹的是切片式提取法與旋轉(zhuǎn)式提取法。因為在提取的時候涉及隧道中軸線法平面方程的計算和曲面擬合算法，所以接下來首先介紹這兩個計算方法。

2.1.1 隧道中心線的法平面方程計算

計算隧道超欠挖面積的過程中，通常會以截取出的隧道橫斷面充當(dāng)參考依據(jù)。截取隧道橫斷面之前必須詳細掌握對應(yīng)截取部位的空間姿態(tài)情況，確保隧道走向和隧道橫斷面截取方向呈正交關(guān)系。隧道中心線法平面方程計算的步驟如下所示。

（1）讀取平曲線、豎曲線、設(shè)計斷面

讀取設(shè)計中心線平曲線文件，把隧道起點、隧道終點等設(shè)計中心線平曲線要素點里程存儲到數(shù)組L1中，將設(shè)計中心線平曲線要素點平面坐標依次存儲到數(shù)組X1、Y1中，把設(shè)計中心線平曲線半徑存儲到數(shù)組中，并將線型存儲到數(shù)組RIH中，各個數(shù)組的下標相互對應(yīng)。

讀取設(shè)計中心線豎曲線文件，把變坡點等設(shè)計中心線豎曲線要素點里程存儲到數(shù)組中，將設(shè)計中心線平曲線要素點高程值存儲到數(shù)組中，將設(shè)計中心線豎曲線半徑存儲到數(shù)組中，并將切線長存儲到數(shù)組RIV中，各個數(shù)組的下標相互對應(yīng)。

讀取設(shè)計斷面文件，把斷面里程依次存儲到數(shù)組Di中。

（2）計算任意設(shè)計斷面與中心線交點的平面坐標

在設(shè)計斷面里程的基礎(chǔ)上確定距離最近的兩個設(shè)計中心線平曲線要素點，它們的里程分別是LI和Li+1。按照SI確定Di的設(shè)計中心線平曲線類型，從Di的設(shè)計中心線平曲線區(qū)間差異來看，主要存在下述這些情況，在此基礎(chǔ)上進行設(shè)計斷面和中軸線交點平面坐標的計算。

2.1.2 曲面擬合算法研究

本文需要進行二次曲面擬合計算，通常使用的計算方法是附有限制性條件的最小二乘曲面擬合算法。通過多次計算來逐步獲取貼近真實情況的數(shù)值，從而實現(xiàn)最優(yōu)化擬合狀態(tài)，也就是說擬合得出的曲面參數(shù)和實際曲面差異不大。

附有限制性條件的最小二乘曲面擬合算法指的是基于曲面方程的點云數(shù)據(jù)擬合模型進行計算，并附加限制性約束條件。

3 三維激光掃描隧道變形分析

本章詳細陳述了各階段點云數(shù)據(jù)相同切割斷面的相互對比，具體來說就是經(jīng)由三維激光掃描點云數(shù)據(jù)預(yù)處理之后，能夠得出隧道橫斷面點云數(shù)據(jù)以及相關(guān)文件，從而針對性剖析隧道橫斷面的變形情況。

3.1 隧道變形分析與工程應(yīng)用

3.1.1 隧道變形分析方法

隧道變形，具體來說是隧道長期承受地面及建筑物荷載、土體擾動、附近施工干擾、隧道結(jié)構(gòu)施工干擾等不利影響而產(chǎn)生的變形。本文借助多個期次隧道點云數(shù)據(jù)來詳細比較和剖析隧道同一部位，從而全方位掌握隧道橫斷面變形情況。本文所使用的分析方法主要有兩種，第一種是整體變形分析法，第二種是局部變形分析法。

3.1.2 工程實例分析

本文將新疆克輪輸氣管線鹽水溝隧道作為分析對象，以該隧道的掃描點云數(shù)據(jù)為基礎(chǔ)進行分析比較。鹽水溝隧道是一條輸氣管線隧道，屬于西氣東輸工程的重要部分，坐落在庫車縣和拜城縣交匯處的山區(qū)之中。

本文以鹽水溝隧道的實測數(shù)據(jù)為基礎(chǔ)詳細分析，首先提取出其中20m的局部點云數(shù)據(jù)，并處理掉噪音，隨后導(dǎo)出三個隧道控制點，這三個控制點充當(dāng)兩期點云數(shù)據(jù)坐標系的共同控制點。

結(jié)論：由于我國沒有統(tǒng)一的隧道收斂變形指標規(guī)范，因此本文選擇按照上海地鐵限制荷載影響下結(jié)構(gòu)變形不超出D/1000的標準，其中D指的是設(shè)計直徑。按照鹽水溝隧道設(shè)計直徑3m來計算，能夠得知管徑測量的中誤差上限是3mm。而上述表格的14個隧道橫斷面半徑擬合中誤差都處于3.94mm到7.24mm范圍內(nèi)，以誤差傳播定律為指導(dǎo)進行計算，可以看出隧道直徑擬合中誤差變化區(qū)間是5.57mm到10.23mm范圍內(nèi)，高于管徑測量中誤差的上限值。由此來看，借助三維激光掃描技術(shù)檢測隧道橫斷面的收斂變形還需繼續(xù)深入研究。

不過從相關(guān)數(shù)據(jù)可以得知，隧道半徑平均變化量是-1.86mm，并且所有隧道斷面中有72%的斷面半徑差值是負數(shù)，即便剩余部分差值為正數(shù)的斷面，其差值也沒有大于2.5mm。倘若調(diào)低標準，那么仍然能夠看出，這些斷面大部分為收斂形態(tài)。

3.2 隧道三維模型重建與可視化分析

三維模型重建技術(shù)指的是以掃描儀器獲取到的點云數(shù)據(jù)為基礎(chǔ)，通過對數(shù)據(jù)進行預(yù)處理，隨后得到建模必需的有效點云數(shù)據(jù)，隨后建立隧道的曲面及實體造型，最終產(chǎn)生三維模型[45，46]。市面上應(yīng)用比較普遍的以點云數(shù)據(jù)為基礎(chǔ)建立三維模型的平臺多種多樣，其中最典型的就是3d max軟件和Cyclone軟件。3d max軟件是Discreet公司推出的多功能3D制作軟件，既可以制作實體造型，又可以渲染與開發(fā)動畫。Cyclone軟件通常使用在處理三維激光掃描系統(tǒng)采集到的數(shù)據(jù)上，可以對點云數(shù)據(jù)進行三維建模等操作。本文將Cyclone軟件切割得到的斷面文件導(dǎo)入至3Dmax軟件，隨后完成建模，最終利用OpenGL平臺制作出3D模型瀏覽程序。

4 隧道工程三維激光掃描數(shù)據(jù)處理原型系統(tǒng)設(shè)計與實現(xiàn)

4.1 需求分析

將點云數(shù)據(jù)以拼接、截取、建模等操作是點云數(shù)據(jù)處理軟件基礎(chǔ)功能，像Leica Cyclone，PointCloud等軟件都是用于數(shù)據(jù)的前期整理，除去精細化的變形監(jiān)測分析應(yīng)用，其他都你能使用與點云數(shù)據(jù)在計算機中自動處理，在坐標系的角度來歸化、數(shù)據(jù)壓縮、斷面切割等許多功能。那么三維激光掃描數(shù)據(jù)處理和模型演示系統(tǒng)就變得相當(dāng)有價值。隧道整體的收斂變形數(shù)據(jù)能夠及時取得，那么檢測效率必然也會得到提升。

4.2 系統(tǒng)的實現(xiàn)

4.2.1 點云數(shù)據(jù)預(yù)處理的實現(xiàn)

以坐標系歸化、斷面切割、壓縮等功能為目的實現(xiàn)點云數(shù)IS，設(shè)置了“PreProces&cpp”和“PreProcess.h”的類文件在系統(tǒng)之中，這些類文件，以變量的手段來對預(yù)處理進行詮釋，此外，還有坐標系歸化類、斷面切割類、點云壓縮類都是設(shè)置而成。

4.2.2 超欠挖面積計算與檢測的實現(xiàn)

對土方量和超欠挖面積進行計算。將“Ca_Scpp”和“Ca_S，h”的類文件設(shè)置在系統(tǒng)中。以變量的方式對超欠挖計算進行詮釋，還有平豎曲線文件讀取、設(shè)計斷面面積計算、隧道斷面截取等功能函數(shù)進行設(shè)置而成。

4.2.3 變形分析和模型演示的實現(xiàn)

本部分是為了將三維模型的讀取和演示還有隧道斷面的變形分析功能具體實現(xiàn)，局部分析法和整體收斂分析法是兩種常用分析方式。“Contactcpp”和 “Contacth”類文件在設(shè)置之后便能實現(xiàn)該功能，變量由該類自身來詮釋，這也包含了模型讀入、數(shù)據(jù)對比分析等操作函數(shù)。附錄3是兩期數(shù)據(jù)對比分析的源碼。

4.3 系統(tǒng)設(shè)計功能特點

VC=++和OPENGL平臺是本系統(tǒng)的開發(fā)依據(jù)，就以當(dāng)前國內(nèi)隧道測量特點和用戶操作習(xí)慣來設(shè)置，可用于協(xié)助三維激光掃描系統(tǒng)對隧道測量整合，主要功能是下列內(nèi)容：

點云三維顯示：使用三維的方式來表示三維激光掃描系統(tǒng)取得的隧道點云測量數(shù)據(jù)，再使用交互瀏覽操作等。

坐標系統(tǒng)一：坐標歸化各期的三維激光掃描點云數(shù)據(jù)，方便日后再進行數(shù)據(jù)分析。

點云數(shù)據(jù)壓縮：壓縮龐大的點云數(shù)據(jù)量，后續(xù)的數(shù)據(jù)分析需要壓縮后的數(shù)據(jù)。

斷面處理：斷面點云通過固定時間點從點云數(shù)據(jù)中截取，斷面進行超欠挖計算和收斂變形分析都要以此設(shè)計。

三維模型演示：交互式瀏覽讀取的3ds模型文件。

結(jié)論

本文的研究的主要內(nèi)容就是三維激光掃描技術(shù)，再結(jié)合內(nèi)外行業(yè)技術(shù)和數(shù)據(jù)的分析，從而獲得了以下結(jié)論：

坐標系在點云數(shù)據(jù)處理中存在著坐標不一致的問題，將公共控制點的坐標系統(tǒng)的歸化直接呈現(xiàn)出，使用控制點的平差計算結(jié)合多個方法，避免出現(xiàn)準度在坐標系歸化中存在偏差。以隧道斷面中圓形的為例，將最小二乘平面擬合方法和穩(wěn)健的空間球殼擬合方法以三個坐標方向上在斷面擬合過程中加入，之后得到了隧道斷面的數(shù)學(xué)方程式。噪聲點會在擬合的過程中被處理掉，避免噪聲點影響到擬合的效果，8毫米成為斷面的半徑擬合的最大誤差。以VC++開發(fā)平臺為研究依據(jù)，從而研發(fā)了三維激光掃描數(shù)據(jù)處理原型系統(tǒng)，該系統(tǒng)集點云數(shù)據(jù)預(yù)處理、后處理和模型演示三項合一，對系統(tǒng)的實用性則適用實測數(shù)據(jù)來證明。

參考文獻

[1] 師海.三維激光掃描技術(shù)在施工隧道監(jiān)測中的應(yīng)用研究[D].北京交通大學(xué)，2017.

[2] 王珂.三維激光掃描技術(shù)在地下洞庫中的應(yīng)用與研究[D].中國地質(zhì)大學(xué)（北京），2019.

[3] 夏金周.移動式三維激光掃描技術(shù)在地鐵隧道變形監(jiān)測中的應(yīng)用研究[D].東南大學(xué)，2019.

[4] 王紅霞.三維激光掃描技術(shù)在橋梁監(jiān)測中的應(yīng)用[D].蘭州理工大學(xué)，2018.

[5] 李宗平，張永濤，楊釗，陳培帥，江鴻.三維激光掃描技術(shù)在隧道變形與斷面檢測中的應(yīng)用研究[J].隧道建設(shè)，2017，37（03）：336-341.