ArcGIS土地利用轉(zhuǎn)移矩陣

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

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

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

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

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

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