?國土空間規(guī)劃怎么會這么忙？？？奉勸看到這里還沒進規(guī)劃院的同志

能騰出點時間去做點自己想做的事的光景一去不返，辭職的念頭一天又一天，我說去練畫都說了有一年連筆桿子都沒摸過。偷得時間寫一篇土地利用轉(zhuǎn)移矩陣。

?文章來源地址http://www.zghlxwxcb.cn/news/detail-431982.html

一、土地利用轉(zhuǎn)移矩陣定義：

?

（太長不看系列）土地利用轉(zhuǎn)移矩陣，就是根據(jù)同一地區(qū)不同時相的土地覆蓋現(xiàn)狀的變化關(guān)系，求得一個二維矩陣。通過對得到的轉(zhuǎn)移矩陣進行分析，能夠得到２個時相，不同的地類之間相互轉(zhuǎn)化的情況，它描述了不同的土地利用的類型在不同年份發(fā)生變化的土地類別以及發(fā)生變化的位置和變化面積。不僅能夠反映上述靜態(tài)的固定區(qū)域固定時間的各個地類面積數(shù)據(jù)，還能夠反映更加豐富的初期各個地類的面積轉(zhuǎn)出以及末期各個地類面積的轉(zhuǎn)入情況。從面積入手，反應(yīng)區(qū)域土地利用變化。面積變化首先反映在不同土地利用類型的總量變化上，通過分析土地利用類型的總量變化，可了解土地利用變化總的變化趨勢和土地利用結(jié)構(gòu)的變化。

?

下面表1便是土地利用轉(zhuǎn)移矩陣，簡單講一句即：行T1表示上一時相，列T2表示下一時相。A1至An表示n種不同的土地覆蓋類型。具體的意思解釋就不寫了，放出來也寫的很繁瑣看的暈，可查看《基于轉(zhuǎn)移矩陣的土地利用變化信息挖掘方法探討》這篇論文。 ? ? ?? ? ??? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?

?

二、開始實驗?

?

數(shù)據(jù)準(zhǔn)備：1、自然資源部global30土地覆蓋數(shù)據(jù)2010/2020。官網(wǎng)可以下載到。2、成都市行政區(qū)劃（我就是想回成都了，我覺得以后的例子都會是成都）

?

將土地覆蓋數(shù)據(jù)加載GIS里后，按研究所需區(qū)域邊界裁剪出來，然后創(chuàng)建為這兩個年份的土地柵格數(shù)據(jù)創(chuàng)建屬性表（他們沒有屬性表，點符號系統(tǒng)里的唯一值就可以創(chuàng)建）?？聪聢D創(chuàng)建后出現(xiàn)了10—100的分類，這些代碼代表著不同的土地利用類型。

下表便是GlobeLand30數(shù)據(jù)的分類系統(tǒng)，需要將每個代碼對應(yīng)的地類連接起來，這一步就是屬性表連接，按照代碼連接即可，不再多說。

?

我這里新建了個字段2010，2020年的就建個2020，連接進來上面的分類后把地類名稱計算過來。

?

2010

2020

?

接下來是柵格轉(zhuǎn)面，將兩個年份的柵格數(shù)據(jù)轉(zhuǎn)為矢量數(shù)據(jù)，記得把簡化面前的√去掉。然后上一步不是新建了字段連接了地類名稱了嗎，字段就不要選value，選2010、2020。

?

?

全都轉(zhuǎn)為矢量數(shù)據(jù)后，按照地類名稱【融合】一下，不然數(shù)據(jù)有點大不太好操作，接下來將兩個年份的矢量數(shù)據(jù)進行【相交】操作。

?

?

?

其實這樣子也初步可以看出來不同時間地類之間的轉(zhuǎn)換關(guān)系了。在相交后的屬性表中添加“面積”字段，計算幾何計算出面積，單位用公頃吧。

?

接下來就要在Excle里操作了，【表轉(zhuǎn)Excle】工具將屬性表轉(zhuǎn)出。Excle插入-數(shù)據(jù)透視表，以2010為行，2020為列，面積求和項為值，得到初步的轉(zhuǎn)移矩陣表。出來之后會發(fā)現(xiàn)有些是空值，這些空值代表地類沒有轉(zhuǎn)化。

?

修改一下表格，將空值用0代替，標(biāo)黃的這條對角線意思是沒有變化的面積，比如最后一行2010年的【水體】—2020年【水體】沒有變化的面積為11985.1公頃，2010年是水體，到2020年還是水體的有11985.1公頃。一般是用百分比來表示，就是占總土地面積的比重，如果有的地類面積很小用面積表示就不怎么明顯。

?

三、怎么理解

?

再說一下這個表怎么看。

2010年是列，2020年是行，以【耕地】為例吧。

橫著看耕地一行，表明2010年是耕地，到2020年轉(zhuǎn)化為冰川、草地、灌木、林地..水體的面積為5635、27、10184……1987公頃。

豎著看耕地一列，表明2020年的耕地由2010年的草地、灌木地、林地…..水體轉(zhuǎn)化來的面積為7982、25.8、9713…..6927公頃。

?

我們也可以計算一下新增和減少，如下表

?

還是以【耕地】為例，同上面的解釋一樣，橫看耕地一行，減少了20243公頃，也就是2010年是耕地，到2020年變?yōu)槠渌仡惖拿娣e共20243公頃。豎看耕地一列，由2010年的其他地類變?yōu)?020年耕地的面積共131533公頃，也就是增加了131533公頃。

?

四、制圖

?

最后就是制圖了，只用文字不顯著，圖文并茂嘛。我們轉(zhuǎn)到GIS里，回到【表轉(zhuǎn)Excle】之前，打開屬性表新建個字段【轉(zhuǎn)換】，用字段計算器寫個這樣的公式（引號是英文狀態(tài)下）

然后這樣就出來啦，可以表示出地類之間的轉(zhuǎn)換關(guān)系。

?

之后再符號系統(tǒng)一系列操作調(diào)一下，就不多說了。

2010年-2020年土地覆蓋變化

2010年-2020年土地覆蓋變化

?

每次作圖就很卡，從來沒好好捯飭過，這次說的比較基礎(chǔ)的東西，還可以計算各時段土地利用覆蓋比例、土地利用變化幅度、土地利用動態(tài)度、凈變化等信息，等下次空下來再寫。如有錯誤多多包涵。

?

?

?

END

如果在實驗過程中有什么地方遇到問題可以加我，竭力解答。需要練習(xí)數(shù)據(jù)也請后臺我【一點規(guī)劃】。

?

?

?

?

?

到了這里，關(guān)于【土地利用變化分析】土地利用轉(zhuǎn)移矩陣的文章就介紹完了。如果您還想了解更多內(nèi)容，請在右上角搜索TOY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章，希望大家以后多多支持TOY模板網(wǎng)！