ENVI遥感实验：农业耕作与城市绿地变化监测-CSDN文库.docx

一、农业耕作变化监测

1.1实验背景

此试验以2008年和2009年获得的LandsatTM数据为数据源，获得两个时相的农业耕作类型。TM数据从美国地质调查局网站下载，经过基于地形的几何校正，经过大气校正的处理，并进行了空间子区裁剪。图像覆盖区域的土地利用类型主要包括水体、沙漠、耕地和城市区域。通过LandsatTM数据提取农业耕作用地分类和农业耕作用地变化信息，变化检测采用分类后比较方法。

1.2实验数据

ag_08_quac.dat

ag_09_quac.dat

1.3技术流程图

图2.1农业耕作变化监测技术流程图

1.4实验步骤

1.4.1LandsatTM农业耕作用地分类

第一步：加载数据了解各时期土地利用情况

打开ENVI，加载数据ag_08_quac.dat、ag_09_quac.dat。在菜单栏中选择Dispaly下的Portal透视工具，移动透视窗口，大致了解两个时期的土地利用情况。

PAGE

图2.2ag_08_quac.dat（左边）、ag_09_quac.dat（右边）

图2.3透视浏览各时期土地利用情况

第二步；绘制掩膜文件

1、在菜单栏中点击文件，选择新建下的VectorLayer，在CreateNewVectorLayer对话框中设置图层名称Mask，矢量数据类型为Polygon，源数据选择ag_08_quac.dat，点击OK。沿着农业用地的边缘绘制多边形矢量作为掩膜文件，绘制好后右击Accept即可，再将其结果另存为Shapefile文件。

图2.4CreateNewVectorLayer参数设置

图2.5Mask.shp数据

第三步：创建训练区样本

1、在图层管理器中，对ag_08_quac.dat数据进行假彩色合成（RGB为432），使植被显示为红色，便于后续解译（这里只介绍ag_08_quac.dat，ag_09_quac.dat同理）。

图2.6假彩色合成显示

2、在右边Toolbox工具箱中选择【ClassificationWorkflow】工具，在弹出的FileSelection对话框中，单击Browse按钮，打开FileSelection对话框，选择“ag_08_quac.dat”文件，回到Classification面板，切换InputMask选项，单击Browse按钮，选择“Mask”文件，单击Next按钮。

图2.7FileSelection参数设置

3、在ClassificationType步骤中，选择UseTrainingData选项，单击Next按钮。

图2.8ClassificationType参数设置

4、在弹出的SupervisedClassification对话框中分别修改以下参数:分类名称(ClassName):绿植耕地、分类颜色(ClassColor):默认、内部填充(FillInterior):Solid。然后在图像上目视解译绘制有绿色农作物的耕地，即显示红色的图像区域。右键选择Accept菜单或者双击鼠标左键结束一个多边形样本的选择。重复操作从而完成训练区样本选取。

图2.9SupervisedClassification参数设置

5、在弹出的Export对话框中设置输出路径及文件名，至此完成ag_08_quac.dat数据的训练样本创建。同理ag_09_quac.Dat。

图2.10设置输出路径

图2.11ag_08_quac.dat分类结果

图2.12SupervisedClassification参数设置

图2.13ag_09_quac.dat分类结果

1.4.2农业耕作用地变化信息提取

第一步：提取农业耕作用地变化信息

1、在Toolbox工具箱中，【ChangeDetection】→【ThematicChangeWorkflow】，选择ThematicChangeWorkflow工具，在FileSelection步骤中，分别为Timel选择“ag_08_classi.dat和为Time2选择ag_09_classi.dat。在InputMask选项卡中，选择“mask.shp”文件作为掩膜文件，单击Next按钮。

图2.14FileSelection参数设置

2、在ThematicChange步骤中，勾选OnlyIncludeAresThatHaveChanged选项，只获得变化的区域。

图2.15ThematicChange参数设置

3、在Cleanup步骤中，设置Smoothing参数去除小斑块，单击Next按钮。最后在ThematicChangeExp