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

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

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

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

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

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

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

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

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

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

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

11、示你的空間不足做這一步之前，請(qǐng)做好前期的地理編碼。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空間查詢，輸入條件語(yǔ)句，比如：93Type=1And 00Type=2；查詢結(jié)果即為第一

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

13、果了）。導(dǎo)出來你用excel打開，想插到那里就插到那里。利用ARCGIS生成土地利用轉(zhuǎn)移矩陣上一篇 / 下一篇  2009-02-01 19:42:58 查看( 825 ) / 評(píng)論( 62 ) 轉(zhuǎn)移矩陣在土地利用變更調(diào)查時(shí)經(jīng)常用到。arcgis里沒有直接的工具來計(jì)算轉(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”選項(xiàng)8.jpg處理好數(shù)據(jù)后，用intersect工具求兩期數(shù)據(jù)的交集0.jpg00.jpg交集的結(jié)果：5.jpg6.jpg將屬性表增加面積字段，并計(jì)算面積9.jpg10.jpg11.jpg結(jié)果:12.jpg從上圖中已經(jīng)可以看出土地轉(zhuǎn)移的情況：dlbm2005是2005年時(shí)該塊地的地類代碼，d