何仲太

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

1 引 言

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

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

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

圖1 研究區(qū)IRS-P5衛(wèi)星立體像對(duì)數(shù)據(jù)示意圖(條帶區(qū)域?yàn)檠芯繀^(qū)域)

3 生成DEM和DOM

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

圖2 DEM提取流程圖

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

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

4 精度評(píng)價(jià)

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

4.1 DEM精度評(píng)價(jià)

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

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

(1)

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

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

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

4.2 DOM精度評(píng)價(jià)

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

(2)

(3)

(4)

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

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

5 結(jié)束語(yǔ)

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

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

參考文獻(xiàn)：

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

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

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

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

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

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

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

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

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

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

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