本篇主要是針對矢量數(shù)據(jù)的分析。

一、不同時期的土地利用矢量數(shù)據(jù)，如何分析其圖形及屬性變化？

土地利用圖（左圖為1993年，右圖為2003年）

思路如下：

可以先對2個圖層進行Union操作，然后在結果中新建一個字段，分析屬性的變化。如果屬性沒變，就賦值為1，如果屬性變化了就賦值為0。

操作步驟如下：

1、將兩個圖層的屬性放到一個圖層里。

使用GP工具Union（聯(lián)合）。注意JoinAttribute，選擇All，也就是結果中包含兩個圖層的所有字段。

結果屬性表如下：

Unoin可以將數(shù)據(jù)的屬性合并到一起，同時，圖斑不一樣的區(qū)域也可以拆分開。更多干貨敬請關注：GIS前沿  見下圖：

2、新建一個短整型字段。

3、使用字段計算器，比較兩個字段的內容。如果屬性沒變，就賦值為1，如果屬性變化了就賦值為0。

Python代碼如下：

變化中寫：（比較的是Class_Name和Class_Name1兩個字段）

isSame( !Class_Name! , !Class_Name_1! )

如果想知道變化的圖斑是從什么類型轉為什么類型，可以新建一個文本型字段，將土地利用類型的兩個字段賦值過去。

結果如下：

二、如何統(tǒng)計不同土地利用類型的面積？

使用GP工具Summary Statistics（匯總統(tǒng)計數(shù)據(jù)），統(tǒng)計各類的總面積。

結果是一張統(tǒng)計表，如下：

三、如何統(tǒng)計不同區(qū)域各類土地利用類型的面積？

可以使用GP工具Tabulate Intersection（交集制表）進行統(tǒng)計。

結果如下：

四、如何將行政區(qū)劃的信息追加到土地利用數(shù)據(jù)中？

有一份行政區(qū)劃矢量數(shù)據(jù)，希望知道每個圖斑隸屬于哪個行政區(qū)劃。數(shù)據(jù)可以是行政區(qū)域，也可以是坡度及其他數(shù)據(jù)。

一個圖斑可能會跨區(qū)域，如果圖斑需要切分，一部分屬于行政區(qū)劃A，另一部分屬于B，可以使用GP工具Intersect（相交），參數(shù)中的JoinAttribute，選擇All。結果如下：

如果不希望圖斑被切分，在屬性表中記錄屬于A、B、C區(qū)，可以使用GP工具Spatial Join（空間連接）。

在區(qū)域名字段上右鍵，合并規(guī)則選擇連接，其他設置默認就可以了。

結果如下：

五、如何制作土地利用變化轉移矩陣？

對于分類后的土地利用柵格數(shù)據(jù)，在ENVI里很容易計算土地利用變化矩陣。針對矢量數(shù)據(jù)，我們可以在ArcGIS中統(tǒng)計數(shù)據(jù)，然后在Excel中制作矩陣。

思路如下：

1、對兩期土地利用數(shù)據(jù)分別進行融合Dissolve（融合）操作，使每個土地利用類型僅有一條記錄；

2、再對融合的結果進行Intersect（相交）操作，分析土地利用變化；

3、轉出為Excel表；

4、在Excel中采用透視表進行統(tǒng)計。

模型如下：

針對上面的操作思路，制作了一個Model模型，轉出結果為Excel。模型內容如下：

模型下載網(wǎng)址，鏈接:https://pan.baidu.com/s/1pKNHDeb 密碼:n9hy

雙擊工具箱中此工具，界面如下：

關于土地利用轉移矩陣模型，適用ArcMap 10.2及以上版本。如何打開模型，編輯模型，見下圖操作。大家可以參考模型中的工具，自己來制作這個模型。

結果為一個Excel表格，制作透視圖方法如下：

1、點擊數(shù)據(jù)透視圖，然后選中所有區(qū)域的數(shù)據(jù)；

2、先勾選第一期的土地利用類型字段；

3、再將第二期的土地利用類型字段拖拽到透視圖的第一行位置；

4、最后將面積字段拖拽到透視圖的中間空白位置；

再修改下字頭內容，添加下邊框，土地利用轉移矩陣就做好了。