何 靜，何忠煥

（1.國家測繪地理信息局 重慶測繪院，重慶 400015）

LiDAR技術(shù)已成為大比例尺測繪產(chǎn)品生產(chǎn)環(huán)節(jié)中一項(xiàng)重要技術(shù)手段。它是采用主動遙感方式，集成GPS、慣性導(dǎo)航、激光測距等先進(jìn)技術(shù)，通過發(fā)射激光脈沖，再由LiDAR系統(tǒng)上的接收單元收集地物的海量反射信號，從而實(shí)現(xiàn)三維地面信息的測繪與定位[1]。與傳統(tǒng)遙感技術(shù)相比，具有對控制測量依賴少、受天氣影響小、自動化程度高、成圖周期短等特點(diǎn)。

1 基于TerraSoLid的LiDAR數(shù)據(jù)處理

本次對LiDAR數(shù)據(jù)的處理在TerraSoLid軟件中進(jìn)行，基于MicroStation平臺開發(fā)，該平臺提供了視圖操作、矢量編輯等功能[2]。TerraSoLid主要包 含 了Terrascan、Terramodel、Terraphoto等 模 塊。LiDAR點(diǎn)云的預(yù)處理、坐標(biāo)轉(zhuǎn)換、點(diǎn)云分類主要使用的是Terrascan及Terramodel,正射影像制作主要利用Terrascan及Terraphoto進(jìn)行。

1.1 DEM制作

1.1.1 點(diǎn)云預(yù)處理

點(diǎn)云數(shù)據(jù)預(yù)處理主要包含獲取的不同航帶點(diǎn)云數(shù)據(jù)的航帶間糾正、坐標(biāo)系統(tǒng)轉(zhuǎn)換以及高程糾正。主要方式是通過人機(jī)交互檢查或者程序自動檢查并糾正的方式進(jìn)行。

1）航帶糾正：由于受到飛機(jī)姿態(tài)、飛行高度等因素影響，不同航帶間的點(diǎn)云數(shù)據(jù)存在高程差異，如果高程較差超出成圖比例尺的數(shù)據(jù)生產(chǎn)允許范圍，則應(yīng)進(jìn)行航帶間點(diǎn)云數(shù)據(jù)糾正。

2）坐標(biāo)系統(tǒng)轉(zhuǎn)換：點(diǎn)云獲取采用的坐標(biāo)與GPS一致，為WGS84地心坐標(biāo)。根據(jù)項(xiàng)目要求，應(yīng)該將點(diǎn)云坐標(biāo)向目標(biāo)數(shù)學(xué)基礎(chǔ)進(jìn)行轉(zhuǎn)換。此處的坐標(biāo)轉(zhuǎn)換包括2個方面，一方面是基于框架的二維坐標(biāo)轉(zhuǎn)換，另一方面是基于不同水準(zhǔn)面的高程值擬合。在Terrascan的坐標(biāo)變換模塊，用于坐標(biāo)轉(zhuǎn)換的控制點(diǎn)位應(yīng)均勻分布于測區(qū)的四周及測區(qū)的中部，坐標(biāo)轉(zhuǎn)換控制點(diǎn)數(shù)量不少于4個。首先定義項(xiàng)目的目標(biāo)投影，主要設(shè)置內(nèi)容包括橢球的長半軸、扁率、七參數(shù)（即目標(biāo)橢球相對于WGS的線元素偏移量、角元素偏移量以及縮放因子）、中央經(jīng)線、東偏及長度縮放比例等。通過預(yù)設(shè)參數(shù)實(shí)現(xiàn)不同基準(zhǔn)的數(shù)據(jù)轉(zhuǎn)換?；诓煌疁?zhǔn)面的高程值擬合一般有2種方法：一種是使用覆蓋測區(qū)若干高程控制點(diǎn)對應(yīng)的同名點(diǎn)云建立轉(zhuǎn)換關(guān)系，從而實(shí)現(xiàn)不同高程基準(zhǔn)的高程數(shù)據(jù)擬合；另一種是通過已知似大地水準(zhǔn)面對應(yīng)測區(qū)的高程異常，直接擬合獲取。本項(xiàng)目采用第一種方式實(shí)現(xiàn)不同水準(zhǔn)面的高程值變換。

1.1.2 點(diǎn)云分類

完成預(yù)處理后的LiDAR數(shù)據(jù)通過檢查合格后開展點(diǎn)云分類，點(diǎn)云分類包括初分類和精分類2步[3]。初分類基于宏執(zhí)行自動分類。采用面向?qū)ο蠓诸惖姆椒ń⒎诸愐?guī)則[4]。本項(xiàng)目提出點(diǎn)云高程值“自上而下”的分類方式，根據(jù)點(diǎn)云數(shù)據(jù)的高度、分布的形狀、密度、坡度等特征編寫宏規(guī)則，進(jìn)行自動分類。

初始分類完成后進(jìn)行手工精分類，準(zhǔn)確分離地面點(diǎn)和非地面點(diǎn)。主要通過人機(jī)交互的方式，通過剖面顯示來進(jìn)行分類。精分類前后地面模型對比如圖1所示。

圖1 點(diǎn)云精分類前后模型渲染示意圖

1.1.3 DEM制作

將精分類以后的地面點(diǎn)在Terramodel中構(gòu)三角網(wǎng)生成數(shù)字高程模型，為了達(dá)到測區(qū)DEM無縫，用來生產(chǎn)數(shù)字高程模型的點(diǎn)云數(shù)據(jù)必須經(jīng)過嚴(yán)格接邊處理。

1.2 DOM制作

在TerraSoLid的Terraphoto中正射影像制作包括以下幾步：

1）影像內(nèi)定向。通過導(dǎo)入相機(jī)檢校文件包括相機(jī)主距、鏡頭畸變等參數(shù)，確定每張影像像主點(diǎn)的位置。

2）導(dǎo)入影像列表。該列表包含測區(qū)每張影像的攝影時刻和獲取的外方位元素。

3）相對定向。手工添加連接點(diǎn)在像對重疊處選取同名點(diǎn)，確定整個測區(qū)模型的各影像間的相對位置關(guān)系。

4）區(qū)域網(wǎng)平差。利用精分類抽稀后的地面點(diǎn)作為平高點(diǎn)[5]，將整個模型納入到地理坐標(biāo)中，解算連接點(diǎn)的絕對坐標(biāo)，利用解算的絕對坐標(biāo)更新影像列表中的外方位元素[6,7]。再將新的外方位元素應(yīng)用到工程中，通過平差計算調(diào)整連接點(diǎn)點(diǎn)位直至平差結(jié)果滿足同等比例尺生產(chǎn)精度要求。

5）正射糾正。消除了像片傾斜誤差和地形起伏的影像后，將中心投影改為正射投影，對影像進(jìn)行正射糾正。

6）勻色、鑲嵌及裁切。對測區(qū)的像片按照模板進(jìn)行勻光勻色，使得整個測區(qū)內(nèi)影像色調(diào)、亮度等信息保持一致，并對完成勻色后的影像在Terraphoto中進(jìn)行鑲嵌線編輯，避免自動鑲嵌后存在地物矛盾的現(xiàn)象。鑲嵌線編輯前后成果如圖2所示。

圖2 房屋鑲嵌線編輯前對比后示意圖

1.3 DLG采集

基于點(diǎn)云數(shù)據(jù)的DLG采集目前主要基于Microstation平臺或者CASS平臺，本項(xiàng)目基于CASS平臺完成DLG采集。主要生產(chǎn)步驟包括：

1）基于DEM成果的等高線生成及處理，利用DEM成果生成成圖比例尺要求的等高線，并將等高線進(jìn)行節(jié)點(diǎn)抽稀和圓滑，既保證等高線精度要求也使等高線的表示看上去比較美觀。

2）將精分類后的地面點(diǎn)云按照規(guī)則格網(wǎng)進(jìn)行抽稀，抽稀后的點(diǎn)云作為高程注記點(diǎn)。

3）對照DOM及點(diǎn)云數(shù)據(jù)進(jìn)行地形要素采集，采集過程中DOM主要用于判斷要素的形狀和性質(zhì)，而點(diǎn)云數(shù)據(jù)用于判斷點(diǎn)云的實(shí)際位置。在實(shí)際生產(chǎn)中可引進(jìn)地面浮雕模型，輔助判斷溝渠、坎、坡等不能直觀從DOM或者海量點(diǎn)云數(shù)據(jù)中發(fā)現(xiàn)的要素。地面浮雕模型如圖3所示。

4）DLG制作完成，通過外業(yè)補(bǔ)測、調(diào)繪后即可進(jìn)行成果整理入庫。

2 精度評定

以“某市1︰2 000信息化測圖”項(xiàng)目為例，采用外業(yè)散點(diǎn)法沿道路施測的平面檢查點(diǎn)和各高程檢查點(diǎn)檢核基于LiDAR數(shù)據(jù)生產(chǎn)DEM、DOM以及DLG的精度。

1）高程精度檢核：將外業(yè)實(shí)測點(diǎn)展于DEM上，對比實(shí)測點(diǎn)與DEM上同名點(diǎn)位高程值差值。外業(yè)實(shí)測949個實(shí)測高程檢查點(diǎn)，實(shí)測點(diǎn)類型、數(shù)量及檢核精度如表1所示。

表1 高程精度檢核表

本次檢核點(diǎn)共分為5個類別，其中野外實(shí)測高程點(diǎn)所占比例最高，為82.5%，且誤差分布服從正態(tài)分布。

1︰2 000測圖高程精度是平地小于0.4 m、丘陵地小于0.5 m、山地小于1.2 m及高山地小于1.5 m。本次所檢核的5類點(diǎn)高程中誤差無論地形如何，DEM精度都高于0.4 m，因此基于LiDAR數(shù)據(jù)生產(chǎn)的DEM產(chǎn)品精度較高，能夠滿足大比例尺測圖的精度要求。

2）將外業(yè)實(shí)測點(diǎn)與采集的DLG疊加，對比同名點(diǎn)點(diǎn)位的平面誤差。外業(yè)實(shí)測1 170個高程檢查點(diǎn)，實(shí)測點(diǎn)類型、數(shù)量及檢核精度如表2所示。

表2 平面精度檢核表

1︰2 000測圖平面精度要求平地、丘陵地位置中誤差應(yīng)小于等于1 m，接邊限差小于等于2個像素。經(jīng)檢核，基于LiDAR數(shù)據(jù)的DOM、DLG產(chǎn)品成果的高程誤差及平面誤差分布符合正態(tài)分布規(guī)律，且個別較大誤差值也在限差范圍內(nèi)。因此，基于LiDAR數(shù)據(jù)生產(chǎn)的3D產(chǎn)品完全能夠滿足1︰2 000比例尺成圖的精度要求。

3 結(jié) 語

通過介紹基于LiDAR數(shù)據(jù)的DEM、DOM、DLG產(chǎn)品制作的技術(shù)流程和方法，提出生產(chǎn)中數(shù)據(jù)處理的關(guān)鍵技術(shù)或應(yīng)注意的環(huán)節(jié)以供參考：

1）數(shù)據(jù)預(yù)處理時，及時對處理成果進(jìn)行過程控制檢查，避免后期出現(xiàn)返工情況。

2）為了測區(qū)DEM實(shí)現(xiàn)完全無縫鑲嵌，除了針對分幅的精分類點(diǎn)云數(shù)據(jù)進(jìn)行接邊外，還應(yīng)對分幅范圍進(jìn)行外擴(kuò)后再導(dǎo)出地面模型。

3）DOM生產(chǎn)中，選取的連接點(diǎn)點(diǎn)位分布及數(shù)量對平差精度有較大影響，選取原則盡量滿足分布影像四角和中心。

4）DLG采集時，與地面存在高差的地物、地貌均應(yīng)參考點(diǎn)云數(shù)據(jù)判斷其實(shí)際位置，如溝、斜坡、坎、高架橋梁等地物。

5）房屋采集是DLG采集的重點(diǎn)，應(yīng)注意：①在建筑密集區(qū)，矮建筑受遮擋，DOM上可能完全看不到低矮建筑的輪廓或者只能看見部分輪廓，此時不能簡單地依據(jù)DOM進(jìn)行繪制，應(yīng)依據(jù)點(diǎn)云的形狀判斷實(shí)際房屋輪廓。②從DOM上如果能直接看見建筑墻基，可以墻基進(jìn)行房屋糾正，但是應(yīng)結(jié)合點(diǎn)云判斷基準(zhǔn)邊。

6）等高線制作的數(shù)據(jù)源盡量使用DEM而不用精分類后的地面點(diǎn)點(diǎn)云，既可以保證地形的整體走勢，也減少了后期的編輯工作量。

[1]劉沛,李英成,薛艷麗,等.基于TerraSoLid與Inpho的LiDAR數(shù)據(jù)處理方法分析與研究[J].遙感應(yīng)用，2010(1)：17-21

[2]李英成，文沃根，王偉.快速獲取地面點(diǎn)位數(shù)據(jù)的LiDAR技術(shù)系統(tǒng)[J].測繪科學(xué)，2002,27(4):35-38

[3]楚長春，麻風(fēng)海.點(diǎn)位激光掃描儀點(diǎn)云數(shù)據(jù)在MicroStation下的處理研究[J].測繪與空間地理信息，2007(10)：114-116

[4]劉建國，張敬，高偉，等.LiDAR點(diǎn)云數(shù)據(jù)中建筑物的快速獲取[J].地球科學(xué):中國地質(zhì)大學(xué)學(xué)報，2006,31(5)：615-618

[5]王佩軍，徐亞明.攝影測量學(xué)[M].武漢：武漢大學(xué)出版社，2008

[6]鄧非.LiDAR數(shù)據(jù)與數(shù)字影像的配準(zhǔn)和地物提取研究[D].武漢：武漢大學(xué)，2006

[7]畢凱，李英成，薛艷麗，等.ERDAS LPS與TerraSoLid軟件相結(jié)合快速制作低空數(shù)碼遙感正射影像圖的方法[J].遙感信息，2008(5)：55-58