張 靖,張愛能,劉國棟
(1.長安大學 地質工程與測繪學院,陜西 西安710064;2.重慶交通大學 土木建筑學院,重慶400074)
三維激光掃描技術是上世紀九十年代中期開始出現的一項高新技術,是繼GPS 空間定位系統之后又一項測繪技術新突破。其突破了傳統單點測量方法,通過高速激光掃描測量的方法,大面積高分辨率地快速獲取被測對象表面的三維坐標數據。三維激光掃描儀技術廣泛應用于變形監測、工程測量、地形測量、古建筑和文物保護、斷面和體積測量等領域,具有不需要合作目標、高精度、高密度、高效率、全數字特征等優點[1-5]。
常用的地形圖測量方法有全站儀數字化測圖、GPS-RTK 測圖、數字攝影測量等方法。但這些方法外業工作量大,且在地形險峻,人難以到達的地方往往顯得無能為力。從三維激光掃描儀的測量原理和掃描過程可知,三維激光掃描以格網掃描方式,高精度、高密度、高速度和免棱鏡地測量地表點,能詳細了解地面點的細節變形和整體變化,同時三維激光掃描技術制作的地形圖精度優于傳統方法,且可大大縮短外業工作時間,將大部分時間轉為在軟件中對掃描數據的內業處理[6-8]。基于三維激光掃描的地形測繪成圖技術的應用,改變了傳統測繪的作業流程,使相關外業測繪流程大大簡化,外業工作時間大大縮短,外業人員的勞動強度大大降低,內業處理的自動化程度也顯著提高。因此將三維激光掃描技術引入到地形測量中是具有重要意義的。
文中主要介紹了三維激光掃描儀在地形測量中的應用。首先對實驗采用的儀器、測區概況、控制網布設及掃描儀掃描方法做了簡要的介紹。然后介紹了將基于掃描儀坐標系統的點云數據轉換到基于施工獨立坐標系統的過程,探討了將平面坐標和高程分開轉換的方法,將轉換后的點云數據導入到Geomagic 中進行重采用,同時進行數據格式轉換,方便cass 繪制等高線。通過全站儀測量一定數量檢核點來驗證三維激光掃描儀精度。
本次實驗使用徠卡ScanStation2 三維激光掃描儀掃描地形,使用尼康DTM -332C 全站儀布設控制網,其測角精度為3″,測距精度為±3 mm +2 ppm.測區位于西安市灞橋區白鹿塬北坡,掃描區域的地形起伏較小,結構較復雜,人流量相對較小,地形表面的雜草和植被較少。
在測區內布設獨立工程控制網,K4點為已知點,平面坐標值設為1 000,1 000 m,高程值為500 m.在測區內布設控制點如圖1,通過全站儀測角、測邊,求得控制點的坐標,使用全站儀三角高程測量獲取控制點高程。

圖1 控制點空間分布圖Fig.1 Control point spatial distribution
將掃描儀安置在控制點K1上,在測區合適位置安置3 個不在同一直線上的藍白標靶,將其中一個藍白標靶安置在控制點k5上。連接掃描儀和電源,打開電源開關讓儀器進行自檢;設置筆記本與掃描儀之間的通訊,啟動配套的Cyclone 軟件并建立數據庫與工程文件,連接掃描儀并通過設置角度范圍對掃描區域進行拍照。設置好參數后進行掃描,找到控制點上標靶大致位置,對標靶進行精細掃描。然后將掃描儀架設到k3上,采用相同方法依次掃描。依次進行,直到所有掃描完成。
點云拼接是多個測站點云數據的整合,是基于掃描儀不同坐標系統轉化統一的過程。點云拼接方法有:基于連接點的拼接;基于標靶的拼接;基于控制點的拼接。本實驗采用基于公共標靶的拼接,拼接后點云數據如圖2 所示。

圖2 測區點云拼接后數據Fig.2 Splicing point cloud data
點云數據預處理過程包括:點云去噪和平滑,點云抽稀壓縮,點云空洞修補等。這些過程都可以通過Cyclone 和Geomagic 軟件完成。

圖3 基于掃描儀平面坐標與施工獨立平面坐標關系Fig.3 Relationship between scanner coordinate and construction coordinate
由于獲取的點云數據平面坐標和高程是基于掃描儀坐標系統的,而地形圖采用的則是獨立工程測量坐標系,因此在用點云數據繪制地形圖前需要進行平面坐標轉換和高程轉換。將基于三維掃描儀點云數據坐標系轉換到獨立工程測量坐標系中。
3.1.1 平面坐標轉換
建立基于掃描儀平面坐標系O-X1Y1,基于工程測量坐標系O-X2Y2,兩個坐標系關系如圖3 所示。采用四參數模型進行平面坐標轉換,包括平移參數:X0,Y0,旋轉參數:a,尺度變化:μ[9].
對于任意一點Pi在兩個坐標系中坐標分別為(X1,Y1),(X2,Y2),存在如下關系

為求出式1 中平移、旋轉和尺度變化參數,至少需要2 個已知平面點,,如多于兩個可采用最小二乘法擬合求解。文中已知k3,k4,k5等3 點在兩坐標系中坐標,見表1.

表1 已知點在不同坐標系平面坐標值Tab.1 Known point in different coordinates plane coordinates
采用最小二乘法,通過已知點平面坐標值求解轉換參數,然后將點云基于掃描儀平面坐標值全部轉到基于工程測量坐標系中,完成坐標轉換。

圖4 高程關系Fig.4 Elevation relationship
3.1.2 高程轉換
掃描儀獲取點云數據中,其高程值起算面是基于掃描儀系統的,把這種高程值記為H點云;通過全站儀獲取的控制點高程值其起算是假定的水準面,在小范圍內可以把其當成平面,把這種高程記為H水準[10-13]。這2 種高程值關系如圖4 所示。根據圖可知 H水準=H點云+N.

表2 高程值Tab.2 Elevation value
采用平面擬合法求取任意點云高程差值N.假設任意點云高程差值N 與該點平面坐標,有

其中 a,b,c 為模型參數。為求出公式2 中的模型參數,至少需要3 個已知公共高程點。文中已知k3,k4,k5,3 點在兩高程系統高程值,見表2.
將3 個公共點高程值數據代入公式2 中,求解模型參數。最后將點云數據高程值轉為基于假定水準面的高程值。
點云數據生成等高線流程如圖5 所示。

圖5 等高線生成流程圖Fig.5 Contour generated flowchart
3.2.1 點云采樣
由于點云數據量大,本次實驗采集的點云數據約為14 萬個,如果直接導入到cass 中肯定不利于坐標數據生成等高線,需對點云數據進行采樣,在Geomagic 中,點云采樣方法有:曲率采樣、等距采樣、統一采樣和隨機采樣。本次實驗采樣曲率采樣,即設定某一百分比值,軟件按照最小二乘法計算某曲面上的曲率,依次類推計算出所有的曲率,然后依照曲率大小劃分的區間對點云進行精簡刪除,得到設定的百分比與點云數值乘積的點值。采樣后點云數據約為1.4 萬個。
3.2.2 等高線繪制
將采樣后的點云數據導入到cass 中,建立DTM 模型,設置等高距,cass 軟件自動生成等高線如圖6 所示。

圖6 等高線Fig.6 The contour lines
本次實驗在測區內均勻布設了8 個檢核點,掃描時在檢核點上安置藍白標靶。通過全站儀獲取的這8 個檢核點的平面坐標和高程與檢核點點云轉換后坐標數據進行比較來驗證地形圖的可靠性。表3 給出了兩者比較的結果。
從表3 可以看出,三維激光掃描儀獲取的檢核點平面坐標與全站儀獲取點的平面坐標的差值的絕對值最小為0.005 m,最大值為0.024 m,平均值為0.012 m,高程差值絕對值最小值為0.006 m,最大值為0.029 m.平均值為0.018 m.掃描儀轉換后的數據誤差主要包括控制測量誤差,掃描儀系統誤差,坐標轉換誤差。比較結果說明,三維激光掃描儀的精度滿足地形圖測量需要。

表3 全站儀測量檢核點坐標與激光掃描儀獲取點坐標比較Tab.3 Comparison of the check points’coordinates measured by total station and scanner
通過實驗,將三維激光掃描儀用于地形測量,實驗結果表明三維激光掃描儀的精度滿足地形圖測量需要。三維激光掃描是一種具有高精度、高密度、高速度、非接觸的掃描手段,在地形測量中,特別是在地形險峻、人難以到達的地方有著較大的優勢。
References
[1] 劉昌軍,趙雨,葉長鋒,等.基于三維激光掃描技術的礦山地形快速測量的關鍵技術研究[J].測繪通報,2012(6):43 -46.LIU Chang-jun,ZHAO Yu,YE Chang-feng,et al. Study of the key technologies for mine rapid topographic survey based on 3D laser measurement[J].Bulletin of Surveying and Mapping,2012(6):43 -46.
[2] 彭維吉,李孝雁,黃 颯.基于地面三維激光掃描技術的快速地形圖測繪[J]. 測繪通報,2013(3):70 -72.PENG Wei-ji,LI Xiao-yan,HUANG Sa. Rapid topographic mapping using terrain 3D laser scanning technique[J]. Bulletin of Surveying and Mapping,2013(3):70 -72.
[3] 李 亮,吳 侃,劉 虎,等.地面三維激光掃描地形測量數據粗差剔除算法及實現[J].測繪科學,2010,35(3):187 -189.LI Liang,WU Kan,LIU Hu,et al. Gross error elimination of terrain survey data obtained by terrestrial 3D laser scanning system[J].Science of Surveying and Mapping,2010,35(3):187 -189.
[4] 何秉順,趙進勇,王 力,等.三維激光掃描技術在堰塞湖地形快速測量中的應用[J]. 防災減災工程學報,2008(3):394 -398.HE Bing-shun,ZHAO Jin-yong,WANG Li,et al.Application of 3D LIDAR in landslide lake topographic measurement[J].Journal of Disaster Prevention and Mitigation,2008(3):394 -398.
[5] 呂寶雄,巨天力.三維激光掃描技術在水電大比例尺地形測量中的應用研究[J]. 西北水電,2011(1):14-16.LV Bao-xiong,JU Tian-li.Study on application of 3D laser scanning technology in hydropower large scale topographical survey[J]. Northwest Hydropower,2011(1):14 -16.
[6] 陳 靜,李清泉,李必軍.激光掃描測量系統的應用研究[J].測繪工程,2001,10(1):49 -52.CHEN Jing,LI Qing-quan,LI Bi-jun. Application research on laser-scanning surveying system[J]. Engineering of Surveying and Mapping,2001,10(1):49 -52.
[7] 梅文勝,周燕芳,周 俊.基于地面三維激光掃描的精細地形測繪[J].測繪通報,2010(1):53 -56.MEI Wen-sheng,ZHOU Yan-fang,ZHOU Jun.Fine topographic mapping based on ground three-dimensional laser scanning[J]. Bulletin of Surveying and Mapping,2010(1):53 -56.
[8] 曹 璽,梁 俊,王延洪,等.三維激光掃描技術在地形測繪成圖中的應用[J]. 人民長江,2011,42(21):22 -24.CAO Xi,LIANG Jun,WANG Yan-hong,et al. Application of 3D LIDAR in terrain surveying and mapping[J].Yangtze Rive,2011,42(21):22 -24.
[9] 宗剛軍,姚頑強.工程測量中新舊坐標轉換的一種實現方法[J].西安科技大學學報,2007,27(3):401 -404.ZONG Gang-jun,YAO Wan-qiang. Implementing of coordinate transfering between two survey system in engineering survey[J]. Journal of Xi’an University of Science and Technology,2007,27(3):401 -404.
[10]劉萬林,郭 嵐,王 利.多面函數法與移動法的加權綜合模型在GPS 水準中的應用[J].西安科技大學學報,2004,24(3):310 -312.LIU Wan-lin,GUO Lan,WANG Li.Application of combining multisurface function andmovement method in GPS leveling[J].Journal of Xi’an University of Science and Technology,2004,24(3):310 -312.
[11]馬素文.三維激光掃描在測量中的應用現狀[J]. 山西建筑,2011,37(9):207 -208.MA Su-wen. Application status of 3D laser scanner in the measurement[J]. Shanxi Architecture,2011,37(9):207 -208.
[12]張啟福,孫現申.三維激光掃描儀測量方法與前景展望[J].北京測繪,2011(1):39 -42.ZHANG Qi-fu,SUN Xian-shen.Measuring principle and developmental prospect of 3D laser scanner[J].Beijing Surveying and Mapping,2011(1):39 -42.
[13]吳 侃,黃承亮,陳冉麗.三維激光掃描技術在建筑物變形監測的應用[J].遼寧工程技術大學學報:自然科學版,2011,30(2):205 -208.WU Kan,HUANG Cheng-liang,CIHEN Ran-li.Application of 3D laser scanning technology in monitoring building deformation[J].Journal of Liaoning Technical University:Natural Science Edition,2011,30(2):205 -208.