《土地利用矢量、柵格轉(zhuǎn)移矩陣計(jì)算.doc》由會(huì)員分享，可在線閱讀，更多相關(guān)《土地利用矢量、柵格轉(zhuǎn)移矩陣計(jì)算.doc（30頁珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

最簡(jiǎn)單的土地利用/景觀轉(zhuǎn)移矩陣的制作（轉(zhuǎn)）(2009-05-19 03:36:24) 標(biāo)簽：雜談? 分類：相關(guān)資料 http://user.qzone.qq.com/124793909 ? 寫文章的過程中，同樣也用到了土地利用/景觀轉(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ì)算），采用Arcmap打開上述兩期文件。 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ù)類型為double C、右鍵點(diǎn)擊剛剛產(chǎn)生的“newarea”，并選擇“calculate values” D、然后點(diǎn)擊“field calculator”對(duì)話框里面的“Advanced”后選擇“help”將 ?????? Dim Output as double Dim pArea as Iarea Set pArea = [shape] Output = pArea.area 拷入到“field calculator”對(duì)話框下面的空白處 E、在對(duì)話框“field calculator”最下面的空白處填上“output” F、在屬性表點(diǎn)擊“options”按鈕，→export（導(dǎo)出），屬性以.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.????分別打開剛剛產(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.????點(diǎn)擊“options”按鈕，選擇“add field”，然后給出一個(gè)新name“newarea”，數(shù)據(jù)類型為double F.????右鍵點(diǎn)擊剛剛產(chǎn)生的“newarea”，并選擇“calculate values” G.????然后點(diǎn)擊“field calculator”對(duì)話框里面的“Advanced”后選擇“help”將??????????Dim Output as double Dim pArea as Iarea Set pArea = [shape] Output = 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_1＝1 點(diǎn)擊apply，此時(shí)會(huì)統(tǒng)計(jì)出該地區(qū)其他年份年仍然為水域的多邊形 C.右鍵點(diǎn)擊“newarea”選擇統(tǒng)計(jì)功能（“statistics”），就會(huì)得到從然為水域的面積 D.分別計(jì)算出GRID_COD_1＝2（林地）、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、??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(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ì)提示你的空間不足 做這一步之前，請(qǐng)做好前期的地理編碼。 2 ArcView3.3加載 spatial analysis模塊 把兩時(shí)相的Vector圖轉(zhuǎn)成grid格式（當(dāng)然中間有一些單位的設(shè)置根據(jù)你做的圖的分辨率來設(shè)置即可）analysis---mapcaculate 直接計(jì)算即可。 3 把兩期解譯完的Vector文件在arctoolbox——overlay——union中疊加，注意：兩個(gè)文件不能用同一個(gè)字段名，比如一個(gè)用93Type，另一個(gè)時(shí)相則用00Type 疊加后的文件在Arcmap中打開，選中文件，然后點(diǎn)右鍵——Property——空間查詢，輸入條件語句，比如：93Type=‘1’And 00Type=‘2’；查詢結(jié)果即為第一種類型轉(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 area→ 2.導(dǎo)入兩期的矢量或柵格圖→,你就點(diǎn)OK,讓arcgis運(yùn)算即可。 3.結(jié)果顯示，在arcgis里導(dǎo)出*.dbf→（這就是結(jié)果了）。導(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)移矩陣的方法。 現(xiàn)有2005、2008兩期土地利用現(xiàn)狀數(shù)據(jù)： 1.jpg 2.jpg 3.jpg 4.jpg 在求轉(zhuǎn)移矩陣前先做好下列工作： 1、備份數(shù)據(jù)。 2、分別將兩期數(shù)據(jù)用dissolve工具按各自的地類編碼合并要素，注意選中最下面的“create multipat features”選項(xiàng) 8.jpg 處理好數(shù)據(jù)后，用intersect工具求兩期數(shù)據(jù)的交集 0.jpg 00.jpg 交集的結(jié)果： 5.jpg 6.jpg 將屬性表增加面積字段，并計(jì)算面積 9.jpg 10.jpg 11.jpg 結(jié)果: 12.jpg 從上圖中已經(jīng)可以看出土地轉(zhuǎn)移的情況：dlbm2005是2005年時(shí)該塊地的地類代碼，dlbm2008是2008年時(shí)該地塊的地類代碼，最后一列是相應(yīng)地塊的面積。 接下來將屬性表導(dǎo)出到dbf文件，然后在excel中打開，執(zhí)行數(shù)據(jù)菜單中的“數(shù)據(jù)透視表和數(shù)據(jù)透視圖”命令 13.jpg 14.jpg 15.jpg 16.jpg 點(diǎn)完成，然后將dlbm2005作為行字段，dlbm2008作為列字段，area作為數(shù)據(jù)拖到相應(yīng)的區(qū)域 17.jpg 生成的數(shù)據(jù)透視表 18.jpg 將數(shù)據(jù)復(fù)制到新的工作表中，適當(dāng)?shù)男拚蠹瓷闪送恋乩棉D(zhuǎn)移矩陣。