田正杰,王雪英,渠 濤,李澤暉

(1. 山西能源學(xué)院,山西 太原 030600; 2. 山西泰翔科技有限公司,山西 太原 030006)

近年來,我國(guó)對(duì)新能源發(fā)電領(lǐng)域較為重視,大力支持新能源發(fā)電與儲(chǔ)能項(xiàng)目,其中光伏發(fā)電和風(fēng)力發(fā)電占據(jù)了重要位置。光伏發(fā)電項(xiàng)目具有占地面積大,以及對(duì)地表坡度、坡向要求高的特點(diǎn)。為合理布設(shè)光伏組件位置,需要提供光伏場(chǎng)區(qū)精確的地形、地貌特征信息,以及地物、地類分布情況;對(duì)光伏組件布置起關(guān)鍵作用的微地形變化情況,更需要在地形圖上詳細(xì)表達(dá),因此需要測(cè)繪1∶500比例尺地形圖。目前,測(cè)繪光伏場(chǎng)區(qū)地形圖主要以低空無人機(jī)航空攝影測(cè)量方法為主。該方法在地表具有植被遮擋的區(qū)域(如林地、玉米地、谷子地、經(jīng)濟(jì)作物等),進(jìn)行地物特征點(diǎn)、高程點(diǎn)數(shù)據(jù)采集困難較大,且數(shù)據(jù)采集人員的技術(shù)參差不齊,易造成較大誤差,給光伏電站設(shè)計(jì)和施工帶來諸多不便。

激光雷達(dá)系統(tǒng)具有多次回波功能,能夠在一定程度上穿透地表遮擋物采集地表點(diǎn)云數(shù)據(jù),并通過點(diǎn)云濾波處理,得到真實(shí)地表信息,生成高精度DEM[1],從而如實(shí)表達(dá)地表微地形特征。此外,無人機(jī)載LiDAR系統(tǒng)高度集成的高清數(shù)碼相機(jī),可在掃描過程中拍攝滿足攝影測(cè)量重疊度要求的高分辨率數(shù)碼像片。內(nèi)業(yè)數(shù)據(jù)處理時(shí),空三解算結(jié)果滿足1∶500比例尺地形圖精度要求后生成高精度DOM,并在DOM的基礎(chǔ)上對(duì)地表地物與地類進(jìn)行準(zhǔn)確判讀。

1 技術(shù)路線

無人機(jī)載LiDAR系統(tǒng)集激光雷達(dá)掃描系統(tǒng)、全球衛(wèi)星系導(dǎo)航統(tǒng)(global navigation satellite system, GNSS)、慣性導(dǎo)航系統(tǒng)(inertial navigation system, INS)、高清數(shù)碼成像系統(tǒng)于一體,對(duì)地面進(jìn)行激光掃描與航拍,獲取地表的三維空間信息與屬性信息[2-5]。機(jī)載LiDAR是一種主動(dòng)遙感設(shè)備,受外界環(huán)境影響較小,能全天候、全天時(shí)工作,可部分穿透植被覆蓋,能迅速、準(zhǔn)確地獲取各類區(qū)域的空間掃描數(shù)據(jù)[6-9],并通過數(shù)據(jù)預(yù)處理獲得激光點(diǎn)云的三維數(shù)據(jù)[10]。

無人機(jī)載LiDAR測(cè)繪1∶500地形圖的過程主要包括外業(yè)數(shù)據(jù)采集、內(nèi)業(yè)數(shù)據(jù)處理及繪制成圖3部分。外業(yè)數(shù)據(jù)采集方案要兼顧內(nèi)業(yè)數(shù)據(jù)處理方法,以保證最終成果精度符合要求。內(nèi)業(yè)數(shù)據(jù)處理包括點(diǎn)云數(shù)據(jù)處理和像片空三解算兩部分。點(diǎn)云數(shù)據(jù)處理可生成符合精度要求的DEM,像片空三解算可生成符合精度要求的DOM[11]。繪制成圖的過程則是利用生成的DEM和DOM,在三維測(cè)圖軟件中采集地物特征點(diǎn)、高程點(diǎn)等,按照《國(guó)家基本比例尺地圖圖式 第1部分:1∶500 1∶1000 1∶2000地形圖圖式》生成1∶500比例尺地形圖[12]。技術(shù)路線如圖1所示。

2 案例與分析

2.1 研究區(qū)域

本文研究選擇的光伏發(fā)電場(chǎng)區(qū)位于山西省臨汾市安澤縣東南側(cè),距安澤縣直線距離為15 km,屬于山區(qū)林光互補(bǔ)和農(nóng)光互補(bǔ)項(xiàng)目,海拔高度為890～1200 m,坡度為25°～60°,地質(zhì)情況以黃褐色粉質(zhì)黏土和泥質(zhì)砂巖為主,受雨水沖刷影響,微地形發(fā)育廣泛。光伏發(fā)電場(chǎng)區(qū)的地類以楊樹林、松樹林、低矮灌木林、耕地、草地、裸地為主,經(jīng)現(xiàn)場(chǎng)踏勘并結(jié)合衛(wèi)星遙感影像,該區(qū)域地表植被覆蓋情況見表1。

表1 測(cè)區(qū)地表覆蓋情況

光伏發(fā)電場(chǎng)區(qū)受地形條件,以及不能占用地類和生態(tài)紅線等光伏組件設(shè)計(jì)要求的制約,導(dǎo)致測(cè)區(qū)較分散,平均每塊測(cè)區(qū)面積為2～5 km2。場(chǎng)地交通情況一般,受山區(qū)地形影響,測(cè)區(qū)內(nèi)網(wǎng)絡(luò)信號(hào)分布不均,大部分區(qū)域信號(hào)較差。

2.2 外業(yè)數(shù)據(jù)采集方案及與鍵點(diǎn)設(shè)置

2.2.1 無人機(jī)載LiDAR的選擇

光伏發(fā)電場(chǎng)區(qū)地表多被不同植被所覆蓋,對(duì)測(cè)區(qū)內(nèi)不同植被覆蓋密度(見表1株距)進(jìn)行調(diào)查發(fā)現(xiàn),大部分植被覆蓋留有一定的間隙,但間隙較小。若采用普通無人機(jī)航測(cè)方法較難準(zhǔn)確采集地面點(diǎn)的空間信息(對(duì)其中某一測(cè)區(qū)進(jìn)行試驗(yàn),成圖結(jié)果未能滿足1∶500比例尺地形圖精度要求)。為保證數(shù)據(jù)采集精度,本文采用低空無人機(jī)載LiDAR的方法獲取地面信息,它可通過植被間間隙掃描到地面點(diǎn)。由于測(cè)區(qū)較分散且單個(gè)測(cè)區(qū)面積較小,綜合考慮技術(shù)與經(jīng)濟(jì)因素,本文采用大疆經(jīng)緯M300 RTK旋翼機(jī),搭載禪思L1激光雷達(dá)進(jìn)行數(shù)據(jù)采集。

2.2.2 起降場(chǎng)地的選擇

起降場(chǎng)地的選擇關(guān)系到數(shù)據(jù)采集能否高效完成。為防止激光掃描過程中UAV與控制器失聯(lián),以及保證起降時(shí)的安全性,起降場(chǎng)地盡量選擇測(cè)區(qū)高處較平坦的位置,且確保視野開闊無遮擋;為便于將UAV運(yùn)送到起降位置,起降場(chǎng)地周圍應(yīng)具有良好的交通條件;為保證LiDAR系統(tǒng)GNSS差分信號(hào)穩(wěn)定可靠,起降場(chǎng)地周圍應(yīng)具有良好的通信網(wǎng)絡(luò)信號(hào)。本文進(jìn)行激光掃描時(shí),LiDAR系統(tǒng)GNSS選擇接入千尋位置服務(wù)網(wǎng)絡(luò)接收差分信號(hào),定位精度可達(dá)厘米級(jí)。

2.2.3 像控點(diǎn)、檢查點(diǎn)布設(shè)

由于LiDAR系統(tǒng)接入GNSS差分信號(hào),并具有IMU高精度慣導(dǎo)數(shù)據(jù),點(diǎn)云與像片外方位元素精度可靠,可大幅減輕外業(yè)布設(shè)像控點(diǎn)的工作量。由于測(cè)區(qū)面積較小,每個(gè)測(cè)區(qū)布設(shè)的像控點(diǎn)數(shù)量滿足控制全局即可,約為6～10個(gè),分別布設(shè)于測(cè)區(qū)四角附近、測(cè)區(qū)較長(zhǎng)邊界的中間位置及測(cè)區(qū)中心區(qū)域。為在后期對(duì)成圖精度進(jìn)行檢驗(yàn),需在測(cè)區(qū)內(nèi)均勻布設(shè)充足的檢查點(diǎn)。布設(shè)像控點(diǎn)與檢查點(diǎn)采用拓普康HiPer VR型GNSS設(shè)備,并接入千尋位置服務(wù)網(wǎng)絡(luò)接收差分信號(hào)。由于測(cè)區(qū)內(nèi)網(wǎng)絡(luò)信號(hào)分布不均勻,某些區(qū)域信號(hào)較差,接收差分信號(hào)較困難,該設(shè)備配備的超遠(yuǎn)距離藍(lán)牙技術(shù),解決了網(wǎng)絡(luò)信號(hào)差給布設(shè)像控點(diǎn)和檢查點(diǎn)帶來的困擾。選擇兩個(gè)試驗(yàn)測(cè)區(qū),根據(jù)范圍和需要,在測(cè)區(qū)1布設(shè)7個(gè)像控點(diǎn),測(cè)區(qū)2布設(shè)8個(gè)像控點(diǎn)。為對(duì)成圖精度進(jìn)行檢驗(yàn),分別在每個(gè)測(cè)區(qū)均勻布設(shè)100個(gè)平高檢查點(diǎn)。

2.2.4 外業(yè)航飛關(guān)鍵點(diǎn)設(shè)置

為獲取高質(zhì)量的點(diǎn)云數(shù)據(jù),考慮點(diǎn)云精度、點(diǎn)云密度及點(diǎn)云分布特征等重要指標(biāo),對(duì)禪思L1激光雷達(dá)性能進(jìn)行綜合分析,發(fā)現(xiàn)L1激光雷達(dá)獲取的點(diǎn)云精度主要與表2列舉的因素相關(guān)。

表2 無人機(jī)載LiDAR獲取的原始點(diǎn)云精度影響分析

由于測(cè)區(qū)大部分面積被不同種類植被所覆蓋,且每種植被的地表覆蓋密度不同。經(jīng)前期踏勘,數(shù)據(jù)采集時(shí)可通過大部分植被間間隙掃描到地面點(diǎn)。但有些植被間隙較小,為獲取充足的地面激光腳點(diǎn),需要增大點(diǎn)云密度,并采用圓錐掃描方式進(jìn)行多次回波掃描。點(diǎn)云密度與航高、航速、采樣頻率、激光脈沖回波數(shù)、點(diǎn)云重疊度等設(shè)置有關(guān),根據(jù)光伏場(chǎng)區(qū)測(cè)圖精度要求及L1激光雷達(dá)特性設(shè)置具體的LiDAR航飛參數(shù),見表3。

表3 無人機(jī)載LiDAR航飛技術(shù)參數(shù)設(shè)置

2.3 內(nèi)業(yè)數(shù)據(jù)處理方案及關(guān)鍵技術(shù)分析

內(nèi)業(yè)數(shù)據(jù)處理包括兩個(gè)主要任務(wù):一是點(diǎn)云數(shù)據(jù)處理,獲取高精度DEM;二是對(duì)航空像片進(jìn)行空中三角測(cè)量解算,獲取高精度DOM。

2.3.1 點(diǎn)云數(shù)據(jù)處理

點(diǎn)云數(shù)據(jù)預(yù)處理采用大疆智圖DJI Terra軟件完成,將原始的點(diǎn)云掃描數(shù)據(jù)、GNSS定位數(shù)據(jù)、IMU定姿數(shù)據(jù)、系統(tǒng)參數(shù)及校準(zhǔn)數(shù)據(jù)等導(dǎo)入DJI Terra軟件中進(jìn)行計(jì)算與坐標(biāo)轉(zhuǎn)換,生成LAS標(biāo)準(zhǔn)格式的點(diǎn)云數(shù)據(jù)。

點(diǎn)云后處理采用數(shù)字綠土LiDAR360 UAV軟件完成。將LAS格式的點(diǎn)云數(shù)據(jù)導(dǎo)入軟件進(jìn)行拼接、坐標(biāo)轉(zhuǎn)換、濾波及分類、人工處理、抽稀等處理后得到地面點(diǎn)云數(shù)據(jù)[13-14]。點(diǎn)云后處理階段的主要環(huán)節(jié)是點(diǎn)云濾波,其作用是將地面點(diǎn)與非地面點(diǎn)分離。由于本文的主要目的是生成高精度DEM,因此地面點(diǎn)數(shù)據(jù)是關(guān)注的要點(diǎn)。采用改進(jìn)的漸進(jìn)加密三角網(wǎng)濾波算法對(duì)點(diǎn)云進(jìn)行濾波處理[15],可有效改進(jìn)點(diǎn)云濾波效果,能夠基本過濾掉大部分林地、玉米稈、雜草、建筑物等地物,但對(duì)于仍無法有效濾除的低矮密集灌木林,需要人工處理,對(duì)于少量生長(zhǎng)茂密的成林區(qū)地面點(diǎn)缺失的情況,需要人工補(bǔ)測(cè)。

2.3.2 點(diǎn)云缺失區(qū)域的補(bǔ)測(cè)方案

通過對(duì)濾波處理后的地面點(diǎn)云進(jìn)行檢查,發(fā)現(xiàn)測(cè)區(qū)內(nèi)有少量地面點(diǎn)缺失的區(qū)域,分析原因并針對(duì)不同情況分別制定了補(bǔ)測(cè)方案。

(1)針對(duì)茂密低矮灌木林導(dǎo)致的地面點(diǎn)不足的問題,采用拓普康HiPer VR型GNSS接收機(jī)接入千尋位置服務(wù)網(wǎng)絡(luò)接收差分信號(hào),利用RTK測(cè)量技術(shù)進(jìn)行補(bǔ)測(cè)。該接收機(jī)支持全星座衛(wèi)星信號(hào),內(nèi)置的IMU和電子羅盤支持傾斜測(cè)量技術(shù),可有效提高補(bǔ)測(cè)精度。此外,地面點(diǎn)缺失的區(qū)域多為網(wǎng)絡(luò)信號(hào)較差的區(qū)域,為有效接入差分信號(hào),設(shè)備采用超遠(yuǎn)距離藍(lán)牙技術(shù),將SIM卡插入手簿,在無信號(hào)遮擋的情況下,可在200 m范圍內(nèi)將差分信息傳給GNSS接收機(jī),測(cè)量時(shí)手簿操作員位于通信網(wǎng)絡(luò)信號(hào)好的區(qū)域,GNSS實(shí)地測(cè)量員進(jìn)入補(bǔ)測(cè)現(xiàn)場(chǎng),通過二人協(xié)調(diào)配合,可有效提高作業(yè)效率。

(2)針對(duì)成林區(qū)地面點(diǎn)缺失區(qū)域的補(bǔ)測(cè),由于成林樹冠遮擋了GNSS衛(wèi)星信號(hào),因此采用三維激光掃描儀完成補(bǔ)測(cè)。研究區(qū)少量的楊樹林與松樹林樹冠長(zhǎng)勢(shì)茂盛,缺少植被間隙,機(jī)載激光掃描難以獲取到充足的地面點(diǎn)。但樹干中下部通視條件較好,株距為2～3 m,且無其他植被生長(zhǎng),較適合地面激光掃描儀掃描。通過現(xiàn)場(chǎng)踏勘與詳細(xì)設(shè)計(jì),選擇無信號(hào)遮擋且有利于地面掃描的設(shè)站點(diǎn)與定向點(diǎn)位置,利用RTK采集設(shè)站點(diǎn)坐標(biāo)和定向點(diǎn)坐標(biāo),將三維激光掃描儀架設(shè)在設(shè)站點(diǎn)上定向后進(jìn)行掃描,這樣保證了三維激光掃描點(diǎn)云數(shù)據(jù)與無人機(jī)載LiDAR掃描點(diǎn)云數(shù)據(jù)的坐標(biāo)系統(tǒng)的統(tǒng)一,避免了點(diǎn)云匹配環(huán)節(jié)造成的誤差。三維激光掃描儀具有測(cè)程長(zhǎng)、精度高、定向及轉(zhuǎn)站方便的特點(diǎn),將補(bǔ)測(cè)點(diǎn)云與航飛點(diǎn)云相結(jié)合,有效解決了地面點(diǎn)缺失問題,效果良好。

2.3.3 地面點(diǎn)云精度檢驗(yàn)

將補(bǔ)測(cè)數(shù)據(jù)導(dǎo)入軟件與LiDAR地面點(diǎn)云結(jié)合,重新濾波處理獲得最終完整的地面點(diǎn)云數(shù)據(jù)。為檢驗(yàn)地面點(diǎn)云精度,在測(cè)區(qū)1、測(cè)區(qū)2內(nèi)分別均勻布設(shè)100個(gè)平高檢查點(diǎn)用于精度檢驗(yàn)。地面點(diǎn)云高程精度分析結(jié)果見表4。高程誤差點(diǎn)狀分布如圖2所示。

表4 地面點(diǎn)云高程精度統(tǒng)計(jì)分析結(jié)果 m

圖2 高程誤差分布

高程精度分析結(jié)果滿足《工程測(cè)量標(biāo)準(zhǔn)》中有關(guān)1∶500比例尺數(shù)字線劃圖高程精度的要求[16],從而利用地面點(diǎn)云生成高精度DEM產(chǎn)品。

2.3.4 航片空中三角測(cè)量

航片空中三角測(cè)量采用大疆智圖DJI Terra軟件完成。由于無人機(jī)載LiDAR數(shù)據(jù)采集過程中獲取了高精度的像片外方位元素,空三解算時(shí)選擇控制全局的少量像控點(diǎn)即可保證整體精度(測(cè)區(qū)1布設(shè)7個(gè)像控點(diǎn),測(cè)區(qū)2布設(shè)8個(gè)像控點(diǎn))。為檢驗(yàn)空三處理精度,將每個(gè)測(cè)區(qū)均勻布設(shè)的100個(gè)平高檢查點(diǎn)代入,一并進(jìn)行空三解算。由測(cè)區(qū)1與測(cè)區(qū)2空三解算報(bào)告統(tǒng)計(jì)分析的平面精度結(jié)果見表5。

表5 平面點(diǎn)位誤差統(tǒng)計(jì)分析結(jié)果 m

根據(jù)以上平面精度分析結(jié)果,空三解算精度滿足《工程測(cè)量標(biāo)準(zhǔn)》中有關(guān)1∶500比例尺數(shù)字線劃圖平面精度的要求[16],從而生成高精度的DOM產(chǎn)品。

2.4 繪制成圖

利用生成的DEM和DOM在EPS軟件中建立三維模型,進(jìn)行地物特征點(diǎn)和高程特征點(diǎn)的采集,生成1∶500比例尺地形圖。由于DEM模型中已經(jīng)過濾掉了地表植被的遮擋,因此提高了地面高程特征點(diǎn)采集的精度,避免了常規(guī)無人機(jī)航測(cè)三維模型中地表植被對(duì)真實(shí)地面的干擾,有效提高了成圖精度。圖3為測(cè)區(qū)1生成的DEM、DOM、DLG產(chǎn)品。

圖3 測(cè)區(qū)1產(chǎn)品

3 結(jié) 語

本文根據(jù)光伏場(chǎng)區(qū)現(xiàn)場(chǎng)實(shí)際地形、地貌及植被覆蓋情況,制定合理的測(cè)圖方案,采用無人機(jī)載LiDAR進(jìn)行光伏場(chǎng)區(qū)1∶500地形圖測(cè)繪,避免了測(cè)區(qū)地表被林地、玉米地、雜草等覆蓋對(duì)數(shù)據(jù)采集精度造成的影響,減少了外業(yè)工作量,提高了測(cè)圖精度和工作效率,為光伏電站設(shè)計(jì)提供了有力支撐,在技術(shù)提升和推動(dòng)行業(yè)應(yīng)用方面具有一定的現(xiàn)實(shí)意義。