【ENVI实验课】02利用ENVI进行农业耕作用地变化监测
以下内容仅提供简单的操作指导及大致思路流程，学术大牛可以当段子看一看!!!
01软件与数据
02主要内容
1、数据预处理：影像裁剪；
2、绘制农耕地边界，利用监督分类方法进行农业耕作用地的分类；3、用ENVI自带的动态监测工具提取农业耕作用地的变化信息。
03实验步骤
◐ 1. 影像裁剪
打开ENVI点击File->Open As-> Landsat-> GeoTIFF with Metadata在对话框中选择两景某城市landsat影像数据的MTL.txt文件和行政界线矢量文件。选择Toolbox->Regions of Interest->SubsetData from ROIs分别对两景图像进行裁剪。
在File→New→Vector Layer，在 Create New Vector Layer 面板中设置：图层名称：Mask；矢量类型：Polygon；源数据：QUAC15。在图像上沿着农业用地的边缘绘制多边形矢量作为掩膜文件。绘制好后单击右键选择Accept菜单。File→SaveAs，将绘制的 Mask 文件保存为外部Shapefile文件。
（1）选择Toolbox-> Classification->ClassificationWorkflow，输入栅格文件QUAC15；切换InputMask选项，选择绘制的农用地边界矢量文件作为掩膜，单击Next。
（2）选择Use Training Data进行监督分类，分为五类样本：绿植耕地、留茬耕地、休耕地/ 裸地、水体和居民用地。在 Algorithm选项中默认选择 Maximum Likelihood 方法。
（3）在 Cleanup 步骤中，设置 Smoothing 参数可以去除椒盐等噪声，设置Aggregation参数可以将小斑块合并到大块中。按照默认设置单击 Next 按钮。
（4）在输出结果（Export）步骤中，在 Export Files 面板，选择分类图像输出路径和文件名。同样的方法将另一幅时相影像进行分类得到分类结果。
◐
4. 农业耕作用地变化信息提取（1）在 Toolbox 工具箱中，双击 Change Detection/Thematic Change Workflow 工具，在 File Selection 步骤中，分别为Time1 “2003wuxi_class.dat” 和Time2 “2004wuxi_class.dat”。选择 Input Mask 选项，在theThematic Change Mask 对话框中,选择 mask.shp文件作为掩膜文件，单击 Next 按钮。
（2）勾选OnlyInclude Areas That Have Changed选项，只获得变化的区域。
（3）在Cleanup步骤中，设置smoothing参数去除小斑块，单击Next按钮。
（4）将结果输出图像和统计报表，单击Finish按钮输出结果。
（5）在变化栅格图像结果的 Classes 图层上右键选择Hide All Classes 菜单，勾选：From‘绿植耕地’to ‘留茬耕地’和From‘绿植耕地’to ‘休耕地/裸地’。
（6）分别右键选择Class Statistics 统计，利用像元数*30*30即可计算出变化的面积。其中有134.7505平方千米的绿植耕地变为留茬耕地，有17.4402平方千米的绿植耕地变为休耕地/裸地。同样，有71.283平方千米由留茬耕地变为绿植耕地，有6.9399平方千米由休耕地/裸地到绿植耕地
Tips
1. 本次操练是对农用地变化情况进行研究，由于选择农用地区域时具有主观性，边界不准确以及监督分类时不能准确区分绿植耕地、留茬耕地和裸地导致结果的不准确。
2. 保证结果的精度需要有一定的先验知识，可以根据已有的土地利用类型判断出一个区域的农用地边界，再通过相关资料查询可以得到该时间的农用地类型，使分类结果更精确可信，就可以得到可信的农用地变化情况。
作者｜不许人间见白头
排版｜Moon
校阅｜数读菌
那今天就到这里结束啦，欢迎留言讨论。文中的图片文字未经许可不要随便“引用”。
如果可以的话，希望能够转发分享，点个在看并且点个赞，给个赞赏~~也欢迎规范转载~
也希望大家和我多留言互动啊！（据说这样可以增加我的推送在你的订阅号里出现的概率）
需要你的“分享”和“在看”
END>
如需全文转载文章、投稿或者合作
可添加微信
(回复超慢！！！)
（不要添加我问各种问题，我大概率不会的==）
（入群请一定要备注入群）
（添加后会在晚上非工作时间通过，请稍安勿躁）
公众号
微博