利用ENVI绘制土地利用图

一、内容说明

下载武汉市的遥感影像图，任意选择一种分类方法，把影像上的土地利用/覆盖的类型分出来（要求将土地利用/覆盖的类型分为农用地、林地、草地、水体、建筑用地、未利用地6大类），并对分类结果进行评价，分析精度，最后制成一幅完整的土地利用图（包括指北针、图例、比例尺等要素）。

二、实验数据介绍

从地理空间数据云网站（http://www.gscloud.cn/）下载Landsat 8 OLI_TIRS 卫星数据，选择武汉（东经113°41′-115°05′、北纬29°58′-31°22′之间）地区云量较小时的遥感图像。

将下载到的数据压缩包解压，得到如下不同波段的遥感图像。
（OLI数据产品为LIT级，已进行几何校正。）

三、实验过程

具体步骤：
1、进行假彩色合成

查阅资料可知，OLI的第5、4、3波段组合的地物图像丰富、鲜明、层次好，可用于植被分类和水体识别，故选择这三个波段进行假彩色合成。

将待合成的三个波段的图像加载入ENVI中。

将三个波段图像进行假彩色合成，并输出合成后的图像OLI-543。

2、将全色波段与多光谱波段融合
将第1步得到的假彩色合成影像与全色波段影像融合，使处理后的遥感图像既提高了空间分辨率，同时又保留了图像原本的多光谱特征，从而达到图像增强的目的，增强判读的准确性。
将多光谱图像（上述假彩色合成后的图像 ）和全色图像（Band8对应图像）读入ENVI中。

选择两图像进行融合。

此处尝试几次报错，发现是输出路径有中文造成的，修改路径后成功输出。

将融合后的图像保存为tif格式。


3、选取城市周边的部分区域裁剪，作为之后处理的图像。

将裁剪后的图像以tif格式保存。


4、进行监督分类。
选择训练的样本。

计算样本的可分离性。


利用支持向量机进行分类。

支持向量机分类结果如图。

5、分类后处理。
利用Majority/Minority 分析去除小斑块。


可以看到原始分类结果的碎斑归为了背景类别中，更加平滑。
然后，利用ENVI Classic进行局部手动修改。

对于局部漏分、错分的像元手动修改。

保存修改结果。

6、精度分析。

利用混淆矩阵计算工具进行验证，结果如表，准确率还是比较高的。

7、制作土地利用/覆盖图。
利用ENVI Classic进行操作。

依次添加地图要素底图、地图名称、地图指北针、地图比例尺和地图图例。

预览效果如图。

保存为tif格式。

四、结论

得到武汉市局部地区土地利用图如下，分类准确率为97.962%。

进行分类统计，结果如图。

（小声说：选择分类样本真的好费眼睛，傻傻看不清是啥，自己选的有一定误差。然后，手动调整的时候也是真的累，一个个小块着实让人抓狂。虽然我的成图可能感觉不太好看，但我已经很开心啦。）