題目要求：數據處理（100分）
對遙感影像進行預處理：波段合成，去除黑邊，并提取出與2015年土地利用數據（class2015.tif）相同的影像區域。（15分）
使用最大似然法監督分類對影像分類，分為5種地物類別：林地，耕地，裸地，建筑物等不透水面，水（請使用提供的矢量數據“訓練樣本2010”作為訓練樣本）。（15分）
對分類后的結果的邊界進行一定的平滑處理，并消除如下圖所示像元個數小于5個像元的獨立區域，用周圍的像元值來代替。（15分）
使用GPS采集了一些點并生成了表格（點位置.txt），需要知道這些點對應于class2015.tif數據中的地物類別。（10分）
提取出2010年是林地/耕地，并且2015年變為建筑用地的區域。（10分）
計算2010年到2015年土地利用轉移矩陣。（15分）
在地圖上展現出2010年到2015年用地類型是如何變化的，制圖并保存為jpg格式。（20分）
那么上次我們已經完成了前4個步驟，這一篇我們完成剩余的3個步驟。
還記得上次的步驟嗎？我們先回顧一下，把前4個步驟的結果都再貼出一下： 
那之所以把接下來的567放在一起來做呢，其實是因為大家可以觀察到，這其實相當于一個步驟，那就是土地利用轉移矩陣的制作，那其實這個東西在我們做作業或者做科研中都是經常碰到的比較基礎的工作，這里呢我們就演示一下，也相當于跟大家探討一下土地利用轉移矩陣的制作了，其實這個有好幾種方法，我們先講第一種，那就是在ArcGIS中使用相交工具：
5&6 使用相交前，我們首先需要將柵格轉為矢量，因為矢量工具是針對矢量數據，我們做出來的class2010和初始給我們的class2015都是柵格數據：
依次將class2010和class2015都轉為矢量數據之后，我們就可以進行相交了，但是首先我們需要查看兩個矢量數據的屬性表，確保他們都有一個土地類型的字段，如果沒有，你用12345來代替土地類型也可以，但是一定要自己記清楚哦~
接著我們就來到了土地利用轉移矩陣的制作了，因為當土地利用轉移矩陣制作完成，我們自然可以快速地提取出2010年是林地/耕地，并且2015年變為建筑用地的區域了：
最終我們得到名為LUCC的土地利用轉移矩陣矢量圖層，我們打開屬性表觀察，位于前面的類型是2010類型，位于后面的類型是2015類型：
那么很簡單的，此時我們先提取出第五步要求的2010年是林地/耕地，并且2015年變為建筑用地的區域，按屬性選擇，然后導出為“耕林至建筑”：
如圖所示第五步我們的結果為：
接下來的第六步和第七步其實都是一個步驟，只不過第六步多出一步導出為EXCEL表格然后計算的過程，那么我們首先需要新建一個面積的字段，然后計算幾何：
接著我們就可以轉EXCEL進行計算了：
我們將無用的列刪掉，只剩下三列，然后插入透視表：
然后在右邊進行如下圖左的設置即可，最終就可以得到土地利用的轉移矩陣，由于研究區域較小，所以這里是以畝為單位的哦~~：
7 我們只要制作出一張2010至2015年的土地利用轉移圖就可以了完成任務了：
我們先新建一個字段，然后在字段計算器中我們只需要按下圖右示范即可，這一步呢，是幫助我們在添加圖例的時候更加簡單和方便：
然后我們在符號系統設置里選擇這個字段，添加所有值：
然后我們調至布局試圖，進行出圖即可，不要忘記添加地圖三要素圖例、指北針和比例尺哦~~
最終，我們得到的土地利用轉移矩陣圖就如下所示了如有錯誤請不吝指正
來源：https://blog.csdn.net/qq_43173805/article/details/120925302
成都途遠GIS是一家專業致力于無人機航空攝影測繪、航空數據處理、GIS地理信息系統研發、數字孿生城市制作、數字沙盤模型等業務的創新型科技公司，為您提供一站式地理信息服務。
