ENVI軟件下ASTER遙感立體像對的DEM提取

1、ENVI軟件下ASTER遙感立體像對的DEM提取 BASED ON THE REMOTE SENSING IMAGE CUBE ASTER TODEM EXTRACTION學院（部）： 專業(yè)班級： 學生姓名： 指導教師： 年 月 日基于ASTER遙感立體像對的DEM提取 摘要文中通過實例研究了ASTER遙感立體像對在高山峽谷地區(qū)提取DEM的精度。首先簡述ASTER的立體像對提取DEM的國內(nèi)外發(fā)展現(xiàn)狀及應用，然后介紹了具有同軌立體測圖能力的ASTER傳感器觀測系統(tǒng)，及其立體影像生成DEM 的算法和DEM 的編輯方法，展示了ASTER 立體影像生產(chǎn)DEM 的實驗結(jié)果，并以試驗結(jié)果說明，其精度可以滿

2、足繪制1100 000 1250 000 比例尺地形圖要求。ASTER 遙感數(shù)據(jù)可以在地球表面自然過程與全球變化等多個研究領(lǐng)域得到廣泛應用。本文是通過使用ENVI4.7對ASTER遙感立體像對影像進行數(shù)據(jù)讀入、輸入控制點、提取連接點、指定參數(shù)、DEM編輯等處理后得到DEM三維立體圖像。關(guān)鍵字： ASTER，DEM，遙感，立體像對，ENVIBased on the remote sensing image cube ASTER to DEM extractionAbstractIn this paper, through the case, study the remote sensing im

3、age cube ASTER in alpine canyon area extraction accuracy of the DEM. First, this paper expounds the ASTER stereo image and abroad for the extraction of the DEM situation and application, then introduces observation system in the same rail stereo mapping capability, ASTER sensor and stereo images and

4、 the algorithm generated DEM editorial method ASTER stereo images, showing the production of experimental results and DEM to test results show that the accuracy can satisfy drawing 1:100 000 1:250 000 scale topographic map requirements. ASTER remote sensing data in the earths surface can be natural

5、process and global change research etc widely.Through the use of ENVI4.7 ,load in, input control points, extracting connection, designated parameters, DEM editing ,in order to obtain DEM.Key word: ASTER，DEM, Remote Sensing，Stereo Image，ENVI目錄摘要（中文） 1摘要（外文） 2第一章 前言41.1 ASTER遙感立體像對的DEM提取的發(fā)展介紹41.2 ASTE

6、R遙感立體像對的DEM提取的在現(xiàn)代的應用51.3 論文研究的內(nèi)容及意義 7第二章 基于ASTER遙感立體像對的DEM提取的原理和方法 72.1 ASTER遙感立體像對的DEM提取的原理 72.2 ASTER遙感立體像對的DEM提取方法及過程8第三章ASTER遙感立體像對數(shù)據(jù)的DEM提取的處理93.1 對ASTER遙感立體像對數(shù)據(jù)的基本處理原理及方法93.1.1 連接點的提取103.1.2 核線影像的生成103.1.3 DEM 的生成103.1.4對生成DEM的處理113.3 ASTER遙感立體像對生成DEM的流程123.3.1 輸入立體影像對123.3.2 定義地面控制點123.3.3 定義聯(lián)

7、接點133.3.4 指定參數(shù)153.3.5生成DEM和三維立體圖 17第四章實驗及結(jié)果分析 18參考文獻20致謝22第一章 前言1.1ASTER遙感立體像對的DEM提取的發(fā)展介紹數(shù)字地面高程模型(Digital Elevation Models ,簡稱DEM)是地表形態(tài)的數(shù)字形式,它由規(guī)則水平間隔處地面點的抽樣高程矩陣組成,在生產(chǎn)中具有很高的利用價值。DEM 數(shù)據(jù)具有廣泛的應用潛力。作為國家空間數(shù)據(jù)基礎(chǔ)設(shè)施(NS2DI) 的框架數(shù)據(jù)之一,DEM 生產(chǎn)技術(shù)已經(jīng)比較成熟,如等高線地形圖生產(chǎn)DEM 技術(shù)、利用航攝像對生產(chǎn)DEM 的數(shù)字化攝影測量方法等?，F(xiàn)在,隨著遙感技術(shù)的不斷發(fā)展,利用遙感影像對提取

8、DEM 數(shù)據(jù)已經(jīng)成為可能,如利用SPOT數(shù)據(jù)生產(chǎn)的DEM 精度可達到10m 以內(nèi)，而利用星載IN2SAR 生成DEM 的高程精度可達到米級,并且平坦地區(qū)的精度優(yōu)于山區(qū)的精度。利用衛(wèi)星數(shù)據(jù)提取DEM 的技術(shù)已經(jīng)成為遙感現(xiàn)今發(fā)展的一個研究課題。ASTER(Advanced Space borne Thermal Emission and Re2flection Radiometer ,高級星載熱輻射反射輻射計) 是美國航空航天局(NASA) 與日本國際經(jīng)貿(mào)商業(yè)部(METI) 合作發(fā)布的高分辨率衛(wèi)星成像設(shè)備,于1999 年12 月搭載NASA 的EOS-AM1 ( Terra) 平臺升空。其目標是獲

9、取地球表面溫度、輻射、反射和高程數(shù)據(jù),研究生物圈、水圈、巖石圈和大氣層之間的互動反應,解決土地利用與覆蓋、自然災害(火山噴發(fā)、水災、森林火災、地震和風暴) 、短期天氣變動、水文等方面的問題。ASTER 在可見光和近紅外區(qū)(VNIR) 有三個波段,在短波紅外區(qū)(SWIR) 有六個波段,在熱紅外區(qū)( TIR) 有五個波段,它們的地面分辨率分別是15m、30m、90m。ASTER 觀測系統(tǒng)的VNIR 子系統(tǒng)在可見光與近紅外區(qū)有三個波段, 空間分辨率為15 m。第三波段有兩個通道3N，3B，其中3B 具有后視(27.6)的能力，可以同軌立體觀測, 其沿徑向的基高比B/H 為0.6，每景大小為60 km

10、 60 km。兩個成像儀可以對相同地物以合理的視差和航向重疊近同步（大約相差55 s）成像。對于這兩幅影像，儀器與地物之間有相同的氣象條件，同時地面的光照和輻射條件也是幾乎相同的，這就極大地提高了影像匹配的成功率。Terra 衛(wèi)星的周期是16 d，VNIR 還具有側(cè)視( 24)能力, 這樣最快5 d。ASTER傳感器發(fā)射的目的在測圖領(lǐng)域就是獲取南緯82和北緯82 之間全球80%區(qū)域的云層覆蓋率小于20%的影像，生產(chǎn)DEM。與其他遙感衛(wèi)星傳感器相比, ASTER有其獨特之處在于：較高空間分辨率和立體觀測能力。VNIR的子系統(tǒng)有一個后視波段(即第三波段,波長范圍0. 78m0. 86m) ,沿軌跡

11、方向(along track) 進行立體觀測,即沿飛行軌道前后掃描立體成像,這樣大大提高了立體像對的接收成功率,也保證了像對間的相關(guān)性和相似性。由此產(chǎn)生的立體像對是提取DEM 的基礎(chǔ),它對于地形解譯具有重要的意義。利用遙感數(shù)據(jù)提取DEM 的精度依賴于基高比率(Base to Height Ratio ,BH) 、空間分辨率和視差測量的誤差。VNIR基高比率是0. 6 ,天底交角為27. 6。如果假定0. 51. 0 像素作為總的視差誤差,則ASTER 的DEM 高度精度應該在550m ,從ASTER 立體數(shù)據(jù)中產(chǎn)生的地圖比例尺將在150 0001250 000 之間。由于ASTER 圖像數(shù)據(jù)本

12、身不像SPOT 等衛(wèi)星數(shù)據(jù)帶有飛行參數(shù)，其內(nèi)部、外部定位要通過精確的、數(shù)量足夠的和平均分布的地面控制點來確定。所以,控制點的精度、分布和數(shù)量是決定提取DEM 數(shù)據(jù)的質(zhì)量的關(guān)鍵。ASTER 是極地軌道環(huán)境遙感衛(wèi)星Terra ( EOS- AM1) 上載有的5 種對地觀測儀器之一,它提供了可見光近紅外(VNIR) 、短波紅外(SWIR) 和熱紅外( TIR) 3 個通道的遙感數(shù)據(jù)。目前,盡管ASTER數(shù)據(jù)可以通過網(wǎng)上免費預訂和下載,但讀取和分析這些數(shù)據(jù),仍然需要依賴于專業(yè)遙感圖像處理軟件,而購買這類軟件，一方面需要大量資金,另一方面,利用這些商業(yè)軟件進行遙感參數(shù)反演,也存在諸多不便之處, 如EOS

13、VIEW 是唯一個專門查看ASTER 數(shù)據(jù)格式(HDF - EOS) 的專業(yè)軟件，但它不能進行數(shù)據(jù)處理，而其它專業(yè)軟件(如ENVI) 雖然可以查看HDF - EOS 格式的遙感數(shù)據(jù)，但遙感反演與處理功能有待進一步完善。因此，利用軟件編程實現(xiàn)ASTER 數(shù)據(jù)讀取，是ASTER 數(shù)據(jù)應用的前提，也是開發(fā)具有國產(chǎn)自主版權(quán)的對地觀測數(shù)據(jù)綜合處理系統(tǒng)的一項基礎(chǔ)工作。為此， 我們課題組對ASTER 數(shù)據(jù)格式進行了全面分析與解讀， 并對ASTER 數(shù)據(jù)讀取關(guān)鍵技術(shù)進行了討論與研究，實現(xiàn)了ASTER 數(shù)據(jù)的讀取及其在定量遙感中的應用。衛(wèi)星遙感立體像對提取DEM是地貌信息獲取的一個重要里程碑，ASTER衛(wèi)星傳感

14、器是可以拍攝立體像對傳感器中的代表，具有數(shù)據(jù)質(zhì)量穩(wěn)定、覆蓋廣泛、價格低廉的特點。本文通過實例研究了ASTER立體像對在高山峽谷地區(qū)提取DEM的精度。首先簡述ASTER的立體像對提取DEM的國內(nèi)外發(fā)展現(xiàn)狀，然后針對一處高程變化顯著地區(qū)在1：10萬比例尺地形圖采集地面控制點（GCP），用1：5萬精度的DEM作檢驗，獲得GCP范圍內(nèi)高程誤差為20.4m，GCP范圍外高程誤差為48.2m，平均誤差是34.3m。這證明可以在小區(qū)域內(nèi)選取GCP控制點，由ASTER立體像大范圍外推生成大范圍DEM，而且采用常規(guī)的技術(shù)手段和普通的商業(yè)軟件就可實現(xiàn)。該方法提取DEM對于我國地形資料缺乏的西部地區(qū)有很強的實用性。

15、1.2 ASTER遙感立體像對的DEM提取的在現(xiàn)代的應用ASTER 數(shù)據(jù)已在多個研究領(lǐng)域得到廣泛應用。下面結(jié)合TM 數(shù)據(jù)，給出ASTER 數(shù)據(jù)在幾個主要研究領(lǐng)域的應用。1) 遙感關(guān)鍵參數(shù)反演。與TM 數(shù)據(jù)相比，ASTER 數(shù)據(jù)在熱紅外有更多的波段，波譜分辨率更高，利用ASTER 數(shù)據(jù)這一特點，可以計算全球與區(qū)域尺度上的土壤水分蒸騰與蒸發(fā)量和植被蒸騰與蒸發(fā)量，將ASTER 數(shù)據(jù)與地面風速、空氣溫度和空氣濕度觀測數(shù)據(jù)相結(jié)合，把“地- 氣”作為一個整體，對低層大氣溫度、濕度垂直廓線和陸面溫度、比輻射率進行參數(shù)同步反演，計算“地- 氣”的水汽通量與熱通量，這對于評價植被受旱狀況，進行作物估產(chǎn)，以及推進

16、全球水分平衡和碳循環(huán)研究具有重要意義。2) 建立全球地面高程模型。TM 數(shù)據(jù)只提供底視圖像，而ASTER 數(shù)據(jù)在可見光波段既提供了底視圖像，也提供了后視圖像。利用同一區(qū)域同一時間獲取的底視和后視圖像，可以在基于包含偽GCP(地面控制點) 條件下，生成一個粗略的地面高程模型。由于該模型采用偽GCP ，它們僅僅具有地面x 、y 值， 而沒有絕對高程信息， 因此， 生成的是相對DEM。為了獲取該地區(qū)絕對地面高程，可以結(jié)合偽GCP ，使用世界USGS GTOPO30 1km DEM ，估計它們的絕對高程值。由于精度和數(shù)據(jù)質(zhì)量等原因，其結(jié)果仍然存在高程誤差，但這比采用相對高程值要準確和方便。3) 土地覆

17、蓋調(diào)查與監(jiān)測。和TM 數(shù)據(jù)相比，ASTER 數(shù)據(jù)在可見光和近紅外波段有更高的空間分辨率，利用可見光和近紅外波段的ASTER 數(shù)據(jù)，可以進行區(qū)域性或全球范圍的土地覆蓋調(diào)查，對不同土地類型進行分類與評估，利用ASTER 圖像數(shù)據(jù)編制土地類型、森林分布、草地等級圖等，并利用同一區(qū)域不同時相的遙感數(shù)據(jù)了解與監(jiān)測土地覆蓋變化，為生態(tài)環(huán)境保護提供參考依據(jù)。4) 冰川覆蓋調(diào)查與變化監(jiān)測。冰川主要分布在南極、北極及高海拔的山地，這些地域環(huán)境惡劣，交通不便，利用ASTER 數(shù)據(jù)具有幾乎覆蓋全球的高分辨率圖像進行冰川研究，快捷方便。若把衛(wèi)星定位和衛(wèi)星遙感技術(shù)相結(jié)合，對比不同時期冰川的衛(wèi)星影像，可以從宏觀上進一步了

18、解現(xiàn)代冰川的近期變化，認識冰川變化與海平面升降關(guān)系，繪制冰川分布圖。目前，國際上一個研究組織正在開展一項利用ASTER 數(shù)據(jù)進行全球陸地冰川變化監(jiān)測的研究。5) 火山監(jiān)測與評估。在過去400 多年里，幾乎每年都有火山噴發(fā)，因此，利用ASTER 提供的熱紅外圖像預報和監(jiān)測火山，具有重要意義。例如，利用ASTER 熱紅外圖像，可以對火山口及其周圍熱異?，F(xiàn)象進行研究，如觀測火山附近的溫泉、熱氣口及火山口湖的溫度在噴發(fā)前的變化; 在火山噴發(fā)時， 通過ASTER 圖像可以發(fā)現(xiàn)火山灰和熔巖流的分布區(qū)域，并繪制熔巖流分布圖。此外，也可以利用ASTER 熱紅外圖像進行熱紅外輻射場監(jiān)測，及時了解地震前出現(xiàn)的熱異

19、?，F(xiàn)象。在強震的孕育和發(fā)生過程中，沿地震構(gòu)造區(qū)帶產(chǎn)生熱輻射是巖石物理實驗和地震現(xiàn)場調(diào)查都證實了的重要物理前兆現(xiàn)象，近年來，中國地震局衛(wèi)星遙感熱紅外地震前兆異常監(jiān)測和預報的部分實驗表明，許多中強地震發(fā)生前，震中區(qū)附近會出現(xiàn)短期(13月) 的熱紅外遙感長波輻射場變異，如1997 年發(fā)生在青藏高原北部的瑪尼7. 5 級地震、1998 年河北張北的6. 2 級地震等，均出現(xiàn)了空間上可識別的衛(wèi)星熱紅外區(qū)域增溫現(xiàn)象。上述情況表明， ASTER熱紅外遙感圖像可作為地震預報的重要參考依據(jù)。1.3論文研究的內(nèi)容及意義 此次論文主要內(nèi)容是使用ENVI軟件對ASTER衛(wèi)星遙感立體影像數(shù)據(jù)進行處理， 提取DEM是地貌

20、信息利用雷達信號的相位信息提取地球表面的高精度三維信息，主要應用于測量地面點的高程及其動態(tài)變化，經(jīng)過對ASTER衛(wèi)星遙感立體影像進行干高精度配準、去平地效應、相位噪聲濾除、相位解纏的處理過生成DEM三維立體圖。 ASTER星載衛(wèi)星數(shù)據(jù)中心可以生成全球的DEM，這似乎解決了所有DEM信息獲取上的問題，但在DEM的應用上仍存在著很大的問題： 1）ASTER的幅寬只有60km和有計劃的全球信息獲取任務只是拍攝總量的8，這制約了東西向流域完整的無云覆蓋的影像（DEM）的獲??；2）生產(chǎn)DEM數(shù)量有限，產(chǎn)量每天動態(tài)變化，一般僅是拍攝影像的數(shù)量的10； 3）生成的DEM產(chǎn)品精度有一定損失，ASTER的影像分

21、辨率為15m，而DEM標準產(chǎn)品只有30m。 4）只能整景處理，對地面控制的點數(shù)量和分布要求高。 5）ASTER已經(jīng)超過了七年的設(shè)計壽命，有隨時停止運轉(zhuǎn)的可能?？晒┨娲牧Ⅲw像對商業(yè)衛(wèi)星遙感平臺還有SPOT、IKNONOS、Quickbird和Cartosat等。 所以如果掌握了立體像對獲取DEM技術(shù)可以，可以不變應百變，方便、靈活、及時的獲取DEM。 此次論文主要是對光學雷達圖像ASTER數(shù)據(jù)進行處理，得到三維地理高程信息，生成三維立體DEM圖像，通過對ASTER星載衛(wèi)星數(shù)據(jù)中心進行處理可以生成全球的DEM，可以得到三維立體模型圖，可用于獲取高精度的地形高程數(shù)據(jù)及監(jiān)測地殼形變等等。本文的主要內(nèi)

22、容就是基于ASTER遙感立體像對提取DEM，重點探討提取DEM的原理、方法，并評定DEM的提取精度。第二章 基于ASTER遙感立體像對的DEM提取的原理和方法 2.1 ASTER遙感立體像對的DEM提取的原理從不同方位對同一地區(qū)攝取兩幅影像構(gòu)成為立體影像對，其實質(zhì)是傳感器模擬人的雙眼，構(gòu)成人體立體視覺，即將像對上的視差反應為人眼的生理視差后得到的立體視覺，一般立體像對的重疊度至少在60%以上。構(gòu)建人造立體視覺需具備三個條件：1）由兩個不同位置（一條基線的兩端）拍攝同一景物的兩張相片稱為立體像對；2）兩只眼睛分別觀察像對中的一張像片；3）觀察時像對上同名像點的連線要同人的眼睛基線大致平行，而且同

23、名點間的距離一般小于眼基線。立體像對獲取DEM的原理簡易闡述為：在天空兩點（相當于人的兩只眼睛）拍攝地面同一點時形成一個角，當天空兩點的空間位置確定后，該角度越大地物點越高；反之，角度越小地物點越低。以此類推，將地面所有的高程解算后就得到了數(shù)字地面模型，如下圖：S1、S2為兩個攝影基站，相當于人的兩只眼睛，A、B、C分別為三個物點，a1、b1、c1和a2、b2、c2分別為三個物點A、B、C經(jīng)S1、S2后形成的像點。2.2 ASTER遙感立體像對的DEM提取方法及過程立體像對獲取DEM的原理簡易闡述是：在天空兩點(P1和P2)拍攝地面同一點A時形成P1AP2夾角，當P1和P2空間位置確定后，該角

24、度越大地物點越高，反之，角度越小地物越低。將地面所有點的高程解算后就得到了數(shù)字地面模型。像對DEM獲取方法有主要兩種：一種是通過GCP和影像上對應的像素來計算衛(wèi)星的外方位元素(衛(wèi)星姿態(tài))，通過像對間的匹配點(TP，Tie Point)配準影像，然后進行后方交匯計算每個地面點的高程；另一種是外方位元素的從衛(wèi)星的星歷中解算，其它步驟同l，此方法要求具有精確的衛(wèi)星星歷，而且生成的DEM為相對高程。本文中采用第一種方法。GCP可以通過三種途經(jīng)獲?。翰捎肎Ps、全站僅等地面測繪儀器獲??；在部分地形圖上采集控制點；在影像上采集x、Y平面坐標，在與影像坐標匹配精度較高的DEM上采集相應坐標的高程z值。第3種

25、方式方便快捷，本文從l：10萬的地形上獲取GCP的三維坐標。本文數(shù)據(jù)處理采用了一種基于ENVI的全自動ASTER DEM提取軟件，其無需用戶指定任何參數(shù)即可自動運行ASTER立體像對的DEM提取。以使用ASTER立體像對進行DEM提取。首先從ASTER HDF文件中提取垂視(3N)和后視(3B)圖像，然后定義輸出DEM的投影方式和分辨率、控制點的投影方式。從后視圖像明顯比垂視圖像存在拉伸畸變效應。立體像對的影像必須配準到相同的地面區(qū)域，這一步是通過在兩幅圖像中選擇同名點的TP點(僅包含在圖像中的相對行列值)，TP的選取要求在整個工作區(qū)中均勻分布。因為后視圖像成像時與垂直成像存在很大的畸變，所以

26、點位的選取尤其在圖像的水平方向即與衛(wèi)星飛行路徑垂直的方向要求密集和均勻一些，這樣可以盡可能地糾正圖像的畸變。根據(jù)研究區(qū)域復雜程度和所取子區(qū)的大小不同，約束點的數(shù)量要求也不同，但是為了后續(xù)的多項式擬合需要，一般最少需要9個點。為了使得到的DEM具有地理參考信息，還要在圖像中選取已知大地坐標和高程信息的GCP來對衛(wèi)星的外方位元素進行結(jié)算，及DEM進行地學校正，一般要求至少選擇4個點。將影像中部偏左地區(qū)的1：10萬地形圖作為GCP的數(shù)據(jù)源，共獲取了15個GCP，此外根據(jù)相對地物明顯程度匹配了25個TP。當圖像配準至相同的地面區(qū)域以后，任何沿軌道方向的位置差異就認為是由于視差引起的，而視差來源于地形起

27、伏。這樣，每個像元的視差可以通過使用三角關(guān)系和衛(wèi)星軌道數(shù)據(jù)(星歷、軌道位置、高度和姿態(tài)等)轉(zhuǎn)換為相對的(無GCP)或絕對(有GCP)高程信息。用ASTER數(shù)據(jù)提取DEM流程圖像準備(導入HDF格式)后視(3B)垂視(3N)投影定義(圖像與控制點)約束點、控制點選取擬合、空間關(guān)系解算DEM提取第三章 ASTER遙感立體像對數(shù)據(jù)的DEM提取的處理3.1 ASTER遙感立體像對數(shù)據(jù)的基本處理立體像對生成DEM 的主要包括連接點的提取，控制點的輸入，核線影像的生成和DEM 的生成。3.1.1 連接點的提取要生成DEM 首先要選擇立體像對上的一些連接點。連接點也是同名點，它是用來建立兩張像片之間關(guān)系的。

28、連接點的提取一般先自動提取，再手工交互編輯。連接點的自動提取采用基于灰度的影像相關(guān)的辦法。本文采用的是相關(guān)系數(shù)法。在3N 影像上取一個的邊長為奇數(shù)的窗口，取一個同樣大小的窗口作模板在3B 影像上以一定的步長移動窗口，求取兩個窗口的相關(guān)系數(shù)，取相關(guān)系數(shù)最大值對應的窗口的中心像素為3N 影像上窗口中心像素的對應點，以一定的步長在3N 影像上移動窗口求取各中心點的對應點，以相關(guān)系數(shù)最大的前若干個點作為對應的同名點。同名點的數(shù)目不能少于9 個，25 個以上最好。因為不需要找出每一個像素的同名點，窗口可以取得大一些，移動的步長也可以是若干個像素。自動提取以后再進行人工編輯，剔出錯誤的連接點，如果點數(shù)太少

29、，人工的選取一些連接點，保證連接點分布均勻。提取完連接點，再讀入控制點3.1.2 核線影像的生成如果直接用影像相關(guān)的算法，求出各個像素的同名點，再計算視差，計算量非常大，一般是先生成核線影像，把二維的相關(guān)問題變成一維的相關(guān)問題。核線在像片上是互不平行的，他們交于核點，如果將像片上的核線投影到一對相對水平的像片-平行于基線的像片對，則核線互相平行。根據(jù)這一原理，在水平像片上建立規(guī)則格網(wǎng)，它的行就是核線。把前面提取的同名點和輸入的控制點代入共線投影方程式： 求出相關(guān)參數(shù)，式中(xt， yt， zt) 為物方，(x， y) 為像方坐標，對于水平相片yt 為常數(shù)，利用求出的參數(shù)將3N，3B 影像轉(zhuǎn)換成

30、核線影像，轉(zhuǎn)換后像素灰度值采用雙線性內(nèi)插的辦法重采樣。3.1.3 DEM 的生成對于生成的核線影像，它已經(jīng)不存在上下視差，只有左右視差，利用圖像相關(guān)算法對3N 和3B 影像進行匹配，這時尋找同名點就是一維相關(guān)。根據(jù)所對應地區(qū)的最大高程值，計算一個最大視差，在3B 影像上的搜索范圍是3N 影像上像素坐標加減這個視差值，這樣可以縮小搜索范圍，減小了計算量。找到兩張像片的重疊區(qū)域的同名點后就可以利用輸入的控制點坐標進行整體平差，解求出地面點的坐標，根據(jù)各點高程的不同賦予相應的灰度值，就生成DEM 的圖像。3.1.4對生成DEM的處理在DEM 提取過程中，一些像素的匹配可能會失敗（由于大量的陰影，對比

31、度較低等原因），圖像由于傳感器自身或外界一些因素的影響會帶有許多噪聲，圖像不清晰，這時就要對生成的DEM 進行編輯。DEM的編輯包括DEM 的內(nèi)插，DEM 的濾波，DEM 的平滑。1）DEM 的內(nèi)插。對于初步生成的DEM，圖像本身因為云的遮擋，或是地面的對比度較低，陰影等因素，一些像素的匹配不正確，像素的值不正常，對于這些不正常區(qū)域采用分塊雙三次多項式（2 式）內(nèi)插的辦法，選擇不正常的區(qū)域，對不正常的像素插值，重復幾次，直到滿意為止。上式有16 個未知參數(shù)，除了已知四個格網(wǎng)數(shù)據(jù)點外，還必須知道每個格網(wǎng)點沿方向的斜率，和該點曲面扭曲。內(nèi)插只能解決小區(qū)域的異常，對于云層遮擋較大的區(qū)域內(nèi)插不一定符合

32、實際的地形，一般將這部分數(shù)據(jù)裁剪，再利用另一個時段生成的DEM 進行拼接。2）噪聲的濾波。在ASTER影像生成的DEM中噪聲主要包括兩類：斑點噪聲和隨機噪聲。對于這兩類不同的噪聲，如果分別用不同的方法進行兩次濾波，容易造成信息的損失，本文采用改進的中值濾波算法去噪。這種算法就是先進行噪聲的判斷，再進行噪聲的去除。具體的算法就是選取或的窗口作模板，在DEM 圖像上移動模板，求取窗口中像素值的平均值和標準差，如果窗口中心的像素值與平均值的差值的絕對值超過2倍的標準差，就認為這是噪聲，就用窗口中像素的中值代替窗口中心的像素值，不斷移動窗口，完成對整個DEM的濾波。這種方法可以有效地消除兩種噪聲。3）

33、DEM 的平滑。經(jīng)過內(nèi)插和去噪的DEM 還會有一些噪聲，有的地方仍存在不正常的點，但是地面的起伏一般是連續(xù)的，這時就需要對DEM 進行平滑。平滑采用的窗口為模板，在DEM 圖像上移動，以平均值代替中心像素值，不斷移動窗口，完成對整個DEM 的平滑。3.2 ASTER遙感立體像對生成DEM的流程這次使用ENVI4.7軟件對ASTER遙感立體像對數(shù)據(jù)生成DEM的流程，其過程如下：3.2.1輸入立體影像對從ENVI的主菜單欄，選擇File Open External File EOS ASTER然后從ENVI的主菜單欄，選擇Topographic DEM Extraction DEM Extract

34、ion Wizard New.選擇立體影像對： 3.2.2定義地面控制點 第二步的步驟，為你提供選項，以定義或使用現(xiàn)有的GCPs。當定義GCPs，你可以選擇任何下列選項: 3.2.3定義聯(lián)接點 立體影像之間的關(guān)系必須定義通過選擇或產(chǎn)生節(jié)點。這節(jié)點是用于定義的點，創(chuàng)造極線幾何的圖像，并經(jīng)常使用DEM提取。當定義聯(lián)接點，可以選擇任何下列選項:點擊下一步生成的垂視(3N)和后視(3B)圖像： 3.2.4指定參數(shù):1)投影DEM的參數(shù)指定輸出:2) DEM提取參數(shù)指定:3.2.5生成DEM和三維立體圖 1）生成DEM2)立體三維DEM效果圖：第四章 實驗及結(jié)果分析 不同的地形地區(qū)，采用不同的控制點數(shù)量

35、，使用多種的立體像對，考慮不同處理方法和軟件的差異，來進行DEM提取和的精度評定的研究。如ASTER數(shù)據(jù)在的平坦農(nóng)業(yè)區(qū)和地表裸露地區(qū)垂直精度10m，丘陵地帶為1520m23， 山區(qū)的精度2030m之間2，4， 有森林覆蓋的坡地可以達到3050m2。在理想條件下Spot5立體像對生成的DEM精度可以達到平面1516m，高程67m5。地面控制點（GCP，Ground Control Point）在立體像對提取DEM的過程中起到關(guān)鍵的作用，對于EOS衛(wèi)星ASTER立體像對的直接對地定位模型，利用國內(nèi)外3個地區(qū)的立體像對進行的實驗結(jié)果表明，無控制對地絕對定位的結(jié)果具有明顯的系統(tǒng)誤差。只需加入一個控制點

36、，則平面和高程精度迅速提高到1Om以內(nèi)。采用一個控制點外推580km時，平面定位精度仍優(yōu)于2Om，高程定位精度約10m，增加控制點的數(shù)量對改善定位精度的效果并不明顯。總體上三維控制點一般從幾個到幾十個不等，控制點越多生成的DEM精度越高，但一般不應少于5個。由上面生成的DEM圖（3.2.5）可見其質(zhì)量與GCP分布的具有很強的相關(guān)性，GCP所在范同及其附近的地區(qū)生成的DEM完整光滑，距GCP越遠DEM質(zhì)量越差，并且出現(xiàn)無值的“孔洞”，同時受高差變化的影響，高差變化顯著的東北地區(qū)未能生成有效高程的面積較大，而距GCP區(qū)距離相同的東南部地區(qū)，地形平坦生成的DEM質(zhì)量較高。這是因為距GCP遠，兩幅影像

37、不能精確配準，高差大后視成像存在著較大的疊掩和陰影效應。在生成的有效DEM地區(qū)，隨機選取20個檢驗點，其中10個在GCP范圍內(nèi)，10個均勻分布在GCP范圍外。采用1：5萬比例尺地形圖生成的DEM相應點作參考。在GCP范圍內(nèi)高程誤差為204 m，GCP范圍外高程誤差為482 m，平均誤差是343 m，超過4000 m以上4個點的誤差達到69 m。這從定量上也說明，在一定的GCP數(shù)量下，GCP的分布和研究區(qū)的地形是影響的DEM質(zhì)量的主要因素。為進一步比較，在DEM的線段位置沿水平方向做地形剖面與相應位置1：5萬DEM的剖面進行比較。兩種地形剖面的總體形態(tài)非常一致，但地形圖DEM過渡更平滑，ASTE

38、RDEM地面細節(jié)信息更豐富而且明顯，如，江河谷所在位置，地形圖DEM剖面上是平底“U”型，而ASTERDEM是尖底的“V”型。這岡為地形圖DEM是通過等高線生成而來，受等高線高程間距的限制，因此過渡平滑，細小且變化劇烈的地貌信息不能充分表達，而ASTER立體相對以15 m的平面分辨率來構(gòu)建DEM，自然能更充分細致的描述地貌的變化。另外，增加TP的數(shù)量可以在有限范圍內(nèi)改善DEM的質(zhì)量，但改善效果并不明顯。本文選取某一山地的ASTER L1A數(shù)據(jù)進行DEM的測圖試驗。試驗中，利用3N和3B 波段的數(shù)據(jù)構(gòu)成的立體像對生成DEM，以生成的DEM 糾正ASTER 1，3N，2 波段構(gòu)成的假彩色影像得到正射影像。精度的檢驗采用檢查點的辦法，以一部分已知坐標的點為檢查點，與所得DEM 對應點的坐標進行比對。對于相對DEM 坐標中誤差為10 30m，對于絕對DEM，當控制點精度在15 30 m 時，坐標的中誤差為7 50m，絕對DEM 的坐標中誤差取決于地面控制點的質(zhì)量，兩種DEM 的坐標中誤差同時又取決于連接點的提取精度和影像自動相關(guān)的結(jié)果。通過研究證實ASTER傳感器的數(shù)據(jù)在地面控制點分布均勻，精度較高的情況下可以較高精度的成圖，它可以用于全球無圖區(qū)的測圖，也可以用于一些低分辨率傳感器數(shù)據(jù)（如MODIS 數(shù)據(jù)）的正射糾正。與異軌立體成像的傳