遙感影像數(shù)字處理課程設計報告-徐州市土地利用動態(tài)變化分析

1、<p>　　存檔日期： 存檔編號： </p><p>　　徐州師范大學科文學院</p><p><b>　　本科生課程設計</b></p><p>　　設 計 題 目：徐州市土地利用動態(tài)變化分析</p><p>　　姓

2、 名： </p><p>　　學 號： 10253006 </p><p>　　院 系： 城市與環(huán)境 </p><p>　　年 級 ： 10級 </p

3、><p>　　指 導 教 師： </p><p><b>　　中文摘要 </b></p><p>　　隨著遙感技術的飛速發(fā)展，遙感應用的逐步深入，遙感圖像處理系統(tǒng)如雨后春筍般</p><p>　　不斷涌現(xiàn)。遙感技術作為對地觀測，提取地表最現(xiàn)勢狀況的最有利工具，被廣泛應用在</p><

4、;p>　　各行各業(yè)，包括測繪，自然資源管理，林業(yè)，水利，交通，環(huán)境保護，電力電信、防震</p><p>　　減災、城市規(guī)劃、國防軍事等，ERDAS IMAGINE 作為遙感界的排頭兵，為大家提供了不</p><p>　　僅僅是增強，濾波，糾正，融合等簡單的基本應用，而是提供了強大的工具，使你在定</p><p>　　量化的分析方面，系統(tǒng)功能的可擴充性方面使您更

5、加得心應手，如專家分類，子象元分</p><p>　　類（混合象元），三維可視化分析，數(shù)字攝影測量等，同時，還給您帶來與GIS 一體化</p><p>　　集成的解決方案，如查詢檢索編輯ArcInfo 的地理信息，建立矢量層后的人工解譯，直</p><p>　　接得到目標的矢量數(shù)椐，還可將分好類的專題影像轉換成ArcInfo 的矢量數(shù)據(jù)（Coverage, <

6、/p><p>　　Shape File）,使分析的結果可以直接為地理信息系統(tǒng)管理與應用，從而發(fā)揮更大的作用。 </p><p>　　關鍵詞：徐州市；遙感圖像處理；土地利用</p><p><b>　　英文摘要 </b></p><p>　　Along with the rapid development of remote

7、sensing technology and application of remote sensing image processing, springing up system is emerging. Sensing technology as the earth observation, the condition of surface processing. along tool, the most w

8、idely used in all walks of life, including the surveying and mapping, and natural resource management, forestry, water conservancy,transportation, environment protection, electric power, telecom mitigation, urban plannin

9、g, national defen</p><p>　　Keywords: Xuzhou City; image processing;landuse</p><p><b>　　目錄</b></p><p>　　實習目的 4</p>

10、<p>　　實習內容 4</p><p>　　實習要求 4</p><p>　　實習步驟 4&

11、lt;/p><p>　　（一）明確目標 4</p><p>　?。ǘ┦占Y料 4</p><p>　　1、數(shù)據(jù)源

12、 4</p><p>　　2、研究區(qū)概況 4</p><p>　?。ㄈ?shù)據(jù)處理 5</p><p>　　1、幾何校正

13、 5</p><p>　　2、圖像剪裁 6</p><p>　?。ㄋ模┍O(jiān)督分類 6</p><p>　　1、構建分類系統(tǒng)

14、 6</p><p>　　2、定義分類模板 7</p><p>　　3、評價分類模板 9</p><p>　　4、執(zhí)行監(jiān)督分類

15、 9</p><p>　　5、評價分類結果 10</p><p>　　6、分類后處理 13</p><p>　　（五）矩陣分析

16、 13</p><p>　　總結 14</p><p>　　參考文獻 15</p><p>　　附錄

17、 15</p><p><b>　　一、實習目的</b></p><p>　　通過課程設計，鞏固和深化《遙感圖像處理方法》課程所學的理論知識，掌握典型地物在TM影像上的構像特征，熟悉遙感圖像處理的方法和步驟，學會運用EDRAS遙感圖像處理軟件進行圖像增強、圖像監(jiān)督分類、對比分析，培養(yǎng)和提高分析和解決實際問

18、題的能力。二、實習內容</p><p>　　1、學會數(shù)據(jù)預處理的主要方法并獨立完成數(shù)字圖像分類的前期準備工作；</p><p>　　2、掌握數(shù)字圖像分類的原理與方法以及詳細的工作流程，并按照要求出圖；</p><p>　　3、遙感圖像的對比分析（矩陣分析）及對提取的用地信息進行簡要的變化分析并討論結果，撰寫課程設計報告。</p><p>&l

19、t;b>　　三、實習要求</b></p><p>　　1．了解項目設計的步驟。 </p><p>　　2．學會如何通過軟件實現(xiàn)數(shù)據(jù)的處理。包括：了解數(shù)據(jù)源、數(shù)據(jù)格式轉化、投影轉換、幾何糾正、圖像監(jiān)督分類、精度評估、轉移矩陣的生成等。 </p><p>　　3．掌握信息提取的不同方法。</p><p><b>　　四

20、、實習步驟</b></p><p><b>　?。ㄒ唬┟鞔_目標 </b></p><p>　　該項目的目標是研究某地區(qū)土地利用動態(tài)變化情況。即對該地區(qū)不同時期的土地利用信息進行提取，進而對比不同時期土地利用現(xiàn)狀以獲取變化情況。分析輸出的結果圖應該清楚地表明哪些用地發(fā)生了變化？是怎樣變化的？變化了多少？說明什么問題。 </p><p>

21、;<b>　　（二）收集資料 </b></p><p><b>　　1. 數(shù)據(jù)源 </b></p><p>　　所需數(shù)據(jù)及詳細使用情況如下：</p><p>　　將這校正好的兩幅圖分別裁剪相同大小的部分進行監(jiān)督分類，然后將個要素提取出來分別對比，比較2000年與1987年的變化之處。</p><p>

22、;<b>　　2、研究區(qū)概況</b></p><p>　　徐州市位于江蘇省的西北部，東經(jīng)116°22′～118°40′、北緯33°43′～ 34°58′之間。東西長約210公里，南北寬約140公里，總面積11258平方公里，占江蘇省總面積的11%。地處蘇、魯、豫、皖四省交界，為懂不沿海與中不地帶、上海經(jīng)濟區(qū)與渤海經(jīng)濟圈的結合部?！皷|襟淮海，西接中原，南

23、屏江淮，北扼齊魯”，素有“五省通衢”之稱。京滬、隴海兩大鐵路在此交匯，京杭大運河傍城而過貫穿徐州南北，公路四通八達，北通京津，南達滬寧，西接蘭新，東抵海濱，為全國重要水陸交通樞紐和東西、南北經(jīng)濟聯(lián)系的重要“十字路口”。徐州劇南京、濟南、鄭州、合肥皆為300余公里。徐州屬暖溫帶半濕潤季風氣候，四季分明，夏無酷暑，冬無嚴寒，年平均氣溫，年降雨量886毫米，無霜期200--220天，日照2100--2400小時。熱量多余華北，光能優(yōu)于蘇南，地

24、形以平原為主，低山和丘陵地相間，海拔19，沂河、沭河、京杭大運河、古黃河穿境而過，駱馬湖、微山湖分布南北兩緣。</p><p><b>　　（三）、數(shù)據(jù)預處理</b></p><p><b>　　1、幾何校正</b></p><p>　　第一步：顯示圖像文件</p><p>　　在ERDAS圖標面板

25、中雙擊Viewer圖標，打開兩個窗口（Viewer#1/Viewer#2），并將兩個窗口平鋪設置，接著執(zhí)行以下操作：</p><p>　?。?）、在ERDAS圖標面板菜單條，單擊Session|Tile Viewers命令。</p><p>　?。?）、在Viewer#1窗口中打開需要校正的TM圖像，即1987年的圖像。</p><p>　?。?）、在Viewer#

26、2窗口中打開作為地理參考的校正過的SPOT圖像，即2000年的圖像。</p><p>　　第二步：啟動幾何校正模塊</p><p>　　（1）、在Viewer#1菜單條，單擊Raster|Geometric Correction命令。</p><p>　?。?）、打開Set Geometric Moder對話框。</p><p>　?。?）、

27、選擇多項式幾何校正計算模型為Polynomial。</p><p>　　（4）、單擊OK按鈕。</p><p>　?。?）、同時打開Geo Correction Tools對話框和Polynomial Model Properties窗口。在Polynomial Model Properties窗口中，定義多項式模型參數(shù)及投影參數(shù)。</p><p>　?。?）、定義

28、多項式次方為2。</p><p>　?。?）、定義投影參數(shù)（Projection）.</p><p>　?。?）、單擊Apply按鈕應用或單擊Close按鈕關閉。</p><p>　?。?）、打開GCP Tool Reference Setup對話框。</p><p>　　第三步：啟動控制點工具</p><p>　　（

29、1）、在GCP Tool Reference Setup對話框中選擇采點模式，即選擇Existing Viewer單選按鈕。</p><p>　?。?）、單擊OK按鈕。</p><p>　?。?）、打開Viewer Selection Instructions 指示器。</p><p>　?。?）、在顯示作為地理參考圖像panAtlanta.img的Viewer#2

30、中單擊。</p><p>　　（5）、打開Reference Map Infornation提示框。</p><p>　?。?）、單擊OK按鈕。</p><p>　?。?）、整個屏幕自動變化，其中包含兩個主窗口、兩個放大窗口、兩個關聯(lián)方框（分別位于兩個窗口中，指示放大窗口與主窗口的關系）、控制點工具對話框和幾何校正工具等。表明控制點工具被啟動，進入控制點采集狀態(tài)。&

31、lt;/p><p>　　第四步：采集地面控制點</p><p>　?。?）、在GCP工具對話框中單擊Select GCP圖標，進入GCP選擇狀態(tài)。</p><p>　　（2）、在GCP數(shù)據(jù)表中將輸入GCP的顏色設置為比較明顯的黃色。</p><p>　?。?）、在Viewer#1中移動關聯(lián)方框位置，尋找明顯的地物特征點，作為輸入GCP。</

32、p><p>　　（4）、在GCP工具對話框中單擊Create GCP圖標，并在Viewer#3中單擊定點，GCP數(shù)據(jù)表將記錄一個輸入GCP,包括其編號、標識碼、X坐標、Y坐標。</p><p>　?。?）、在GCP工具對話框中單擊Select GCP圖標，重新進入GCP選擇狀態(tài)。</p><p>　　（6）、在GCP數(shù)據(jù)表中將輸入GCP的顏色設置為比較明顯的紅色。<

33、;/p><p>　　（7）、在Viewer#2中移動關聯(lián)方框位置，尋找明顯的地物特征點，作為輸入GCP。</p><p>　?。?）、在GCP工具對話框中單擊Create GCP圖標，并在Viewer#4中單擊定點，GCP數(shù)據(jù)表將記錄一個輸入GCP,包括其編號、標識碼、X坐標、Y坐標。</p><p>　?。?）、在GCP工具對話框中單擊Create GCP圖標，重新進

34、入GCP選擇狀態(tài)；并將光標移回到Viewer#1，準備采集另一個輸入控制點。</p><p>　?。?0）、不斷重復上述步驟，采集若干GCP，直到滿足所選定的幾何校正模型為止。而后，每采集一個Input GCP,系統(tǒng)就自動產(chǎn)生一個Ref.GCP,通過移動Ref.GCP可以逐步優(yōu)化校正模型。</p><p>　　第五步：采集地面檢查點</p><p>　　以上采集的G

35、CP的類型均為控制點，用于控制計算、建立轉換模型及多項式方程。下面要采集的GCP的類型均為檢查點。如果控制點的誤差比較小的話，也可以不采集地面檢查點。</p><p><b>　　第六步：圖像重采樣</b></p><p>　　在Geo Correction Tools對話框中單擊Image Resample圖標，打開Resample對話框，在Resample對話框中

36、，定義重采樣參數(shù)。</p><p>　?。?）、輸出圖像文件名（Output File）為1987_rectify.img.</p><p>　?。?）、選擇重采樣方法（Resample Method）為Nearest Neighbor.</p><p>　　（3）、定義輸出圖像范圍（Output Corners）,在ULX、ULY、URX、LRY微調框中分別輸入需

37、要的數(shù)值。</p><p>　?。?）、定義輸出象元大?。∣utput Cell Sizes）,X值30/Y值30.</p><p>　?。?）、設置輸出統(tǒng)計中忽略零值，即選中Ignore Zero in Stats 復選框。</p><p>　?。?）、設置重新計算輸出默認值（Recalcuate Output Dafaults），設Skip Factor 為10

38、。</p><p>　　、單擊OK按鈕(關閉Resaample 對話框，啟動重采樣過程)。</p><p><b>　　2、圖像裁剪</b></p><p>　　在ERDAS圖標面板菜單條單擊Main|Data Preparation|Subset Image命令，打開Subset對話框或在ERDAS圖標面板工具條單擊Data Prep 圖標|

39、Subset命令，打開Subset對話框</p><p>　　在Subset對話框中需要設置下列參數(shù)：</p><p>　?。?）、輸入文件名稱（Input File）為1987_rectify.img。</p><p>　?。?）、輸出文件名稱（Output File）為1987_sub.img。</p><p>　?。?）、坐標類型（Co

40、ordinate Type）為File。</p><p>　?。?）、確定裁剪范圍（Subset Definition）,在ULX、ULY、LRX、LRY微調框中分別輸入需要的數(shù)值，使得圖像大小為512*512。</p><p>　?。?）、輸出數(shù)據(jù)類型（Output Data Type）為Unsigned 8 Bit。</p><p>　　（6）、輸出文件類型（O

41、utput Layers Type）為Continuous。</p><p>　?。?）、輸出統(tǒng)計忽略零值，選中Ignore Zero In Output Stats復選框。</p><p>　　（8）、輸出象元波段（Select Layers）為2:5（表示選則2、3、4、5這四個波段）。</p><p>　　單擊OK按鈕（關閉Subset Image對話框，執(zhí)行

42、圖像裁剪）。</p><p><b>　　、圖像分類</b></p><p><b>　　1、構建分類系統(tǒng)</b></p><p>　　監(jiān)督分類比非監(jiān)督分類更多地要用戶來控制，常用于對研究區(qū)域比較了解的情況。在監(jiān)督分類過程中，首先選擇可以識別或者借助其他信息可以斷定其類型的象元建立模板，然后基于該模板使計算機系統(tǒng)自動識別具

43、有相同特性的象元。對分類結果進行評價后再對模板進行修改，多次反復后建立一個比較標準的模板，并在此基礎上最終進行分類。監(jiān)督分類一般要經(jīng)過一下幾個步驟：建立模板（訓練樣本）、評價模板、確定初步分類結果、檢驗分類結果、分類后處理、分類特征統(tǒng)計、柵格矢量轉換。</p><p><b>　　2、定義分類模板</b></p><p>　　第一步： 顯示需要分類的圖像</p&

44、gt;<p>　　在窗口中顯示1987_sub.img, 選擇 Fit to Frame，其他使用默認設置）。</p><p>　　第二步：打開分類模板編輯器</p><p>　　在ERDAS 圖標面板菜單條，單擊Main|Image Classification|Classification|Signature Editor 命令，打開Signature Editor 窗口

45、。</p><p>　　在ERDAS圖標面板工具條，單擊Classifier 圖標| Classification|Signature Editor 命令，打開Signature Editor 窗口。</p><p>　　第三步：調節(jié)分類屬性字段</p><p>　　在Signature Editor 窗口菜單條，單擊View|Column 命令，打開View Si

46、gnature Columns 對話框。</p><p>　　在View Signature Columns 對話框，</p><p>　　1.單擊第一個字段的Column 列并向下拖動鼠標直到最后一個字段，此時，所有字段都被選擇上，并用黃色（默認色）表示出來。</p><p>　　2.按住Shift鍵的同事分別單擊Red、Green、Blue三個字段，Red、Gr

47、een、Blue三個字段將分別從選擇幾種被清除。</p><p>　　3.單擊Apply 按鈕，關閉View Signature Columns 對話框 。</p><p>　　從View Signature Columns 對話框可以看到，Red、Green、Blue三個字段將不再顯示。</p><p>　　第四步：獲取分類模板信息</p><

48、p>　　應用AOI擴展工具在原始圖像獲取分類模板信息</p><p>　　在1987_sub.img圖像的菜單條件單擊AOI|Seed Properties 命令，打開Region Growing Properties 對話框。</p><p>　　在Region Growing Properties 對話框中，設置下列參數(shù)：</p><p>　　（1）選

49、擇相鄰象元擴展方式（Neighhborhood）：選擇按4個相鄰象元擴展，包括兩種擴展方式：表示以上、下、左、右4個象元作為相鄰進行擴展和表示種子象元周圍的9個象元都是擴展的相鄰象元。</p><p>　?。?）選擇區(qū)域擴展的地理約束條件（Geographic Constrains），包括約束條件：</p><p>　　A.面積及約束Area確定每個AOI所包含的最多象元數(shù)（或者面積）。&

50、lt;/p><p>　　B.距離約束Distance 確定AOI所包含象元距種子象元的最大距離。</p><p>　　C.這兩個約束可以只設置一個，也可以設置兩個或者一個也不設。</p><p>　?。?）確定區(qū)域擴展的地理約束參數(shù)，約束面積大小為300/采用單位象元Pixels。</p><p>　　（4）設置波普歐氏距離為10（這是判斷相鄰象

51、元與種子象元平均值之間的最大波普歐氏距離，大于該距離的相鄰象元將不被接受）。</p><p>　　（5）單擊Options按鈕，打開Region Grow Options 對話框，設置區(qū)域擴展過程中的算法。</p><p>　　（6）選擇Include Island Polygons和 Update Region Mean 算法。</p><p>　　下面的操作將在

52、Raster工具面板、圖像窗口、Signature Editor對話框三者之間交替進行。</p><p>　?。?）在Raster工具面板繪制單擊，進入AOI繪制狀態(tài)。</p><p>　?。?）在圖像窗口中選擇綠色區(qū)域（農田），繪制一個多邊形AOI。</p><p>　?。?）在Signature Editor 對話框，單擊Create New Signature

53、 圖標，將多邊形AOI區(qū)域架子到Signature Editor分類模板屬性表中。</p><p>　?。?）在圖像窗口中選擇另一個綠色區(qū)域，再繪制一個多邊形AOI。</p><p>　　（5）同樣在Signature Editor 對話框，單擊Create New Signature 圖標，將多邊形AOI區(qū)域架子到Signature Editor分類模板屬性表中。</p>

54、<p>　?。?）重復上述兩部操作過程，選擇圖像中您認為屬相相同的多個綠色區(qū)域單擊生成若干擴展AOI，并將其作為模板一次加入到Signature Editor分類模板屬性表中。</p><p>　?。?）按住Shift 鍵，同時在Signature Editor 分類模板屬性表中依次單擊選擇Class#字段下面的分類編碼，將您上面加入的多個綠色區(qū)域擴展AOI模板全部選定。</p><

55、p>　?。?）在Signature Editor 工具條，單擊Merge Signatures 圖標，將多個綠色區(qū)域AOI模板合并，生成一個綜合的新模板，其中包含了合并前的多個模板。</p><p>　　（9）在Signature Editor 菜單條，單擊Edit|Delete 命令，刪除合并前的多個模板。</p><p>　?。?0）在Signature Editor 屬性表，

56、改變合并生成的分類模板的屬性：包括名稱與顏色分類名稱（Signature Name ）：Forest/顏色（Color）：綠色。</p><p>　?。?1）重復上述所有操作過程，在圖像窗口選擇多個黑色區(qū)域定義種子象元，自動生成擴展AOI（水體），一次加載到Signature ediotor 分類屬性表，并執(zhí)行合并生成綜合的水體分類模板，然后確定分類模板名稱和顏色。</p><p>　?。?/p>

57、12）同樣重復上述所有操作過程，繪制多個藍色區(qū)域AOI（建筑）、多個綠色區(qū)域AOI（農田）等，加載、合并、命名，建立新的模板。</p><p>　?。?3）如果將多有的類型都建立了分離模板的話，就可以保存分類模板。</p><p>　　第五步：保存分類模板</p><p>　　第（4）步分別用不同方法產(chǎn)生了分類模板，現(xiàn)在需要將分類模板保存起來，以便隨后依據(jù)分類模板進

58、行監(jiān)督分類。</p><p>　　在Signature ediotor窗口菜單條，</p><p>　?。?）單擊File|Save 命令。</p><p>　?。?）打開Save Signature File As 對話框。</p><p>　　（3）確定是保存多有模板（All）或只保存被選中的模板（Selected）。</p>

59、<p>　?。?）確定保存分類模板文件的目錄和文件名（*.sig）。</p><p>　?。?）單擊OK按鈕（關閉 Save Signature File As 對話框，保存模板）。</p><p><b>　　評價分類模板</b></p><p><b>　　直方圖繪制</b></p><

60、;p>　　直方圖繪制工具通過分類分析類別的直方圖模板進行評價和比較，本功能可以同時對一個或多個類別制作直方圖，如果處理對象是耽擱類別（選擇Single Signature），那就是當前活動類別（位于“>”符號旁邊的那個類別），如果是多個類別的直方圖，那就是處于選擇幾種的類別。</p><p>　　在Signature Editor 窗口菜單條，單擊View|Histograms 命令，打開Histog

61、rams Control Panel對話框中，需要設置下列參數(shù)：</p><p>　　確定分類模板數(shù)量（Signature）為Single Signature。</p><p>　　確定分類波段數(shù)量（Bands）為Single Signature。</p><p>　　確定應用哪個波段（Band No.）為4。</p><p>　　單擊Pio

62、t按鈕（繪制分類直方圖）。</p><p>　　顯示繪制分類的分類直方圖。</p><p><b>　　4、執(zhí)行監(jiān)督分類</b></p><p>　　在ERDAS圖標面板菜單條，單擊Main|Image Classification 命令，打開Classification 對話框。</p><p>　　在ERDAS圖標

63、面板工具條，單擊Classifier 圖標|Classification |Supervised Classification 命令，打開Supervised Classification對話框。</p><p>　　在Supervised Classification對話框中，需要確定一下參數(shù)：</p><p>　?。?）確定輸入原始文件（Input Raster File）為1987_

64、sub.img。</p><p>　?。?）定義輸出分類文件（Classified File）為1987_sub_sup.img。</p><p>　?。?）確定分類模板文件（Input Signature File）為1987_sub.img。</p><p>　?。?）選擇分類輸出距離文件為Distance File(用于分類結果進行閾值處理)。</p&g

65、t;<p>　?。?）定義分類距離文件(Filename)為1987_sub_dist.img。</p><p>　　（6）選擇非參數(shù)規(guī)則（Non-Parametric Rule）為Feature Space。</p><p>　　（7）選擇疊加規(guī)則（Overlay Rule）為Parametric Rule。</p><p>　　（8）選擇未分類規(guī)則

66、（Unclassified Rule）為Parametric Rule。</p><p>　　(9)選擇參數(shù)規(guī)則（Parametric Rule）為Maximum Likelihood。</p><p>　　(10)取消選中Classify zeros 復選框（分類過程中是否包括0值）。</p><p>　　(11)單擊OK按鈕（執(zhí)行監(jiān)督分類，關閉Supervise

67、d Classification 對話框）。</p><p>　　在Supervised Classification 對話框中，還可以定義分類圖的屬性表項目：</p><p>　　單擊Attribue Options按鈕。</p><p>　　打開Attribue Options對話框。</p><p>　　在Attribue Option

68、s對話框上進行選擇。</p><p>　　單擊OK（關閉Attribue Options對話框）。</p><p>　　返回Supervised Classification 對話框。</p><p><b>　　5、評價分類結果</b></p><p><b>　　閾值處理</b></p&g

69、t;<p>　　第1步：顯示分類圖像并啟動閾值處理</p><p>　　首先要在窗口中打開分類后的專題圖像，然后啟動閾值處理功能：</p><p>　　在ERDAS圖標面板菜單條，單擊Main|Image Classification|Classification命令。</p><p>　　在ERDAS圖標面板工具條，單擊Classifier 圖標|

70、Classification |Threshold命令，打開Threshold窗口。</p><p>　　第2步：確定分類圖像和距離圖像</p><p>　　在Threshold窗口菜單條，單擊File|Open命令，打開Open Files對話框。</p><p>　　在Open Files對話框：</p><p>　　確定專題分類圖像（C

71、lassified Image）為1987_sub_sup.img。</p><p>　　確定分類距離圖像（Distance Image）為1987_sub_dist.img。</p><p>　　單擊OK按鈕（關閉Open Files對話框）。</p><p>　　返回Threshold窗口</p><p>　　第3步：視圖選擇及直方圖計算

72、</p><p>　　在Threshold窗口菜單條，單擊View|Select Viewer命令，單擊顯示分類專題圖像的窗口。</p><p>　　在Threshold窗口菜單條，單擊Histogram|Compute命令（計算各個類別的距離直方圖，如果需要的話該直方圖可通過Threshold窗口菜單單擊Histogram|Save而保存為一個模板文件*.sig文件）。</p>

73、;<p>　　第4步：選擇類別并確定閾值</p><p>　　在Threshold窗口的分類屬性表中，選擇專題類別。</p><p>　　移動“>”符號到指定的專題類別旁邊。</p><p>　　在菜單條單擊Histogram|View命令。</p><p>　　選定類別的Distance Histogram被現(xiàn)實出來。&

74、lt;/p><p>　　拖動Histogram X軸上的箭頭到想設置閾值的位置，Threshold窗口中的Chi-square值自動發(fā)生變化，表明該類別的閾值設定完畢。</p><p>　　重復上述步驟，依次設定每一個類別的閾值。</p><p>　　第5步：顯示閾值處理圖像</p><p>　　在Threshold窗口菜單條，單擊View|Vi

75、ew Colors|Default Colors命令（進行環(huán)境設置。選擇默認色彩是將閾值以外的象元顯示成黑色，而將屬于分類閾值之內的象元以類別顏色顯示）。</p><p>　　單擊Process|To Viewer命令。</p><p>　　閾值處理圖像將顯示在分類圖像之上，及形成一個閾值掩膜。</p><p>　　第6步：觀察閾值處理圖像</p>&

76、lt;p>　　在菜單條單擊Utility|Flicker對話框，可以進行如下設置：</p><p>　　設置閃爍速度（Speed）為500.</p><p>　　設置自動閃爍狀態(tài)，選擇Auto Mode(觀察類別與原圖像之間的對于關系)。</p><p>　　單擊Cancel按鈕（關閉Viewer Flicker對話框）。</p><p&g

77、t;　　第7步：保存閾值處理圖像</p><p>　　在Threshold窗口菜單條，</p><p>　　單擊Process|To Viewer命令。</p><p>　　打開Threshold to File對話框。</p><p>　　在Output Image 中確定要產(chǎn)生的文件的名字和目錄。</p><p>

78、<b>　　單擊OK按鈕。</b></p><p><b>　　分類精度評估</b></p><p>　　第1步：打開分類前原始圖像</p><p>　　在Viewer中打開分類前的原始圖像，以便進行精度評估。</p><p>　　第2步：打開精度評估對話框</p><p>

79、　　在ERDAS圖標面板菜單條，單擊Main|Image Classification |Classification|Accuracy Assessment命令。</p><p>　　在ERDAS圖標面板工具條，單擊Classifier 圖標|Classification |Accuracy Assessment命令，打開Accuracy Assessment窗口。</p><p>　　

80、第3步：打開分類專題圖像</p><p>　　在Accuracy Assessment窗口工具條，單擊Open File圖標或在Accuracy Assessment窗口菜單條,</p><p>　　(1)單擊 File|Open命令。</p><p>　　(2)打開Classified Image對話框。</p><p>　　(3)在Cla

81、ssified Image對話框中確定與窗口中對應的分類專題圖像。</p><p>　　(4)單擊OK按鈕（關閉Classified Image對話框）。</p><p>　　(5)返回Accuracy Assessment窗口。</p><p>　　第4步：鏈接原始圖像與精度評估窗口</p><p>　　在Accuracy Assessme

82、nt窗口工具條，單擊Select Viewer圖標。或在Accuracy Assessment窗口菜單條，</p><p>　　(1)單擊View|Select Viewer命令。</p><p>　　(2)將光標移在顯示原始圖像的窗口中單擊一下。</p><p>　　(3)原始圖像窗口與精度評估窗口相連接。</p><p>　　第5步：設置

83、隨機會點的色彩</p><p>　　在Accuracy Assessment窗口菜單條，</p><p>　　(1)單擊View|Change Colors命令。</p><p>　　(2)打開Change Colors對話框。</p><p>　　(3)在Points with no reference文本框確定沒有真實參考值的點的顏色。&

84、lt;/p><p>　　(4)在Points with no reference文本框確定有真實參考值的點的顏色。</p><p>　　(5)單擊OK按鈕（執(zhí)行參數(shù)設置）。</p><p>　　(6)返回Accuracy Assessment對話框。</p><p><b>　　第6步：產(chǎn)生隨機點</b></p>

85、<p>　　在Accuracy Assessment菜單條，</p><p>　　(1)單擊Edit|Create/Add Random Points命令。</p><p>　　(2)打開Add Random Points對話框。</p><p>　　(3)在Search Count微調框中輸入“1024”。 </p><p&g

86、t;　　(5)在Number of Points微調框中輸入“10”。</p><p>　　(6)在Distribution Parameters選項組中選擇Random單選按鈕。</p><p>　　(7)單擊OK按鈕（按照參數(shù)設置產(chǎn)生隨機點）。</p><p>　　(8)返回Accuracy Assessment窗口。</p><p>　

87、　第7步：顯示隨機點及其類別</p><p>　　在Accuracy Assessment窗口菜單條，</p><p>　　單擊View|Show All命令（所以隨機點均以第5步所設置的顏色顯示在窗口中）。</p><p>　　單擊Edit|Show Class Value命令（各點的類別號出現(xiàn)在數(shù)據(jù)表的class字段中）。</p><p&g

88、t;　　第8步：輸入?yún)⒖键c實際類別</p><p>　　在Accuracy Assessment對話框精度評估矩陣。</p><p>　　在Reference字段輸入各個隨機點的實際類值（只要輸入?yún)⒖键c的實際分類值，它在窗口中的色彩就變?yōu)榈?步設置的顏色）。</p><p>　　第9步：輸出分類評價報告</p><p>　　單擊Report|

89、Options命令。</p><p>　　分類評價報告輸出內容選項，單擊選擇參數(shù)。</p><p>　　單擊Report|Accuracy Report命令。</p><p>　　單擊Report|Cell Report命令。</p><p>　　所有報告將顯示在ERDAS文本編輯器窗口，可以保存為文本文件。</p><p

90、>　　單擊File|Save Table命令。</p><p>　　單擊File|Close命令。</p><p><b>　　6、分類后處理</b></p><p>　　在ERDAS圖標面板菜單條，單擊Main|Image Interpreter|GIS Analysis|Recode命令，打開Recode對話框。</p>

91、<p>　　在ERDAS圖標面工具條，單擊Interpreter 圖標|GIS Analysis|Recode命令，打開Recode對話框。</p><p>　　在Recode對話框，需要確定下列參數(shù)：(1)確定輸入文件（Input File）為1987_sub_sup.img 。</p><p>　　定義輸出文件（Output File）為1987_sub_recode.i

92、mg。</p><p>　　設置新的分類編碼（Setup Recode），單擊Setup Recode 按鈕。</p><p>　　打開Themaic Recode 表格。</p><p>　　根據(jù)需要改變New Value 字段的取值（直接輸入）。</p><p>　　單擊OK按鈕（關閉Themaic Recode 表格，完成新編碼輸入）。

93、</p><p>　　確定輸出數(shù)據(jù)類型（Output）5Unsigned 4 Bit。</p><p>　　單擊OK按鈕（關閉Recode對話框，執(zhí)行圖像重編碼，輸出圖像將按照New Value變換專題分類圖像屬性，產(chǎn)生新的專題分類圖像）。</p><p><b>　　、矩陣分析</b></p><p>　　矩陣分析功能

94、（Matrix）可以將兩個輸入分類專題圖像或矢量地圖數(shù)據(jù)，按照其專題屬性在空間上的重疊性產(chǎn)生一個新的圖像文件，新文件包含兩個輸入文件中重疊的分類專題屬性。</p><p>　　在ERDAS圖標面板菜單條，單擊Main|Image Interpreter |GIS Analysis |Matrix 命令，打開Matrix 對話框。</p><p>　　在Matrix 對話框中，需要確定下列參

95、數(shù)：</p><p>　　確定第一個輸入文件（Image or Vector File#1）為1987_sub_sup.img。</p><p>　　重編碼設置（Setup Recode）為打開重編碼表格進行編碼。</p><p>　　確定第二個輸入文件（Image or Vector File#2）為2000_sub_sup.img。</p><

96、;p>　　重編碼設置（Setup Recode）為打開重編碼表格進行編碼。</p><p>　　確定輸出文件（Output Recode）為1987_2000.img。</p><p>　　確定輸出統(tǒng)計時忽略零值，選中Ignore Zeor in Stats復選框。</p><p>　　設置兩幅圖像運算規(guī)則（Select Area by），選擇Union（并集

97、）單選按鈕。</p><p>　　去頂射出數(shù)據(jù)類型（Output）為Unsigned 16 bit。</p><p>　　單擊OK按鈕（關閉Matrix對話框，執(zhí)行矩陣分析）。</p><p><b>　　結果分析：</b></p><p>　　從兩幅圖對比來看，從1987到2000年居民點的分布不斷擴大，特別是沿云龍湖

98、地區(qū)。2000年裸地相對減少，有些山地的面積相對減少，林地面積也相對減少，耕地面積相對有所增加，河流面積也增多，擴大了灌溉的能力。</p><p>　　從整體上看，居民點空間擴散，變大，集中分布在沿河、沿湖地區(qū)，裸地顯著減少農田稍微比以前增多一點，變動不大。綠地減少，湖泊減少。</p><p><b>　　總結</b></p><p>　　遙感

99、導論課程設計結束了,在老師的耐心的講解和自主學習下,本次課程設計還算順利,每一個步驟都未出現(xiàn)大的問題,一些小問題也都能夠通過分析原因自己解決。 在本次實習中我收獲頗多。首先,在本次實習中我復習了該學期《遙感導論》課程中學過的理論知識,并且將本來分散的知識點融會貫通并付諸于實踐,加深了對知識的理解,鍛煉了理論聯(lián)系實際的能力,其次,我初步學會了圖像處理軟件ERDAS IMAGING的基本操作,學習遙感難免要處理圖像,正是這次課程設計讓我全

100、面接觸了圖像處理軟件,了解了一些基本功能,對未來的工作和進一步深造打下了基礎。最后,由于是使用ERDAS軟件時間不長,實習過程中難免遇到問題，例如對它全英文的界面很不熟悉,稍微沒注意老師的講解,就不知所措但最終都通過各種辦法解決了，這使我分析問題和解決問題的能力得以提高，也促進了同學之間的交流與合作。 本次課程設計的前三天是在學校的機房里度過的，在此期間進行了對具體工作的模擬，我們第一次完成了圖像預處理、幾何糾正、圖像裁剪、圖像監(jiān)督分類

101、、專題信息提取和最后的專題圖制作等操作。該過程讓我們了解了用遙感影像制作專題圖的作業(yè)流程及各步驟的方法與可能出現(xiàn)的問</p><p><b>　　六、參考文獻 </b></p><p>　　[1] 《ERDAS IMAGINE 遙感圖像處理方法》 .黨安榮、王曉棟、陳曉峰、張建寶. 《清華大學出版社》 2003年4月。 </p><p>