遥感ENVI试验报告材料

1、文档目录前言 3一、实验目的 3二、实验容 3三、实验时间 3四、组织人员 31. 专题概述 42. 处理流程介绍 42.1 图像获取 42.2 数据读取和定标 42.3 图像配准 52.4 大气校正 52.5 反演模型构建及模型应用 52.6 植被变化 63. 详细处理过程 73.1 数据预处理 73.1.1 安装环境小卫星数据处理补丁73.1.2 数据处理和定标73.1.3 工程区裁剪93.1.4 图像配准143.1.5 大气校正173.1.6 裁剪浑善达克区233.2 植被覆盖度反演 273.2.1 计算归一化植被指数273.2.2 计算植被覆盖度283.3 植被变化监测 29文档3.3

2、.1植被覆盖区提取293.3.2 植被变化检测313.4 成果后期处理与应用 323.4.1植被变化区域图的背景值处理323.4.2 植被变化区域制图33实验心得 36文档前言、 实验目的1、掌握 ENVI软件的基本操作。2、掌握卫星影像的预处理的基本流程。3、通过实习，学会自己去处理一些问题。4、进一步提高学生分析问题、解决问题的能力，增强实践技能，并培养学生勇于 动手、勤于动手、热爱本专业的思想。5、深刻地理解和巩固基本理论知识 , 掌握基本技能和动手操作能力 , 提高综合 观察分析问题的能力、实习容1、了解 ENVI的基本操作。2、实现影像图像的几何校正、融合、镶嵌及剪裁。3、掌握 EN

3、VI对影像信息的提取4、了解 ENVI的一些应用分析文档专题：基于环境小卫星的草原荒漠化监测一、专题概述浑善达克地区位于草原阴山北麓锡林郭勒高原中部，是亚洲草原荒漠化土地东部 边 缘区的重要组成部分，经纬度在东经 114 5511638，北纬 41 464307之 间，平均海拔 高度在 1100 米左右。退化区属温带半干旱区、中温带干旱大陆性季风气 候。浑善达克退化 土地多为沙地，植被稀少，特别是春季地面回暖解冻，地表裸露，多 细沙土，狂风起时沙尘 弥漫，形成沙尘天气。近年来频频发生在京津地区春季的沙尘天 气与该地区生态环境恶化， 人地关系严重失调相关。据统计，京津地区沙尘暴70%的沙源来自于

4、这个区域。根据该区域 不同时期的植被覆盖数据可以实现该区域的植被变化监 测。 本专题详细介绍了利用环境小卫星 CCD-1A 图像反演植被覆盖图的完整流程，专题 涉及 环境小卫星的数据读取、辐射定标、大气校正、植被覆盖反演模型的建立、遥感反 演过程、 植被覆盖变化监测等容。 除了使用 ENVI 主模块功能外， 还需要用到 FLAASH 大 气校正扩 展模块、 IDL 开发的环境小卫星数据读取补丁、波段运算等功能。二、处理流程介绍根据环境小卫星 CCD 数据特点及草原植被变化监测的要求，采用的技术路线为：先 对 环境小卫星 CCD 数据进行数据预处理： 数据读取、 辐射定标、 大气校正、 研究区裁

5、剪， 建 立反演模型，利用波段运算工具，反演出整个浑善达克地区的归一化植被指数、植被 覆盖度， 根据不同时期的植被覆盖数据，实现草原土地退化的遥感监测。 流程说明：2.1 图像获取本文的基础数据包括 2006 年土地覆盖图数据、 环境小卫星 CCD-1A 数据。环境小卫 星 CCD-1A 数据可以在环保部卫星环境应用中心免费下载获取， 土地覆盖图数据的制作时 间是 2006 年 8 月份，本专题选择了影像质量良好的 2009 年 8 月 11 日的环境星数据 作为数据源。 环境小卫星 CCD 数据下载地址： .secmep.2.2 数据读取和定标网上免费获取的 HJ-1A 卫星 CCD的分发格

6、式为 Geotiff ，每一个波段为一个 Geotiff 文件， 并提供一个元数据说明（ .XML）。可以依次打开每个波段，并用元数据说明中提供 的定标参 数用 Band Math 工具进行手动定标，再用 Layer stcaking 功能将 4 个波段 的图像合成为一个 多波段的数据文件。本专题中用的是环境小卫星的读取补丁，直接读 取 CCD 数据，输出结 果为一个多波段的 ENVI 标准栅格文件，并带有中心波长等信息， 可直接输出定标结果（辐 射亮度）。补丁下载地址为：bbs.esrichina-bj./ESRI/viewthread.php?tid=75575文档2.3 图像配准获取的

7、2006 年的土地利用分类图已经过精确的地理定位， 以该景作为基准影像， 对 环境小卫星数据进行配准，便于后续的变化监测。2.4 大气校正对于环境小卫星数据，提供光谱响应函数， 可用 FLAASH 大气校正模块中进行大气 校正， 波谱响应函数下载地址：.cresda./n16/n1115/n1522/n2118/index.html 。2.5 、反演模型构建及模型应用文档2.6 植被变化监测根据 2006 年 8 月和 2009 年 8 月的植被覆盖数据， 采用波段运算， 实现草原植被变化 的遥感监测文档三、详细处理过程本专题的数据存放在“ 19- 草原植被变化遥感监测”文件夹。3.1 数据预

8、处理3.1.1 安装环境小卫星数据处理补丁将 ENVI_HJ1A1B_Tools.sav 补丁放在： homeITTIDLIDL80productsenvi48save_add 目录下。3.1.2 数据读取和定标主菜单 -File-Open External File-HJ-1A/1B Tools，打开环境小卫星数据读取补丁。在 HJ-1A/1B Tools V3.0 面板中，选择 CCD，点击 Input File 输入“ 1- 环境小卫 星数据 HJ1A-CCD2-2-64-20090811-L20000154793154793 ”文件夹中的 .xml 文件， 点击 Output Path

9、 设置数据的输出路径，勾选“ Calibration ”“ Layer Stacking ”两个选项， 单击 Apply 按钮。文档3.1.3 工程区裁剪 由于整景影像数据围非常大， 本专题工程区只是其中一小部分， 在进行大气校正之前， 先将浑善达克以及周边区域裁剪出来。（1）打开上一步生成的文件： HJ1A-CCD2-2-64-20090811-L20000154793_Calbrated_LayerStacking.img ，文档2）主菜单 -File-Save File As-ENVI Standard，弹出 New File Builder 面板；3）在 New File Builde

10、r 面板中，单击 Import File ，弹出的 Create New File Input File 面板；4）在 Create New File Input File面板中，选中 Select Input File 列表中的裁剪数据，单击 Spatial Subset 按钮；文档5）在 Select Spatial Subset面板中，单击 Image ，弹出 Subset by Image 对话框，6）在 Subset by Image 对话框中，按住鼠标左键拖动图像中的红色矩形框确定裁剪区域， 裁剪出包括浑善达克区域的一部分，单击 OK；7）在 Select Spatial Subs

11、et面板中，可以看到裁剪区域信息，单击 OK；文档8）在 Create New File Input File对话框中，单击 OK；及路9）在 New File Builder ，单击 Choose 设置输出文件名 20090811-Cal-sub.img 径, 单击 OK。文档3.1.4 图像配准下面以土地利用图作为基准影像对环境小卫星图像进行图像配准。1）分别打开和显示基准影像“ 19- 基于环境小卫星的草原植被变化监测 浑善达克 2006 年 8 月土地利用分类图 .img ”。文档2）主菜单 -Map-Registration-Select GCPs:Image to Image，打开

12、几何校正模块。3) 选择显示 2006 年土地利用分类图文件的 Display 为基准影像( Base Image)，显示环境星文件的 Display 为待校正影像( Warp Image)，点击 OK 进入采集地面控制点4) 打开 Tools-Link-Geografic link，将两个窗口都选择为 on ，单击确定，找到定位的大 致区域后，再 Tools-Link-Geografic link ，改为 off ，关闭。5) 在两个 Display 中找到相同区域，在 Zoom 窗口中，点击左小下角第三个按钮， 打开定位十字光标，将十字光标到相同点上，点击 Ground Control P

13、oints selection上的 Add Point 按钮，将当前找到的点加入控制点列表。文档6）用同样的方法继续寻找其余的点，当选择控制点的数量达到3 时， RMS被 自动计算。 Ground Control Points Selection上的 Predict 按钮可用，选择Options-Auto Predict ， 打开自动预测功能。这时在 Base Image 上面定位点， Warp Image 上会自动预测区域。7）完成控制点的选择， RMS值 小于 1 个像素，点击 Ground Control Points，将控制点保存Selection 上的 File-Save Coeff

14、icients to ASCII8）在 Ground Control Points Selection Map） ，选 择校正文件 （HJ 数据文件 ）上，选择 Options- Warp File （as Image文档30 米。(9) 在校正参数面板中(图专 5-7 )，默认投影参数和像元大小与基准影像一致，(10) 重采样选择 Nearest Neighor ，背景值( Background)为 0.(11) Output Image Extent ：默认是根据基准图像大小计算，可以做适当的调整。(12) 选择输出路径和文件名，单击 Ok 按钮。文档3.1.5 大气校正 环境小卫星提供了

15、波谱响应函数，以文本形式提供，第一列表示波长（ nm），后面四 列 分别表示 4 个波段对应波长的波谱响应值。需要制作波谱曲线来描述波谱响应函数， 用于大气校正。制作波谱曲线 :（1）主菜单 Window-Start New Plot Window ，打开 ENVI Plot Window 面板，在 波谱绘制 窗口中,选择 File-Input Data-ASCII ，导入“ 681_HJ1ACCD2.txt” 文本文件，如图，自动将第一列作为 x 轴，后面 4 列作为 y 轴，波长单位 选择 Nanometers ，单击 OK。文档2) 如图，在绘制窗口生成了 4 条曲线，选择 Edit-D

16、ata Parameters ，编辑每条线的名称 为 b1 ，b2，b3，b4，便于区分文档3）选择 File-Save Plot As-Spectral Library，在 Output Plots to SpectralLibrary 面板中，单击 Select All Items，单击 OK。4）在 Output Spectral Library面板中，有输出曲线相关参数设置，按默认，选择输出路 径和文件名，单击 OK，将波谱曲线保存为波谱库文件：环境 1A 星 CCD2 光谱响应 .sli 。文档FLAASH大 气校正 :第一步：数据准备FLAASH 对图像文件有以下几个要求：(1)

17、数据是经过定标后的辐射亮度(辐射率)数据，单位是： (W)/ (cm2*nm*sr)。(2) 数据带有中心波长( wavelenth )值，如果是高光谱还必须有波段宽度( FWH)M, 这两 个参数都可以通过编辑头文件信息输入( Edit Header )。(3) 数据类型：支持四种数据类型：浮点型( floating )、长整型 (long integer ) 、整型 (integer )和无符号整型 (unsigned int) 。数据存储类型： ENVI 标准栅格格式文 件， 且是 BIP 或者 BIL 。 4) 波谱围： 4002500nm。本次用的环境小卫星经过以上处理，已经定标为

18、W*m(-2)*sr(-1)*um(-1) 单位、 浮点 型的辐射率数据，有中心波长信息，下面将 BSQ 格式转成 BIL 格式。选择主菜单 Basic Tools-Convert Data(BSQ、BIL、 BIP)，选择已经经过定标和配 准的数 据 20090811-cal-jz.img ，在 Convert File Parameters 中， Output Interleave 选择 BIL ，选择 Convert In Place ：yes ，单击 OK。第二步：设置参数进行 FLAASH 大气校正(1)主菜单 Spectral-FLAASH 打开 FLAASH 大气校正模块；(2)

19、 点击 Input Radiance Image ，选择 BIL 格式的环境小卫星数据20090811-cal-jz.img ，在 Radiance Scale Factors 面板中选择 Use single scale factor for all bands ，由于定标的辐 射量数据与 FLAASH 的辐射亮度的单位相差 10 倍，所以在此 Single scale factor 选择 默认： 10，单击 OK；(3) 设置输出文件及路径设置(4) 传感器基本信息设置： 成像中心点经纬度、传感器高度、成像区域平均高度、成像 时间设置，这些都可以 从数据头文件中读取：HJ1A-CCD2-2

20、-64-20090811-L20000154793.XML 。(5) 大气模型，选择 MLS，气溶胶模型，选择 Rural ，气溶胶反演方法选择 None，能见 度设置为 40km。文档6）单击 Multispectral Setting按钮，在 Filter Function File 导入之前做好的光谱响应曲线 “环境 1A 星 CCD2 光谱响应.sli ”，单击 OK;7）单击 Advanced Settings ，在高级设置中， Tile Size 默认的是 Cash size 的大小， 手动 改为 100Mb，单击 OK；文档8）设置好后，在大气校正模块面板中，单击 Apply 。

21、（9）大气校正完成后，检查大气校正的结果， 分别加载校正前后的图像， 将两幅影像进 行 地理，移动到植被区域（植被的波谱曲线比较特殊） ，在影像上右键，选择 Z Profile （Spectrum）打开光谱曲线窗口，显示两幅图像同一位置的光谱曲线图。文档3.1.6 裁剪浑善达克区（1）显示大气校正后的图像 20090811-cal-jz-FLAASH.img ，在 image 窗口选择Overlay-Vectors ，打开 Vector Parameters 面板，选择 File-Open Vector File ， 打开 hunshandake.evf ；文档2）在 Available Ve

22、ctors List面板中选择该矢量文件，点击 Load Selected ，选择显示图像 的 Display ，单击 OK，矢量叠加在影像上；3）在 Available Vectors List 面板中，选择 File-Export Layers to ROI，在 SelectData File to Associate with new ROI 面板中，选择 20090811-cal-jz-FLAASH.img ，文档单击 OK，在 Export EVF Layers to ROI 中，选择 Convert all records of an EVFl ayerto one ROI ，单击

23、 OK，将 矢量转为一个 ROI ；4）在图像窗口，选择 Overlay-Region of Interest ，打开 ROI 面板，浑善达克 ROI 显示在 图像上，在 ROI Tool 面板中，选择 File-Subset Data via ROIs，在 Select Input File to Subset via ROI面板中，选择 20090811-cal-jz-FLAASH.img ，单击 OK；文档5）在 Spatial Subset via ROI Parameters面板中选择浑善达克 ROI ，Mask pixelsoutside of ROI选择 Yes ；6）设置输出路

24、径及文件名 20090811-yanjiuqu.img ，单击 OK。文档到此， 数据预处理工作已经全部完成， 下面介绍如何应用 ENVI 进行植被覆盖度反演和 植被变化监测。3.2 植被覆盖度反演3.2.1 计算归一化植被指数（1） 选择主菜单 -File- Transform -NDVI ，打开 NDVI 计算模块 （2） 选择裁剪后的数据： 20090811-yanjiuqu.img 。3） 在 NDVI Calculation Parameters 面板中进行相应设置，环境星波段参数与 TM 数 据相似，应用被植被强吸收的红光波段（环境星第 3 波段）和被植被强反射的近 红外波 段（环

25、境星第四波段）计算归一化植被指数。数据类型默认为浮点型。4） 选择文件保存名（ 20090811-yanjiuqu-NDVI ）和路径，单击 OK 执行文档3.2.2 计算植被覆盖度计算植被覆盖度 FC 采用的是混合像元分解法， 将整景影像的地类大致分为水体、 植 被 和建筑，具体的计算公式如下：FC = (NDVI- NDVISoil) (NDVIVeg - NDVISoil) (1)其中， NDVI 为归一化植被指数， NDVISoil 为完全是裸土或无植被覆盖区域的 NDVI 值， NDVIVeg 则代表完全被植被所覆盖的像元的 NDVI 值，即纯植被像元的 NDVI 值。 取经验值 N

26、DVIVeg = 0.70 和 NDVISoil = 0.00 ，且有，当某个像元的 NDVI 大于 0.70 时， FC 取值为 1； 当 NDVI 小于 0.00 ，FC 取值为 0。利用 ENVI主菜单-Basic Tools-Band Math ，在公式输入栏中输入：(b1 gt 0.7)*1+(b1 lt 0.)*0+(b1 ge 0 and b1 le 0.7)*(b1-0.0)/(0.7-0.0)b1: 选择 NDVI 图像文档3.3 植被变化监测3.3.1 植被覆盖区提取(1)2009年 8月植被覆盖区提取利用 ENVI主菜单-Basic Tools-Band Math ，在公

27、式输入栏中输入：(b1 le 0.3)*0 +(b1 gt 0.3)*1b1: 选择 HJ1A-CCD2-20090811-cal-sub-flaash-ref-hunshandake-VCI 图像， 0.3 为经验值。选择文件保存名( 2009 年 8 月植被覆盖区)和路径，单击 OK 执行，得到 2009 年 8 月 植被盖度区。文档（2）2006年 8月植被覆盖区提取获取的 2006年 8 月标准的土地利用分类图， 分类代码与地物类型如下表所示， 耕地、 林 地、草地，属于植被覆盖区域。ENVI 主菜单 -Basic Tools-Band Math ，在公式输入栏中输入：(b1 ge 1 and b1 le 3)*1+(b1 lt 1)*0+(b1 gt 3)*0 b1: 选择“浑善达克 2006 年 8 月土地利用分类图” 。选择文件保存名 （200608 植被覆盖图 .img ）和路径，单击 OK 执