土地利用矢量、柵格轉(zhuǎn)移矩陣計算

1、最簡單的土地利用/景觀轉(zhuǎn)移矩陣的制作（轉(zhuǎn)）(2009-05-19 03:36:24)標簽：雜談  分類：相關(guān)資料 寫文章的過程中，同樣也用到了土地利用/景觀轉(zhuǎn)移矩陣。轉(zhuǎn)移矩陣的制作很多軟件可以實現(xiàn)，有Arcview、Arcmap、還有人用matlab。但是網(wǎng)上都寫的都不詳細，在這里，本文才用arcmap的疊加操作excel的數(shù)據(jù)統(tǒng)計功能，可以輕松完成這一過程。相信看完了本文，初學者或者對GIS了解不多的人，也能輕松的制作出轉(zhuǎn)移矩陣。方法一：1 數(shù)據(jù)準備   準備好兩期的coverage或shape文件。（注意：拓撲關(guān)系要建好，不然無法計算），采用Arcm

2、ap打開上述兩期文件。2 疊加操作選擇Arcmap里面的“ArcToolbox”按鈕下面的“Analysis Tools”工具下面的疊加分析模塊（Overlay）下面的交叉分析功能（Intersect）。選擇已經(jīng)打開的兩期數(shù)據(jù)進行疊加分析。疊加的結(jié)果存在一個可以找到的地方。同時，把疊加后的結(jié)果添加在Arcmap里。3 輸出疊加文件的屬性數(shù)據(jù)   A、右鍵打開intersect產(chǎn)生的矢量文件的數(shù)據(jù)屬性表（open attribute table）。B、點擊“options”按鈕，選擇“add field”，然后給出一個新name“newarea”，數(shù)據(jù)類型為doubleC、右

3、鍵點擊剛剛產(chǎn)生的“newarea”，并選擇“calculate values”D、然后點擊“field calculator”對話框里面的“Advanced”后選擇“help”將       Dim Output as doubleDim pArea as IareaSet pArea = shapeOutput = pArea.area拷入到“field calculator”對話框下面的空白處E、在對話框“field calculator”最下面的空白處填上“output”F、在屬性表點擊“options”按鈕，export（導出

4、），屬性以.dbf格式結(jié)果存儲。4  Excel進行轉(zhuǎn)移矩陣制作A、剛才存儲的.dbf文件可以使用excel打開，打開的結(jié)果重新保存為.xls（excel）文件。（dbf文件不能保存）  B、重新打開剛存儲的excel文件。  C、選中所有數(shù)據(jù)，選擇數(shù)據(jù)數(shù)據(jù)透視表和數(shù)據(jù)透視圖。將兩期數(shù)據(jù)的id值分別拖入行列字段，然后以newarea字段作為數(shù)據(jù)項拖入計算區(qū)域。得到的便是轉(zhuǎn)移矩陣。如果發(fā)現(xiàn)矩陣的形式不美觀，可以對單位進行修改。方法二：作某一地物與其他不同年份地類例如：1、做水體和其它年份分類地物的交集A.  

5、0; 分別打開剛剛產(chǎn)生的水域多邊形和該地區(qū)其他年份的分類矢量文件B.    選擇gis里面的“AacToolbox”按鈕下面的“Analysis Tools”工具下面的疊加分析模塊（Overlay）下面的交叉分析功能（Intersect）C.    分別輸入上面抽出的水域多邊形和該地區(qū)其他年份區(qū)的分類后的多邊形進行intersect分析D.    右鍵打開intersect產(chǎn)生的矢量文件的數(shù)據(jù)屬性表（open attribute table）E. 

6、0;  點擊“options”按鈕，選擇“add field”，然后給出一個新name“newarea”，數(shù)據(jù)類型為doubleF.    右鍵點擊剛剛產(chǎn)生的“newarea”，并選擇“calculate values”G.    然后點擊“field calculator”對話框里面的“Advanced”后選擇“help”將          Dim Output as doubleDim pA

7、rea as IareaSet pArea = shapeOutput = pArea.area拷入到“field calculator”對話框下面的空白處H.在對話框“field calculator”最下面填上“output”2、然后分別計算水體在該地區(qū)其他年份年轉(zhuǎn)入和轉(zhuǎn)出情況A.點擊“potions”下面的“select by attribute”按鈕，并選擇“GRID_COD_1”B.令GRID_COD_11 點擊apply，此時會統(tǒng)計出該地區(qū)其他年份年仍然為水域的多邊形C.右鍵點擊“newarea”選擇統(tǒng)計功能（“statistics”），就會得到從然為水域的面積D.分別計算出GRI

8、D_COD_12（林地）、3（草地）、4（城鎮(zhèn)點）、5（耕地）、6（灘涂地）的面積就得到從由水域轉(zhuǎn)入的面積。以此類推網(wǎng)上有人還采用的方法：1、  Erdas：erdas-main-interpreter-gis analysis-matrix輸入兩個年份的分類圖，然后進行重編碼即可。2、  Arcview：可以作個運算，其中一期的圖*10或者100或者1000（土地利用類型取1位數(shù)時候×10，兩位100),然后再加上另一期的土地利用圖得出一個圖的屬性，肯定12，15，或者1221，或者1531之類的數(shù)據(jù)，應該就可以實現(xiàn)了。3、 

9、60;Matlab：聽說，還沒有親眼見過，據(jù)說可以成功4、ArcGIS：在一期數(shù)據(jù)的基礎上，繪制出其變化的部分（很復雜，工作了量很大）。根據(jù)你的數(shù)據(jù)類型選用不同的數(shù)據(jù)生成方法若你的數(shù)據(jù)是Raster格式：則有如下方法1 Erdas Imagine-Interpreter-Gis Analysis-Matrix，輸入兩個時相的Raster數(shù)據(jù)即可做這一步之前記得先對兩時相的數(shù)據(jù)進行重編碼（nterpreter-Gis Analysis-Recode）一般運行如果出現(xiàn)錯誤肯定是重編碼沒做好，請繼續(xù)查證。2 先在Erdas中利用 Modeler 計算如下公式NC(I,J)=NC(I)*10+NC(J

10、),(J>I)其中：NC(I,J)表示i，j 兩年份的土地利用變化圖；NC（i）表示i年份遙感分類影像；NC（j）表示j年份的遙感分類影像。在此計算的基礎上，將以上變化影像圖轉(zhuǎn)化為BIL格式，再利用ARC/INFO GRID模塊將影像轉(zhuǎn)為GRID格式，然后利用GRID模塊中的屬性表（vat）查看命令對影像灰度值進行統(tǒng)計，最后得出土地利用轉(zhuǎn)化舉證。（注:此方法本人尚未實現(xiàn)過，不知可行否）。若數(shù)據(jù)是Vector格式1 Erdas Imagine-Interpreter-Gis Analysis-Matrix，輸入兩個時相的Vector數(shù)據(jù)即可此時注意輸出柵格大小不應設的太小要不一運行就會提示

11、你的空間不足做這一步之前，請做好前期的地理編碼。2 ArcView3.3加載 spatial analysis模塊把兩時相的Vector圖轉(zhuǎn)成grid格式（當然中間有一些單位的設置根據(jù)你做的圖的分辨率來設置即可）analysis-mapcaculate直接計算即可。3 把兩期解譯完的Vector文件在arctoolboxoverlayunion中疊加，注意：兩個文件不能用同一個字段名，比如一個用93Type，另一個時相則用00Type疊加后的文件在Arcmap中打開，選中文件，然后點右鍵Property空間查詢，輸入條件語句，比如：93Type=1And 00Type=2；查詢結(jié)果即為第一種類

12、型轉(zhuǎn)化為第二種類型的圖形，可以另建一圖層比如：12，把查詢結(jié)果復制到12圖層上。統(tǒng)計出面積，依進行，就可以得到土地利用類型轉(zhuǎn)移矩陣。剛才看了緣老大發(fā)的帖子，題目是：利用arcigis生成土地利用轉(zhuǎn)移矩陣。方法流程講的很詳細。最近我也用了好多方法來做，其中個人認為最簡便的方法是利用arcgis的spatial analyst tools-zonal-tabulate area詳細步驟請見-方法很簡單。另提醒（矢量，柵格都可以哦）1.找到zonal-tebulate area2.導入兩期的矢量或柵格圖,你就點OK,讓arcgis運算即可。3.結(jié)果顯示，在arcgis里導出*.dbf（這就是結(jié)果了）

13、。導出來你用excel打開，想插到那里就插到那里。利用ARCGIS生成土地利用轉(zhuǎn)移矩陣上一篇 / 下一篇  2009-02-01 19:42:58 查看( 825 ) / 評論( 62 ) 轉(zhuǎn)移矩陣在土地利用變更調(diào)查時經(jīng)常用到。arcgis里沒有直接的工具來計算轉(zhuǎn)移矩陣，但我們可以利用已有的工具間接來求。下面介紹一下利用arcgis求轉(zhuǎn)移矩陣的方法?，F(xiàn)有2005、2008兩期土地利用現(xiàn)狀數(shù)據(jù)：1.jpg2.jpg3.jpg4.jpg在求轉(zhuǎn)移矩陣前先做好下列工作：1、備份數(shù)據(jù)。2、分別將兩期數(shù)據(jù)用dissolve工具按各自的地類編碼合并要素，注意選中最下面的“create multipat features”選項8.jpg處理好數(shù)據(jù)后，用intersect工具求兩期數(shù)據(jù)的交集0.jpg00.jpg交集的結(jié)果：5.jpg6.jpg將屬性表增加面積字段，并計算面積9.jpg10.jpg11.jpg結(jié)果:12.jpg從上圖中已經(jīng)可以看出土地轉(zhuǎn)移的情況：dlbm2005是2005年時該塊地的地類代碼，d