鄒 馨, 劉 健, 張 杰, 陳曉勇

(1.東華理工大學(xué)測(cè)繪工程學(xué)院，330013，南昌；2.中國科學(xué)院空天信息創(chuàng)新研究院，100094，北京)

0 引言

數(shù)字高程模型(Digital Elevation Model, DEM)作為測(cè)繪科學(xué)的數(shù)字化產(chǎn)品，廣泛應(yīng)用于各個(gè)領(lǐng)域。DEM實(shí)現(xiàn)了區(qū)域地形表面的數(shù)字化表達(dá)，是對(duì)二維地理空間上具有連續(xù)變化特征地理現(xiàn)象通過有限的地形高程數(shù)據(jù)實(shí)現(xiàn)對(duì)地形曲面的數(shù)字化模擬[1-2]。建立DEM的方法有很多種，除了傳統(tǒng)的航空攝影測(cè)量外，生成DEM的途徑還有通過野外測(cè)量或者從現(xiàn)有地形圖上采集高程點(diǎn)或者等高線內(nèi)插生成DEM；通過雷達(dá)技術(shù)獲取DEM等方法。

近年來隨著立體衛(wèi)星的發(fā)展，通過衛(wèi)星影像立體像對(duì)可以快速地獲取大范圍的DEM數(shù)據(jù)。牛云璞等基于IRS-P5立體影像數(shù)據(jù)來提取DEM，證明了基于IRS-P5衛(wèi)星影像提取DEM可以滿足一些水利工程規(guī)劃和設(shè)計(jì)階段的需要[3]。丁輝等利用IKONOS立體像對(duì)提取翠華山地區(qū)的DEM，采用高精度GPS儀野外采集控制點(diǎn)，對(duì)RPC模型進(jìn)行修正，結(jié)果表明DEM精度明顯提高，DEM和正射影像均滿足1:5 000比例尺地形圖和遙感平面制圖要求[4]。王秋玲等利用ZY-3立體像對(duì)，在無地面控制點(diǎn)情況下DEM精度在20 m以內(nèi)，在有控情況下滿足山區(qū)地形1:50 000 DEM測(cè)繪的精度要求[5]。陶蘭花等分別采用PCI GXL軟件和INPHO軟件對(duì)GF-7立體像對(duì)進(jìn)行數(shù)字表面模型生產(chǎn)實(shí)驗(yàn)，結(jié)果表明這2款軟件生產(chǎn)的產(chǎn)品精度均達(dá)到了平地1:10 000數(shù)字表面模型精度要求，但是并未討論在有控制點(diǎn)和無控制點(diǎn)情況下DEM的提取精度[6]。胡柳茹等采用高分7號(hào)衛(wèi)星獲取的多源遙感數(shù)據(jù)進(jìn)行平面和高程精度優(yōu)化，結(jié)果表明利用激光測(cè)高數(shù)據(jù)經(jīng)過中值模型和線性模型明顯提升DSM高程精度，利用足印影像可改善DOM的平面精度[7]。

高分衛(wèi)星作為遙感衛(wèi)星的代表，影像含有豐富大量有效的信息，其中高分7號(hào)衛(wèi)星是我國首顆民用亞米級(jí)高分辨率的立體測(cè)繪衛(wèi)星，搭載了前、后視2臺(tái)高分辨率遙感相機(jī)分別從前、后2個(gè)方向?qū)Φ孛嫱痪拔镞M(jìn)行不同角度的觀測(cè)[8]。本文以龍巖市為研究區(qū)域，基于高分7號(hào)立體像對(duì)進(jìn)行DEM的提取。

1 高分7號(hào)影像概況

高分7號(hào)衛(wèi)星(簡(jiǎn)稱GF-7)，是高分專項(xiàng)工程中的重要衛(wèi)星之一，是推動(dòng)空間地理信息化的主要力量，正在對(duì)中國測(cè)繪產(chǎn)業(yè)的發(fā)展產(chǎn)生深遠(yuǎn)影響[9]。高分7號(hào)影像，數(shù)據(jù)質(zhì)量不僅穩(wěn)定且可靠，且其立體測(cè)繪有利于減小大高差地區(qū)的遮擋，后視影像可以制作良好的正射影像。全色相機(jī)的前視分辨率優(yōu)于0.8 m，后視分辨率優(yōu)于0.65 m。多光譜相機(jī)有藍(lán)、綠、紅、近紅外4個(gè)波段，分辨率優(yōu)于為2.6 m[10]。該衛(wèi)星為自然資源的管理提供了良好服務(wù)，在高精度立體測(cè)繪、調(diào)查資源等方面提供了自然資源現(xiàn)狀數(shù)據(jù)、立體影像數(shù)據(jù)和輔助數(shù)據(jù)。并且在各個(gè)領(lǐng)域發(fā)揮著重要的作用，服務(wù)于1:10 000比例尺全國基礎(chǔ)地理信息測(cè)繪產(chǎn)品，進(jìn)一步提高了自主高分衛(wèi)星對(duì)自然資源主體業(yè)務(wù)的支撐服務(wù)能力。表1為高分7號(hào)衛(wèi)星的主要參數(shù)。

表1 高分7號(hào)衛(wèi)星參數(shù)

龍巖市地勢(shì)東高西低，北高南低，全市山區(qū)居多，同時(shí)有丘陵和平原，全市平均海拔652 m，具有一定的代表性。因此本文所使用的影像是福建省龍巖市的2021年2月高分7號(hào)立體像對(duì)，為1對(duì)2景GF-7立體全色影像，經(jīng)度范圍為116.916°～117.189°，緯度范圍為24.990°～25.259°，產(chǎn)品級(jí)別為1A級(jí)影像(經(jīng)過預(yù)處理的輻射校正產(chǎn)品)，云量為5%，太陽高度角為46°。2幅影像的覆蓋范圍如圖1所示，從圖1中可見立體像對(duì)的重疊度較高，前后視重疊度約為96%，大部分為山區(qū)，其中影像的左上部分由于有云存在遮擋，總體來說，影像質(zhì)量良好，特征點(diǎn)明顯可以進(jìn)行DEM提取實(shí)驗(yàn)。

圖1 GF-7衛(wèi)星立體像對(duì)

2 基于立體像對(duì)的DEM提取的原理與方法

立體像對(duì)是指按照人造立體效應(yīng)的原理，從攝影基線兩端攝取的具有重疊影像的一對(duì)像片，一般像片的重疊度大于60%，其目的是為了增強(qiáng)解譯。在攝影測(cè)量中，要構(gòu)成立體模型則必須從不同方位拍攝同一地物的2張像片[11]。

立體像對(duì)的DEM提取方法及過程：利用處于不同位置的兩個(gè)傳感器對(duì)地面上同一個(gè)點(diǎn)拍攝，形成一個(gè)夾角，當(dāng)確定了天空中傳感器的位置之后，夾角越小說明拍攝的地物點(diǎn)越低，反之，夾角越大地物就越高[12]。當(dāng)?shù)孛嫠悬c(diǎn)的高程值解算出后，就得到了數(shù)字地面模型。其原理如圖2所示。

3 DEM提取流程

立體像對(duì)提取DEM技術(shù)路線如圖3所示，主要包括讀取前視和后視影像，選擇是否定義控制點(diǎn)以及定義控制點(diǎn)的個(gè)數(shù)，在立體像對(duì)中選擇同名點(diǎn)的連接點(diǎn)，生成核線影像，提取DEM等步驟。

圖2 基于立體像對(duì)提取DEM原理示意圖

圖3 基于立體像對(duì)提取DEM技術(shù)路線

3.1 輸入立體像對(duì)

在ENVI軟件主界面下，讀取高分七號(hào)衛(wèi)星影像信息，如圖4所示。

圖4 前視影像(左)和后視影像(右)

3.2 定義地面控制點(diǎn)

本文的地面控制點(diǎn)均為實(shí)地野外測(cè)量，基于JSCORS網(wǎng)絡(luò)RTK技術(shù)施測(cè)，控制點(diǎn)的布設(shè)方式采取區(qū)域網(wǎng)布點(diǎn)法進(jìn)行布設(shè)，控制點(diǎn)的選擇和布設(shè)不僅和布設(shè)方案有關(guān)，還要考慮誤差改正等對(duì)控制點(diǎn)的具體點(diǎn)位要求，所以控制點(diǎn)點(diǎn)位選取應(yīng)遵循分布均勻、特征明顯、易于判別、交通便利、便于測(cè)量等原則。由于部分地方難以測(cè)量，所以控制點(diǎn)布設(shè)在平地以及交通便利的山上。

高分7號(hào)影像提供了RPC參數(shù)文件來解算傳感器模型，本文在不定義控制點(diǎn)、1個(gè)地面控制點(diǎn)、5個(gè)地面控制點(diǎn)這3種情況下提取DEM。在無控的情況下得到的DEM是相對(duì)高程，即以衛(wèi)星地勢(shì)面為基準(zhǔn)。在有控的情況下得到的DEM是絕對(duì)高程。本文選擇的1個(gè)控制點(diǎn)和5個(gè)控制點(diǎn)在影像上的分布如圖5所示。

圖5 1個(gè)控制點(diǎn)分布(左)和5個(gè)控制點(diǎn)分布(右)

3.3 連接點(diǎn)提取與匹配

通過連接點(diǎn)建立立體像對(duì)之間的關(guān)系，先進(jìn)行連接點(diǎn)的自動(dòng)提取，對(duì)連接點(diǎn)進(jìn)行編輯，再進(jìn)行連接點(diǎn)的匹配。本文中連接點(diǎn)的自動(dòng)提取采用基于區(qū)域灰度匹配法自動(dòng)尋找重疊區(qū)的連接點(diǎn)，連接點(diǎn)數(shù)目初始設(shè)置為100個(gè)。在立體像對(duì)(帶有粗略地理坐標(biāo))2個(gè)圖像上找到一個(gè)同名點(diǎn)，尋找地形特征匹配的小區(qū)域。通過手工交互編輯保證每個(gè)連接點(diǎn)匹配的準(zhǔn)確性，只有當(dāng)最大Y方向視差的值小于10時(shí)，才可以利用連接點(diǎn)計(jì)算生成核線圖。編輯完成后Y方向視差如圖6所示，連接點(diǎn)分布如圖7所示。

圖6 Y方向視差改正前后值

圖7 連接點(diǎn)選擇結(jié)果

3.4 核線影像的生成

生成核線影像，將二維的相關(guān)問題變成一維的相關(guān)問題。在水平像片上建立規(guī)則格網(wǎng)，它的行就是核線。把提取的同名點(diǎn)和輸入的控制點(diǎn)代入共線投影方程式：

(1)

求出相關(guān)參數(shù)，式中(xt，yt，zt)為物方，(x，y)為像方坐標(biāo)，對(duì)于水平相片yt為常數(shù)，利用求出的參數(shù)將影像轉(zhuǎn)換成核線影像，轉(zhuǎn)換后像素灰度值采用雙線性內(nèi)插的辦法重采樣。ENVI得到的前后視核線影像如圖8所示。

圖8 前視(左)和后視(右)核線影像

4 DEM提取結(jié)果及精度評(píng)價(jià)

4.1 DEM提取結(jié)果

設(shè)置提取DEM的相關(guān)參數(shù)，得到了在無控制點(diǎn)、1個(gè)控制點(diǎn)、5個(gè)控制點(diǎn)情況下生成的DEM，如圖9所示。

(a)無控 (b)1控 (c)5控

4.2 DEM精度評(píng)定

高程信息作為DEM中主要的價(jià)值體現(xiàn)，DEM精度非常重要，在DEM精度評(píng)價(jià)中，高程誤差的評(píng)價(jià)是最重要的方面。DEM的精度與地形因素、控制點(diǎn)的數(shù)量、精度、同名點(diǎn)的個(gè)數(shù)等因素有著很大的關(guān)系。因此，在DEM提取完成后，要對(duì)提取的DEM的精度進(jìn)行評(píng)價(jià)，即DEM精度評(píng)定。DEM精度評(píng)定主要有2種方式，分別為對(duì)DEM的平面精度和高程精度一起評(píng)定、對(duì)DEM的平面精度和高程精度分開評(píng)定。本文對(duì)DEM提取的高程精度進(jìn)行評(píng)價(jià)分析。

目前，評(píng)定DEM精度評(píng)定的主要方法有檢查點(diǎn)法、等高線法、剖面法以及影像分析法對(duì)DEM進(jìn)行精度評(píng)價(jià)[13]。由于實(shí)地測(cè)量了影像部分區(qū)域內(nèi)的高程值，故本文采用檢查點(diǎn)法進(jìn)行DEM精度評(píng)定。

檢查點(diǎn)法指將檢查點(diǎn)隨機(jī)分布在影像上，對(duì)生成的DEM在檢查點(diǎn)上取值進(jìn)行DEM提取值與DEM真值進(jìn)行比較，得到各個(gè)檢查點(diǎn)的提取高程值與實(shí)際高程值的誤差表。本文基于實(shí)地測(cè)量的控制點(diǎn)，抽取部分控制點(diǎn)進(jìn)行DEM高程精度評(píng)定。導(dǎo)入實(shí)地測(cè)量的控制點(diǎn)到提取的DEM中，通過ENVI中的Geographic Link得到控制點(diǎn)處的高程值，提取的DEM對(duì)應(yīng)的高程值、高程值與真實(shí)值之間的差值見表2。

表2 檢查點(diǎn)與DEM提取高程對(duì)比/m

根據(jù)表2中可見提取的DEM高程與真實(shí)高程，及在無控制點(diǎn)和有控制點(diǎn)的情況下提取高程值與真實(shí)值之間的差值。在無控制點(diǎn)情況下，GF-7立體像對(duì)生成的高程最大誤差為21.4 m，最小誤差為7.2 m；在1個(gè)控制點(diǎn)情況下，GF-7立體像對(duì)生成的高程最大誤差為5.2 m，最小誤差為0.018 m近似為0；在5個(gè)控制點(diǎn)情況下，GF-7立體像對(duì)生成的高程最大誤差為6.9 m，最小誤差為0.2 m。根據(jù)表2得到不同控制點(diǎn)個(gè)數(shù)下，檢查點(diǎn)真實(shí)高程與提取高程值圖(圖10)、檢查點(diǎn)高程誤差圖(圖11)。

圖10 GF-7 DEM檢查點(diǎn)真實(shí)高程與提取高程值

圖11 檢查點(diǎn)高程誤差

根據(jù)公式計(jì)算表中的中誤差。

(2)

式中：Hi為實(shí)地測(cè)量的DEM真值，Zi為基于高分7號(hào)影像提取DEM的高程值，n為控制點(diǎn)的個(gè)數(shù)。

從表2中可知一共有18個(gè)檢查點(diǎn)，根據(jù)均方誤差公式，分別計(jì)算無控制點(diǎn)、1個(gè)控制點(diǎn)、5個(gè)控制點(diǎn)情況下的高程中誤差。RMSE分別為14.6 m、2.4 m、2.6 m。在有控制點(diǎn)情況下，立體像對(duì)提取DEM的誤差滿足地形1:10 000 DEM的精度[14]。

5 結(jié)束語

本文利用GF-7立體像對(duì)福建龍巖地區(qū)提取DEM，在無控制點(diǎn)和有控制點(diǎn)的情況下，基于同一立體像對(duì)、相同的連接點(diǎn)進(jìn)行DEM的提取，并對(duì)其進(jìn)行精度評(píng)定。通過DEM精度評(píng)定發(fā)現(xiàn)，在無控制點(diǎn)的情況下，DEM的中誤差為14.6 m；在有控的情況下DEM高程精度比無控情況下精度明顯高，中誤差降到2.4 m；隨著控制點(diǎn)個(gè)數(shù)增加，高程精度變化不大。在有控的情況下DEM精度比無控情況下精度高。在有控制點(diǎn)的情況下，DEM高程精度可達(dá)到山地地區(qū)1:10 000的精度要求。

產(chǎn)生DEM誤差的原因可能有以下幾點(diǎn)：1)檢查點(diǎn)的選取與分布對(duì)DEM評(píng)定有一定的影響。在山體等高程起伏較大的地方，連接點(diǎn)布設(shè)應(yīng)該按照山體的起伏情況進(jìn)行相應(yīng)的改變；2)選擇一個(gè)控制點(diǎn)進(jìn)行DEM提取時(shí)，部分檢查點(diǎn)的高程值不減反增，這是因?yàn)槠鸱^大，該控制點(diǎn)在平地，導(dǎo)致提取的DEM值與真實(shí)值相差較大；當(dāng)在高山處和平地處增加控制點(diǎn)后，DEM提取效果比1個(gè)控制點(diǎn)時(shí)提取效果好，誤差基本上在4 m以內(nèi)。