接着上一次的前四个步骤→波段合成，去除黑边并提取土地利用数据（第七届GIS技能大赛解题）上
我们这次来做后面的三个步骤（文末可获取练习数据）
首先我们回顾一下题目要求：数据处理（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年变为建筑用地的区域，按属性选择，然后导出为“耕林至建筑”：
- END -
波段合成，去除黑边并提取土地利用数据....下