楊青山，湯學(xué)才，范彬彬，許正鵬

(1. 巴音郭楞蒙古自治州國土資源勘測規(guī)劃設(shè)計(jì)院，新疆 庫爾勒 841000； 2. 新疆巴州國源測繪規(guī)劃中心，新疆 庫爾勒 841000)

輕型機(jī)載LiDAR在農(nóng)村宅基地房屋測量中的應(yīng)用

楊青山1，湯學(xué)才1，范彬彬1，許正鵬2

(1. 巴音郭楞蒙古自治州國土資源勘測規(guī)劃設(shè)計(jì)院，新疆 庫爾勒 841000； 2. 新疆巴州國源測繪規(guī)劃中心，新疆 庫爾勒 841000)

選擇新疆焉耆縣北大渠鄉(xiāng)北大渠村房屋密集區(qū)作為研究區(qū)，通過三角翼搭載輕型機(jī)載LiDAR進(jìn)行了低空(170 m)交叉航線高密度點(diǎn)云數(shù)據(jù)采集，對點(diǎn)云數(shù)據(jù)進(jìn)行拼接、校正、紋理信息增強(qiáng)，在點(diǎn)云上對房屋進(jìn)行矢量化，并對結(jié)果進(jìn)行精度檢查。點(diǎn)云數(shù)據(jù)房角點(diǎn)采集率為83.3%，中誤差為4.8 cm。采用機(jī)載LiDAR測量房角點(diǎn)能夠大大減少外業(yè)房角點(diǎn)測量的工作量。機(jī)載LiDAR的應(yīng)用鮮有紋理信息的提取，本次測試通過對高密度點(diǎn)云數(shù)據(jù)進(jìn)行有效的處理，提取了點(diǎn)云數(shù)據(jù)紋理信息，為機(jī)載LiDAR數(shù)據(jù)紋理信息的應(yīng)用提供了參考，并對其精度有了明確的認(rèn)識(shí)，可為后期相關(guān)工作的開展提供借鑒。

機(jī)載LiDAR；房屋測量；紋理提?。痪确治?/p>

Abstract： Taking a built-up area in Beidaqu village, Beidaqu township, Yanqi County, Xinjiang Uygur Autonomous Region as the study area, we acquired high density point cloud data in a cross-route way at low altitude (170 m) with the delta wing equipped with light airborne LiDAR. Then, we spliced, corrected and enhanced texture information of the point cloud data. Finally, we vectorized the housing information with the processed point cloud data and checked the vectorization results’ accuracy. The results show that the collection rate of housing corner points is 83.3% with the point cloud data, and the median error is 4.8 cm. Using the airborne LiDAR, we can greatly reduce the workload of measuring the housing corner points. In view of rarely studies of texture information extracting with airborne LiDAR, this paper offers a reference for the application of it, and gives a clear understanding of its accuracy. In addition, the results can provide reference to the further related work.

Keywords： airborne LiDAR；building surveying；texture；precision analysis

2016年農(nóng)村地籍調(diào)查加入了房屋測量的內(nèi)容，要求采用解析法測量房屋角點(diǎn)坐標(biāo)，實(shí)地丈量房屋邊長，并采用幾何要素計(jì)算房屋面積，確保不動(dòng)產(chǎn)單元的界址清楚、面積準(zhǔn)確。農(nóng)村宅基地調(diào)查房角點(diǎn)測量精度按照《地籍調(diào)查規(guī)程》中界址點(diǎn)測量精度執(zhí)行，即房角點(diǎn)測量精度應(yīng)滿足中誤差≤5 cm[1]。采用航空攝影測量的技術(shù)手段獲取房角點(diǎn)坐標(biāo)，精度無法滿足要求；而采用GPS RTK或全站儀等測量手段完全實(shí)地測量，測量效率低，外業(yè)工作量巨大[2]。

機(jī)載LiDAR是一種集激光測距、全球定位系統(tǒng)(GPS)、慣性導(dǎo)航系統(tǒng)(INS)和CCD相機(jī)為一體的搭載在固定翼飛機(jī)或直升機(jī)上的激光探測和測距系統(tǒng)，精度高，能夠?qū)崟r(shí)獲取地形表面的三維空間信息[3-7]。然而受航空管制、飛行條件等限制，以飛機(jī)搭載LiDAR設(shè)備進(jìn)行測繪作業(yè)，在時(shí)間和成本方面遇到極大挑戰(zhàn)。輕型LiDAR具有體積小、重量輕的特點(diǎn)，適合搭載在三角翼或無人機(jī)等低空飛行、起飛條件簡易、航空管制少的航空器上[8-9]，在節(jié)省數(shù)據(jù)采集成本的同時(shí)，降低了機(jī)載LiDAR的技術(shù)門檻，受到越來越多的關(guān)注。

鑒于三角翼搭載輕型LiDAR具備精度高、飛行條件簡易等特點(diǎn)[10-12]，本文考慮將該技術(shù)應(yīng)用于農(nóng)村地籍調(diào)查房角點(diǎn)測量，以期減少外業(yè)房角點(diǎn)實(shí)地測量的工作量，提高工作效率。為明確輕型機(jī)載LiDAR應(yīng)用于農(nóng)村宅基地房屋調(diào)查的技術(shù)流程、效率和精度，選擇焉耆縣境內(nèi)北大渠鄉(xiāng)北大渠村為測試對象，通過三角翼搭載輕型LiDAR進(jìn)行超低空、高密度和高重疊采集數(shù)據(jù)，充分掌握輕型機(jī)載LiDAR的作業(yè)流程，明確其平面精度及限制因素，確定輕型機(jī)載LiDAR測量房角點(diǎn)精度能否滿足不動(dòng)產(chǎn)測量規(guī)范的精度要求，為相關(guān)的測繪工作提供借鑒。

1 測區(qū)狀況及儀器設(shè)備

1.1 測區(qū)概況

本次測試選擇焉耆縣境內(nèi)北大渠鄉(xiāng)北大渠村中央密集居民點(diǎn)作為測試區(qū)，北大渠村位于平原地帶，南北貫通，地勢平坦。居民點(diǎn)內(nèi)房屋密集，土房、老磚房、新建混房、多層樓房等多種房屋共存，為不同結(jié)構(gòu)類型的房屋測量提供了樣本，具有較好的農(nóng)村宅基地房屋調(diào)查代表性。航攝范圍如圖1所示，黑色選框?yàn)闇y試區(qū) ，測區(qū)總面積約5 km2。

圖1 航攝范圍

1.2 儀器設(shè)備

本次測試采用73.5 kW、最大載重250 kg的動(dòng)力三角翼作為航飛平臺(tái)，搭載集成了機(jī)載LiDAR激光發(fā)射系統(tǒng)、航空慣性導(dǎo)航系統(tǒng)、相機(jī)，總重量為13 kg的輕型機(jī)載LiDAR,平面定位精度為10 mm+1×10-6D，高程定位精度為20 mm+1×10-6D。

2 數(shù)據(jù)獲取

2.1 機(jī)載LiDAR數(shù)據(jù)采集方法

2.1.1 航飛高度的確定

機(jī)載LiDAR的量測精度主要由激光測距(S)、GPS定位(Xs、Ys、Zs)，以及IMU姿態(tài)(Ω、ω、κ)和掃描角(θ)等8個(gè)測量值的精度所決定。在各項(xiàng)參量量測誤差一定的情況下，各項(xiàng)系數(shù)是隨掃描角θ及測距值S而變化的[13]，掃描角誤差是固定的，可在出廠時(shí)測定，因此各項(xiàng)系數(shù)隨測距值S變化而變化，S值與航飛高度直接相關(guān)，航飛高度越低，S值越小，各項(xiàng)系數(shù)誤差越小。綜合考慮飛行安全的因素，將航飛高度設(shè)定為170 m。

2.1.2 航線設(shè)計(jì)

采用順著房屋排列方向和垂直房屋排列方向交叉飛行的方式進(jìn)行點(diǎn)云數(shù)據(jù)采集，通過低空高重疊度交叉航線飛行，盡可能提高了點(diǎn)云密度(機(jī)載三維激光掃描數(shù)據(jù)采集技術(shù)參數(shù)見表1)，充分獲取了房屋的紋理信息。

測區(qū)面積約5 km2，面積較小，只需在測區(qū)中心架設(shè)一個(gè)基站。在進(jìn)行飛行作業(yè)的同時(shí)，基站開機(jī)與三角翼上搭載的GPS同步觀測，且需在三角翼起飛前5 min打開GPS接收機(jī)，在三角翼停機(jī)后5 min關(guān)閉GPS接收機(jī)。

表1 數(shù)據(jù)采集技術(shù)參數(shù)

2.2 地面控制布設(shè)方法

點(diǎn)云數(shù)據(jù)無法具體標(biāo)記出特征點(diǎn)，由因此需要提取具有特征的斜面來控制點(diǎn)云的平面精度[14-15]?？刂茦?biāo)志結(jié)構(gòu)如圖2所示，由背靠背的2個(gè)正方形板構(gòu)成，板的尺寸為1 m×1 m，板與地面的角度可以調(diào)整為30°、45°和60°。測區(qū)內(nèi)一共布設(shè)控制標(biāo)志物10處，共布設(shè)20個(gè)斜面均勻分布于整個(gè)測區(qū)。

2.3 地面測量

新疆巴州CORS已覆蓋測區(qū)，本次地面控制測量采用網(wǎng)絡(luò)RTK方式，在測量地面控制點(diǎn)時(shí)，同時(shí)記錄WGS-84坐標(biāo)系和CGCS2000坐標(biāo)系成果，用于測量坐標(biāo)系轉(zhuǎn)換。

2.3.1 控制標(biāo)志的測量

在布設(shè)的控制面板上，挑選4個(gè)測量點(diǎn)的三維坐標(biāo)，保證4個(gè)點(diǎn)連線構(gòu)成的閉合4邊形面積大于0.5 m2，其分布形式如圖2所示。本次測試共測量斜面20個(gè)，用于點(diǎn)云的絕對精度校正。

圖2 斜面選點(diǎn)分布形式

2.3.2 檢查點(diǎn)的測量

挑選沒有房檐的房屋作為檢查點(diǎn)，測量時(shí)攀爬到房頂，將GPS立于房屋的最外檐進(jìn)行測量，共計(jì)測量均勻分布在測區(qū)的房角點(diǎn)36個(gè)。

3 數(shù)據(jù)處理

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

3.1.1 點(diǎn)云相對精度拼接

選取架構(gòu)航向，以架構(gòu)航線為基準(zhǔn)調(diào)整縱向航線的位置和姿態(tài)，使縱向航線點(diǎn)云和架構(gòu)航線點(diǎn)云完全融合到一起，進(jìn)一步調(diào)整其他的橫向航線，將全部點(diǎn)云拼接到一起。點(diǎn)云拼接調(diào)整順序和策略見表2。其中，北、東、高、滾轉(zhuǎn)、俯仰、方位分別簡稱N、E、H、R、P、Y。

表2 點(diǎn)云數(shù)據(jù)參數(shù)調(diào)整策略

3.1.2 點(diǎn)云絕對精度校正

用控制數(shù)據(jù)校正已拼接的點(diǎn)云，測量的20個(gè)斜面全部參與平差，殘差小于0.03 m，完成校正后點(diǎn)云數(shù)據(jù)完全融合，無分層現(xiàn)象，融合后多航帶點(diǎn)云與單一航帶點(diǎn)云厚度一致。

3.1.3 去除冗余點(diǎn)云數(shù)據(jù)

低空高重疊度交叉航線的數(shù)據(jù)采集方式獲取的點(diǎn)云數(shù)據(jù)密度大，大量點(diǎn)云對于本次測試為冗余數(shù)據(jù)。圖3(a)為點(diǎn)云投影到二維平面的頂視圖。冗余的點(diǎn)云將會(huì)嚴(yán)重影響矢量化數(shù)據(jù)加載速度，且不利于判斷建筑物的邊緣，因此將平面(X、Y方向)的點(diǎn)云數(shù)據(jù)進(jìn)行抽稀。X、Y方向抽稀為0.5 m/點(diǎn)，抽稀后的點(diǎn)云數(shù)據(jù)如圖3(b)所示。進(jìn)一步將高度低于1 m的地面點(diǎn)云數(shù)據(jù)刪除，利于后期建筑物邊緣的采集。剔除地面點(diǎn)云后數(shù)據(jù)頂視圖如圖3(c)所示。

圖3 點(diǎn)云投影到二維平面的頂視圖

3.2 建筑物界址點(diǎn)采集

將點(diǎn)云數(shù)據(jù)導(dǎo)入南方CASS軟件，房屋邊緣因立面(Z方向)點(diǎn)云疊加，形成了較為明顯的紋理特征，對這些紋理進(jìn)行矢量化，即得到了房屋的線劃圖，矢量化時(shí)取(Z方向)點(diǎn)云形成紋理中心。

3.3 坐標(biāo)轉(zhuǎn)換

機(jī)載LiDAR采用WGS-84坐標(biāo)系，而農(nóng)村宅基地房屋調(diào)查采用CGCS2000坐標(biāo)系，因此，實(shí)際應(yīng)用過程中，涉及WGS-84坐標(biāo)系與CGCS2000坐標(biāo)系轉(zhuǎn)換。本次測試項(xiàng)目區(qū)范圍較小，在測區(qū)的周邊和中心選取5組測量點(diǎn)，同時(shí)測量WGS-84坐標(biāo)和CGCS2000坐標(biāo)，用這5組點(diǎn)計(jì)算適用于本測區(qū)的小范圍四參數(shù)，然后將成果統(tǒng)一轉(zhuǎn)換為CGCS2000坐標(biāo)成果。

4 精度分析

本次測試共計(jì)測量36個(gè)房角點(diǎn)作為檢查點(diǎn)，均勻地分布于測試區(qū)范圍內(nèi)。因遮擋和遮蔽導(dǎo)致墻面點(diǎn)云稀疏，部分紋理信息丟失，房角點(diǎn)未能在點(diǎn)云上捕捉。36個(gè)檢查點(diǎn)中，共捕捉到30個(gè)房角點(diǎn)，采集率83.3%，中誤差0.048 m。檢查統(tǒng)計(jì)結(jié)果見表3。

表3 機(jī)載LiDAR測量精度統(tǒng)計(jì)

根據(jù)國土資發(fā)〔2014〕101號(hào)文件，農(nóng)村宅基地調(diào)查房角點(diǎn)測量精度按照《地籍調(diào)查規(guī)程》(TD/T 1001—2012)中界址點(diǎn)測量精度執(zhí)行。

參照《地籍調(diào)查規(guī)程》中界址點(diǎn)精度指標(biāo)，對機(jī)載LiDAR點(diǎn)云矢量化房角點(diǎn)結(jié)果進(jìn)行分析，統(tǒng)計(jì)結(jié)果見表4。本次測試房角點(diǎn)采集率為83.3%。其中41.7%(14個(gè))的房角點(diǎn)點(diǎn)位誤差≤5 cm，能夠滿足《地籍調(diào)查規(guī)程》一級(jí)界址點(diǎn)中誤差精度要求；77.8%(28個(gè))的房角點(diǎn)點(diǎn)位中誤差≤10 cm，能夠滿足《地籍調(diào)查規(guī)程》一級(jí)界址點(diǎn)允許誤差精度要求房；所有矢量化數(shù)據(jù)83.3%(30個(gè))房角點(diǎn)點(diǎn)位誤差均能夠滿足《地籍調(diào)查規(guī)程》二級(jí)界址點(diǎn)允許誤差精度要求。

表4 北大渠村房角點(diǎn)允許誤差精度統(tǒng)計(jì) (%)

5 討 論

上述測試方法中，測量房角點(diǎn)中誤差為4.8 cm，捕捉率為83.3%，應(yīng)用于農(nóng)村宅基地房屋調(diào)查能夠大量減少外業(yè)房角點(diǎn)的測量。采用低空、大重疊度交叉航線飛行，能夠提高點(diǎn)云精度，增強(qiáng)房屋的紋理信息，但航飛效率較低，將該測試方法應(yīng)用于大面積的農(nóng)村地籍調(diào)查項(xiàng)目需飛行員在低空環(huán)境下長時(shí)間作業(yè)，風(fēng)險(xiǎn)極高。機(jī)載LiDAR的量測精度影響因子中， IMU姿態(tài)(Ω、ω、κ)對機(jī)載LiDAR精度影響最大，制約著航飛高度，隨著輕小型慣導(dǎo)技術(shù)的進(jìn)一步發(fā)展，能夠在滿足精度的前提下增加航飛高度，實(shí)現(xiàn)大面積安全作業(yè)。采用地面布設(shè)人工控制標(biāo)志物可提高點(diǎn)云拼接的精度，但布設(shè)過程需花費(fèi)大量的時(shí)間和人力，并且標(biāo)志物需要在航飛之前完成，制約著航飛的進(jìn)度。測試中發(fā)現(xiàn)測區(qū)中分布著自然的標(biāo)志物，如路基斜坡、渠面等，可通過后期在測區(qū)挑選自然地物作為控制標(biāo)識(shí)，完成點(diǎn)云數(shù)據(jù)的控制拼接。房角點(diǎn)精度滿足1∶500大比例尺地形圖測圖要求,將機(jī)載LiDAR應(yīng)用于大比例地形圖測量，能夠有效減少界址點(diǎn)測量的工作量，提高工作效率。隨著測繪地理信息技術(shù)的發(fā)展，掃描的三維點(diǎn)云信息可在今后的不動(dòng)產(chǎn)測繪、地籍三維模型建立中發(fā)揮積極作用，實(shí)現(xiàn)一次數(shù)據(jù)采集，重復(fù)數(shù)據(jù)利用。

6 結(jié) 論

本次測試選擇焉耆縣北大渠鄉(xiāng)北大渠村房屋密集區(qū)作為研究區(qū)，采用三角翼搭載輕型機(jī)載LiDAR低空(170 m)交叉航線高密度點(diǎn)云數(shù)據(jù)采集，地面合理布設(shè)人工控制標(biāo)志，對點(diǎn)云數(shù)據(jù)拼接校正，抽稀平面(X、Y)方向的點(diǎn)云數(shù)據(jù)(0.5 m/點(diǎn))，并將相對高度低于1 m的點(diǎn)云數(shù)據(jù)刪除，將處理好的點(diǎn)云數(shù)據(jù)導(dǎo)入南方CASS軟件通過點(diǎn)云數(shù)據(jù)的紋理特征進(jìn)行矢量化，利用測區(qū)內(nèi)均勻分布的36個(gè)檢查點(diǎn)，參照《地籍調(diào)查規(guī)程》對結(jié)果進(jìn)行了比對分析，得到以下結(jié)論：

(1) 點(diǎn)云數(shù)據(jù)紋理信息弱，通過低空、高重疊、交叉航線的方式增加點(diǎn)云密度，保證點(diǎn)云的精度，立面(Z)方向的點(diǎn)云可形成能夠識(shí)別的紋理。因遮擋和遮蔽導(dǎo)致墻面點(diǎn)云稀疏，部分紋理信息丟失，部分房角點(diǎn)未能在點(diǎn)云上捕捉，房角點(diǎn)采集率為83.3%。

(2) 機(jī)載LiDAR測量房角點(diǎn)中誤差為4.8 cm，根據(jù)《地籍調(diào)查規(guī)程》， 41.7%的房角點(diǎn)滿足一級(jí)界址點(diǎn)中誤差精度要求，77.8%滿足一級(jí)界址點(diǎn)允許誤差精度要求，所有矢量化數(shù)據(jù)中83.3%均能夠滿足二級(jí)界址點(diǎn)允許誤差精度要求。

(3) 上述測試方法雖然夠大量減少外業(yè)房角點(diǎn)測量的工作量，降低人員攀爬房屋的風(fēng)險(xiǎn)，但航飛效率較低，不具有大規(guī)模作業(yè)的可行性。機(jī)載LiDAR的應(yīng)用鮮有紋理信息的提取，本次測試通過獲取高密度點(diǎn)云數(shù)據(jù)進(jìn)行有效的處理，提取點(diǎn)云數(shù)據(jù)紋理信息矢量化，為機(jī)載LiDAR數(shù)據(jù)紋理信息的應(yīng)用提供了參考，并對其精度有了明確的認(rèn)識(shí)，可為后期相關(guān)工作的開展提供借鑒。

[1] 李君，李貴中，李樹云.應(yīng)用無人機(jī)數(shù)字影像進(jìn)行農(nóng)村宅基地確權(quán)技術(shù)研究[J].測繪通報(bào)，2016(4)：83-84.

[2] 范彬彬，楊青山，包琦，等.農(nóng)村地籍調(diào)查中輪臺(tái)1∶1000真正射影像生產(chǎn)及其應(yīng)用評估[J].測繪通報(bào)，2016(10)：76-80.

[3] 覃昕垚，張建軍，王勇，等.機(jī)載liDAR定位精度分析[J].測繪工程，2016，25(5)：32-35.

[4] 張衡，李通，肖鵬，等.機(jī)載LiDAR技術(shù)在山區(qū)高速公路勘測設(shè)計(jì)中的研究與應(yīng)用[J].北京測繪，2016(6)：95-98.

[5] 黃克城，宋時(shí)文，閻鳳霞.機(jī)載LiDAR技術(shù)在地形圖測繪中的應(yīng)用[J].地理空間信息，2016，14(4)：99-101.

[6] 陳良良，隋立春，蔣濤，等.地面三維激光掃描數(shù)據(jù)配準(zhǔn)方法[J].測繪通報(bào)，2014(5)：80-82.

[7] AXELSSON P.Processing of Laser Scanner Data——algorithms and Applications[J]. Journal of Photogrammetry and Romote Sensing,1999(54):138-147.

[8] 陳君朝，劉新舟，陳楚江.高速公路改擴(kuò)建精密機(jī)載三維激光掃描測量方法[J].公路交通科技,2015(1)：4-7.

[9] 石德斌.輕型LiDAR設(shè)備在鐵路工程勘測中的應(yīng)用分析[J].鐵道工程學(xué)報(bào)，2016，33(8):30-33.

[10] 張玉方，程新文，歐陽平，等.機(jī)載LiDar數(shù)據(jù)處理及其應(yīng)用綜述[J].工程地球物理學(xué)報(bào)，2008，5(1)：119-124.

[11] 袁豹.基于平面約束的機(jī)載LiDAR航帶平差方法研究[J].測繪工程，2015，24(10)：27-30.

[12] 史照良，曹敏.基于LiDAR技術(shù)的海島礁、灘涂測繪研究[J]. 測繪通報(bào)，2007(5):49-52.

[13] 賴旭東.機(jī)載激光雷達(dá)基礎(chǔ)原理與應(yīng)用[M].北京：電子工業(yè)出版社，2010.

[14] 徐辛超，徐愛功，于丹.地面三維激光掃描點(diǎn)云拼接影像因素分析[J].測繪通報(bào)，2017(2)：14-18.

[15] 徐壽志，程鵬飛，張玉，等.地面三維激光掃描儀的檢校與測量精度評定[J].測繪通報(bào)，2016(2)：79-83.

TheApplicationofLightAirborneLiDARonRuralHomesteadBuildingSurveying

YANG Qingshan1，TANG Xuecai1，F(xiàn)AN Binbin1，XU Zhengpeng2

(1. Bazhou Land Resources Surveying and Planning Institute,Korla 841000，China; 2. Xinjiang Bazhou State Source Surveying and Mapping Planning Center,Korla 841000，China)

P27

A

0494-0911(2017)09-0088-04

2017-07-06

楊青山(1970—)，男，高級(jí)工程師，主要從事不動(dòng)產(chǎn)測繪、工程測量和攝影測量方面的工作。E-mail：137268099@qq.com

楊青山，湯學(xué)才，范彬彬，等.輕型機(jī)載LiDAR在農(nóng)村宅基地房屋測量中的應(yīng)用[J].測繪通報(bào)，2017(9)：88-91.

10.13474/j.cnki.11-2246.2017.0294.