《土地利用矢量、柵格轉移矩陣計算.doc》由會員分享，可在線閱讀，更多相關《土地利用矢量、柵格轉移矩陣計算.doc（30頁珍藏版）》請在裝配圖網(wǎng)上搜索。

最簡單的土地利用/景觀轉移矩陣的制作（轉）(2009-05-19 03:36:24) 標簽：雜談? 分類：相關資料 http://user.qzone.qq.com/124793909 ? 寫文章的過程中，同樣也用到了土地利用/景觀轉移矩陣。轉移矩陣的制作很多軟件可以實現(xiàn)，有Arcview、Arcmap、還有人用matlab。但是網(wǎng)上都寫的都不詳細，在這里，本文才用arcmap的疊加操作＋excel的數(shù)據(jù)統(tǒng)計功能，可以輕松完成這一過程。相信看完了本文，初學者或者對GIS了解不多的人，也能輕松的制作出轉移矩陣。 方法一： 1 數(shù)據(jù)準備 ?? 準備好兩期的coverage或shape文件。（注意：拓撲關系要建好，不然無法計算），采用Arcmap打開上述兩期文件。 2 疊加操作 選擇Arcmap里面的“ArcToolbox”按鈕下面的“Analysis Tools”工具下面的疊加分析模塊（Overlay）下面的交叉分析功能（Intersect）。選擇已經打開的兩期數(shù)據(jù)進行疊加分析。疊加的結果存在一個可以找到的地方。同時，把疊加后的結果添加在Arcmap里。 3 輸出疊加文件的屬性數(shù)據(jù) ?? A、右鍵打開intersect產生的矢量文件的數(shù)據(jù)屬性表（open attribute table）。 B、點擊“options”按鈕，選擇“add field”，然后給出一個新name“newarea”，數(shù)據(jù)類型為double C、右鍵點擊剛剛產生的“newarea”，并選擇“calculate values” D、然后點擊“field calculator”對話框里面的“Advanced”后選擇“help”將 ?????? Dim Output as double Dim pArea as Iarea Set pArea = [shape] Output = pArea.area 拷入到“field calculator”對話框下面的空白處 E、在對話框“field calculator”最下面的空白處填上“output” F、在屬性表點擊“options”按鈕，→export（導出），屬性以.dbf格式結果存儲。 4??Excel進行轉移矩陣制作 A、剛才存儲的.dbf文件可以使用excel打開，打開的結果重新保存為.xls（excel）文件。 （dbf文件不能保存） ??B、重新打開剛存儲的excel文件。 ??C、選中所有數(shù)據(jù)，選擇數(shù)據(jù)\數(shù)據(jù)透視表和數(shù)據(jù)透視圖。將兩期數(shù)據(jù)的id值分別拖入行列字段，然后以‘newarea’字段作為數(shù)據(jù)項拖入計算區(qū)域。得到的便是轉移矩陣。如果發(fā)現(xiàn)矩陣的形式不美觀，可以對單位進行修改。 方法二：作某一地物與其他不同年份地類 例如： 1、做水體和其它年份分類地物的交集 A.????分別打開剛剛產生的水域多邊形和該地區(qū)其他年份的分類矢量文件 B.????選擇gis里面的“AacToolbox”按鈕下面的“Analysis Tools”工具下面的疊加分析模塊（Overlay）下面的交叉分析功能（Intersect） C.????分別輸入上面抽出的水域多邊形和該地區(qū)其他年份區(qū)的分類后的多邊形進行intersect分析 D.????右鍵打開intersect產生的矢量文件的數(shù)據(jù)屬性表（open attribute table） E.????點擊“options”按鈕，選擇“add field”，然后給出一個新name“newarea”，數(shù)據(jù)類型為double F.????右鍵點擊剛剛產生的“newarea”，并選擇“calculate values” G.????然后點擊“field calculator”對話框里面的“Advanced”后選擇“help”將??????????Dim Output as double Dim pArea as Iarea Set pArea = [shape] Output = pArea.area拷入到“field calculator”對話框下面的空白處 H.在對話框“field calculator”最下面填上“output” 2、然后分別計算水體在該地區(qū)其他年份年轉入和轉出情況 A.點擊“potions”下面的“select by attribute”按鈕，并選擇“GRID_COD_1” B.令GRID_COD_1＝1 點擊apply，此時會統(tǒng)計出該地區(qū)其他年份年仍然為水域的多邊形 C.右鍵點擊“newarea”選擇統(tǒng)計功能（“statistics”），就會得到從然為水域的面積 D.分別計算出GRID_COD_1＝2（林地）、3（草地）、4（城鎮(zhèn)點）、5（耕地）、6（灘涂地）的面積就得到從由水域轉入的面積。 以此類推 網(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、??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),(J>I) 其中：NC(I,J)表示i，j 兩年份的土地利用變化圖；NC（i）表示i年份遙感分類影像；NC（j）表示j年份的遙感分類影像。 在此計算的基礎上，將以上變化影像圖轉化為BIL格式，再利用ARC/INFO GRID模塊將影像轉為GRID格式，然后利用GRID模塊中的屬性表（vat）查看命令對影像灰度值進行統(tǒng)計，最后得出土地利用轉化舉證。（注:此方法本人尚未實現(xiàn)過，不知可行否）。 若數(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文件在arctoolbox——overlay——union中疊加，注意：兩個文件不能用同一個字段名，比如一個用93Type，另一個時相則用00Type 疊加后的文件在Arcmap中打開，選中文件，然后點右鍵——Property——空間查詢，輸入條件語句，比如：93Type=‘1’And 00Type=‘2’；查詢結果即為第一種類型轉化為第二種類型的圖形，可以另建一圖層比如：12，把查詢結果復制到12圖層上。統(tǒng)計出面積，依進行，就可以得到土地利用類型轉移矩陣。  剛才看了緣老大發(fā)的帖子，題目是：利用arcigis生成土地利用轉移矩陣。方法流程講的很詳細。最近我也用了好多方法來做，其中個人認為最簡便的方法是利用arcgis的spatial analyst tools---zonal---tabulate area 詳細步驟請見------方法很簡單。另提醒（矢量，柵格都可以哦 ） 1.找到zonal--tebulate area→ 2.導入兩期的矢量或柵格圖→,你就點OK,讓arcgis運算即可。 3.結果顯示，在arcgis里導出*.dbf→（這就是結果了）。導出來你用excel打開，想插到那里就插到那里。  利用ARCGIS生成土地利用轉移矩陣 上一篇 / 下一篇 ?2009-02-01 19:42:58 查看( 825 ) / 評論( 62 ) 轉移矩陣在土地利用變更調查時經常用到。arcgis里沒有直接的工具來計算轉移矩陣，但我們可以利用已有的工具間接來求。下面介紹一下利用arcgis求轉移矩陣的方法。 現(xiàn)有2005、2008兩期土地利用現(xiàn)狀數(shù)據(jù)： 1.jpg 2.jpg 3.jpg 4.jpg 在求轉移矩陣前先做好下列工作： 1、備份數(shù)據(jù)。 2、分別將兩期數(shù)據(jù)用dissolve工具按各自的地類編碼合并要素，注意選中最下面的“create multipat features”選項 8.jpg 處理好數(shù)據(jù)后，用intersect工具求兩期數(shù)據(jù)的交集 0.jpg 00.jpg 交集的結果： 5.jpg 6.jpg 將屬性表增加面積字段，并計算面積 9.jpg 10.jpg 11.jpg 結果: 12.jpg 從上圖中已經可以看出土地轉移的情況：dlbm2005是2005年時該塊地的地類代碼，dlbm2008是2008年時該地塊的地類代碼，最后一列是相應地塊的面積。 接下來將屬性表導出到dbf文件，然后在excel中打開，執(zhí)行數(shù)據(jù)菜單中的“數(shù)據(jù)透視表和數(shù)據(jù)透視圖”命令 13.jpg 14.jpg 15.jpg 16.jpg 點完成，然后將dlbm2005作為行字段，dlbm2008作為列字段，area作為數(shù)據(jù)拖到相應的區(qū)域 17.jpg 生成的數(shù)據(jù)透視表 18.jpg 將數(shù)據(jù)復制到新的工作表中，適當?shù)男拚蠹瓷闪送恋乩棉D移矩陣。