“All roads lead to Rome”
本文以中國2000年和2010年的歐空局300米分辨率LUCC數據為例，運用 ArcGIS + Excel 獲取柵格影響的土地利用轉移矩陣。
分三步(如果已分類完成，則可直接跳過第一步；如果僅想得到轉移矩陣結果，可以跳過第二步，直接利用Tabulate Area工具，依次添加初始年份、終止年份信息，導出的結果中，左側為轉出結果，即初始年份情況，右側為轉入結果，為變化情況)。
1. 重分類土地利用類型。
2. 一個像元的屬性值表征多期土地利用類型信息。
3. 結果導出。
1. 重分類土地利用類型
首先加載兩期柵格數據：
根據需要，將若干類別，划分為幾大類別：
在Excel中新建表格，將原地類編碼和重分類后編碼一一對應，這里分被設置屬性名為【Value】和【RecCode】,保存為【.csv】或者【.txt】后加載進當前工作空間。
利用 【Reclass by Table】工具進行重分類
得到處理后的影像。
2. 利用柵格計算器合並多時段屬性信息。
此處利用柵格計算器合並屬性信息。
思路是：單個像元從2000年的耕地【Value=1】轉變為2010年的建設用地【Value=5】，
則可以用Value2000*10 +Value2010 得到新值為15，（當然，地類數較多時乘100也是沒問題的）
即公式為【"Reclass2000.tif"*10+"Reclass2010.tif"】
打開柵格計算器（可以批處理）
得到結果的屬性即包含了土地利用的變化情況。
打開屬性表，所見如下：
第二步即完成。
3.結果導出。
打開【Tabulate Area】面積制表工具，批處理，得到屬性表。
可以看出，有187.891651的耕地未發生變化，1.546543的1（耕地）轉變為5（建設用地）
然后用【Table to Excel】導出至EXCEL表中即可。
由於導出結果是單列數據，並且其中可能有部分地類轉移未出現，需要首先檢核，並手動補齊空行。（如上圖僅出現55，需要手動補齊51、52、53、54、56、57）。
在此之后可以在Excel中利用【數據透視表】或者Offset函數等調整格式，轉為轉移矩陣，這里不做具體講解。
希望這篇文章能對你有所幫助~