石光 胥洪川 彭暢 龍洪 馮國正
摘要:高原湖區(qū)氣候惡劣、交通不便、風(fēng)力較大,采用傳統(tǒng)測量方式或航測方式較為困難。根據(jù)高原湖區(qū)植被稀疏、地物較少、地勢較陡、湖泊航行條件良好的特點,可開展船載三維激光掃描系統(tǒng)湖區(qū)陸上地形測量工作。介紹了船載三維激光系統(tǒng)坐標轉(zhuǎn)換原理,分析各傳感器及其集成誤差,采用誤差傳播定律進行系統(tǒng)誤差估計。經(jīng)實例數(shù)據(jù)進行精度評定,證實了船載三維激光掃描系統(tǒng)在高原測量可靠性好、精度高,可為高原測湖陸上地形測量提供新的解決方案。
關(guān)鍵詞:船載三維激光掃描系統(tǒng);系統(tǒng)誤差;誤差傳播定律;精度分析;高原湖泊測量
中圖法分類號:P343.3
文獻標志碼:A
DOI: 10.15974/j.cnki.slsdkb.2020.05.001
船載三維激光掃描系統(tǒng)構(gòu)建技術(shù)成熟,并已在內(nèi)陸河道地形測量、海岸線測繪、島礁測繪和港口測量中得到應(yīng)用[1-3]。我國青藏高原地區(qū)湖泊眾多,卻因技術(shù)和安全原因,大多并未開展過系統(tǒng)科學(xué)的測量工作。因此,對這些高原湖泊開展測量工作,對填補國情資料空白、湖泊保護研究等具有重要意義。
高原湖區(qū)氣候惡劣、交通不便,難以采用傳統(tǒng)陸地測量方式。因高原氣壓、風(fēng)力等氣候因素影響,且申請空域、任務(wù)報批難度大,作業(yè)成本高,故也不宜采用航空攝影或機載雷達作業(yè)。青藏高原植被稀疏、地物較少、地形陡峭、湖泊航行條件良好,適宜船載三維激光掃描系統(tǒng)進行陸上地形測繪。在闡述船載三維掃描系統(tǒng)誤差來源、各傳感器精度估計的基礎(chǔ)上,結(jié)合船載三維激光掃描儀系統(tǒng)時空配準原理,利用誤差傳播定律對船載三維激光系統(tǒng)集成精度進行評估與分析,并采用實例數(shù)據(jù)評定系統(tǒng)精度。
1 系統(tǒng)坐標轉(zhuǎn)換原理
船載三維激光掃描系統(tǒng)中涉及到幾個坐標系之間的轉(zhuǎn)換,如圖1所示。
假設(shè)在KL時刻,其對應(yīng)地面點的WGS-84坐
2 系統(tǒng)誤差來源
船載三維激光誤差主要為各傳感器誤差和系統(tǒng)集成誤差。
2.1 傳感器相關(guān)誤差
2.1.1 GNSS定位誤差
GNSS主要誤差源見圖2。GNSS技術(shù)出現(xiàn)較早,相應(yīng)的誤差處理技術(shù)較為成熟,主要通過模型法、求差法、平差法予以消除或削弱[4-5]。
2.1.2 激光掃描儀誤差
按照誤差來源,激光掃描儀誤差分為儀器本身誤差、與信號傳播相關(guān)的外界、環(huán)境引起的誤差和反射目標特性引起的誤差,見圖3。
(1)儀器本身誤差。主要是指儀器制造加工過程中的機械誤差。這些誤差難以消除或減弱,但可通過儀器檢校、檢定等方式予以減弱或抵償。
(2)外界環(huán)境引起的誤差。主要是由于測距、測角環(huán)境引起的三維激光掃描儀誤差。其誤差影響同全站儀,注意選擇良好的作業(yè)環(huán)境。
(3)反射目標特性引起的誤差。三維激光掃描儀的激光反射目標特性,如大小、形狀、材質(zhì)等對激光掃描儀測程和測距造成不同程度的影響。由于反射面的種類千差萬別,盡管國內(nèi)外學(xué)者進行了相關(guān)研究,但難以消除其誤差影響,只有研究其變化規(guī)律進行改正,才能盡量減少反射目標特性引起的誤差[6]。
2.1.3 0ctans誤差
Octans是采用“Sagnac”效應(yīng)實現(xiàn)為載體提供六自由度運動姿態(tài)和真北方向,其誤差來源主要來自測區(qū)緯度與船速。
(1)緯度。Octans獲取的真北方向與測區(qū)緯度的正割值成正比[7-8]。其誤差影響見圖4。
(2)船速。羅經(jīng)定向、姿態(tài)誤差受靈敏度影響,船速均會對其產(chǎn)生誤差影響,在IS0 8728國際標準中規(guī)定:船速對航向角的影響可表述為V/5π(V為航速北向分量,單位為節(jié))。以O(shè)ctans的標稱精度±0.10為例,則對應(yīng)的速度誤差V=0.1×5π≈1.57節(jié)。而當Octans不輸入航速時,以實際測量中航速約6節(jié)為例,最大航向角誤差為6/5π≈ 0.38°。由此看出,有必要在測量過程中保持勻速測量,在船體調(diào)頭時北向速度變化最大,所以在測量過程中應(yīng)盡量保證沿直線行駛,船舶調(diào)頭時應(yīng)暫停記錄數(shù)據(jù)[6]。
削弱緯度和船速誤差簡單有效的方法為將GNSS定位數(shù)據(jù)(GPGGA)與對地速度(VTG)數(shù)據(jù)輸入Octans,通過對Octans輸入高精度船速和緯度數(shù)據(jù)達到提升其精度的目的。
(3)計算誤差。主要包括截斷誤差和舍人誤差。在姿態(tài)、航向解算過程中,會引入高階項的截斷誤差和解算傳播方程舍人誤差,可通過提高更新率削弱誤差。
2.2 與系統(tǒng)集成相關(guān)的誤差
2.2.1 安置誤差
船載三維激光掃描系統(tǒng)傳感器安置過程中,無法做到傳感器同心同軸安裝,安裝的偏心距和偏心角帶來安置誤差。通常Octans的中心為船載移動掃描系統(tǒng)中心,所以安置誤差主要指GNSS相位中心與Octans中心的偏心距,三維激光掃描儀與Octans中心偏心距和偏心角。安置誤差會使船載移動掃描系統(tǒng)數(shù)據(jù)結(jié)果出現(xiàn)系統(tǒng)性偏差,因此,系統(tǒng)作業(yè)前需要確定安置誤差。本文通過安裝平臺CAD設(shè)計圖及各傳感器中心與外型尺寸說明精確確定安置位置,可有效削弱安置誤差。
2.2.2 時間同步誤差
由于傳感器之間的差異,各傳感器的時間基準也不盡相同,需要將各自的時間系統(tǒng)統(tǒng)一到標準UTC時間系統(tǒng)上。在多源數(shù)據(jù)融合時,如果各傳感器時間系統(tǒng)不一致或存在偏差,就會對系統(tǒng)的解算結(jié)果產(chǎn)生影響。船載移動掃描系統(tǒng)集成的關(guān)鍵因素之一就是精密的時間同步。目前移動掃描系統(tǒng)的時間同步誤差控制在1ms范圍內(nèi)[9-10]。
2.2.3 坐標轉(zhuǎn)換誤差
坐標轉(zhuǎn)換誤差為船載移動掃描系統(tǒng)進行多源數(shù)據(jù)融合時,由于各傳感器坐標系轉(zhuǎn)換引起的誤差。此誤差一般影響很小,可忽略[11-12]。
2.2.4 內(nèi)插誤差
內(nèi)插誤差是由于船載移動掃描系統(tǒng)的各傳感器的采樣頻率不同,需要對原始觀測數(shù)據(jù)進行內(nèi)插,由數(shù)據(jù)內(nèi)插產(chǎn)生的誤差。目前,激光掃描儀的數(shù)據(jù)采樣頻率可達240 Hz,Octans的數(shù)據(jù)采樣頻率最高為200 Hz,而GNSS的數(shù)據(jù)采樣頻率最高為20Hz。因此,為了獲得每個激光腳點的位置信息和姿態(tài)信息,就需要對GNSS數(shù)據(jù)內(nèi)插。插值的算法很多,一般采用多項式插值計算。
3 系統(tǒng)誤差分析
3.1 GNSS誤差影響
船載三維激光掃描系統(tǒng)定位主要采用差分GPS技術(shù),其對船載移動掃描系統(tǒng)的誤差影響直接移植到系統(tǒng)的點云數(shù)據(jù)成果中。
3.2 激光掃描儀誤差影響
激光掃描儀誤差主要來自測距、測角誤差,因此激光掃描儀誤差對激光點云位置誤差表達式為
3.4 安置誤差影響
三維激光掃描儀設(shè)備安置過程中,由于激光掃描儀坐標軸系和船體坐標軸系不平行,導(dǎo)致兩坐標軸系存在3個偏心角,通過安置角校正可解算出偏心角。假設(shè)確定3個偏心角的誤差分別為Aa、△β、△γ(小角度值),可得到偏心角誤差旋轉(zhuǎn)矩陣。4系統(tǒng)精度估計
船載三維激光掃描系統(tǒng)最終點云數(shù)據(jù)成果精度主要與GNSS定位精度、Octans測角精度、三維激光掃描儀的測距、測角精度有關(guān),根據(jù)誤差傳播定律及各傳感器的標稱精度,對系統(tǒng)進行理論精度估計。
其中,GNSS系統(tǒng)和安置向量誤差是線元素,其精度直接傳遞到數(shù)據(jù)成果中。運動傳感器和安置角是角元素,對數(shù)據(jù)成果精度影響類似,其誤差傳播式如下:
4.2 系統(tǒng)精度估算
三維激光掃描儀測距精度取8mm@150m。由于船載三維激光為2D線掃描,故不含其水平角誤差,水平角φ取值為90°,動態(tài)測角精度取2.7”。GNSS定位水平精度取0.10 m,高程精度取0.15 m。同步誤差取1 ms。設(shè)備安置向量誤差2 mm。測船盡量保持勻速直線運動,內(nèi)插誤差按船速每節(jié)取1 cm。分別計算當掃描儀豎直方向夾角為450時,船速為2,4,6節(jié)對應(yīng)不同得掃描距離三維坐標精度見表1。
5 精度評定
以青藏高原格仁錯湖泊為例進行精度評定。
5.1 測區(qū)概況
本次測量高原湖泊為格仁錯湖區(qū)。格仁錯湖泊位于西藏自治區(qū)那曲地區(qū)申扎縣境內(nèi),北緯30°57- 31°19,平均海拔4650m。湖面長約60.0 km,最大寬14.0 km,平均寬7.93 km,面積475.9 km2。
5.2 數(shù)據(jù)采集及處理
船載三維激光掃描系統(tǒng)集成各硬件應(yīng)相對穩(wěn)定,且設(shè)備安裝平臺與船體應(yīng)保證牢固安裝。設(shè)備安裝平臺強度高,且不應(yīng)造成由船體震動而引發(fā)永久變形。通過平臺設(shè)計圖與各硬件的幾何中心確定,其相互位置關(guān)系,系統(tǒng)集成見圖5。
數(shù)據(jù)采集主要采用Qinsy軟件控制三維激光掃描儀V2-2000、運動傳感器OCTANS和相機的工作。V2-2000采用2D線掃描模式,量程設(shè)置為750m。船沿著水岸線前進過程中,距離岸邊約100 m,V2-2000采集點云數(shù)據(jù),Octans采集姿態(tài)、航向數(shù)據(jù),GNSS接收機提供導(dǎo)航定位信息及時間同步數(shù)據(jù)(GPGGA、ZDA、IPPS),同軸相機按照一定的時間間隔拍照記錄岸上數(shù)據(jù)特征;內(nèi)業(yè)數(shù)據(jù)處理時,首先由Qinsy軟件輸出經(jīng)校準后的點云數(shù)據(jù),導(dǎo)入三維點云數(shù)據(jù)后處理軟件測量RISCAN Pro,后用EPS三維模塊生成數(shù)字地形圖(DLG)成果,其流程見圖6。
5.3 精度評定
為評定系統(tǒng)精度,在點云中提取特征點,利用RTK測取點云“同名點”坐標,利用二者較差計算中誤差,計算式為式中,n為觀測個數(shù);Vi為觀測值較差。船載激光掃描系統(tǒng)與RTK測得坐標較差見表2。
根據(jù)式(28)及表2統(tǒng)計結(jié)果計算得出,X中誤差為+0.105 m,y中誤差為+0.079 m,平面位置S中誤差為+0.131 m,高程Z中誤差為+0.092 m。
6 結(jié)論
利用船載三維激光掃描系統(tǒng)進行高原湖泊岸上地形測量,闡述了船載三維激光掃描系統(tǒng)各傳感器誤差來源及削弱措施。分析了系統(tǒng)集成坐標轉(zhuǎn)換原理,并推到坐標轉(zhuǎn)換實現(xiàn)過程。利用各傳感器精度及誤差傳播定律,評估船載三維激光掃描系統(tǒng)誤差,并用實例數(shù)據(jù)進行精度評定,得出如下結(jié)論。
(1)提出的船載三維激光掃描系統(tǒng)各傳感器及其集成誤差來源及削弱措施可為系統(tǒng)集成和數(shù)據(jù)采集精度控制提供參考。
(2)根據(jù)船載三維激光掃描系統(tǒng)集成坐標轉(zhuǎn)換原理,得出系統(tǒng)坐標轉(zhuǎn)換式,利用誤差傳播定律得出系統(tǒng)誤差方程,對系統(tǒng)集成誤差評估及傳感器選擇提供依據(jù)。
(3)經(jīng)實例證實,船載三維激光掃描系統(tǒng)在高原地區(qū)適用性強,測量數(shù)據(jù)可靠且精度高,能夠滿足高原測湖精度要求。
參考文獻:
[1]周建紅,馬耀昌,劉世振,等,水陸地形三維一體化測量系統(tǒng)關(guān)鍵技術(shù)研究[J]人民長江,2017,48(24):61-65.
[2] 鄧神寶,沈清華,王小剛,船載激光三維掃描系統(tǒng)構(gòu)建與應(yīng)用[J].人民珠江,2016,37(10):23-26.
[3] 邊志剛,王東.船載水上水下一體化綜合測量系統(tǒng)技術(shù)與應(yīng)用[J].港工技術(shù),2017,54(1):109-112.
[4]劉占超,龐建成,劉百奇,一種改進的高精度POS時間同步方法[J].儀器儀表學(xué)報,2011,32 (10):2198-2203.
[5] 徐壽志.車載移動測量系統(tǒng)檢校技術(shù)及其評定方法[D].武漢:武漢大學(xué),2016.
[6]黃亞洲,曹忠良,航速及緯度誤差對Octans航向值的影響[J].中國水運,2018,18(6):5-6.
[7] 汪連賀.三維激光移動測量系統(tǒng)在海島礁測量中的應(yīng)用[J].海洋測繪,2015,35(5):79-82.
[8]陳尚登,楊俊凱,楊振林,等,基于船載三維激光掃描系統(tǒng)的航道沿岸大比例尺測圖[J].中國水運,2017,36(11):197-200.
[9]趙珞成.GPS輸出的NMEA0183信號的時間特性分析[J].海洋測繪,2002,22(3):19-21.
[10]趙珞成.利用CPS導(dǎo)航信號進行水下地形測量的時間同步[J].海洋測繪,2002,22(5):11-14.
[11]曹璽,梁俊,王延洪,等.三維激光掃描技術(shù)在地形測繪成圖中的應(yīng)用[J].人民長江,2011,41(21):22-24.
[12]李志鵬,張辛,喻守剛,等,基于三維激光掃描的大比例尺地形測繪方法研究[J]人民長江,2014.44(7): 70-73.
(編輯:唐湘茜)
作者簡介:石光,男,工程師,主要從事河道勘測工作。E-mail: 57371054l@qq.com