赣南地区1991年-2017年土地利用变化检测分析

第一部分:研究区概况

  • 地理位置:江西省赣州市东南部、广东省梅州市北部以及福建省龙岩市西部地区;
  • 气候条件:赣南地处中亚热带,呈典型的亚热带季风性湿润气候,此地四季分明;
  • 地形地貌特点:地形复杂,以山地、丘陵、盆地为主,地域差异大。

【erdas做土地利用变化检测分析】_变化检测

第二部分:数据选择与下载

我选择的是Landsat的数据的条带号是【121043】,在“地理空间数据云”网站上下载了1991年3月和2017年11月的数据

第三部分:数据处理

1.数据预处理

1.1合成波段

【Raster】→【Layer Selection and Stacking】→选则输入影像及其波段(不选BQA波段),然后定义输出位置以及文件名。Landsat 8 只要添加前七个波段即可

【erdas做土地利用变化检测分析】_变化检测_02

1.2裁剪影像

由于合成后的影像存在瑕疵,且两个年份的影像并不是完全统一的,所以需要对它们进行裁剪,使其保持统一;

下图为裁剪之前的影像:

【erdas做土地利用变化检测分析】_变化检测_03

使用【Drawing】工具,勾画出需要用到的影像,然后使用“subset”工具导出这部分的影像;

【erdas做土地利用变化检测分析】_gis_04


【erdas做土地利用变化检测分析】_gis_05

下图为裁剪之后的影像:

【erdas做土地利用变化检测分析】_聚类_06

1.3统一影像的空间分辨率和坐标系

由于裁剪出来的影像分辨率不一致,所以要统一影像分辨率,使用“Resample”工具,导入需要更改分辨率的影像,定义好输出的像素大小,忽略零值。

【erdas做土地利用变化检测分析】_gis_07

2.监督分类

2.1监督分类规则:

1.水体/养殖区(蓝色)
2.林地(墨绿)
3.水田(浅蓝)
4.旱地(黄色)
5.城区(红色)
6.草洲/湿地(Caozhou (Wetland))
7.湖滩(Beach)(橙色)
8.开发区/裸地(棕色)
9.草甸(meadow)
10.不可利用(Not available)
11.云 (cloud)(白色)

监督分类需要三个统一:

  • a)控制范围统一(用一样的AOI文件裁剪)
  • b)坐标系和分辨率统一
  • c)编码统一(分类模板内容一致)

2.2监督分类步骤

2.2.1添加裁剪后的遥感影像
2.2.2对添加后的遥感影像进行监督分类

首先用【Signature Editor】选取控制点,然后对同类的点进行融合

【erdas做土地利用变化检测分析】_变化检测_08

【erdas做土地利用变化检测分析】_gis_09

【erdas做土地利用变化检测分析】_数据_10

融合结果如下:

【erdas做土地利用变化检测分析】_变化检测_11


【erdas做土地利用变化检测分析】_聚类_12

2.2.3监督分类

【erdas做土地利用变化检测分析】_变化检测_13

2.2.3.1重编码

在找完点之后,右键图像,查看Attribute Table,会看到数百个编码,因此要对影像进行重编码,使用【Recode】工具,修改相应的New Value,然后点击“Apply”,关闭窗口,重新添加影像,使用模型构建器,输入上一步重编码之后的影像,然后定义输出(输出的时候选择“Thematic”),Run模型,完成之后再放进来,这时候打开影像的Attribute Table,就发现编码已经重编好了。

【erdas做土地利用变化检测分析】_变化检测_14


【erdas做土地利用变化检测分析】_变化检测_15


【erdas做土地利用变化检测分析】_gis_16


【erdas做土地利用变化检测分析】_变化检测_17


【erdas做土地利用变化检测分析】_聚类_18


【erdas做土地利用变化检测分析】_变化检测_19

输出影像:

【erdas做土地利用变化检测分析】_数据_20

2.2.3.2 聚类和剔除

由于监督分类后的影像中,很多细碎的地类得以保留,影响后期处理的精度,将8个像元相邻的聚类到一起,小于8个像元的自动被归为其附近最大的地类中。然后使用剔除,设置想要保留的像元大小的阈值,我设置的阈值为10,即小于10个像元的地类将被剔除,并入旁边最大的地类中去。

【erdas做土地利用变化检测分析】_数据_21


【erdas做土地利用变化检测分析】_聚类_22

2.2.3.3提取地物并建立缓冲区

使用模型提取地物:

【erdas做土地利用变化检测分析】_聚类_23


【erdas做土地利用变化检测分析】_聚类_24

使用模型建立缓冲区:

【erdas做土地利用变化检测分析】_数据_25

【erdas做土地利用变化检测分析】_变化检测_26

水体缓冲结果:

【erdas做土地利用变化检测分析】_gis_27

2.2.3.4使用AOI修改某些地类

新建AOI

【erdas做土地利用变化检测分析】_数据_28


【erdas做土地利用变化检测分析】_变化检测_29


【erdas做土地利用变化检测分析】_变化检测_30

【erdas做土地利用变化检测分析】_聚类_31


【erdas做土地利用变化检测分析】_聚类_32

【erdas做土地利用变化检测分析】_数据_33

输出的时候选择专题图形式

【erdas做土地利用变化检测分析】_变化检测_34

第四部分:精度评估

1.精度评估步骤:

选择【Accuracy Assessment】,选择要进行评价的点,每类最少30个点。然后选择需要评估的类别,没有分的类别就不用选择。之后在Reference里输入类别,链接裁剪后的影像和监督分类之后的影像,选择Show All,点变黄了即验证成功。

【erdas做土地利用变化检测分析】_变化检测_35

【erdas做土地利用变化检测分析】_gis_36


【erdas做土地利用变化检测分析】_数据_37


【erdas做土地利用变化检测分析】_gis_38


【erdas做土地利用变化检测分析】_数据_39


【erdas做土地利用变化检测分析】_变化检测_40


【erdas做土地利用变化检测分析】_数据_41

2.精度评估结果

【erdas做土地利用变化检测分析】_聚类_42

【erdas做土地利用变化检测分析】_gis_43

第五部分:计算变化检测结果

【erdas做土地利用变化检测分析】_数据_44

【erdas做土地利用变化检测分析】_聚类_45