本篇主要是針對(duì)矢量數(shù)據(jù)的分析。
一、不同時(shí)期的土地利用矢量數(shù)據(jù)，如何分析其圖形及屬性變化？
土地利用圖（左圖為1993年，右圖為2003年）
思路如下：
可以先對(duì)2個(gè)圖層進(jìn)行Union操作，然后在結(jié)果中新建一個(gè)字段，分析屬性的變化。如果屬性沒變，就賦值為1，如果屬性變化了就賦值為0。
操作步驟如下：
1、將兩個(gè)圖層的屬性放到一個(gè)圖層里。
使用GP工具Union（聯(lián)合）。注意JoinAttribute，選擇All，也就是結(jié)果中包含兩個(gè)圖層的所有字段。
結(jié)果屬性表如下：
Unoin可以將數(shù)據(jù)的屬性合并到一起，同時(shí)，圖斑不一樣的區(qū)域也可以拆分開。更多干貨敬請(qǐng)關(guān)注：GIS前沿 見下圖：
2、新建一個(gè)短整型字段。
3、使用字段計(jì)算器，比較兩個(gè)字段的內(nèi)容。如果屬性沒變，就賦值為1，如果屬性變化了就賦值為0。
Python代碼如下：
1 def isSame(x, y):
2 fieldA = str(x)
3 fieldB = str(y)
4 if fieldA == fieldB:
5 return 1
6 else:
7 return 0變化中寫：（比較的是Class_Name和Class_Name1兩個(gè)字段）
isSame( !Class_Name! , !Class_Name_1! )
如果想知道變化的圖斑是從什么類型轉(zhuǎn)為什么類型，可以新建一個(gè)文本型字段，將土地利用類型的兩個(gè)字段賦值過去。
結(jié)果如下：
二、如何統(tǒng)計(jì)不同土地利用類型的面積？
使用GP工具Summary Statistics（匯總統(tǒng)計(jì)數(shù)據(jù)），統(tǒng)計(jì)各類的總面積。
結(jié)果是一張統(tǒng)計(jì)表，如下：
三、如何統(tǒng)計(jì)不同區(qū)域各類土地利用類型的面積？
可以使用GP工具Tabulate Intersection（交集制表）進(jìn)行統(tǒng)計(jì)。
結(jié)果如下：
四、如何將行政區(qū)劃的信息追加到土地利用數(shù)據(jù)中？
有一份行政區(qū)劃矢量數(shù)據(jù)，希望知道每個(gè)圖斑隸屬于哪個(gè)行政區(qū)劃。數(shù)據(jù)可以是行政區(qū)域，也可以是坡度及其他數(shù)據(jù)。
一個(gè)圖斑可能會(huì)跨區(qū)域，如果圖斑需要切分，一部分屬于行政區(qū)劃A，另一部分屬于B，可以使用GP工具Intersect（相交），參數(shù)中的JoinAttribute，選擇All。結(jié)果如下：
如果不希望圖斑被切分，在屬性表中記錄屬于A、B、C區(qū)，可以使用GP工具Spatial Join（空間連接）。
在區(qū)域名字段上右鍵，合并規(guī)則選擇連接，其他設(shè)置默認(rèn)就可以了。
結(jié)果如下：
五、如何制作土地利用變化轉(zhuǎn)移矩陣？
對(duì)于分類后的土地利用柵格數(shù)據(jù)，在ENVI里很容易計(jì)算土地利用變化矩陣。針對(duì)矢量數(shù)據(jù)，我們可以在ArcGIS中統(tǒng)計(jì)數(shù)據(jù)，然后在Excel中制作矩陣。
思路如下：
1、對(duì)兩期土地利用數(shù)據(jù)分別進(jìn)行融合Dissolve（融合）操作，使每個(gè)土地利用類型僅有一條記錄；
2、再對(duì)融合的結(jié)果進(jìn)行Intersect（相交）操作，分析土地利用變化；
3、轉(zhuǎn)出為Excel表；
4、在Excel中采用透視表進(jìn)行統(tǒng)計(jì)。
模型如下：
針對(duì)上面的操作思路，制作了一個(gè)Model模型，轉(zhuǎn)出結(jié)果為Excel。模型內(nèi)容如下：
模型下載網(wǎng)址，鏈接:https://pan.baidu.com/s/1pKNHDeb 密碼:n9hy
雙擊工具箱中此工具，界面如下：
關(guān)于土地利用轉(zhuǎn)移矩陣模型，適用ArcMap 10.2及以上版本。如何打開模型，編輯模型，見下圖操作。大家可以參考模型中的工具，自己來制作這個(gè)模型。
結(jié)果為一個(gè)Excel表格，制作透視圖方法如下：
1、點(diǎn)擊數(shù)據(jù)透視圖，然后選中所有區(qū)域的數(shù)據(jù)；
2、先勾選第一期的土地利用類型字段；
3、再將第二期的土地利用類型字段拖拽到透視圖的第一行位置；
4、最后將面積字段拖拽到透視圖的中間空白位置；
再修改下字頭內(nèi)容，添加下邊框，土地利用轉(zhuǎn)移矩陣就做好了。
聯(lián)系客服
