遥感提取土地利用（从下载到ENVI到ACRGIS 【简单粗暴版】）

 

一、影像下载

下载地址：地理空间数据云http://www.gscloud.cn/search，可以用邮箱注册一个账户。

1. 选择影像：数据集根据自己所需要的影像精度选择

   影像选择注意点：

• 可选择的数据集包括：Landsat8或landsat7（SLE-ON）或landsat4/5TM ，这些卫星的服务时间不一样，根据你的需要选择。
• 要优先选择云量小的图幅，尽量小于5，不然云量太大不好解译。
• 多年解译最后每一年的时间接近，比如，我们都选广州的夏秋的影像（因为冬天和春天落叶影响植被的解译）

（4）选中影像，点击下载

2. 打开影像

（1）打开ENVI

打开envi-file-open-后缀为.MTL 的文件

（2）如果打开以后是这样子的，那么在 调整一下拉伸度就会变亮。

3.监督分类

（1）右边Toolbox-classification-classification workflow

（2）Browse（选择要分类影像）

在选择分类方法选择面板中（Classification Type），选择Use Training Data单击Next按钮；

注：No Training Data选项是非监督分类方法，不需要选择分类样本。

在左边Layer manager 中的影像图层上，右键选择Change RGB Bands，选择波段组合信息，突出选择地物的信息。解译城市信息依次序选741波段.

  

（4）在Supervised Classification面板中，添加类别，在右边修改类别名称、颜色：

选择样本，选中urban类，用鼠标在图上目视判断建设用地的区域，即显示紫色的图像区域选择作为样本，右键选择Accept结束一个多边形样本的选择（或双击左键）。重复操作选择几个多边形样本。同理，选择水体、林地、草地、农田的样本。

（6）勾选Preview可以预览分类结果；

r

4.分类后处理

（10）在Cleanup面板中，可以设置平滑阈值和聚类阈值，以去除分析结果中的"小碎块"；

         Enable Smoothing: 主要去除椒盐噪声，如默认3x3，就是在3x3范围内的中心点的像素值会被9个像素内最多像元数的类别代替。

         Enable Aggregation: 主要去除小区域，如默认9，就是小于等于9个像素的区域重新合并到临近的、更大的区域。

（11）这里选择默认参数，勾选Preview预览结果；

（12）单击Next按钮。

  

5.输出结果

（13）export files中只要选择储存为image（栅格图像）；additional export中选择储存输出分类结果的统计文件

（15）单击Finish，输出结果。

 

6.导入arcgis的数据处理

1. 打开分类后结果，将数据转为TIFF格式（File-save file as-TIFF/GEO TIFF）
2. 在arcgis里打开tiff文件
3. 在右边table of contents中选中图层双击，在右边show中选择unique values（唯一值），单击ok

 

PS:若输出的数据有黑色边框，如：则新建一个面状shp，画一个想要保留下的范围，用剪裁工具剪裁。

 

7.图像剪裁

1. 添加广州市的行政边界面，在acrgis中裁剪分类结果（arctoolbox-spatial analyst tools-extract by mask），输出数据

 

下图是之前做完15年的广州的土地利用分类。

 

8.图像镶嵌（研究范围覆盖多张影像,以此功能拼接）

若有多幅影像需要拼接则在arcgis中（arctoolbox-data management tools-raster-raster dataset-mosaic to new raster Mosaic to new raster界面中将需要拼接的tiff文件（上一步剪裁过的）全部导入

Output Location—输出位置

Raster Dataset Name with Extension—xxx.tif(后缀为.tif)

Number Of Bands—1

ok

完成恭喜你完成解译！

 

BY 婷仔&橄榄

EMMMMMM