何仲太

(1.北京大學(xué) 地球與空間科學(xué)學(xué)院，北京 100871；2.中國地震局地殼應(yīng)力研究所 地殼動力學(xué)重點實驗室，北京 100085)

1 引 言

高分辨率遙感影像的快速發(fā)展極大地促進(jìn)了利用高分辨率立體像對提取DEM的應(yīng)用。市場上目前有立體測圖能力的衛(wèi)星主要有SPOT-5、IKONOS、IRS-P5、WorldView-I/Ⅱ、ALOS、QuickBird、GeoEye和國產(chǎn)立體測圖衛(wèi)星ZY-3、TH-1等?；诟叻直媛市l(wèi)星數(shù)據(jù)進(jìn)行DEM提取的研究目前國內(nèi)外已有很多，但基于IRS-P5立體像對提取DEM的研究應(yīng)用目前還較少。IRS-P5是印度于2005年5月5日發(fā)射的一顆高分辨率制圖衛(wèi)星。衛(wèi)星的軌道高度為618km，搭載有2個分辨率約為2.5m的全色傳感器，沿軌道方向連續(xù)推掃，形成同軌立體像對[1]。IRS-P5數(shù)據(jù)具有成本低、覆蓋廣泛、質(zhì)量穩(wěn)定、分辨率適中的特點，基于IRS-P5立體像對提取DEM的研究具有實際應(yīng)用價值和意義[2]。近年來數(shù)字?jǐn)z影測量技術(shù)取得了巨大進(jìn)展，目前常用來從IRS-P5立體像對提取DEM的軟件主要有我國自主研發(fā)的全數(shù)字?jǐn)z影工作站(VirtuoZo)、數(shù)字?jǐn)z影測量工作站(系統(tǒng))(JX-4 DPS)、一體化測圖系統(tǒng)PixelGrid以及遙感軟件PCI(GeomatieaorthoEngine模塊)、ERDAS(LPS模塊)、ENVI(DEM Extraction模塊)。基于數(shù)字?jǐn)z影測量軟件從IRS-P5立體像對提取DEM步驟相對比較復(fù)雜，但生成的DEM具有更好的質(zhì)量。近年來DEM和DOM數(shù)據(jù)在活斷層定量研究中得到了廣泛應(yīng)用，不同的活斷層定量參數(shù)研究對DEM和DOM數(shù)據(jù)精度要求不同。為了通過大青山山前斷裂上升盤基巖山區(qū)中沖溝的裂點分布特征來獲取斷裂晚第四紀(jì)古地震序列，因此需要獲得能夠識別高差約為3m～7m沖溝裂點的DEM數(shù)據(jù)。本文選取內(nèi)蒙古大青山山前斷裂地區(qū)的IRS-P5立體像對為例，采用在有理函數(shù)模型的幾何定位基礎(chǔ)上，利用PixelGrid軟件生成DEM、DOM數(shù)據(jù)，并進(jìn)行精度評價。

2 研究區(qū)與數(shù)據(jù)

大青山山前斷裂位于內(nèi)蒙古境內(nèi)，西起包頭市，東至呼和浩特，全長約200km，總體走向NEE，傾向S，第四紀(jì)以來活動非常強烈[3]。斷裂晚更新世以來的活動在空間上的分布是不均勻的，前人根據(jù)古地震時空分布以及斷裂的幾何特點、活動時間、活動速率及構(gòu)造地貌等方面的實際資料將其劃5段[4-5]。斷裂兩側(cè)地表保留了大量斷裂晚更新世和全新世強烈活動的信息，本文選取斷裂及兩側(cè)范圍的IRS-P5像對數(shù)據(jù)進(jìn)行研究，為斷裂微地貌定量研究提供基礎(chǔ)數(shù)據(jù)。本文所選用的數(shù)據(jù)是覆蓋了整個研究區(qū)的5個IRS-P5立體像對，數(shù)據(jù)均為標(biāo)準(zhǔn)3級產(chǎn)品，且均沒有云、雪覆蓋，地物的幾何結(jié)構(gòu)和紋理信息清晰(圖1)?？刂瀑Y料主要選取已有的1∶10000紙質(zhì)地形圖資料。依據(jù)在地形圖中讀取的控制點坐標(biāo)及高程，利用PixelGrid軟件對P5立體像對提取5m格網(wǎng)間距的DEM，并生成2.5m分辨率的正射影像，然后對正射影像進(jìn)行調(diào)色、鑲嵌等，完成DOM的制作。

圖1 研究區(qū)IRS-P5衛(wèi)星立體像對數(shù)據(jù)示意圖(條帶區(qū)域為研究區(qū)域)

3 生成DEM和DOM

IRS-P5向用戶提供了有理函數(shù)模型(RFM)，一般情況下數(shù)據(jù)供應(yīng)商提供的有理函數(shù)模型參數(shù)沒有地面控制點的參與，所能達(dá)到的精度受制于星歷數(shù)據(jù)與姿態(tài)參數(shù)的精度。大量實驗證明，在無地面控制點的情況下IRS-P5遙感影像數(shù)據(jù)平面和高程的定位精度約為70m[6]，因此在實際應(yīng)用中需要使用地面控制點來提高影像數(shù)據(jù)精度。按照圖2所示流程進(jìn)行DEM提取，首先將掃描后的地形圖進(jìn)行模型糾正，并在糾正后的地形圖上選取合理有效的控制點。控制點應(yīng)立體、均勻地分布在影像區(qū)內(nèi)，尤其應(yīng)加強測區(qū)邊沿的控制點布設(shè)，并且控制點盡量避免布設(shè)在近似一條直線或近似一個平面上。將選取的控制點以Excel表格的形式記錄其點號坐標(biāo)和點位描述等內(nèi)容，形成一個PixelGrid規(guī)定的*.gcp格式的文件，對照刺點片在影像上量測控制點。

圖2 DEM提取流程圖

量測控制點后進(jìn)行平差定向，先進(jìn)行“RPC+平移”的粗差剔除，后進(jìn)行“RPC+二維仿射變換”的最終解算。剔除粗差是用“RPC+平移”方法進(jìn)行平差，誤差大于20的點視為粗差，將該點刪除，在刪除點的區(qū)域要重新選擇控制點并進(jìn)行量測。重復(fù)平差、刪點，直到剔除全部粗差。最終解算使平差結(jié)果滿足誤差要求，即定向誤差小于3個像素，每個點的x、y、z誤差均不大于2倍中誤差。平差定向結(jié)果滿足要求后生成近似核線，并進(jìn)行影像匹配與編輯。通過匹配種子點線量測模塊，檢查模型是否存在視差，如果出現(xiàn)視差需要修改點位、添加連接點并重新定向。初始影像匹配僅匹配出稀疏的點，可在匹配種子點線量測模塊中進(jìn)行編輯，在山脊、山谷、影像質(zhì)量不好的地方添加特征線能提高“高精度匹配”的精度。根據(jù)初始匹配及編輯的結(jié)果，匹配大量高精度的點，并生成后綴名為dtm的文件。選擇模型及對應(yīng)的dtm文件進(jìn)行立體編輯，對匹配不好的地方在立體編輯下進(jìn)行修改。每個模型的點被分成兩行多列的矩形塊，編輯時從左到右、從上到下逐塊進(jìn)行。將編輯之后的成果轉(zhuǎn)換為tif格式，加載正確的投影和坐標(biāo)系，按照作業(yè)范圍對成果進(jìn)行拼接和裁切，形成最終的DEM成果(圖3(a))。利用PixelGrid軟件將P5數(shù)據(jù)的后視tif文件在影像格式轉(zhuǎn)換及預(yù)處理模塊進(jìn)行格式轉(zhuǎn)換，再將原文件包中與文件名對應(yīng)的.raw.pxy，.raw.rpc，.raw.spt 3個文件拷貝到新建文件包中。使用拼接一起后的DEM數(shù)據(jù)在正射影像生成模塊最終生成DOM(圖3(b))。

圖3 DEM(a)和DOM(b)成果圖

4 精度評價

實際應(yīng)用中常用的DEM精度評價方法主要有：檢查點法、剖面線法及等高線套合法。其中檢查點法最為簡單易行，是在研究區(qū)范圍內(nèi)選取檢查點并獲得這些點的實際高程值，同時在生成的DEM上找到相應(yīng)點的高程值，計算二者之間的差值，檢查其是否滿足國家關(guān)于DEM精度標(biāo)準(zhǔn)的相應(yīng)規(guī)定。目前研究普遍認(rèn)為在良好的地面控制點的支持下，對于IRS-P5立體像對提取的DEM的平面和高程精度均可達(dá)到5m，滿足國家1∶50000基本比例尺地形圖測圖規(guī)范要求。還有一些學(xué)者試驗研究認(rèn)為，IRS-P5立體像對可以完成常規(guī)航空攝影測量無法進(jìn)行的困難區(qū)域山地1∶10000地形圖立體測圖，其平面、高程精度能達(dá)到山地和高山地測圖中誤差要求[1，7-9]。關(guān)于地面控制點選取對于精度的影響認(rèn)識有所不同，有的學(xué)者研究認(rèn)為在至少采用4個地面控制點的情況下，平面和高程方向影像定位精度能夠達(dá)到2m的精度，6～8個控制點即可達(dá)到最優(yōu)精度，平面和高程精度在1.2m～1.7m之間[10]。還有學(xué)者通過實驗分析得出每個立體像對只需要5 個地面控制點即可達(dá)到最好的平面和高程精度，過多的地面控制點對于精度的提高沒有實質(zhì)影響[11]。

4.1 DEM精度評價

DEM的精度評價在內(nèi)容上可通過兩種不同的方式進(jìn)行，一種是平面精度和高程精度分開評價，另一種是兩種精度同時評價。在實際應(yīng)用中由于平面精度比較容易控制，一般只討論DEM的高程精度評價問題并用高程精度表達(dá)DEM精度。本次針對DEM的精度檢測，也是只針對DEM的高程精度的檢測。

采用攝影測量方法和野外實測方法生成的數(shù)字高程模型，其格網(wǎng)點高程中誤差應(yīng)不大于相應(yīng)比例尺地形圖測圖規(guī)范或編繪規(guī)范中規(guī)定的等高線高程中誤差。本文對于DEM的精度檢測，采用的是實際生產(chǎn)中最常用的檢查點法，在研究區(qū)均勻選取18個檢查點，事先將檢查點按格網(wǎng)或任意形式進(jìn)行分布，對生成的DEM在這些點處進(jìn)行檢查。將這些點處的內(nèi)插高程和實際高程逐一比較得到的誤差，然后計算中誤差。設(shè)檢查點的高程為Zk(k=1，2，3，…，n)，在建立DEM之后，由DEM內(nèi)插出的這些點的高程為Rk，則DEM的精度即高程中誤差為：

(1)

表1 DEM精度檢查結(jié)果

結(jié)果(表1)顯示檢查點的高程精度均可達(dá)到5m，滿足國家1∶50000基本比例尺地形圖測圖規(guī)范要求。中誤差的計算結(jié)果也顯示IRS-P5立體像對生成的DEM高程精度能達(dá)到1∶10000地形圖山地和高山地測圖中誤差要求。

本文在DEM生成過程中誤差主要可歸結(jié)為像點誤差與定向誤差。像點誤差主要表現(xiàn)在進(jìn)行絕對定向前各種系統(tǒng)誤差與偶然誤差的綜合誤差，受制于衛(wèi)星影像本身質(zhì)量。定向誤差包括相對定向和絕對定向的誤差，相對定向誤差主要取決于同名像點量測的誤差，其定向精度可根據(jù)平差結(jié)果進(jìn)行估計。絕對定向誤差主要與采用的模型、控制點的誤差、分布等有關(guān)。其中控制點應(yīng)立體、均勻分布，避免在近似一條直線的似一個平面布設(shè)。另外，在建立DEM的過程中，由于需要經(jīng)過內(nèi)插計算和建模處理，內(nèi)插點的計算高程總是與實際量測高程之間存在差異，這種誤差一方面與內(nèi)插算法有關(guān)，另一方面也與原始數(shù)據(jù)的分布和密度有關(guān)。

4.2 DOM精度評價

將控制資料上的控制點的坐標(biāo)值與在生成的DOM上量取的坐標(biāo)值相比較，按下式統(tǒng)計中誤差。結(jié)果(表2)顯示18個檢查點的水平向中誤差均小于理論中誤差限差7.5m，說明DOM平面精度符合1∶10000地形圖精度要求。利用DEM實現(xiàn)影像的微分糾正是生產(chǎn)DOM的關(guān)鍵技術(shù)，因此DOM的成圖精度與所用DEM的精度具有很大的關(guān)系。另外本試驗中DOM的精度還與控制點精度密切相關(guān)，因此應(yīng)盡量避免前期地形圖掃描和模型糾正過程誤差的產(chǎn)生，并合理有效地選取控制點。

(2)

(3)

(4)

其中，Xi、Yi為GPS控制點坐標(biāo)值；xi、yi為糾正后DOM上量取的坐標(biāo)值；Mx、My分別為檢查點在X、Y方向的點位中誤差；n為檢測點個數(shù)；M為檢測點點位中誤差。

表2 DOM精度檢查結(jié)果

5 結(jié)束語

本文以內(nèi)蒙古大青山地區(qū)的IRS-P5全色立體像對為例，在1∶10000地形圖上選取控制點，在PixelGrid軟件平臺上進(jìn)行試驗。結(jié)果顯示試驗區(qū)IRS-P5全色立體像對可以生成5m格網(wǎng)間距的DEM和2.5m分辨率的DOM。以1∶10000地形圖為參照，用檢查點法對DEM、DOM的精度進(jìn)行評價。結(jié)果顯示IRS-P5立體像對生成的DEM滿足國家1∶50000基本比例尺地形圖測圖規(guī)范要求，高程精度能達(dá)到1∶10000地形圖山地和高山地測圖中誤差要求，DOM平面精度符合1∶10000地形圖精度要求。本文研究顯示，利用IRS-P5立體像對生成的DEM數(shù)據(jù)滿足活斷層定量研究中對于微地貌沖溝裂點的識別和提取，可應(yīng)用于活斷層定量研究中的古地震研究。另外，IRS-P5全色立體像對生成高分辨率DEM和DOM相對傳統(tǒng)數(shù)據(jù)源具有覆蓋廣泛、成本低、時效高、更新快速的優(yōu)勢，在地形圖的數(shù)據(jù)獲取和更新、災(zāi)后快速評估以及其他定量地貌研究中具有廣泛的應(yīng)用前景。

致謝：感謝北京東方道邇信息技術(shù)股份有限公司。

參考文獻(xiàn)：

[1] 趙巍，楊靜，王海燕，等.基于IRS-P5衛(wèi)星立體影像的地形圖測繪試驗[J].測繪標(biāo)準(zhǔn)化，2008，24(2)：34-38.

[2] 黨軍宏，孫美萍.基于IRS-P5衛(wèi)星影像的DEM制作與方法研究[J].山西建筑，2012，38(3)：218-219.

[3] 國家地震局鄂爾多斯周緣活動斷裂系課題組.鄂爾多斯周緣活動斷裂系[M].北京：地震出版社，1988.

[4] 冉勇康，張培震，陳立春.河套斷陷帶大青山山前斷裂晚第四紀(jì)古地震完整性研究[J].地學(xué)前緣，2003，10：207-216.

[5] 何仲太，馬保起，盧海峰.大青山山前活動斷裂帶分段與潛在震源區(qū)劃分[J].地震地質(zhì)，2007，29(4)：765-775.

[6] 趙利平，劉鳳德，李健，等.印度測圖衛(wèi)星IRS-P5定位精度初步研究[J].遙感信息，2007，(2)：28-32.

[7] 王國昌.衛(wèi)星影像制圖在鐵路勘測設(shè)計中的應(yīng)用研究[J].鐵道工程學(xué)報，2010(2)：14-17，31.

[8] 邸國輝，程崇木，劉傳生，等.基于IRS-P5衛(wèi)星影像測圖的關(guān)鍵技術(shù)研究[J].地理空間信息，2012，10(3)：27-30，53.

[9] 范蘭，趙巍，陳杰，等.P5衛(wèi)星影像立體測圖的試驗和分析[J].測繪通報，2013(2)：54-58.

[10] 范興旺，翁永玲，胡伍生，等.IRS-P5立體像對提取DEM及精度評價[J].遙感技術(shù)與應(yīng)用，2010，25(4)：547-551.

[11] 黃華平.IRS-P5衛(wèi)星影像立體制圖精度的研究[J].測繪，2012，35(5)：235-237.