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

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

2、 名： </p><p>　　學(xué) 號(hào)： 10253006 </p><p>　　院 系： 城市與環(huán)境 </p><p>　　年 級(jí) ： 10級(jí) </p

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

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

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

6、/p><p>　　Shape File）,使分析的結(jié)果可以直接為地理信息系統(tǒng)管理與應(yīng)用，從而發(fā)揮更大的作用。 </p><p>　　關(guān)鍵詞：徐州市；遙感圖像處理；土地利用</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>　　實(shí)習(xí)目的 4</p>

10、<p>　　實(shí)習(xí)內(nèi)容 4</p><p>　　實(shí)習(xí)要求 4</p><p>　　實(shí)習(xí)步驟 4&

11、lt;/p><p>　　（一）明確目標(biāo) 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>　　（四）監(jiān)督分類 6</p><p>　　1、構(gòu)建分類系統(tǒng)

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

15、 9</p><p>　　5、評(píng)價(jià)分類結(jié)果 10</p><p>　　6、分類后處理 13</p><p>　?。ㄎ澹┚仃嚪治?

16、 13</p><p>　　總結(jié) 14</p><p>　　參考文獻(xiàn) 15</p><p>　　附錄

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

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

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

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

21、;<b>　?。ǘ┦占Y料 </b></p><p><b>　　1. 數(shù)據(jù)源 </b></p><p>　　所需數(shù)據(jù)及詳細(xì)使用情況如下：</p><p>　　將這校正好的兩幅圖分別裁剪相同大小的部分進(jìn)行監(jiān)督分類，然后將個(gè)要素提取出來(lái)分別對(duì)比，比較2000年與1987年的變化之處。</p><p>

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

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

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

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

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

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

28、多項(xiàng)式次方為2。</p><p>　　（2）、定義投影參數(shù)（Projection）.</p><p>　?。?）、單擊Apply按鈕應(yīng)用或單擊Close按鈕關(guān)閉。</p><p>　?。?）、打開GCP Tool Reference Setup對(duì)話框。</p><p>　　第三步：?jiǎn)?dòng)控制點(diǎn)工具</p><p>　?。?/p>

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

46、gnature Columns 對(duì)話框。</p><p>　　在View Signature Columns 對(duì)話框，</p><p>　　1.單擊第一個(gè)字段的Column 列并向下拖動(dòng)鼠標(biāo)直到最后一個(gè)字段，此時(shí)，所有字段都被選擇上，并用黃色（默認(rèn)色）表示出來(lái)。</p><p>　　2.按住Shift鍵的同事分別單擊Red、Green、Blue三個(gè)字段，Red、Gr

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

68、s對(duì)話框上進(jìn)行選擇。</p><p>　　單擊OK（關(guān)閉Attribue Options對(duì)話框）。</p><p>　　返回Supervised Classification 對(duì)話框。</p><p><b>　　5、評(píng)價(jià)分類結(jié)果</b></p><p><b>　　閾值處理</b></p&g

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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