利用GIS软件通过遥感波段文件提取用地类型

一、遥感技术

遥感技术是指从高空或外层空间接收来自地球表层各类地理的电磁波信息，并通过这些信息进行扫描、摄影、传输和处理，对地表各类地物和现象进行远距离控测和识别的综合技术。从高空或外层空间接收来自地球表层各类地理的电磁波信息，并通过这些信息进行扫描、摄影、传输和处理，对地表各类地物和现象进行远距离控测和识别的综合技术。

二、数据获取

城市规划中，我们多采用卫星遥感数据。目前国内权威的卫星遥感数据可以在地理空间数据云平台获得。根据数据质量，分为免费和付费，本教程中采用免费数据进行叙述。

地理空间数据云网站：http://www.gscloud.cn/

完成登录注册后，可在数据资源中查找需要的遥感数据。

三、软件平台

完整的遥感图像处理软件ENVI拥有更全面的数据处理工具，但在支援ICM使用时，我们只希望获取用地类型的数据，所以本教程中采用了ArcGIS进行处理。步骤解释中，截图为英文版，请使用中文版的同学们自行寻找对应的工具。

四、用地解译

1. 数据确认

在地理空间数据云上，选择公开数据。

并选择landset8数据，搜索文件编号。本教程以LC81200382016088LGN00文件为例。

在右侧点击下载。

下载后的数据包括了12个TIF文件和1个TXT文件，其中我们使用后缀为1-7的前7个波段文件。

2. 波段合成
   

打开ArcGIS，并在ArcToolBox中找到Data Management Tools-栅格-栅格处理-合成波段。

找到我们需要合成的前7个TIF文件。

合成完成后，得到完整的栅格文件。

3. 划定区域
   

新建多边形图层，并缩放到规划区域，在研究范围内绘制一个多边形。

在ArcToolBox中找到Spatial Analyst Tools-提取分析-以长方形提取。

在提取工具中选择要提取的栅格文件和提取的多边形区域，并设置保存路径。

提取完成后得到区域的栅格文件。

4. 波段显示
   

我们在开始时说明了选用了前7个波段文件，对于这些波段，我们可以通过选择显示栅格文件的波段来形成不同的样式以方便我们找到对应的用地类型。

在这里我们推荐使用真彩色波段。右键栅格文件，打开属性。

在符号系统中可以通过选择每个频道所对应的波段。真彩色为red-Band4；green-Band3；blue-Band2。

常用的波段显示排序有：

Band_543：标准假彩色合成，用于植被相关的监测，植被显示为红色、越健康越亮
Band_764：假彩色合成，用于城市监测
Band_564：假彩色合成，有效区分陆地和水体
Band_753：假彩色合成，具有良好的大气透射
Band_652：假彩色合成，用于农作物监测，农作物显示为高亮绿、裸地为品红
Band_632：假彩色合成，突出裸露地表的一些景观
Band_571：假彩色合成，有效监测植被和水体

5. 创建特征文件
   

我们在教程中，先设定4中用地类型，分别是“建筑、水体、绿地”。首先，新建一个多边形shp文件。在shp图层的表格窗口中新建class列，并打开ArcMap的编辑模式。

Shp文件准备完成后，我们开始在Map窗口中绘制个类型的点。

选取具有代表性的区域，并在class列中标注用地类型。

绘制完成后，保存编辑内容，推出编辑模式。在ArcToolBox中找到Spatial Analyst Tools-多源分析-创建特征文件

“栅格波段”选择区域选择后的文件，“输入栅格数据或要素样本数据”选择创建的shp文件，并选择class列作为样本列。

完成后，将会在设置的保存路径中创建一个GSG文件。

6. 用地分类识别
   

在ArcToolBox中找到Spatial Analyst Tools-多源分析-最大似然法分类工具。

“栅格波段”选择区域选择后的文件，特征文件选择刚才生成的GSG文件，其余参数设为默认。

完成后，我们能够得到如下图像，各颜色对应用地类型。

7. 转换输出

上述步骤完成后，我们希望导出shp文件提供给ICM使用。

在ArcToolBox中找到Conversion Tools-From Raster-Raster to Polygon

在输入栅格中选择生成的栅格文件，Field选择Classname列，设置输出目录后，完成导出。

导出后我们可以在ICM或ArcMap中导入shp文件查看。