山西省高精度數字高程模型建設方案設計與實現

楊庚印

(山西省測繪工程院，山西太原 030002)

山西省高精度數字高程模型建設方案設計與實現

楊庚印

(山西省測繪工程院，山西太原 030002)

結合山西省高精度數字高程模型建設項目，詳細介紹利用機載激光雷達進行大面積數據獲取的技術方案和數據處理關鍵技術，為進行大面積機載激光雷達數據獲取和處理積累了豐富的經驗。該項目的建成，滿足了省內各類不同行業用戶以及科學研究機構對本省精密地表數據的需求。

機載激光雷達；高精度；數字高程模型；方案設計

一、引 言

隨著時代的發展，建立一個高精度、高時空分辨率、高效率、高覆蓋率的全省數字高程模型數據庫，成為國土測繪、城市規劃、環境監測、防災減災、交通建設、電力建設、水利建設、地面沉降測量、工程測量以及城市地理信息系統等領域的迫切需求。為及時滿足山西省內各類不同行業用戶以及科學研究機構對本省精密地表數據的需求，實施了山西省高精度數字高程模型建設工程。

數據獲取方法選擇。傳統數字高程模型建立的方法有多種，從數據源及采集方式講有：①直接從地面測量，如用GPS、全站儀、野外測量等；②根據航空或航天影像，通過攝影測量途徑獲取，如立體坐標儀觀測及空三加密法、解析測圖、數字攝影測量等；③從現有地形圖上采集，如格網讀點法、數字化儀手扶跟蹤及掃描儀半自動采集，然后通過內插生成DEM等方法。

最新方法：①合成孔徑雷達干涉測量(InSAR)，它是使用雷達信號的相位信息提取地球表面三維信息，分為機載和星載兩種；②機載激光雷達(Li-DAR)，它是一種安裝在飛機上的機載激光探測和測距系統。

根據以上所述，傳統方法已不適合建立高精度的數字高程模型的要求。在選擇新技術時，考慮到機載激光雷達是采用主動對地觀測技術，并具有全天候、作業周期短、精度高、費用低等優點，并具有攜帶中畫幅可量測相機，同時獲得了影像數據，故在項目建設中選擇了機載激光雷達。

二、激光雷達系統測量原理及系統簡介

1.激光雷達系統測量原理

激光掃描測量是通過激光掃描器和距離傳感器來獲取被測目標的表面形態。激光掃描器一般由激光發射器、接收器、時間計數器、微電腦等組成。激光脈沖發射器周期地驅動一激光二極管發射激光脈沖，然后由接收透鏡接收目標表面后向反射信號，產生一接收信號，利用一穩定的石英時鐘對發射與接收時間差做計數，經由微電腦對測量資料進行內部微處理，顯示或存儲輸出距離和角度資料，并與距離傳感器獲取的數據相匹配，最后經過相應系統軟件進行一系列處理，獲取目標表面三維坐標數據，從而進行各種量算或建立立體模型。

2.系統簡介

根據市場現有機載激光雷達設備情況和項目需求，采用了Leica公司的ALS60機載激光雷達測量系統。

(1)系統主要組成

系統包含激光掃描儀、激光控制器、系統控制器、相機系統、POS系統(GPS＋IMU)、飛行控制系統。

(2)主要技術參數

1)飛行相對高度200～6000 m范圍。

2)最大脈沖頻率200 000 Hz。

3)最大掃描頻率100 Hz。

4)回波/強度信號次數4/3。

5)最大掃描角75°。

6)自適應滾動補償最大75°。

7)激光發散角0.22mr＠1/e2(0.15mr＠1/e)。

8)IPAS型號CUS6。

(3)測量精度

表1為該系統在掃描角為40°時的平面和高程測量精度。

表1

三、方案設計

1.航線設計

1)航測基準面高程要略低于航測分區的平均高程。

2)對于大面積、高航高測量，當地面地形起伏的相對高差較小，即一般情況下小于1000 m時，可采用MPiA技術(空中內插多脈沖)；否則應采用SPiA技術(單脈沖)。

3)如果設備飛行中采用了MPiA技術，應該務必保證相對航高的范圍在設備回波有效范圍之內，否則會出現山頭或谷底沒有回波的情況。

4)對于特殊高差較大地區，可以采用分區飛行。

5)航線航向的布設按攝區的地形來定，盡量減少同一條航線中地形起伏變化太大，如果地面起伏變化大于航高的1/3時，應分區進行設計。

6)航線設計時向外擴的安全距離為相對航高× 0.7。

7)每條航線的設計長度要小于飛機25 min的航程，設計時要考慮飛機的類型和順風與逆風的速度。

8)航線重合度：旁向15%以上。

9)激光設備的技術參數文件在設計時需要唯一確定，并在AeroPlan中指定。

10)編制航線示意圖(?.jpg)，編制飛行員GPS領航表(?.xls)。

11)提供航線飛行數據(?.fpb)。

12)航線設計結束之后，要在Google Earth中檢查設計的結果。

2.激光雷達檢校場選擇

激光掃描儀校驗場應盡量選擇太原市、長治市、大同市、晉中市等具有1∶500或1∶2000比例尺地形圖的地區，其條件包括以下3個方面：

1)校準控制場。選擇類似機場跑道或公路、長至少2000 m、寬20～50 m的較大平坦區域。

2)校準建筑物。校準區內必須有一個較大“人”字形尖頂建筑物。建議建筑物長度大于100 m、寬度大于80 m、高度大于5 m。用于校準側滾和俯仰姿態。

3)盡量遠離水面(如湖、江)等低反射率的地區。

3.激光雷達檢校場布點

布設激光檢校控制點要求：

1)直線控制點：沿著道路中心，每隔5 m一個控制點，長度大約2 km。

2)零散控制點：在中心區域均勻布設10～15個控制點。

3)所有控制點都布設在路面上，且地物材料均勻。避免高低反射率交接地區，避免周圍地物遮擋，避免在陡坎和地物過度邊界、便道邊緣布設。

4.飛行要求

1)飛行前，嚴格按照飛行檢查單的要求進行飛行前檢查，確保設備安裝和各項設置正確無誤。

2)飛行轉彎角不大于22°。

3)飛行時，飛行Roll角值小于5°，Pitch角值小于6°。

4)在入第一條航線前和出最后一條航線后，都必須進行IMU初始化飛行，即：進入第一條航線前先平飛5 min，然后作“8”字轉彎飛行，再平飛2 min進入航線；出最后一條航線后，直飛2 min，再作“8”字轉彎飛行。

5)在設備出現故障時，重新啟動設備，在進線之前需要進行“8”字重新飛行。

6)如果下一個測區距現有測區超過25 min的直飛距離，需要在進下一個測區之前飛“8”字。

7)飛機上同步觀測數據在進入首航線的IMU初始化前5 min至出最后航線IMU初始化后5 min的時間段全程記錄。

8)在規定的航測期限內，選擇地表植被及其他覆蓋物(如積雪)對LiDAR測量影響較小、云霧少、無揚沙的季節進行攝影。

9)要根據GPS衛星預報，選擇PDOP值小于3的時間段進行航攝；若3＜PDOP＜6的時間段在20min之內時，該時間段也容許進行航攝，但這種時間段一天的航攝過程中不能多于3個。

10)在每個架次飛行時，應保證相鄰兩條航線有一次同向飛行或一條航線往返飛行。

5.地面基站

1)基準站采用山西省連續運行基準站網(SXCORS)，為滿足測量精度要求，最好選擇雙基站，基站距離測區的邊緣不應超過50 km。

2)測繪基準：國家CGCS2000。

3)觀測要求：

a.與機載GPS同步觀測。

b.GPS衛星高度角≥15°。

c.有效觀測衛星個數≥5。

d.衛星觀測值象限分布(25±10)%。

e.GPS接收機采樣間隔為0.5″。

f.提供基站觀測數據(?.Dat、?.T01)。

6.檢查點布點方案及要求

依據全省1∶10 000圖幅分幅范圍，根據地形類別，山地、高山地四幅圖布設一個平面點和一個高程點，平地、丘陵地每幅圖布設一個平面點和一個高程點，在全省呈均勻分布形態。

檢查點平面點點位應選擇在明顯地物的角點上，并盡量在平坦地區；高程點點位必須選擇在面積不小于10 m2的平坦地面上。

所有點位的選擇必需滿足RTK觀測要求。

四、數據處理

1.工作流程

機載激光雷達數據處理工作流程如圖1所示。

圖1 機載激光雷達數據處理工作流程

2.數據釋放及處理

(1)激光點云數據定位定向

對原始IPAS數據進行處理，得到精確的導航文件。一般分為3步：

1)解壓原始IPAS數據，得到GPS數據和IMU數據。

2)進行GPS差分解算，提高GPS的定位精度。

3)融合GPS和IMU數據，進行內插處理，得到精確的導航文件。

POS系統是機載激光探測和測距系統的關鍵，其核心是采用動態差分技術和慣性導航測量裝置，在航攝飛行過程中測量傳感器的位置和姿態，經嚴格的聯合數據處理，獲得高精度的傳感器外方位元素，從而實現沒有或極少地面控制的傳感器定位和定向。

在該項目中GPS動態差分采用SXCORS系統站點進行解算。

(2)生成標準激光數據LAS文件

利用原始激光點云文件、導航文件和激光檢校文件，在ALS Post Processing軟件中生成標準激光點LAS文件。此文件可以在任何激光數據處理軟件中打開。如TerraSolid、Mars View等軟件。

3.激光點云數據的檢校

在飛行過程中，IMU和激光掃描儀的相對姿態可能會發生微小的變化，從而對激光數據產生影響，為消除這種影響，通常要對定位定向后的激光數據進行檢校。

檢校分為視準軸檢校(boresight calibration)和距離檢校(range calibration)，通常需要檢校的參數為偏心分量：側滾角(roll)，俯仰角(pitch)和航偏角(heading)。roll視準軸偏差定義了IMU和激光發射器之間在X軸方向上的偏差，單位是弧度。任何掃描鏡編碼器的誤差也包含在roll誤差中。pitch視準軸偏差定義了IMU和激光發射器之間在Y軸方向上的偏差，單位是弧度。heading視準軸偏差定義了IMU和激光發射器之間在Z軸方向上的偏差，單位是弧度。range糾正可以補償激光掃描儀中電子器件延遲所產生的誤差。這個誤差一般是一個常數。由于大量數據同時運行速度較慢，實際生產中，為較快地得到較好的檢校參數，首先在檢校場數據中選擇一塊典型地形進行自動檢校，讓系統自動生成全區最優的檢校參數。如果在實際測區發現有細微的差別，這時可以進行手動微調，即可滿足對測區激光數據的檢校。

針對rroll、pitch、heading模型檢校結果應滿足表2要求。

表2m

4.數據分類及成果制作

在機載LiDAR預處理數據基礎上進行手動精細分類濾波，將LiDAR點云數據地表點(含植被、水系、構筑物等自然人工地物表面點)精確分類出來，對手動精細分類后的地表點云數據進行最終DSM、DEM產品的構建與輸出。點云分類要求如下。

(1)點云數據自動濾波、分類

利用TerraSolid公司提供的專用軟件，在Terra-Scan和TerraModel模塊中對激光條帶數據進行分塊處理，并利用編寫的宏命令逐塊進行濾波分類，區分出地面點、植被、水域等類別。

(2)點云精細分類

對自動分類的成果參照獲取的影像數據進行手動精細分類濾波，以獲得高精度的數字高程模型。結合項目要求，將LiDAR點云數據分為water、ground、noise、default、temp、lowpoint六類。

(3)點云分類的原則

1)地貌、土質以及與地面相連而成的道路、水系、高臺等地物應表示。

2)建筑物、高塔等人工地物(如架空跨越的橋梁)、植被、地物應與地面點進行區分。

3)臨時地物(如臨時土堆等靜地物，車輛、行人、飛鳥等動地物)、粗差點等應剔除。

4)地物尺寸小于產品格網間距的可作忽略處理。

5)地物高程差異小于相應比例尺DEM中誤差的可做忽略處理。

6)對于制作DEM需要濾除的地物應放到相應的類別，不能刪除。

7)手動分類前對于靜止水域如池塘、湖泊應進行置平處理。

8)河流水涯線上的高程應從上游到下游逐漸降低，水涯線同一平面位置水涯線高程值置平。

(4)水域置平處理

在點云分類前應先進行水域置平處理，步驟如下：

1)在Microstation軟件中畫水域邊界范圍線，統一歸層。

2)將范圍線內ground類點云分到water類中。

3)將水域所有范圍線畫完后，參考范圍線內的已有water層點云，統一修改閉合多邊形內符合實際情況的指定線高度。

4)關閉所有點云層，選中全部所畫水域范圍線，在Terramodel中進行Triangulate Survey(壓平、置平)，然后使用Insert Breakline Element工具參考閉合多邊形內插特征點。

5)將生成的水域邊界點導出生成XYZ文本格式文件，在Terrascan中讀入，將該文件和原始點云文件一起保存，并將線上特征點分到water類。

(5)點云分類方法

1)對高程突變的區域，調整參數或算法，重新進行小面積的自動分類。

2)對分類錯誤的點重新進行分類。分類時，以點云切剖面為主要依據，DOM僅作參考；在比較平滑、直線區域切剖面時，剖面寬度可適當放大，在拐角尤其是立交橋、高架公路等接地與架空的臨界區域，切剖面一定要盡量窄，務求精確；溝渠寬度5 m以內的可不作水域點的精細分類；寬度大于5 m和面狀水域在制作DEM時均需要置平。生成的水域邊界breakline特征點在讀入后分類為temp。

(6)DSM成果制作

利用TerraSolid公司提供的專用軟件，在Terra-Scan模塊下，利用分類點云數據建立數字表面模型后，再輸出符合設計要求的數據格式。

(7)DEM成果制作

對經過自動濾波、手動精細分類后的地面點云數據作為特征點以及帶有高程信息的特征線進行數字高程模型的構建。

五、精度檢測

通過對布設的1865個野外檢查點檢查，高程方向上的差值基本呈正態分布，h方向最大較差為＋0.49 m；高程中誤差為±0.18 m。整體精度優于項目設計中規定的精度指標。

六、結束語

本文結合山西省高精度數字高程模型建設項目，詳細論述了高精度數字模型數據獲取、數據處理的關鍵技術及技術方案，為進行大面積機載激光雷達數據獲取積累了豐富的經驗。同時，該項目的完成，將用來全面更新山西省數字高程模型數據庫，為全省的基礎測繪、應急保障、災害分析、研究機構及各行各業應用提供了精確地表數據，將在全省的轉型跨越發展中發揮重要作用。

[1] 張小紅.機載激光雷達測量技術理論與方法[M].武漢：武漢大學出版社，2007.

[2] 鞏淑楠，陳云.徐敏，等.機載激光雷達數據處理方法的研究與應用[J].測繪與空間地理信息，2010，33(5)：165-168.

[3] 鐘良，支曉棟，劉鵬飛，等.Leica機載激光雷達應用研究[J].測繪通報，2009(6)：72-73.

[4] 李清泉，李必軍.陳靜，等.激光雷達測量技術及其應用研究[J].武漢測繪科技大學學報，2000，25(5)：387-391.

[5] 賴旭東.機載激光雷達基礎原理與應用[M].北京：電子工業出版社，2010.

[6] 徐祖艦，王滋政，陽鋒.機載激光雷達測量技術及工程應用實踐[M].武漢：武漢大學出版社，2009.

[7] 周淑芳，李增元，范文義，等.基于機載激光雷達數據的DEM獲取及應用[J].遙感技術與應用，2007，22(3)：356-360.

[8] 喻雄.機載激光雷達在山區高速公路勘測中的應用[J].測繪通報，2011(2)：31-34.

[9] 徐景中，萬幼川，張圣望.LiDAR地面點云的簡化方法研究[J].測繪信息與工程，2008(1)：32-34.

[10] 吳華意，宋愛紅，李新科.機載激光雷達系統的應用與數據后處理技術[J].測繪與空間地理信息，2006(3)：58-63.

Design and Implementation of Shanxi High-precision Digital Elevation Models

YANG Gengyin

P237

B

0494-0911(2014)09-0076-05

2013-07-17

楊庚印(1974—)，男，山西萬榮人，高級工程師，主要研究方向為數字航空攝影及數據處理。

楊庚印.山西省高精度數字高程模型建設方案設計與實現[J].測繪通報，2014(9)：76-80.

10.13474/j.cnki.11-2246.2014.0299