1、一、數(shù)據(jù)準備（圖1）準備兩幅不同時相的土地利用現(xiàn)狀圖（shp格式），每幅圖的屬性表都要有一個表示土地利用類型的字段，并且要使用不同的名稱加以區(qū)分，如Type1995，Type2000。土地利用類型名稱必須統(tǒng)一，并且完整，如都使用“城鎮(zhèn)用地”、“有林地”等。二、數(shù)據(jù)融合（圖2）在ArcMap里分別打開兩個時相的圖層，打開ArcToolbox，選擇Data Management Tools | Generalization | Dissolve工具。Input Feature選擇要融合的圖層，Output Feature Class選擇輸出結果存儲的位置及名稱，Dissolve Field(s)選

2、擇土地利用類型字段（如Type1995），然后勾選Creat multipart features選項，點擊OK完成。重復此過程，對另一時相數(shù)據(jù)進行融合。此步驟使相同利用類型的記錄融合為一個記錄，以提高后面步驟的計算速度。三、疊置分析（圖3）在ArcMap中打開兩個時相融合后的數(shù)據(jù)，在ArcToolbox中選擇Analysis Tools | Overlay | Intersect工具，Input Features選擇兩個時相的圖層，Output Feature Class選擇疊加結果存儲的位置及名稱，其余選項可以忽略，單擊【OK】完成。四、計算面積并導出屬性表（圖4-6）在ArcMap中打開

3、疊加后的圖層數(shù)據(jù)，在該圖層上右鍵打開屬性表，選擇Option | Add field 新建一個字段，命名為NewArea。在Editer工具條中選擇Editer | Start Editing，然后在屬性表中NewArea字段上單擊右鍵選擇Calculate Geometry ，在打開的Calculate Geometry對話框中，Property選擇Area，Units選擇要使用的面積單位，單擊【OK】完成圖斑面積計算。依次選擇Editer | Save Edits / End Editing保存和退出編輯狀態(tài)。在屬性表中選擇Option | Export 將屬性表保存為dbf文件。五、制作

4、轉移矩陣（圖7-10）（以Excel2007為例）在Excel中打開上一步保存的dbf，另存為Excel格式并打開。在Excel中選中所有數(shù)據(jù)（不要點左上角，只選擇有效數(shù)據(jù)），點擊【插入】選項卡，選擇【數(shù)據(jù)透視表】|【數(shù)據(jù)透視表】，點擊【確定】。在打開的數(shù)據(jù)透視表中按圖示將字段拖入相應區(qū)域。Excel自動計算矩陣，將該表稍事整飾就得到美觀的土地利用轉移矩陣。矩陣中r(I, j)就表示i類型向j類型轉移的土地面積，空值表示i類型向j類型沒有轉移。閱讀全文(131) / 評論 / 扔小紙條 / 給曾經(jīng)善良過留言 收藏: QQ書簽 del.icio.us / 訂閱: Google 抓蝦 【分享】利用

5、ARCGIS求算土地轉移矩陣作者: gaochangju (站內(nèi)聯(lián)系TA) 收錄: 2009-06-19 發(fā)布: 2009-06-19-利用ARCGIS生成土地利用轉移矩陣轉移矩陣在土地利用變更調(diào)查時經(jīng)常用到。arcgis里沒有直接的工具來計算轉移矩陣，但我們可以利用已有的工具間接來求。下面介紹一下利用arcgis求轉移矩陣的方法。現(xiàn)有2005、2008兩期土地利用現(xiàn)狀數(shù)據(jù)：在求轉移矩陣前先做好下列工作：1、備份數(shù)據(jù)。2、分別將兩期數(shù)據(jù)用dissolve工具按各自的地類編碼合并要素，注意選中最下面的“create multipat features”選項處理好數(shù)據(jù)后，用intersect工具求

6、兩期數(shù)據(jù)的交集交集的結果：將屬性表增加面積字段，并計算面積結果:從上圖中已經(jīng)可以看出土地轉移的情況：dlbm2005是2005年時該塊地的地類代碼，dlbm2008是2008年時該地塊的地類代碼，最后一列是相應地塊的面積。接下來將屬性表導出到dbf文件，然后在excel中打開，執(zhí)行數(shù)據(jù)菜單中的“數(shù)據(jù)透視表和數(shù)據(jù)透視圖”命令點完成，然后將dlbm2005作為行字段，dlbm2008作為列字段，area作為數(shù)據(jù)拖到相應的區(qū)域生成的數(shù)據(jù)透視表若數(shù)據(jù)是Vector格式1 Erdas Imagine-Interpreter-Gis Analysis-Matrix，輸入兩個時相的Vector數(shù)據(jù)即可 此時注意 輸出柵格大小不應設的太小 要不一運行就會提示你的空間不足做這一步之前，請做好前期的地理編碼。2 ArcView3.3加載 spatial analysis模塊把兩時相的Vector圖轉成grid格式（當然中間有一些單位的設置根據(jù)你做的圖的分辨率來設置即可）analysis-mapcaculate 直接計算即可。3 把兩期解譯完的Vector文件在arctoolboxoverlayunion中疊加，注意：兩個文件不能用同一個字段名，比如一個用93Type，另一個時相則用 00Type疊加后的文件在Arcmap中打開，選中文件，然后點右鍵Property空間查詢