劉昌軍，丁留謙，張順福，廖井霞

中國水利水電科學(xué)研究院，北京 100038

0 引言

在隧道（洞）施工過程中，巖體結(jié)構(gòu)面的存在對隧洞安全有著重要影響?？焖?、準確地獲取巖體結(jié)構(gòu)面的數(shù)量、幾何信息和優(yōu)勢產(chǎn)狀等參數(shù)是隧洞工程巖體穩(wěn)定性評價和地質(zhì)災(zāi)害防治等工作的前提和基礎(chǔ)[1-2]。目前隧洞巖體結(jié)構(gòu)面信息獲取主要以人工現(xiàn)場接觸測量為主，勞動強度大、效率低下，且難以滿足快速施工要求。隨著激光測量和攝影測量的快速發(fā)展，非接觸測量已經(jīng)應(yīng)用在各行各業(yè)中。近幾年來在巖土行業(yè)中得到快速發(fā)展，如高陡邊坡地質(zhì)勘察和測量、危巖體調(diào)查、巖體結(jié)構(gòu)面統(tǒng)計、邊坡變形監(jiān)測等[3-6]。

近年來國內(nèi)外學(xué)者也開始嘗試應(yīng)用攝影測量技術(shù)[3]和激光測量技術(shù)等非接觸測量方法進行高陡邊坡的巖體結(jié)構(gòu)面的調(diào)查，取得了一定成果[4-10]。然而利用非接觸測量方法研究隧洞巖體結(jié)構(gòu)面的研究還較少，一般情況下，隧洞內(nèi)可見光不足，采用攝影測量方法進行調(diào)查往往受到一定限制，而激光測量方法不受可見光影響，在隧洞調(diào)查中具有較強的應(yīng)用優(yōu)勢。

針對隧洞巖體結(jié)構(gòu)面人工現(xiàn)場量測和統(tǒng)計分析的難點，結(jié)合廣泛使用的激光測量技術(shù)和激光點云數(shù)據(jù)處理的難點，提出了一整套應(yīng)用激光測量技術(shù)進行隧洞巖體結(jié)構(gòu)面的全自動統(tǒng)計分析方法，特別是針對隧洞表面激光點云數(shù)據(jù)的特點，應(yīng)用柱面投影的Delaunay生長法對激光點云數(shù)據(jù)進行了三角網(wǎng)重構(gòu)?；诩す鉁y量技術(shù)和FKM（Fuzzy KMean）群聚方法進行隧洞巖體結(jié)構(gòu)面優(yōu)勢產(chǎn)狀統(tǒng)計分析的流程與方法如下：1）首先利用激光掃描儀獲取隧洞表面巖體的激光點云數(shù)據(jù)，結(jié)合全站儀測定的反射片的坐標數(shù)據(jù)將各站激光點云數(shù)據(jù)進行坐標轉(zhuǎn)換和拼接，使得拼接后的點云數(shù)據(jù)的y坐標與北方向重合。2）對處理后的激光點云數(shù)據(jù)進行柱面坐標投影，采用Delaunay生長算法對投影后的激光點云數(shù)據(jù)構(gòu)建三角網(wǎng)，將三維點坐標數(shù)據(jù)按照平面建立的點云拓撲關(guān)系進行點號連接，得到隧洞巖體表面的空間TIN網(wǎng)格模型。3）利用TIN網(wǎng)格模型，求得的每個三角網(wǎng)平面方程的法向向量，以三角網(wǎng)平面方程的法向向量計算其產(chǎn)狀。4）采用FKM模糊群聚方法和改進的方位矩陣方法對每個小三角網(wǎng)的產(chǎn)狀進行全自動群聚分類，用三角網(wǎng)的產(chǎn)狀聚類結(jié)果代表不同組的優(yōu)勢產(chǎn)狀，并用不同顏色表示結(jié)構(gòu)面的群聚分類結(jié)果。

基于以上原理，采用自主研發(fā)的海量激光點云數(shù)據(jù)的處理和巖體結(jié)構(gòu)面統(tǒng)計分析軟件FSS[4]對杭州某隧洞的巖體結(jié)構(gòu)面進行了統(tǒng)計分析，得到了各組結(jié)構(gòu)面的優(yōu)勢產(chǎn)狀及其空間分布。工程實例的研究結(jié)果表明，該方法在隧洞內(nèi)巖體結(jié)構(gòu)面的調(diào)查中具有廣闊的應(yīng)用前景，為隧洞巖體結(jié)構(gòu)面優(yōu)勢產(chǎn)狀的快速調(diào)查和自動統(tǒng)計分析開辟了新的途徑。

1 基于柱面投影的點云重構(gòu)方法的實現(xiàn)

三維激光點云構(gòu)網(wǎng)方法主要采用基于面剖分、體剖分和面投影的3種方法建立[11-14]。基于面投影是先將三維空間數(shù)據(jù)投影到平面上進行平面構(gòu)網(wǎng)，再映射到三維空間最終達到三維構(gòu)網(wǎng)目的。而基于面投影法可分為基于平面投影、基于臨近區(qū)域切平面投影、圓柱面投影、球面投影和橢球面投影等方法。

結(jié)合隧洞激光點云數(shù)據(jù)的特點，采用基于柱面投影的Delaunay構(gòu)網(wǎng)方法對隧洞表面激光點云進行了構(gòu)網(wǎng)?；舅悸啡缦拢簩⑥D(zhuǎn)換成統(tǒng)一空間坐標后的激光點云數(shù)據(jù)以柱面為媒介投影到平面上，得到平面點坐標，在平面內(nèi)構(gòu)造平面三角網(wǎng)，最后將平面三角網(wǎng)映射到空間完成三維構(gòu)網(wǎng)。

1.1 柱面坐標系的建立及投影

1.1.1 柱面坐標系

利用柱面投影坐標系將三維激光點云投影在二維平面內(nèi)構(gòu)三角網(wǎng)是一種快捷有效的構(gòu)網(wǎng)方法。鄭德華等[15]和徐秀川等[11]采用柱面投影法研究了三維激光點云的三角網(wǎng)構(gòu)建方法，其核心是在保證封閉物體表面點之間的連接關(guān)系的前提下建立柱面坐標系統(tǒng)。

如圖1a所示建立柱面坐標系，三維直角坐標系為左手坐標系，圓柱面底面圓心O為坐標中心，圓柱底面半徑為R，圓柱面切割母線AA′。在上述坐標系中，設(shè)P（X，Y，Z）為任一點激光點的坐標數(shù)據(jù)，在圓柱面上投影坐標為P′（r，ω，h）。圓柱面長度為h。

圖1 柱面坐標系和投影后的平面坐標系Fig.1 Cylindrical coordinate system and plane coordinate system after the projection

1.1.2 柱面坐標系統(tǒng)的轉(zhuǎn)換

三維激光掃描獲取的點云坐標系一般是以測站為原點的左手坐標系，為了建立三角形格網(wǎng)首先應(yīng)將其轉(zhuǎn)換到柱面坐標系統(tǒng)中。具體方法如下：三維激光掃描點云坐標中心點與柱面坐標系統(tǒng)的底面圓心點重合，將三維點坐標沿柱面半徑向圓柱基準面投影[16]。三維空間直角坐標系和柱面坐標系的轉(zhuǎn)換公式為

1.1.3 投影后的平面坐標轉(zhuǎn)換為空間坐標

將柱面坐標系轉(zhuǎn)換為平面二維坐標系，具體方法為文獻[16]提出的沿圓柱基準面并沿切割母線AA′切開并展開的方法。展開后的平面坐標原點為圖1b中的A′點，其中左側(cè)的AA′軸為x軸，向上為正向，A′A′為y軸，向右為正向，見圖1b。其轉(zhuǎn)換坐標系矩陣為

式中，θ為切割母線在圓柱基準面上順時針到XZ平面的角度。

1.2 基于柱面投影的構(gòu)網(wǎng)方法的實現(xiàn)

激光點云數(shù)據(jù)柱面投影構(gòu)網(wǎng)的基本思路是，首先確定圓柱面底面中心及半徑以確定進行投影的圓柱面，然后將點云數(shù)據(jù)投影到圓柱面上，將直角坐標變換為柱面坐標。再將柱面坐標展開到平面坐標，在平面坐標系中進行構(gòu)網(wǎng)。在IDL具體實現(xiàn)過程中，需構(gòu)網(wǎng)的點云數(shù)據(jù)的數(shù)組為pointData（一個3×N的2維數(shù)組）。完成柱面投影構(gòu)網(wǎng)算法包含以下幾個步驟：

1）選擇柱面底圓中心orign：根據(jù)圖建立的柱面坐標系，確定激光點云數(shù)據(jù)坐標中心平移到柱面坐標系的底圓圓心坐標。

2）確定圓柱底圓半徑radius：柱面底圓半徑的賦值原則是讓所確定的球與點云中每個點都盡可能的接近，這樣能保證點云投影到柱面上盡可能地減少變形。這需要制定擬合優(yōu)化策略，但在本算法中這并不是重點，因而在此作簡單處理。在pointData中隨機抽取sampleNum個點，將其投影到柱面底面上，將各點坐標平均即可得到柱面底圓半徑，經(jīng)實驗這種方法效果也較好。

3）點云投影到柱面上：柱面確定后將點云數(shù)據(jù)以底圓圓心為原點按照公式（1）進行坐標變換，利用IDL中的函數(shù)CV＿Coord將激光點云從直角坐標變換為柱面坐標pointData1。

4）柱面坐標系轉(zhuǎn)換為平面坐標：將柱面坐標按照圖1b方法展開為平面二維坐標，在平面坐標內(nèi)利用IDL函數(shù)QHull構(gòu)建Delaunay三角網(wǎng)，得到點云數(shù)據(jù)的拓撲關(guān)系，并存放至triangle數(shù)組中。

5）平面坐標轉(zhuǎn)換為空間坐標：利用公式（2）將平面坐標轉(zhuǎn)換為空間坐標，將三維空間點按照平面三角構(gòu)網(wǎng)確定的拓撲關(guān)系根據(jù)點號進行連接，形成空間Delaunay三角網(wǎng)。

2 FKM方法模糊聚類的基本理論

2.1 FKM模糊聚類方法

本研究將采用FKM計算方法進行結(jié)構(gòu)面產(chǎn)狀分類的分析[17]。FKM算法通過引入隸屬度函數(shù)表示屬于不同類別的程度對數(shù)據(jù)進行劃分：首先估計聚類的中心，其次反復(fù)調(diào)整聚類中心，使數(shù)據(jù)集中的每個點距每個中心的距離之和達到最小或滿足終止條件。

給定n個數(shù)據(jù)樣本集合X＝（x1，x2，…，xn），觀測樣本xi的特性指標為xi＝（xi1，xi2，…，xin），也即每個樣本具有m個屬性，n個數(shù)據(jù)對象分成k（2≤k≤n）個不同類。數(shù)據(jù)集合X中的對象xi以一定的隸屬度隸屬某一類，每一類對應(yīng)一個模糊矩陣U。如果X中含有k個類別，則有這k個聚類中心構(gòu)成的向量矩陣為V，由此可建立FKM聚類算法的目標函數(shù)：

式中：u為法向向量對各群心法向量的隸屬度矩陣；Xj≡（xj，yj，zj）為平面的單位法向向量，xj，yj，zj為Xj的分量，計算見式（4）；Vi≡（xi，yi，zi）為各群聚的群心法向向量；m為模糊度，由于m控制隸屬度在各類之間共享的程度，所以m越大，模糊性就越大，d（Xj，Vi）為Xj與Vi的距離。

式中：β為巖體結(jié)構(gòu)面的傾角；α為結(jié)構(gòu)面的傾向。

R.E.Hammah和J.H.Curran[8]提出將結(jié)構(gòu)面產(chǎn)狀轉(zhuǎn)換成在球體（半徑為1）上的單位法向向量，然后計算單位法向向量與群聚的群心法向向量的正弦值，作為此兩向量間的距離值：

FKM聚類算法是在滿足隸屬度uij值為

的條件下，求解目標函數(shù)J的最小值，即

且元素uij滿足以下條件：

利用 R.E.Hammah和J.H.Curran[8]提出改進的方位矩陣S＊計算新的群心法向量[2，10]：

利用改進的方位矩陣可以改進新群心法向向量在橢圓形中不均勻分布情況。

取m＝2，利用反復(fù)迭代運算，尋找最終群心法向向量及相應(yīng)隸屬度矩陣。迭代收斂以兩次隸屬度矩陣的各元素差異值的絕對值小于允許值ε為準，即

式中：u為第一次隸屬度矩陣；u′為第二次隸屬度矩陣。

2.2 激光掃描和FKM聚類耦合方法的技術(shù)流程

基于激光掃描和FKM聚類耦合方法的巖體結(jié)構(gòu)面優(yōu)勢產(chǎn)狀統(tǒng)計的技術(shù)流程如圖2所示。

2.3 激光掃描和FKM聚類耦合方法的適用條件

基于激光掃描技術(shù)和FKM聚類耦合方法研究巖體結(jié)構(gòu)面優(yōu)勢產(chǎn)狀比傳統(tǒng)方法具有一定的優(yōu)勢，主要體現(xiàn)在以下幾個方面：

1）采用激光掃描技術(shù)可以進行遠距離巖體結(jié)構(gòu)面快速掃描和統(tǒng)計，特別適合現(xiàn)場無法到達和存在危險區(qū)域的巖體結(jié)構(gòu)面統(tǒng)計分析。

2）采用改進的方位矩陣和FKM模糊聚類方法可以模擬巖體結(jié)構(gòu)面的發(fā)育分布規(guī)律，彌補了等密度圖、玫瑰圖等傳統(tǒng)統(tǒng)計方法的不足，量化的給出了結(jié)構(gòu)面自動分組及各組的優(yōu)勢產(chǎn)狀，且可以分組進行三維空間顯示。

3）由于三維激光掃描是自發(fā)光源，在一些無可見光或可見光不足情況下（如地下洞室、隧洞和巷道等）仍然可以使用，可以實時獲取和分析施工過程巖體裂隙發(fā)育情況，用于指導(dǎo)施工和設(shè)計。

4）該方法是利用巖體表面三角網(wǎng)的群集特性進行的巖體結(jié)構(gòu)面統(tǒng)計，巖體表面所有的結(jié)構(gòu)面都參與了群集分類，三角網(wǎng)平面產(chǎn)狀的群聚結(jié)果可以作為巖體結(jié)構(gòu)面各組的優(yōu)勢產(chǎn)狀。而傳統(tǒng)的結(jié)構(gòu)面優(yōu)勢產(chǎn)狀統(tǒng)計都是只有局部或部分結(jié)構(gòu)面產(chǎn)狀參與了統(tǒng)計分析。

圖2 基于激光點云的巖體結(jié)構(gòu)面FKM聚類算法流程圖Fig.2 Flow diagram of FKM clustering algorithm based on laser scanner

雖然該方法在巖體結(jié)構(gòu)面自動統(tǒng)計中具有較大優(yōu)勢，但也存在一定的局限性：

1）該算法特別適合于裸露巖體，在一些植被茂密的巖體上，該方法有一定的局限性。

2）由于該方法是用三角網(wǎng)的產(chǎn)狀替代裂隙產(chǎn)狀，因此在一些巖體表面較為光滑，且不能區(qū)分結(jié)構(gòu)面時，就會大大降低其統(tǒng)計結(jié)果的精度。

3）對于一些特別細小的結(jié)構(gòu)面，在遠距離掃描時，由于掃描精度問題，細小的結(jié)構(gòu)面不能被識別出來，將影響統(tǒng)計結(jié)果的精度。

4）如果巖體表面特別潮濕或有地下水滲出的情況，該方法的使用會受到一定的限制。

2.4 激光測量和FKM聚類耦合結(jié)果的分析方法

使用基于激光測量和FKM聚類的耦合分析方法研究巖土結(jié)構(gòu)面的優(yōu)勢產(chǎn)狀與傳統(tǒng)的統(tǒng)計方法不同。該方法得到的結(jié)果是巖體表面三角網(wǎng)的優(yōu)勢產(chǎn)狀的聚類結(jié)果，用三角網(wǎng)產(chǎn)狀的聚類結(jié)果來代替巖體結(jié)構(gòu)面分組優(yōu)勢產(chǎn)狀，而且隨著不同群聚數(shù)的增加，三角網(wǎng)產(chǎn)狀的聚類組數(shù)也相應(yīng)增加。將該方法應(yīng)用到實際工程時，需要利用人工現(xiàn)場判別結(jié)構(gòu)面的分組將不同群聚數(shù)的分組結(jié)果加以區(qū)分，對人工開鑿和非結(jié)構(gòu)面的聚類結(jié)果進行剔除，進而確定合適的群聚數(shù)和分組結(jié)果。

圖3 FSS軟件系統(tǒng)框架和工作流程Fig.3 The system framework and workflow of FSS software

3 軟件編制

采用IDL語言編制了基于激光點云數(shù)據(jù)的巖體結(jié)構(gòu)面處理軟件FSS，該軟件主要包括數(shù)據(jù)處理模塊、結(jié)構(gòu)面統(tǒng)計分析模塊和制圖輸出模塊3個主要應(yīng)用模塊。其中：數(shù)據(jù)處理模塊主要包括激光點云的處理功能（點云的局部刪除、抽稀和坐標轉(zhuǎn)換等）、點云的構(gòu)網(wǎng)（基于平面、球面和柱面等投影算法的生長法構(gòu)網(wǎng)）等功能；結(jié)構(gòu)面統(tǒng)計分析模塊包括單個結(jié)構(gòu)面的空間擬合、繪制跡線和巖體結(jié)構(gòu)面的模糊群聚分析等功能[4]。FSS軟件的系統(tǒng)框架和工作流程見圖3，軟件界面見圖4。

4 應(yīng)用案例

4.1 研究區(qū)基本情況

以杭州將臺山某人工開挖隧洞的巖體結(jié)構(gòu)面統(tǒng)計分析為例。杭州將臺山隧洞位于將臺山東南麓，毗鄰八卦田，遙望錢塘江。隧洞位于將臺山東南麓邊坡的西側(cè)，為丁字形交叉洞室，洞口呈半圓形。該洞總長約200m，裂隙較為發(fā)育，多閉合，無充填，延伸長度較短。

4.2 現(xiàn)場掃描及數(shù)據(jù)處理

首先利用地面激光掃描儀對巖洞內(nèi)的表面進行測量，獲取巖洞內(nèi)巖體的激光點云坐標約70萬個，同時利用全站儀獲取現(xiàn)場布置的15個圓形反射片中心點大地坐標。利用15個反射片中心點大地坐標將掃描獲取的激光點云數(shù)據(jù)進行坐標轉(zhuǎn)換，將轉(zhuǎn)換成大地坐標的激光點云數(shù)據(jù)導(dǎo)入FSS軟件，有關(guān)FSS軟件功能及界面見文獻[4]。本次調(diào)查選擇了東西向的裂隙分布較為密集的一段山體隧洞，從掃描數(shù)據(jù)可查詢該巖洞長約15m，高約5m，直徑約6.5m。圖5為洞口附近現(xiàn)場照片和洞口附近激光點云數(shù)據(jù)（點云賦值照片上的RGB色彩）。圖6為隧洞研究段的點云圖。

根據(jù)獲取的激光點云數(shù)據(jù)，為了構(gòu)網(wǎng)方便，將掃描獲取的激光點云數(shù)據(jù)抽稀一半，對抽稀后的激光點云采用第2節(jié)介紹的方法進行空間三角網(wǎng)重構(gòu)，構(gòu)網(wǎng)后的隧洞三維網(wǎng)格圖見圖7。

圖4 FSS軟件界面圖Fig.4 The interface diagram of FSS software

圖5 洞口局部的現(xiàn)場照片和激光點云數(shù)據(jù)Fig.5 Local scene photos and laser point cloud data of the entrance of cave

圖6 隧洞激光點云圖Fig.6 Diagram of laser point cloud of the tunnel

圖7 隧洞三維TIN網(wǎng)圖Fig.7 3-D mesh figure of the tunnel

4.3 結(jié)構(gòu)面模糊統(tǒng)計分析

從現(xiàn)場觀察可知，該隧洞主要有5組巖體結(jié)構(gòu)面較為發(fā)育。利用隧洞表面的空間TIN網(wǎng)，采用FKM的模糊群聚方法對隧洞存在的結(jié)構(gòu)面優(yōu)勢產(chǎn)狀進行統(tǒng)計分析。圖8給出了不同群聚數(shù)下（k＝4～9）的隧洞巖體結(jié)構(gòu)面的空間分布和優(yōu)勢產(chǎn)狀，表1給出了相應(yīng)群聚數(shù)下（k＝4～9）的結(jié)構(gòu)面優(yōu)勢產(chǎn)狀統(tǒng)計值和現(xiàn)場實測值。

由表1可知：第6組結(jié)構(gòu)面在k＝6以后和第3組結(jié)構(gòu)面傾向和傾角較為接近，第7組結(jié)構(gòu)面在k＝7以后和第1組結(jié)構(gòu)面的傾向和傾角較為接近；2組結(jié)構(gòu)面在k＝6以后其優(yōu)勢產(chǎn)狀的變化量變小，且在k＝6以后才能區(qū)分2組結(jié)構(gòu)面。而第9組結(jié)構(gòu)面在k＝9以后才出現(xiàn)，且均不和前8組結(jié)構(gòu)面相近。因此，對于該隧洞的巖體結(jié)構(gòu)面取k＝6的分類結(jié)果較為合適。從分類的結(jié)果圖（圖8）上可看出，第8組結(jié)構(gòu)面為隧洞的地面和頂面，且不同群聚數(shù)下的該組結(jié)構(gòu)面的分類結(jié)果呈現(xiàn)離散性，不作為單獨一組。因此從k＝6～9的分類結(jié)果看，隧洞內(nèi)的巖體結(jié)構(gòu)面的優(yōu)勢產(chǎn)狀共5組，且和實際測量值接近。另外，由以上模糊群聚的自動分析結(jié)果可知，計算的群心法向向量（即優(yōu)勢產(chǎn)狀）隨群聚數(shù)的增加，均出現(xiàn)變小的結(jié)果，即優(yōu)勢產(chǎn)狀隨著群聚數(shù)增加，其值變化不大，且能將產(chǎn)狀接近的不同結(jié)構(gòu)面加以辨識。

由上面的應(yīng)用案例的統(tǒng)計結(jié)果可看出，本文提出的基于激光點云數(shù)據(jù)和FKM的巖體結(jié)構(gòu)面群聚分析方法可用于復(fù)雜洞室的結(jié)構(gòu)面優(yōu)勢產(chǎn)狀統(tǒng)計和分類，計算結(jié)果可靠，計算速度快，能夠適應(yīng)于隧洞內(nèi)特殊條件的巖體穩(wěn)定性分析。另外，由激光點云數(shù)據(jù)可以對隧洞裂隙的幾何信息（如跡長、間距和縫寬等）進行分析，也可給出其地質(zhì)編錄圖，詳細研究成果見文獻[4]。

圖8 不同群聚數(shù)下的巖體結(jié)構(gòu)面模糊群聚結(jié)果Fig.8 Clustering results of rock mass discontinuity occurrence under different clustering numbers

表1 不同群聚數(shù)下的巖體結(jié)構(gòu)面優(yōu)勢產(chǎn)狀統(tǒng)計Table1 Statistics of different number of the attitude of rock mass discontinuity in advantage

5 結(jié)論

1）將激光掃描技術(shù)和FKM的模糊群聚方法應(yīng)用于隧洞等復(fù)雜條件下的巖體結(jié)構(gòu)面統(tǒng)計分析中，該方法可以實現(xiàn)非接觸測量的巖體結(jié)構(gòu)面優(yōu)勢產(chǎn)狀的自動聚類分析。群聚結(jié)果表明，部分群聚的幾何中心（群心值）隨著群聚數(shù)目增加，其變化量將會減小，該群心值類別的顏色可反映其結(jié)構(gòu)面的空間分布，且群心值可作為該結(jié)構(gòu)面優(yōu)勢產(chǎn)狀的代表值。工程實例應(yīng)用表明，該方法計算結(jié)果可靠、合理，具有較強的通用性和推廣性。

2）采用的基于柱面投影的Delaunay生長法進行復(fù)雜激光點云的構(gòu)網(wǎng)方法速度快、構(gòu)網(wǎng)質(zhì)量高、重構(gòu)得到的三角網(wǎng)表面與被采樣的巖體表面拓撲差別最小。同時還提供了基于球面投影和平面投影的激光點云的生長法的構(gòu)網(wǎng)算法，解決了復(fù)雜條件下的激光點云的快速構(gòu)網(wǎng)問題。

3）利用自主研發(fā)的基于海量激光點云數(shù)據(jù)的巖體結(jié)構(gòu)面統(tǒng)計分析軟件FSS，可以快速對掃描獲取的巖體表面的激光點云數(shù)據(jù)進行巖體結(jié)構(gòu)面統(tǒng)計分析和編錄工作，減少了高陡邊坡人為攀爬的危險，減少了現(xiàn)場工作和數(shù)據(jù)后處理工作。

（References）：

[1]王述紅，張靖杰，李云龍，等.工程巖體開挖過程全空間塊體搜索及系統(tǒng)研制[J].東北大學(xué)學(xué)報：自然科學(xué)版，2010，31（7）：1026-1029.Wang Shuhong，Zhang Jingjie，Li Yunlong，et al.Quickly Spatial Identification Process of Blocks During Rock Excavation and Its System Modeling[J].Journal of Northeastern University：Natural Science，2010，31（7）：1026-1029.

[2]王述紅，楊勇，王洋，等.基于數(shù)字攝影測量的開挖空間模型及不穩(wěn)定塊體的快速識別[J].巖石力學(xué)與工程學(xué)報，2010，29（增刊1）：3432-3437.Wang Shuhong，Yang Yong，Wang Yang，et al.Spatial Modeling and Quickly Identification of Unsteady Rock Blocks Based on Digital Photogrammetry [J].Chinese Journal of Rock Mechanics and Engineering，2010，29（Sup.1）：3432-3437.

[3]王鳳艷，陳劍平，楊國東，等.基于數(shù)字近景攝影測量的巖體結(jié)構(gòu)面幾何信息解算模型[J].吉林大學(xué)學(xué)報：地球科學(xué)版，2012，42（6）：1839-1846.Wang Fengyan，Chen Jianping，Yang Guodong，et al.Solution Models of Geometrical Information of Rock Mass Discontinuity Based on Digital Close Range Photogrammetry[J].Journal of Jilin University：Earth Science Edition，2012，42（6）：1839-1846.

[4]劉昌軍，丁留謙，孫東亞.基于激光點云數(shù)據(jù)的巖體結(jié)構(gòu)面全自動模糊群聚分析及幾何信息獲取[J].巖石力學(xué)與工程學(xué)報，2011，30（2）：358-365.Liu Changjun，Ding Liuqian，Sun Dongya.Automatic Fuzzy Clustering Analysis and Geometric Information Acquisition of Rock Mass Discontinuity Based on Laser Point Cloud Data[J].Chinese Journal of Rock Mechanics，2011，30（2）：358-365.

[5]劉昌軍，高立東，丁留謙，等.應(yīng)用激光掃描技術(shù)進行巖體結(jié)構(gòu)面的半自動統(tǒng)計研究[J].水文地質(zhì)工程地質(zhì)，2011（2）：52-57.Liu Changjun，Gao Lidong，Ding Liuqian，et al.Research of Semi-Automatic Statistics of Rock Mass Discontinuity Applying Laser Scanning Technology[J].Hydrogeology and Engineering Geology，2011（2）：52-57.

[6]劉昌軍，丁留謙，孫東亞.三星堆月亮灣城墻遺址覆土方量計算[J].水利水電科技進展，2011（2）：81-84.Liu Changjun，Ding Liuqian，Sun Dongya.Calculation of Soil Quantity in Moon Bay Wall Site of Sanxingdui[J].Advances in Science and Technology of Water Resources，2011（2）：81-84.

[7]Feng Q，Sj?gren P，Stephansson O，et al.Measuring Fracture Orientation at Exposed Rock Faces by Using a Non-Reflector Total Station [J].Engineering Geology，2001，59（1／2）：133-146.

[8]Hammah R E，Curran J H.Fuzzy Cluster Algorithm for the Automatic Identification of Joints Sets[J].International Journal of Rock Mechanics and Mining Sciences，1998，35（7）：889-905.

[9]王鳳艷，陳劍平，付學(xué)惠，等.基于VirtuoZo的巖體結(jié)構(gòu)面幾何信息獲取研究[J].巖石力學(xué)與工程學(xué)報，2008，27（1）：169-175.Wang Fengyan，Chen Jianping，F(xiàn)u Xuehui，et al.Study on Geometrical Information of Obtaining Rock Mass Discontinuites Based on VirtuoZo[J].Chinese Journal of Rock Mechanics and Engineering，2008，27（1）：169-175.

[10]洪子恩，馮正一，吳宗江.應(yīng)用三維激光掃描于巖坡露頭位能的量測[J].水土保持學(xué)報，2007，39（3）：247-267.Hong Zien， Feng Zhengyi， Wu Zongjiang.Application of 3DLaser Scanning on Measuring Attitudes of Rock Slop Outcrops[J].Journal of Soil and Water Conversion，2007，39（3）：247-267.

[11]徐秀川，花向紅，龔子楨，等.基于球面投影的散亂點云三維建模算法實現(xiàn)與效果分析[J].測繪工程，2011，20（3）：5-8.Xu Xiuchuan，Hua Xianghong，Gong Zizhen，et al.The Implementation Algorithm of Three-Dimensional Modeling to Point Cloud Based on Spherical Surface Projection and Effect Analysis[J].Engineering of Survey and Mapping，2011，20（3）：5-8.

[12]武曉波，王世新，肖春生.Delaunay三角網(wǎng)的生成算法研究[J].測繪學(xué)報，1999，28（1）：28-35.Wu Xiaobo，Wang Shixin，Xiao Chunsheng.A New Study of Delaunay Triangulation Creation[J].Acta Geodaetica et Cartographic Sinica，1999，28（1）：28-35.

[13]張帆，黃先鋒，李德仁.基于球面投影的單站地面激光掃描點云構(gòu)網(wǎng)方法[J].測繪學(xué)報，2009，38（1）：48-54.Zhang Fan，Huang Xianfeng，Li Deren.Spherical Projection Based Triangulation for One Station Terrestrial Laser Scanning Point Cloud[J].Acta Geodaetica et Cartographic Sinica，2009，38（1）：48-54.

[14]李逢春，龔俊，王青.基于三維TIN的精細表面建模方法[J].計算機應(yīng)用研究，2006，23（8）：159-161.Li Fengchun，Gong Jun，Wang Qing.Highly Refined Surface Reconstruction Method Based on 3DTIN[J].Application Research of Computers，2006，23（8）：159-161.

[15]鄭德華，張云濤.基于物體表面散亂三維激光掃描點的三角形格網(wǎng)建立[J].測繪工程，2004，13（4）：62-65.Zheng Dehua，Zhang Yuntao.Triangular Establishment Based on Disordered 3DLaser Scanning Point of Object Surface[J].Engineering of Surveying and Mapping，2004，13（4）：62-65.

[16]潘華偉，鄒北驥.一種圓柱形全景圖生成新算法及其實現(xiàn)[J].計算工程與科學(xué)，2003，25（6）：13-16.Pan Huawei，Zou Beiji.A New Algorithm for Generating Cylinder Panorama and Its Implementation[J].Computer Engineering and Science，2003，25（6）：13-16.

[17]田赤英，張旭男，宋士吉，等.改進的模糊C均值聚類算法及其在海底熱液硫化物組分分析中的應(yīng)用[J].海洋學(xué)研究，2010，28（4）：22-28.Tian Chiying，Zhang Xunan，Song Shiji，et al.Algorithm of Improved FCM Clustering and Its Application to Seafloor Hydrothermal Sulfide Composition Analysis[J].Journal of Marine Science，2010，28（4）：22-28.