波段合成，去除黑边并制作土地利用转移矩阵（附练习数据下载）

这次呢，我们来做一做第七届全国大学生GIS技能大赛试题A（上午），因为这次的大赛其实有蛮多套试题，所以我们可能会多写一点类似的文章哦~~试题A告诉我们：土地利用每年都在发生变化，在统计一些大区域的土地利用情况时，我们可以使用遥感数据来进行粗略的估计。根据以下要求，研究2010年到2015年土地利用的变化情况。

首先我们看到题目要求：数据处理（100分）

1.对遥感影像进行预处理：波段合成，去除黑边，并提取出与2015年土地利用数据（class2015.tif）相同的影像区域。（15分）

2.使用最大似然法监督分类对影像分类，分为5种地物类别：林地，耕地，裸地，建筑物等不透水面，水（请使用提供的矢量数据“训练样本2010”作为训练样本）。（15分）

3.对分类后的结果的边界进行一定的平滑处理，并消除如下图所示像元个数小于5个像元的独立区域，用周围的像元值来代替。（15分）
4.使用GPS采集了一些点并生成了表格（点位置.txt），需要知道这些点对应于class2015.tif数据中的地物类别。（10分）

5.提取出2010年是林地/耕地，并且2015年变为建筑用地的区域。（10分）

6.计算2010年到2015年土地利用转移矩阵。（15分）

7.在地图上展现出2010年到2015年用地类型是如何变化的，制图并保存为jpg格式。（20分）

那么由于步骤较多，这次我们先完成前4个步骤，下一篇我们完成生于3个步骤。

首先给我们的数据有：（文末有练习数据下载链接）

1. lt51190382010144bjc00文件夹：2010年的影像数据存放在此文件夹中。

2. class2015.tif：2015年的土地利用结果数据。

3. 训练样本2010.shp：对2010年影像执行最大似然分类法所使用的训练样本数据。

4. 点位置.txt：GPS采集的点生成的表格。

开始操作

01

首先我们把数据都丢进ArcGIS中查看，第一步，我们需要在ArcGIS中进行波段合成：