实验910遥感技术应用.docVIP

PAGE  PAGE 22 实验9 遥感技术应用 练习1植物遥感判读与植被指数制图 1.实验目的与任务： 利用遥感影像求植被指数，学习用波段运算和了解用IDL两种方法 2.实验设备与数据： 设备：遥感图像处理系统ENVI 数据：can-tmr.img或者bhtmr.img 植被指数 NDVI的计算 NDVI—是Normalized Difference Vegetation Index的缩写，即归一化植被指数：NDVI=(NIR-R)/(NIR+R)，或两个波段反射率的计算。 其取值范围 -1~1. NDVI的说明： 1）、NDVI的应用：检测植被生长状态、植被覆盖度和消除部分辐射误差等； 2）、-1=NDVI=1，负值表示地面覆盖为云、水、雪等，对可见光高反射；0表示有岩石或裸土等，NIR和R近似相等；正值，表示有植被覆盖，且随覆盖度增大而增大； 3）、NDVI的局限性表现在，用非线性拉伸的方式增强了NIR和R的反射率的对比度。对于同一幅图象，分别求RVI和NDVI时会发现，RVI值增加的速度高于NDVI增加速度，即NDVI对高植被区具有较低的灵敏度； 4）、NDVI能反映出植物冠层的背景影响，如土壤、潮湿地面、学、枯叶、粗超度等，且与植被覆盖有关； 计算步骤：方法一：利用波段运算 打开ENVI自带的数据中的遥感影像：bhtmref.img 使用BAND MATH 将经过了辐射定标(或者是定标+大气纠正)的图像（此实验仅仅用自带的数据，数值范围0-100）,运用工具中的band math 采用公式：NDVI=（b4-b3）/（b4+b3） 注意：在envi中计算的时候用 float（b4-b3）/（b4+b3） 其中b4近红外的反射率 b3红光的反射率 例如: 第二种方法：利用envi-transform下的NDVI直接自动计算。 第三种方法使用IDL编程计算NDVI（了解就可以，感兴趣的同学可以尝试做一下） IDL简介 IDL(interactive data language)交互式数据语言是进行数据分析,可视化表现和应用开发的软件工具.其特性包括: 高级图像处理能力,交互式二维和三维图形技术,面对对象的编程方式,opengl图形加速,量化可视化表现,集成数学与统计学算法,灵活的数据输入输出方式,跨平台图形用户界面工具包,连接ODBC兼容数据库存及多种程序连接工具等. 以下是idl的主界面: IDL程序的特点: 分隔符为”,”,而非空格; (b) 不分大小写; (c)变量无需事先申明 IDL的符号(需注意的) ?$ –作为一行的第一个字符时，返回到操作系统下，如：$ dir –作为一行的最后一个字符时，相当于一行未写完，换行。 ?；后面是注释 ?@ 批作业 如：@test IDL程序的主要结构: Pro 程序名 (程序体) END 另外我们可以在命令行里直接输入命令进行运行. 关于IDL的知识,大家可以找相应的资料看 下面是简单利用一下IDL做波段运算 A 运用IDL编写波段运算的公式,再在BAND MATH中调用 例如: 1）对此程序做调试,CTRL+F5 2）选择ENVI中的BAND MATH, 3）选择相应的波段进行运算. B 直接运用IDL编程,计算NDVI 例如: 直接在IDL的输入行中输入ndvi-math，回车；选择使用的影像就可以了 练习2 地表水信息提取与制图 1.实验目的与任务： （1）利用NDVI或者其它指数进行水体信息的提取； （2）熟悉遥感制图方法。 2.实验设备与数据： 设备：遥感图像处理系统ENVI 数据：lyg_water.img 或者焦作数据，或者软件自带数据。只要有水体就可以。 3 采用方法：利用植被指数NDVI，进行水体信息的提取。 从上面的说明可以看出，NDVI的负值为云、水和雪等。如果所用影像没有云和雪，那么NDVI小于0的区域就应该是水体。本实验就是据此进行水体提取（注意：有的时候水体比较浅，NDVI也会大于0）。 4 水体提取的具体步骤： 1）在ENVI中打开lyg_water.img影像； 2）计算此影像的NDVI (方法可以采用ENVI—basic tools- band math方法；也采用ENVI-Transform- NDVI直接计算得到影像的NDVI数值)； 3） 从计算的NDVI影像中提取数值小于0的区域； 方法：点击ENVI—basic tools- band math； 然后在出现的波段运算的对话框中输入b1 LT 0，然后点击OK;(意思是选择波段1中数值小于0的像元)； 然后在出现的要求选择b1波