1、土地利用解譯操作流程準備工作：在空間數(shù)據(jù)云平臺（）網(wǎng)站上下載1992年-2012年遙感影像數(shù)據(jù)（注意：為確保鄭州市在圖幅范圍內(nèi)，檢索時候輸入與鄭州市相近的經(jīng)緯度坐標，控制云量大小，確保數(shù)據(jù)時間在5-10月份之間）考慮到下載的數(shù)據(jù)質(zhì)量比較好，數(shù)據(jù)都自帶坐標且云量基本都為0，省去配準和大氣校正步驟，直接加載已知1995年的土地分類矢量圖，并對遙感圖像進行裁剪。1. 由于下載的數(shù)據(jù)都是分波段的，因此先進行波段合成，具體操作如下：打開ENVI4.8，打開File->Open external file->Landsat->GeoTIFF，選擇7個TIF格式文件打開，選擇RGB顯示，勾

2、選742波段，Load RGB，打開Basic Tools->Layer Stacking點擊Import File，輸入要合成的波段，選擇742波段進行合成并輸出一個文件（名為1992）2加載1995年鄭州市土地利用矢量圖，具體操作如下：打開File->Open Vetor Filel瀏覽到要打開的SHP格式文件，點擊打開，點OK。選擇剛剛打開的SHP圖層并加載到之前打開的RGB圖層，如下圖選擇Display #1:得到如下結(jié)果：3.根據(jù)加載的矢量圖對遙感影像進行不規(guī)則裁剪，具體操作如下：打開圖像感興趣區(qū)域（ROI）如下所圖示：點擊File->Subset Data via

3、 ROIs，對圖像進行裁剪輸出文件名為1992裁剪，如圖所示：加載裁剪圖像，顯示裁剪結(jié)果：4.按照同樣步驟，將矢量圖加載到裁剪后圖像上（參考第二步）略5將剪裁圖像另存為.img格式文件，命名為1992裁剪.img。6.將1992裁剪.img和1992土地利用矢量文件加載到ArcMap中，如下圖所示：7.新建Shpfile圖層，定義坐標系（與遙感影像坐標系保持一致），打開編輯和捕捉工具，選擇editor->snapping->snapping toolbar如下圖所示：8.比對1994年與1995年的土地利用變化，將增加變化類型提取矢量化，并編輯屬性表，打開新建圖層屬性表，并添加字段

4、信息，如下圖所示1-3代表該增加部分是由類型1用地轉(zhuǎn)化成類型3用地的。9.增加部分數(shù)字化提取結(jié)果如下圖所示：注意：（1） 在提取增加信息的過程中，由于變化部分圖斑面積普遍較小，矢量化過程中記得打開捕捉工具，避免兩個斑塊離得太近造成拓撲錯誤；（2） 數(shù)字化過程中可以利用查詢工具快速查找原本地類編碼，便于數(shù)字化過程中屬性表變化信息編輯；（3） 將1995年土地利用屬性表中的地類編碼進行合并，如將123、113等全部合并為第1類，將51、52、53等全部合并為第5類，整理屬性表信息，刪除冗余部分，主要屬性信息應包括ID、地類編碼、土地面積、變化信息、變化面積等基本信息。說明：對增加部分進行提取時，是

5、以1995年為基準進行推算提取的，對于1995年以前的年份所提取出來的增加信息實際上來講是返減信息，如1-3類型變化，按照時間順序來看，應該是3-1的變化，而對于1995年后的數(shù)據(jù)則無此問題。10.對圖層zhengzhou1995和B-1994年變化信息圖層做合并操作，得到1994年土地利用圖層打開ArcToolbox->Analysis Tools->Overlay->Union分析工具如上圖所示，將B-1994年變化圖層和圖層zhengzhou1995添加到Input Features中，輸出結(jié)果命名為1994，得到結(jié)果如下圖所示：說明：對于合并后屬性信息，刪除矢量化過程

6、中造成的明顯拓撲錯誤信息，刪除無用字段并添加1994_class新字段信息來表明1994年新增加的地類信息，保留1995-1994年變化信息和1995_class字段信息，方便變化前信息與變化后信息對比，查找變化信息。（注意：前期可以跳過1113步的獲取土地轉(zhuǎn)移矩陣部分，直接做14步）11. 分別對1994年、1995年土地利用變化圖層做數(shù)據(jù)融合操作打開ArcToolbox，選擇Data Management Tools ->Generalization ->Dissolve工具。Input Feature選擇要融合的圖層，Output Feature Class選擇輸出結(jié)果存儲的

7、位置及名稱，Dissolve Field(s)選擇土地利用類型字段（如Type1995），然后勾選Creat multipart features選項，點擊OK完成，如下圖所示：重復此過程，對1995年數(shù)據(jù)也進行融合，查看屬性表和融合后圖層變化如下圖：12. 對融合后數(shù)據(jù)進行疊置分析在ArcMap中打開兩個時相融合后的數(shù)據(jù)，在ArcToolbox中選擇Analysis Tools ->Overlay ->Intersect工具，Input Features選擇兩個時相的圖層，Output Feature Class選擇疊加結(jié)果存儲的位置，命名為A，其余選項可以忽略，單擊【OK】完成

8、，如圖所示：查看A圖層屬性表信息，刪除無用字段，添加面積字段Area，并在屬性表中Area字段上單擊右鍵選擇Calculate Geometry ，在打開的Calculate Geometry對話框中，Property選擇Area，Units選擇要使用的面積單位，單擊【OK】完成圖斑面積計算，屬性信息如圖所示：13. 制作轉(zhuǎn)移矩陣保存ARCGIS中所得結(jié)果，用Excel打開A.dbf文件，并另存為.xls格式文件。在Excel中選中所有數(shù)據(jù)（不要點左上角，只選擇有效數(shù)據(jù)），點擊【插入】選項卡，選擇【數(shù)據(jù)透視表】|【數(shù)據(jù)透視表】，點擊【確定】。在打開的數(shù)據(jù)透視表中按圖示將字段拖入相應區(qū)域。得到結(jié)

9、果如下圖所示：14. 遙感解譯結(jié)果精度驗證（1）對1999_B數(shù)據(jù)進行融合處理打開ArcToolbox，選擇Data Management Tools ->Generalization ->Dissolve工具。Input Feature選擇要融合的圖層1999_B，Output Feature Class選擇輸出結(jié)果存儲的位置及名稱，Dissolve Field(s)選擇土地利用類型字段（如B_class），然后勾選Creat multipart features選項，點擊OK完成，如下圖所示：完成后查看新生成文件1999_B_Dis的屬性表，發(fā)現(xiàn)全部融合為5類，如下圖所示：（2

10、） 在ArcGIS中生成隨機點如下圖所示在工具箱里打開數(shù)據(jù)管理工具->要素類->創(chuàng)建隨機點如下圖填寫輸出路徑、輸出名稱、約束要素類選擇上一步融合處理結(jié)果，長整型改成50，下邊的線性單位可設為0m。（3） 添加XY坐標如下圖所示打開工具箱中的數(shù)據(jù)管理工具->要素->添加XY坐標選擇上一步生成的圖層1999_B_R，并點擊確定。結(jié)果如下圖所示：打開1999_B_R屬性表可以發(fā)現(xiàn)增加了X、Y坐標值如下圖：（4） 結(jié)合谷歌地球進行校驗并調(diào)整地類將上一步得到的1999_B_R.dbf格式文件用Excel方式打開，可以得到每種地類20個隨機點的X、Y坐標值。打開谷歌地球，打開工具-

11、>選項->3D視圖，勾選“通用橫軸墨卡托投影”（5）對比校驗，調(diào)整地類點擊添加地標按鈕分別將不同地類的X、Y坐標如下圖所示輸入Google地球中，并放大每一個標記后點，查看標記后點所處地類是否與之前隨機點地類相一致，如果有變化則在Excel表格中對地類信息進行調(diào)整編輯（無變化不用調(diào)整，且只調(diào)整地類，如1調(diào)為2，坐標值無需調(diào)整）全部調(diào)整完畢，將新變化地類全部進行重新歸類，如之前2地類變?yōu)?地類，則將其直接移至1地類部分，每個地類單獨按如下圖格式粘貼到記事本中，存為.TXT格式：（6）建立感興趣區(qū)域打開前期處理好的遙感影像，通過Overlay->Region of Insterest.打開感興趣區(qū)，并建立5個新的地類，如