基于ERDAS遙感制圖課程設(shè)計(jì)實(shí)驗(yàn)報(bào)告

1、基于ERDAS遙感制圖課程設(shè)計(jì) 一、課程設(shè)計(jì)目的 經(jīng)過一次系統(tǒng)的遙感圖像數(shù)據(jù)的課程設(shè)計(jì)，使同學(xué)們掌握ERDAS遙感圖像處理軟件的圖像處理技巧，通過數(shù)據(jù)的下載，融合，校正，裁剪，拼接、以及遙感專題圖的制作等過程，提高同學(xué)們的實(shí)際動(dòng)手能力。二、課程設(shè)計(jì)內(nèi)容1、在馬里蘭大學(xué)地球科學(xué)數(shù)據(jù)下載數(shù)據(jù)。第一步：打開馬里蘭大學(xué)網(wǎng)站圖1 馬里蘭大學(xué)衛(wèi)星影像下載頁面第二步：在網(wǎng)頁上選擇上圖紅色箭頭所指MapSearch進(jìn)入，選擇TM數(shù)據(jù)，經(jīng)緯度Lat/Long標(biāo)簽，初步選定山東省的經(jīng)緯度區(qū)間為精度114122，緯度3438圖2需要的圖像范圍第三步：點(diǎn)擊Draw標(biāo)簽，畫出所要的地圖區(qū)域；第四步：點(diǎn)擊Preview&

2、amp;Download進(jìn)入下載頁面圖3 下載所需要的數(shù)據(jù)第五步：下載數(shù)據(jù)。（注：下載使用的數(shù)據(jù)和姚曉玲的數(shù)據(jù)為同一數(shù)據(jù)）2、在國際科學(xué)數(shù)據(jù)平臺(tái)上下載遙感數(shù)據(jù)圖4 國際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)下載方法同上面基本相同 ，先注冊(cè)自己的用戶名，讓后登陸根據(jù)代號(hào)下載自己所需的影像數(shù)據(jù)即可，此處不再重復(fù)下載。3、遙感影像數(shù)據(jù)幾何校正第一步：打開要處理的影像。第二步：?jiǎn)?dòng)幾何校正模塊，在View#視窗里選擇Raster|Geometric Correction.打開Set Geometric Model 對(duì)話框。選擇多項(xiàng)式幾何校正模型 PolynomialOK程序自動(dòng)打開Geo Correction Tools

3、對(duì)話框和Polynomial Model Properties 對(duì)話框.先選擇Close 關(guān)閉Polynomial Model Properties 對(duì)話框.程序自動(dòng)打開GCP Tool Reference Setup 對(duì)話框.選擇Keyboard Only.OK.程序自動(dòng)打開 Reference Map Information 提示框。.選擇Map Units: Meters.添加地圖投影參數(shù)，如下圖：選擇OK 確定地圖投影參數(shù)，并關(guān)閉上圖。.選擇OK，確定 Reference Map Information，并關(guān)閉提示框。.并自動(dòng)打開采集控制點(diǎn)對(duì)話框。答案：圖象幾何校正第一步：顯示圖象文件

4、在視窗中打開需要校正的Landsat TM 圖象：lanzhoucity.img.第二步：?jiǎn)?dòng)幾何校正模塊在Viewer#1 的菜單條中，選擇Raster|Geometric Correction打開Set Geometric Model 對(duì)話框.選擇多項(xiàng)式幾何校正模型 PolynomialOK.程序自動(dòng)打開Geo Correction Tools 對(duì)話框和Polynomial Model Properties 對(duì)話框. .先選擇Close 關(guān)閉Polynomial Model Properties 對(duì)話框.程序自動(dòng)打開GCP Tool Reference Setup 對(duì)話框.選擇Keyboa

5、rd Only.OK.程序自動(dòng)打開 Reference Map Information 提示框。.選擇Map Units: Meters.添加地圖投影參數(shù)，如下圖：.選擇OK 確定地圖投影參數(shù)，并關(guān)閉上圖。.選擇OK，確定 Reference Map Information，并關(guān)閉提示框。.并自動(dòng)打開采集控制點(diǎn)對(duì)話框。GCP 的具體采集過程:在圖像幾何校正過程中，采集控制點(diǎn)是一項(xiàng)非常重要和相當(dāng)重要的工作，具體過程如下：.在GCP 工具對(duì)話框中點(diǎn)select GCP 圖標(biāo)，進(jìn)行GCP 選擇狀態(tài)。.在view#1 中移動(dòng)關(guān)聯(lián)方框位置，尋找明顯地物特征點(diǎn)，作為輸入GCP。.在GCP 工具對(duì)話框中點(diǎn)擊

6、Great GCP 圖標(biāo)，并在view#2 中點(diǎn)擊左鍵定點(diǎn)，GCP 數(shù)據(jù)表將記錄一個(gè)輸入的GCP，包括編號(hào)、標(biāo)識(shí)碼、X、Y 坐標(biāo)。.在GCP 工具對(duì)話框中輸入地圖參數(shù)坐標(biāo)X、Y。.不斷重復(fù)上述步驟，采集若干GCP，直到滿足所選是的幾何校正模型為止。采集地面檢查點(diǎn)以上所采集的GCP 為控制點(diǎn)，用于建立轉(zhuǎn)換方模型及 多項(xiàng)式方程，地面檢查點(diǎn)，則用于檢驗(yàn)所建立的轉(zhuǎn)換方程的精度和實(shí)用性，具體過程如下：.在GCP TOOL 菜單條中選擇GCP 類型：Edit/Set Point Type check。.在GCP TOOL 菜單條中確定GCP 匹配參數(shù)：Edit/Point matching 打開GCP M

7、atching 對(duì)話框，并確定參數(shù)。.確定地面檢查點(diǎn)，其操作與選擇控制點(diǎn)完全一樣。.計(jì)算檢查點(diǎn)誤差：在GCP TOOL 工具條中，點(diǎn)擊Compute Error 圖標(biāo)，檢查點(diǎn)的誤差就會(huì)顯示在GCP TOOL 的上方(如下圖)，只有所有檢查的誤差小于一個(gè)像元時(shí)，才能進(jìn)行以下的步驟。.在Geo Correction Tools 對(duì)話框中選擇Model Properties 圖標(biāo)打開選擇或檢查參數(shù)，然后選擇close 關(guān)閉。圖象重采樣.在Geo Correction Tools 對(duì)話框中選擇Image Resample 圖標(biāo)打開Image Resample 對(duì)話框，并定義重采樣參數(shù)。.輸出圖像文件名

8、(output .選擇重采樣方法(Resample Method):Nearest Neighbor.定義輸出圖像在圖與像元大小。.設(shè)查輸出統(tǒng)計(jì)中忽略零值。.選擇OK 啟動(dòng)重采樣進(jìn)程，并關(guān)閉Image Resample 對(duì)話框。4、遙感影像數(shù)據(jù)融合第一步：在ERDAS 圖標(biāo)面板工具條中，點(diǎn)擊 Interpreter-Utilities-Layer Stack，出現(xiàn) Layer Selection and Stacking 的對(duì)話框。第二步：在 Layer Selection and Stacking 對(duì)話框中，依次選擇并加載（Add）單波段img 圖像： 輸入單波段圖像文件（Input Fil

9、e: *.img）: p122r33_5t930615_nn1.imgAdd 輸入單波段圖像文件（Input File: *.img）: p122r33_5t930615_nn2.imgAdd 輸入單波段圖像文件（Input File: *.img）: p122r33_5t930615_nn3.imgAdd輸入單波段圖像文件（Input File: *.img）: p122r33_5t930615_nn4.imgAdd輸入單波段圖像文件（Input File: *.img）: p122r33_5t930615_nn5.imgAdd輸入單波段圖像文件（Input File: *.img）: p1

10、22r33_5t930615_nn6.imgAdd.輸入單波段圖像文件（Input File: *.img）: p122r33_5t930615_nn7.imgAdd選擇組合多波段圖像文件路徑和名字（Output File:*.img）: p122r33_5t930615_Result.img 第三步：選擇OK ，執(zhí)行并完成波段組合。5、遙感圖像的拼接第一步：在 ERDAS 圖標(biāo)面板工具條中，點(diǎn)擊Datapreparation-Mosaicc lmages，出現(xiàn)Mosaic Tool 視窗 第二步：加載 Mosaic 圖像。在Mosaic Tool視窗菜單條中，Edit-Add images，

11、出現(xiàn) Add Images for Mosaic 對(duì)話框。依次加載窗拼接的圖像。 Image Area Option選Crop Area。第三步：加載圖像完成后，在 Mosaic Tool 視窗工具條中，點(diǎn)擊 set Input Mode 圖標(biāo) ，進(jìn)入設(shè)置圖象模式的狀態(tài)，利用所提供的編輯工具，進(jìn)行圖象疊置組合調(diào)查。 第四步：圖象匹配設(shè)置，點(diǎn)擊 Edit -Image Matching ，出現(xiàn) Matching options 對(duì)話框，設(shè)置匹配，方法：Overlap Areas。 第五步：在 Mosaic Tool 視窗菜單條中，點(diǎn)擊 Edit/set Overlap Function打開 se

12、t Overlap Function 對(duì)話框 設(shè)置以下參數(shù)： .設(shè)置相交關(guān)系(Intersection Method):No Cutline Exists。 .設(shè)置重疊圖像元灰度計(jì)算(select Function):Average。 Apply close 完成。 第六步：運(yùn)行 Mosaic 工具 在 Mosaic Tool 視窗菜單條中,點(diǎn)擊 Process-Run Mosaic打開 Run Mosaic 對(duì)話框。 設(shè)置下列參數(shù)：確定輸出文件名：shandong.img 確定輸出圖像區(qū)域：ALL OK ，進(jìn)行圖像拼接。6、遙感影像數(shù)據(jù)裁剪不規(guī)則裁剪步驟：第一步：采用網(wǎng)上下載的MapInfo

13、格式的數(shù)據(jù)導(dǎo)入SuperMap Deskpro 6,然后導(dǎo)出為.shp文件。第二步：打開圖標(biāo)面板，點(diǎn)擊InterPreter|Utilities|Vector To Raster，打開Vector To Raster設(shè)置參數(shù)如下，Input Vector File選擇shandong.shp Output Image File取名mask.img Cell Size 為 X:0.01 Y:0.01OK，轉(zhuǎn)換成功。Mask.img圖像如圖：第三步：InterPreter|Utilities|Mask菜單，選擇Mask選項(xiàng)打開Mask對(duì)話框，在Mask對(duì)話框中設(shè)置以下參數(shù)：輸入文件名稱：shand

14、ong.img輸入掩膜文件：mask.img輸出文件名稱：shandong_result.img;OK，開始裁剪，等待完成。完成后裁剪結(jié)果在View窗口中效果如圖：7、遙感影像數(shù)據(jù)增強(qiáng)。7.1 卷積增強(qiáng)處理ERDAS 圖標(biāo)面板工具條，點(diǎn)擊Interpreter/spatial Enhancement / convolution打開convolution 對(duì)話框，并設(shè)置如下參數(shù)： 圖 卷積增強(qiáng)處理前和卷積增強(qiáng)處理7.2 紋理分析紋理分析通過在一定的窗口內(nèi)進(jìn)行二次變異分析或三次對(duì)稱分析，使圖像的紋理結(jié)構(gòu)得到增強(qiáng)，具體過程如下：ERDAS 圖板面板工具條中，點(diǎn)擊Interpreter|spatial Enhancement|Texture，打開Te