地表反射率,温度,植被指数
- 格式:doc
- 大小:1.48 MB
- 文档页数:31
基本原理一)地表反射率是指地表物体向各个方向上反射的太阳总辐射通量与到达该物体表面上的总辐射通量之比。
反照率可以通过遥感成像提供的辐射亮度值L 或反照率p ,二向性反射率分布函数BRDF 来获得:地物反射率的光谱特征差异是从遥感影像中识别地表不同类型地物的基本依据,也是地表其他各种物理、生物物理参数反演的依据地表。
地表反射率的计算步骤:1、辐射定标:根据遥感影像DN 值计算到达传感器的各波段辐射亮度也就是将传感器记录的辐射量化值(Digital Number ,DN )转换成绝对辐射亮度值、表观反射率,或者表观温度的过程。
绝对定标:通过各种标准辐射源,建立辐射亮度值与辐射量化值(DN )之间的定量关系式中,辐射亮度值L 的常用单位为W/(m2.μm.sr),或者μW/(cm2.nm.sr) 。
1W/(m2.μm.sr)=0.1 μW/(cm2.nm.sr)2、各波段表观反射率计算3、大气辐射校正(ENVI FLAASH/QUAC )绝对大气辐射校正:消除大气辐射衰减效应,将遥感影像的DN 值转换为地表反射率、辐亮度、地表温度等的方法,此过程包含了辐射定标。
相对大气辐射校正:将遥感影像的DN 值转换为类似的整型数,同时消除大气辐射衰减效应。
FLAASH 是用数学建模辐射的物理行为,纠正波长在可见光至近红外和短波红外区域,最多3微米。
(对于热地区,使用基本工具>预处理>校准工具>热大气压校正菜单选项。
)不同于预先计算模拟结果的数据库内插辐射传输特性许多其他大气校正程序, FLAASH 采用了MODTRAN4辐射传输代码。
MODTRAN4并入ENVI FLAASH 的版本被修改,以校正在HITRAN -96水行参数的误差。
可以选择任何一种标准MODTRAN 大气模型和气溶胶类型,FLAASH 还包括以下功能:校正邻近效应(像素混合是由于表面反射辐射的散射) 计算场景的平均能见度(气溶胶/雾量)。
地表反射率、温度、植被指数、几何精纠正和Landsat影像影像几何精纠正1.深入理解影像几何精纠正的原理2.学会使用影像对影像的几何精纠正方法和具体操作步骤1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
2. 利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正。
1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
2. 利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正Landsat影像数据下载、导入、目视解译与分析1.学会在美国NASA/USGS网上进行Landsat卫星影像的检索和下载,认识影像名编号意义。
2.初步掌握ENVI/IDL影像处理软件的使用方法,熟悉软件的用户界面、功能模块,掌握基本功能的使用。
复习遥感导论课程中的遥感影像目视解译环节,选择自己熟悉的地区,进行Landsat TM/ETM+影像的目视解译。
实习内容1、课前准备:根据自己感兴趣的地区,下载一景Landsat TM/ETM+影像。
2、将单波段分别添加在ENVI中(采用file| open image file工具,注意在添加波段数据时,热红外波段影像应另存为一个文件)3、利用Basic tools | layer stscking视窗选择多波段影像进行添加,在available中选择刚才添加的影像,选择RGB color分别选取4、3、2和7、4、2波段组合进行彩色合成,此时,可以再打开一个视窗(new display)观察地物的色调变化。
列表说明上述地物分别在两种波段组合下的颜色特征。
4、查询并记录影像文件的基本信息、投影信息,以及各个波段直方图信息。
5、将影像缩小、放大、漫游工具识别影像中的土地利用/土地覆盖类型,可能的土地利用/土地覆盖类型包括:(1) 耕地farmland(8) 公路/铁路road or railway(2) 草地grassland(9) 河流stream(3) 裸地barren land(10) 水库reservoir(4) 森林forest(11) 冰雪ice and snow(5) 城镇居民地town(12) 云cloud(6) 农村居民地village(13) 阴影shadow(7) 沙漠desert6、利用load RGB将图像显示,后用视窗中的光谱剖面工具,提取上述地物在不同波段的数值(Digital Number,DN);要求针对影像中的6种地物至少各采集10个样本,取平均值,做光谱剖面图,分析不同地物的灰度值随波段变化的特点。
操作方法及过程1、使用ENVI对landsat 7 ETM+原始数据进行辐射定标:①对1、2、3、4、5、7波段进行辐射定标。
利用ENVI中的File |Open External File |Landsat Geo TIFF with MetaData加载威武市Landsat ETM+原始影像数据中的_MTL文件,再利用Basic Tools |Preprocessing |Calibration Utilities |Landsat Calibration 在弹出的对话对话框中选择包含1、2、3、4、5、7波段的_MTL文件,将Calibration Type选为Radiance,然后选择输出路径保存为radiance。
②对61和62波段进行辐射定标。
步骤和上面的一样,只是选择输入文件时为包含61和62波段的_MTL文件,将结果保存为radiance_band6。
2、将BSQ格式的影像数据转化为BIL:利用Basic Tools |Convert Data,弹出对话框中选择Radiance,Output Interleave中选择BIL,选择输出路径保存为radiance_BIL。
3、使用FLAASH大气辐射校正模型进行地表反射率的计算:①利用Spectral |FLASSH弹出大气校正模型参数设置窗口如下:分别按照以上所示的内容进行参数设置,将输入文件设为radiance_BIL,输出文件设为flassh,设置Scene Center Location时,打开原始影像在头文件中找到行和列,算出中心行和列,利用Pixel Locator工具找到中心点的经纬度。
将Sensor Type设为Landsat TM7。
设置Ground Elevation时,利用裁剪工具在亚洲幅SRTM DEM影像数据中裁剪该地区的DEM数据,再用统计功能算出高程的平均值为2058m。
在头文件中找到Flight Data:1999年8月10日,Flight Time GTM:3时36分39秒。
MODIS数据介绍数据概况1999年2月18日,美国成功地发射了地球观测系统(EOS)的第一颗先进的极地轨道环境遥感卫星Terra。
它的主要目标是实现从单系列极轨空间平台上对太阳辐射、大气、海洋和陆地进行综合观测,获取有关海洋、陆地、冰雪圈和太阳动力系统等信息,进行土地利用和土地覆盖研究、气候季节和年际变化研究、自然灾害监测和分析研究、长期气候变率的变化以及大气臭氧变化研究等,进而实现对大气和地球环境变化的长期观测和研究的总体(战略)目标。
2002年5月4日成功发射Aqua星后,每天可以接收两颗星的资料。
搭载在Terra和Aqua两颗卫星上的中分辨率成像光谱仪(MODIS)是美国地球观测系统(EOS)计划中用于观测全球生物和物理过程的重要仪器。
它具有36个中等分辨率水平(0.25um~1um)的光谱波段,每1-2天对地球表面观测一次。
获取陆地和海洋温度、初级生产率、陆地表面覆盖、云、汽溶胶、水汽和火情等目标的图像。
本网站提供的MODIS陆地标准产品来自NASA的陆地过程分布式数据档案中心(The Land Processes Distributed Active Archive Center,LP DAAC/NASA)。
包括:基于Terra星和Aqua星数据的地表反射率(250m,daily;500m,daily;250m,8days;500m,8day)、地表温度(1000m,daily;1000m,8days;5600m,daily)、地表覆盖(500m,96days;1000m,yearly)、植被指数NDVI&EVI (250m,16daily;500m,16days;1000m,16days;1000m,monthly;、温度异常/火产品(1000m,daily;1000m,8days)、叶面积指数LAI/光合有效辐射分量FPAR(1000m,8days)、总初级生产力GPP(1000m,8days)。
地表反射率、温度、植被指数、几何精纠正和Landsat影像一、DNVI建模【地表反射率】第3波段第4波段【DNVI】【3、4波段表观反射率和地表反射率的线性关系】【表观反射率和地表反射率的线性关系数学表达式】波段关系式波段关系式1波段y=0.8933*x+0.0473 4波段y=0.9401*x+0.0065 2波段y=0.8801*x+0.0242 5波段y=0.9399*x+0.001 3波段y=0.9161*x+0.0143 7波段y=0.9584*x+0.0004【部分地物的DNVI值】地物DNVI值min max mean stdev Reservior 0.057713 0.338587 0.145087 0.038598 Snow -0.12395 0.152669 0.025088 0.031572Urban -0.356923 0.038094 -0.273288 0.045284Plant 0.333387 0.786695 0.656094 0.081619Desert 0.071897 0.155663 0.100783 0.014291River 0.043469 0.429917 0.127503 0.08131【结果与分析】:通过对提取地物的DNVI值的可以发现,绿色的DNVI值比较高,原因是绿色植物叶绿素引起的红光吸收和叶肉组织引起的近红外光反射使得植被在近红外波段和红光波段有很大的差异;水体和裸地在红光波段和近红外波段反射率相当,因此水库和裸地的NDVI值接近0;雪地NDVI最低值中出现负值,是由于在近红外波段比可见光波段有较低的反射率;沙漠中植被很少,因此其近红外波段和红光波段的反射情况和裸地类似,因此其NDVI值接近于0;河流的NDVI值稍大于由于河流中存在一定的含沙量,使得地物在近红外波段的反射率大于近红外波段。
二、温度反演【温度反演】低增益温度反演高增益温度反演【第1波段部分地物低增益温度反演数据】开尔文温度摄氏温度反演温度地物min max mean stdev min max mean Reservior 287.47641 289.289886 288.13127 0.388036 14.32641 16.13989 14.98127 Snow 273.154785 293.990417 278.177771 3.788266 0.004785 20.84042 5.027771 Bare Land 295.989319 310.676086 303.445647 2.819391 22.83932 37.52609 30.29565 Urban 300.165253 310.928528 307.469228 1.530421 27.01525 37.77853 34.31923 Plant 294.278015 305.525879 298.698402 2.333251 21.12802 32.37588 25.5484 Desert 302.605286 309.915955 306.491575 1.39902 29.45529 36.76596 33.34158影像几何精纠正1.深入理解影像几何精纠正的原理2.学会使用影像对影像的几何精纠正方法和具体操作步骤1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
资料范本本资料为word版本,可以直接编辑和打印,感谢您的下载植被指数整理介绍地点:__________________时间:__________________说明:本资料适用于约定双方经过谈判,协商而共同承认,共同遵守的责任与义务,仅供参考,文档可直接下载或修改,不需要的部分可直接删除,使用时请详细阅读内容植被指数介绍目录TOC \o "1-3" \h \z \u HYPERLINK \l "_Toc408667889" 1. 植被指数概述 PAGEREF _Toc408667889 \h 3HYPERLINK \l "_Toc408667890" 2. 植被指数的分类 PAGEREF _Toc408667890 \h 3HYPERLINK \l "_Toc408667891" 2.1不考虑影响因子 PAGEREF_Toc408667891 \h 3HYPERLINK \l "_Toc408667892" 2.2考虑影响因子 PAGEREF_Toc408667892 \h 4HYPERLINK \l "_Toc408667893" 2.2.1 消除土壤因子 PAGEREF _Toc408667893 \h 4HYPERLINK \l "_Toc408667894" 2.2.2 消除大气因子 PAGEREF _Toc408667894 \h 4HYPERLINK \l "_Toc408667895" 2.2.3 消除综合因子 PAGEREF _Toc408667895 \h 5HYPERLINK \l "_Toc408667896" 3. 植被指数的应用 PAGEREF _Toc408667896 \h 5HYPERLINK \l "_Toc408667897" 3.1生态 PAGEREF_Toc408667897 \h 5HYPERLINK \l "_Toc408667898" 3.2林业 PAGEREF_Toc408667898 \h 7HYPERLINK \l "_Toc408667899" 3.3农业 PAGEREF_Toc408667899 \h 9HYPERLINK \l "_Toc408667900" 3.4环境 PAGEREF_Toc408667900 \h 10HYPERLINK \l "_Toc408667901" 3.5海洋 PAGEREF_Toc408667901 \h 11HYPERLINK \l "_Toc408667902" 参考文献 PAGEREF_Toc408667902 \h 12植被指数概述植被指数是用不同波段的植被-土壤系统的反射率因子以一定形式组合成的参数,它与植被特征参数间的函数联系比单一波段值更稳定、可靠[1]。
操作方法及过程1、使用ENVI对landsat 7 ETM+原始数据进行辐射定标:①对1、2、3、4、5、7波段进行辐射定标。
利用ENVI中的File |Open External File |Landsat Geo TIFF with MetaData加载威武市Landsat ETM+原始影像数据中的_MTL文件,再利用Basic Tools |Preprocessing |Calibration Utilities |Landsat Calibration 在弹出的对话对话框中选择包含1、2、3、4、5、7波段的_MTL文件,将Calibration Type选为Radiance,然后选择输出路径保存为radiance。
②对61和62波段进行辐射定标。
步骤和上面的一样,只是选择输入文件时为包含61和62波段的_MTL文件,将结果保存为radiance_band6。
2、将BSQ格式的影像数据转化为BIL:利用Basic Tools |Convert Data,弹出对话框中选择Radiance,Output Interleave中选择BIL,选择输出路径保存为radiance_BIL。
3、使用FLAASH大气辐射校正模型进行地表反射率的计算:①利用Spectral |FLASSH弹出大气校正模型参数设置窗口如下:分别按照以上所示的内容进行参数设置,将输入文件设为radiance_BIL,输出文件设为flassh,设置Scene Center Location时,打开原始影像在头文件中找到行和列,算出中心行和列,利用Pixel Locator工具找到中心点的经纬度。
将Sensor Type设为Landsat TM7。
设置Ground Elevation时,利用裁剪工具在亚洲幅SRTM DEM影像数据中裁剪该地区的DEM数据,再用统计功能算出高程的平均值为2058m。
在头文件中找到Flight Data:1999年8月10日,Flight Time GTM:3时36分39秒。
植被的反射特性植被不仅仅是我们家园中的装饰,它还在有效地调节着地表的热量,保护我们的环境免受极端的气温,而最显著的功能之一就是反射。
反射是植物态植物与地表表面之间的相互作用,它主要是由入射太阳光反射离地表表面而形成的。
通过反射,可以将太阳光中的热量减少。
即使在非常高温的天气,植物也可以帮助控制地表温度在一定的范围内,降低地表表面的温度,进而有利于我们身体健康。
反射率是表示植物表面反射太阳光到地表的热量的指标,它在不同季节、不同的植物和不同的地表表面都会有差别。
植物反射率的比较因子有叶片的大小、叶片的排列方式、叶色、叶片表面的反光性、植物的类型、地表表面的坡度、植被覆盖率和地表表面形态等。
根据不同的季节,反射率也会不同,经过改变季节,反射率也随之改变,从而对地表表面的温度也会产生影响。
植被为我们提供反射效果,不仅可以降低地表表面的温度,还可以减少热量的蒸发。
它可以阻止太阳光直接照射地表表面,阻挡蒸发的热量,减缓地表温度的升高。
而且在炎热的盛夏,植物反射的太阳光可以遮挡太阳直射,对空气湿度的影响也有一定的减小作用。
另外,植物的反射特性还可以使周围的空气净化,它可以完成光谱转换,将太阳辐射转换成低能量的辐射,既保护环境,又不影响人体正常生理过程。
同时,植物反射的太阳光可以把紫外线、可见光以及近红外辐射转换成低能量的可见光,以及对地球的温度和亮度有一定的影响。
植物的反射特性,不仅能够帮助我们调节温度,也可以有效地净化我们的空气,很好地保护我们的环境,使我们的家园更加绿色、健康和美丽。
它不仅仅是一种装饰,更是一种社会责任,是我们应该去把握的每一个可见的机会,以保护我们的环境。
只有这样,我们才能把家园变成一个更加绿色、健康和美丽的地方,给子孙后代留下更好的环境。
植被指数(Vegetable Index)植被指数是不同遥感光谱波段间的线性或非线性组合,被认为能作为反映绿色植被的相对丰度和活性的辐射量值(无量纲)的标志,是绿色植被的叶面积指数(LAI)、盖度、叶绿素含量、绿色生物量以及被吸收的光合有效辐射(APAR)的综合体现。
目前,在科学文献中发布了超过150种植被指数模型,这些植被指数中只有极少数是经过系统的实践检验。
植被指数主要反映植被在可见光、近红外波段反射与土壤背景之间差异的指标,各个植被指数在一定条件下能用来定量说明植被的生长状况。
1、健康的绿色植被在NIR和R的反射差异比较大,原因在于R对于绿色植物来说是强吸收的,NIR则是高反射高透射的;2、建立植被指数的目的是有效地综合各有关的光谱信号,增强植被信息,减少非植被信息3、植被指数有明显的地域性和时效性,受植被本身、环境、大气等条件的影响几种常用的植被指数及其应用(一)比值植被指数(RVI)公式:RVI=ρNIR/ρRED(近红外波段反射率/红光波段反射率)特征:植被覆盖度影响RVI,当植被覆盖度较高时,RVI对植被十分敏感;当植被覆盖度<50%时,这种敏感性显著降低;值的范围是0-30+,一般绿色植被区的范围是2-8。
RVI受大气条件影响,大气效应大大降低对植被检测的灵敏度,所以在计算前需要进行大气校正,或用反射率计算RVI。
应用:①利用比值植被指数研究城市建设用地扩张速率,预测或规划城市未来今年的发展前景。
不同用地的地表温度由高到低排序是城镇用地、工矿与交通用地、农村宅基地、林地、旱地,说明建设用地的地表温度较高,其比值植被指数较非建设用地小。
RVI的平均值M和标准差D可以作为定量指标来提取建设用地:RVI ≤M-D/2为建设用地;RVI>M-D/2为非建设用地。
②可用于实时、快速、无损监测作物氮素状况,这对于精确氮肥管理有重要意义。
利用高光谱比值指数RSI(990,720)来估算小麦叶片氮积累量为便携式小麦氮素监测仪的研制开发及遥感信息的快速提取提供了适用可行的波段选择与技术依据。
常用的植被指数,土壤指数,水体指数有哪些?植被指数与土壤指数一、RVI——比值植被指数:RVI=NIR/R,或两个波段反射率的比值。
1、绿色健康植被覆盖地区的RVI远大于1,而无植被覆盖的地面(裸土、人工建筑、水体、植被枯死或严重虫害)的RVI在1附近。
植被的RVI通常大于2;2、RVI是绿色植物的灵敏指示参数,与LAI、叶干生物量(DM)、叶绿素含量相关性高,可用于检测和估算植物生物量;3、植被覆盖度影响RVI,当植被覆盖度较高时,RVI对植被十分敏感;当植被覆盖度<50%时,这种敏感性显著降低;4、RVI受大气条件影响,大气效应大大降低对植被检测的灵敏度,所以在计算前需要进行大气校正,或用反射率计算RVI。
二、NDVI——归一化植被指数:NDVI=(NIR-R)/(NIR+R),或两个波段反射率的计算。
1、NDVI的应用:检测植被生长状态、植被覆盖度和消除部分辐射误差等;2、-1<=NDVI<=1,负值表示地面覆盖为云、水、雪等,对可见光高反射;0表示有岩石或裸土等,NIR和R近似相等;正值,表示有植被覆盖,且随覆盖度增大而增大;3、NDVI的局限性表现在,用非线性拉伸的方式增强了NIR和R的反射率的对比度。
对于同一幅图象,分别求RVI和NDVI时会发现,RVI值增加的速度高于NDVI增加速度,即NDVI对高植被区具有较低的灵敏度;4、NDVI能反映出植物冠层的背景影响,如土壤、潮湿地面、学、枯叶、粗超度等,且与植被覆盖有关;三、DVI\EVI——差值\环境植被指数:DVI=NIR-R,或两个波段反射率的计算。
1、对土壤背景的变化极为敏感;四、SAVI\TSAVI\MSAVI——调整土壤亮度的植被指数:SAVI=((NIR-R)/(NIR+R+L))(1+L),或两个波段反射率的计算。
1、目的是解释背景的光学特征变化并修正NDVI对土壤背景的敏感。
与NDVI 相比,增加了根据实际情况确定的土壤调节系数L,取值范围0~1。
在遥感应用领域,植被指数已广泛用来定性和定量评价植被覆盖及其生长活力。
由于植被光谱表现为植被、土壤亮度、环境影响、阴影、土壤颜色和湿度复杂混合反应,而且受大气空间—时相变化的影响,因此植被指数没有一个普遍的值,其研究经常表明不同的结果。
研究结果表明,利用在轨卫星的红光和红外波段的不同组合进行植被研究非常好,这些波段在气象卫星和地球观测卫星上都普遍存在,并包含90%以上的植被信息,这些波段间的不同组合方式统被称为植被指数。
植被指数有助于增强遥感影像的解译力,并已作为一种遥感手段广泛应用于土地利用覆盖探测、植被覆盖密度评价、作物识别和作物预报等方面,并在专题制图方面增强了分类能力。
植被指数还可用来诊断植被一系列生物物理参量:叶面积指数(LAI)、植被覆盖率、生物量、光合有效辐射吸收系数(APAR)等;反过来又可用来分析植被生长过程:净初级生产力(NPP)和蒸散(蒸腾)等。
为了估算和监测植被覆盖,最早发展了比值植被指数(RVI)。
但RVI对大气影响敏感,而且当植被覆盖不够浓密时(小于50%),它的分辨能力也很弱,只有在植被覆盖浓密的情况下效果最好。
归一化差异植被指数(NDVI)对绿色植被表现敏感,它可以对农作物和半干旱地区降水量进行预测,该指数常被用来进行区域和全球的植被状态研究。
对低密度植被覆盖,NDVI对于观测和照明几何非常敏感。
但在农作物生长的初始季节,将过高估计植被覆盖的百分比;在农作物生长的结束季节,将产生估计低值。
继之,将各波段反射率以不同形式进行组合来消除外在的影响因素,如遥感器定标、大气、观测和照明几何条件等。
这些线性组合或波段比值的指数发展满足特定的遥感应用,如作物产量、森林开发、植被管理和探测等。
农业植被指数(A VI)针对作物生长阶段测量绿色植被;多时相植被指数(MTVI),将两个不同日期的数值简单相减,是为了观测两个日期植被覆盖条件的变化和作物类型的分类,并用来探测由于火灾和土地流失造成的森林覆盖变化。
高光谱遥感技术在农林植被调查方面的应用高光谱分辨率遥感(简称高光谱遥感),是20世纪末迅速发展起来的一项集探测器技术、精密光学机械、微弱信号检测、计算机和信息处理技术于一体的全新遥感技术。
它能够获得地物的连续光谱信息,实现地物图像信息与光谱信息的同步获取,因而在地质、林业、农业、生态环境、海洋、军事等领域具有巨大的应用价值和广阔的发展前景。
植被作为遥感观测和记录的第一表层,是遥感数据反映的最直接信息。
目前,多光谱遥感已被广泛应用于植被的长势研究、沙漠化研究、气候演变规律分析等方面,但传统的多波段遥感数据对于植被的研究和应用仍仅限于一般性的红光吸收特征与近红外的反射特征及中红外的水吸收特征波段,由于受波段宽度和波段数以及波长位置的限制,往往对植被类型不敏感,对植被长势反映不理想,而高光谱遥感在对目标的空间特征成像的同时,对每个像元可在更宽范围上,形成几十个乃至几百个窄波段连续的光谱覆盖,使更深入地考察植被光谱的响应机制和物理机制成为可能,因此成为植被和林业方面监测的强有力工具。
1. 高光谱遥感在植被调查方面的优势高光谱遥感在光谱分辨率方面的提高,使地物目标的属性信息探测能力有所增强,因此,较之全色和多光谱遥感,高光谱遥感有以下显著优势:(1)成像光谱仪所获取的地物连续光谱比较真实,能全面反映自然界各种植被所固有的光谱特征以及其间的细节差异性,从而大大提高地物分类的精细程度和准确性,使得高光谱图像数据与光谱仪地面实测光谱曲线数据之间的直接匹配成为可能。
(2)高光谱图像数据提高了根据混合光谱模型进行混合像元分解的能力,减少了土壤等植被生长背景地物的影响,从而能够获取最终光谱端元的真实光谱特性曲线数据。
(3)高光谱分辨率的植被图像数据将对传统的植被指数运算予以改进,大大提高了植被指数所能反演的信息量,使人们可以直接收获诸如植被叶面积指数、生物量、光合有效吸收系数等植被生物物理参量。
(4)提高遥感高定量分析的精度和可靠性,基于高光谱分辨率的光谱吸收特征信息提取可以完成部分植被生物化学成分(如植被干物质和水分含量等)定量填图。
几种常见植被指数标准化管理部编码-[99968T-6889628-J68568-1689N]植被指数主要反映植被在可见光、近红外波段反射与土壤背景之间差异的指标,各个植被指数在一定条件下能用来定量说明植被的生长状况。
在学习和使用植被指数时必须由一些基本的认识:1、健康的绿色植被在NIR和R的反射差异比较大,原因在于R对于绿色植物来说是强吸收的,NIR则是高反射高透射的;2、建立植被指数的目的是有效地综合各有关的光谱信号,增强植被信息,减少非植被信息3、植被指数有明显的地域性和时效性,受植被本身、环境、大气等条件的影响一、RVI——比值植被指数:RVI=NIR/R,或两个波段反射率的比值。
1、绿色健康植被覆盖地区的RVI远大于1,而无植被覆盖的地面(裸土、人工建筑、水体、植被枯死或严重虫害)的RVI在1附近。
植被的RVI通常大于2;2、RVI是绿色植物的灵敏指示参数,与LAI、叶干生物量(DM)、叶绿素含量相关性高,可用于检测和估算植物生物量;3、植被覆盖度影响RVI,当植被覆盖度较高时,RVI对植被十分敏感;当植被覆盖度<50%时,这种敏感性显着降低;4、RVI受大气条件影响,大气效应大大降低对植被检测的灵敏度,所以在计算前需要进行大气校正,或用反射率计算RVI。
二、NDVI——归一化植被指数:NDVI=(NIR-R)/(NIR+R),或两个波段反射率的计算。
1、NDVI的应用:检测植被生长状态、植被覆盖度和消除部分辐射误差等;2、-1<=NDVI<=1,负值表示地面覆盖为云、水、雪等,对可见光高反射;0表示有岩石或裸土等,NIR和R近似相等;正值,表示有植被覆盖,且随覆盖度增大而增大;3、NDVI的局限性表现在,用非线性拉伸的方式增强了NIR和R的反射率的对比度。
对于同一幅图象,分别求RVI和NDVI时会发现,RVI值增加的速度高于NDVI增加速度,即NDVI对高植被区具有较低的灵敏度;4、NDVI能反映出植物冠层的背景影响,如土壤、潮湿地面、学、枯叶、粗超度等,且与植被覆盖有关;三、DVI\EVI——差值\环境植被指数:DVI=NIR-R,或两个波段反射率的计算。
植被光谱分析与植被指数计算解读植被光谱分析与植被指数计算在遥感中,常常结合不同波长范围的反射率来增强植被特征,如植被指数(vegetation indices ——VI)的计算,植被指数(VI)是两个或多个波长范围内的地物反射率组合运算,以增强植被某⼀特性或者细节。
⽬前,在科学⽂献中发布了超过150种植被指数模型,这些植被指数中只有极少数是经过系统的实践检验。
本⽂总结现有植被指数,根据对植被波谱特征产⽣重要影响的主要化学成份:⾊素(Pigments)、⽔分(Water)、碳(Carbon)、氮(Nitrogen),总结了7⼤类实⽤性较强的植被指数,即:宽带绿度、窄带绿度、光利⽤率、冠层氮、⼲旱或碳衰减、叶⾊素、冠层⽔分含量。
这些植被指数可以简单度量绿⾊植被的数量和⽣长状况、叶绿素含量、叶⼦表⾯冠层、叶聚丛、冠层结构、植被在光合作⽤中对⼊射光的利⽤效率、测量植被冠层中氮的相对含量、估算纤维素和⽊质素⼲燥状态的碳含量、度量植被中与胁迫性相关的⾊素、植被冠层中⽔分含量等。
包括以下内容:●植被光谱特征●植被指数●HJ-1-HSI植被指数计算1.植被光谱特征植被跟太阳辐射的相互关系有别于其他物质,如裸⼟、⽔体等,⽐如植被的“红边”现象,即在<700nm附近强吸收,>700nm⾼反射。
很多因素影响植被对太阳辐射的吸收和反射,包括波长、⽔分含量、⾊素、养分、碳等。
研究植被的波长范围⼀般为400 nm to 2500 nm,这也是传感器设计选择的波长范围。
这个波长范围可范围以下四个部分:●可见光(Visible):400 nm to 700 nm●近红外(Near-infrared——NIR):700 nm to 1300 nm●短波红外1(Shortwave infrared 1—— SWIR-1):1300 nm to 1900 nm●短波红外2(Shortwave infrared 2——SWIR-2):1900 nm to 2500 nm其中NIR和SWIR-1的过渡区(1400nm附近)是⼤⽓⽔的强吸收范围,卫星或者航空传感器⼀般不获取这范围的反射值。
地表反射率、温度、植被指数、几何精纠正和Landsat影像影像几何精纠正1.深入理解影像几何精纠正的原理2.学会使用影像对影像的几何精纠正方法和具体操作步骤1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
2. 利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正。
1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
2. 利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正Landsat影像数据下载、导入、目视解译与分析1.学会在美国NASA/USGS网上进行Landsat卫星影像的检索和下载,认识影像名编号意义。
2.初步掌握ENVI/IDL影像处理软件的使用方法,熟悉软件的用户界面、功能模块,掌握基本功能的使用。
复习遥感导论课程中的遥感影像目视解译环节,选择自己熟悉的地区,进行Landsat TM/ETM+影像的目视解译。
实习内容1、课前准备:根据自己感兴趣的地区,下载一景Landsat TM/ETM+影像。
2、将单波段分别添加在ENVI中(采用file| open image file工具,注意在添加波段数据时,热红外波段影像应另存为一个文件)3、利用Basic tools | layer stscking视窗选择多波段影像进行添加,在available中选择刚才添加的影像,选择RGB color分别选取4、3、2和7、4、2波段组合进行彩色合成,此时,可以再打开一个视窗(new display)观察地物的色调变化。
列表说明上述地物分别在两种波段组合下的颜色特征。
4、查询并记录影像文件的基本信息、投影信息,以及各个波段直方图信息。
5、将影像缩小、放大、漫游工具识别影像中的土地利用/土地覆盖类型,可能的土地利用/土地覆盖类型包括:(1) 耕地farmland(8) 公路/铁路road or railway(2) 草地grassland(9) 河流stream(3) 裸地barren land(10) 水库reservoir(4) 森林forest(11) 冰雪ice and snow(5) 城镇居民地town(12) 云cloud(6) 农村居民地village(13) 阴影shadow(7) 沙漠desert6、利用load RGB将图像显示,后用视窗中的光谱剖面工具,提取上述地物在不同波段的数值(Digital Number,DN);要求针对影像中的6种地物至少各采集10个样本,取平均值,做光谱剖面图,分析不同地物的灰度值随波段变化的特点。
对几种常用植被指数的认识植被指数主要反映植被在可见光、近红外波段反射与土壤背景之间差异的指标,各个植被指数在一定条件下能用来定量说明植被的生长状况。
在学习和使用植被指数时必须由一些基本的认识:1、健康的绿色植被在NIR和R的反射差异比较大,原因在于R对于绿色植物来说是强吸收的,NIR则是高反射高透射的;2、建立植被指数的目的是有效地综合各有关的光谱信号,增强植被信息,减少非植被信息3、植被指数有明显的地域性和时效性,受植被本身、环境、大气等条件的影响一、RVI——比值植被指数:RVI=NIR/R,或两个波段反射率的比值。
1、绿色健康植被覆盖地区的RVI远大于1,而无植被覆盖的地面(裸土、人工建筑、水体、植被枯死或严重虫害)的RVI在1附近。
植被的RVI通常大于2;2、RVI是绿色植物的灵敏指示参数,与LAI、叶干生物量(DM)、叶绿素含量相关性高,可用于检测和估算植物生物量;3、植被覆盖度影响RVI,当植被覆盖度较高时,RVI对植被十分敏感;当植被覆盖度<50%时,这种敏感性显著降低;~4、RVI受大气条件影响,大气效应大大降低对植被检测的灵敏度,所以在计算前需要进行大气校正,或用反射率计算RVI。
二、NDVI——归一化植被指数:NDVI=(NIR-R)/(NIR+R),或两个波段反射率的计算。
1、NDVI的应用:检测植被生长状态、植被覆盖度和消除部分辐射误差等;2、-1<=NDVI<=1,负值表示地面覆盖为云、水、雪等,对可见光高反射;0表示有岩石或裸土等,NIR和R近似相等;正值,表示有植被覆盖,且随覆盖度增大而增大;3、NDVI的局限性表现在,用非线性拉伸的方式增强了NIR和R的反射率的对比度。
对于同一幅图象,分别求RVI和NDVI时会发现,RVI值增加的速度高于NDVI增加速度,即NDVI对高植被区具有较低的灵敏度;4、NDVI能反映出植物冠层的背景影响,如土壤、潮湿地面、学、枯叶、粗超度等,且与植被覆盖有关;三、DVI\EVI——差值\环境植被指数:DVI=NIR-R,或两个波段反射率的计算。
地表反射率、温度、植被指数、几何精纠正和Landsat影像Basic Tools|Band Math,在Band Math对话框中输入公式,公式中的b3和b4分别选取第3和第4波段的地表反射率。
然后导出结果。
二、地表温度反演1、计算辐射亮度。
加载htm影像,根据头文件中的数据,得到1、2波段的辐射亮度的计算公式0.067086617777667001*b1+(-0.067086617777667001)和0.037204722719868001*b2+(3.1627953249638470),步骤同上,得出辐射量度的计算结果。
2、辐射反演。
利用公式T=k2/ln(k1/Lλ+1)算地物的辐射反演,其中T为开尔温度;查找参数值:k1=666.09; k2=1282.71;Lλ分别利用步骤1中的波段1和波段2的辐射量度。
3、统计反演后的地物的温度值,并比较其差异。
打开反演后的温度影像,右击影像选择ROI Tool,统计各种地物值的最大值,最小值,均值,标准差,将其统计到Excel中,比较其差异。
结果与分析一、DNVI建模【地表反射率】第3波段第4波段【DNVI】【3、4波段表观反射率和地表反射率的线性关系】【表观反射率和地表反射率的线性关系数学表达式】波段关系式波段关系式1波段y=0.8933*x+0.0473 4波段y=0.9401*x+0.00652波段y=0.8801*x+0.0242 5波段y=0.9399*x+0.0013波段y=0.9161*x+0.0143 7波段y=0.9584*x+0.0004【部分地物的DNVI值】地物DNVI值min max mean stdevReservior 0.057713 0.338587 0.145087 0.038598Snow -0.12395 0.152669 0.025088 0.031572Bare Land 0.105628 0.374843 0.192701 0.043621Urban -0.356923 0.038094 -0.273288 0.045284Plant 0.333387 0.786695 0.656094 0.081619Desert 0.071897 0.155663 0.100783 0.014291River 0.043469 0.429917 0.127503 0.08131【结果与分析】:通过对提取地物的DNVI值的可以发现,绿色的DNVI值比较高,原因是绿色植物叶绿素引起的红光吸收和叶肉组织引起的近红外光反射使得植被在近红外波段和红光波段有很大的差异;水体和裸地在红光波段和近红外波段反射率相当,因此水库和裸地的NDVI值接近0;雪地NDVI最低值中出现负值,是由于在近红外波段比可见光波段有较低的反射率;沙漠中植被很少,因此其近红外波段和红光波段的反射情况和裸地类似,因此其NDVI值接近于0;河流的NDVI值稍大于由于河流中存在一定的含沙量,使得地物在近红外波段的反射率大于近红外波段。
二、温度反演【温度反演】低增益温度反演高增益温度反演【第1波段部分地物低增益温度反演数据】开尔文温度摄氏温度反演温度地物min max mean stdev min max mean影像几何精纠正1.深入理解影像几何精纠正的原理2.学会使用影像对影像的几何精纠正方法和具体操作步骤1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
2. 利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正。
1. 扫描地形图(宁夏中卫地区1:25万104811.img )进行几何精纠正(包括投影参数、单位的调整将选定参考点保存)。
2. 利用纠正好的全色波段高分辨率影像完成同景多光谱影像的几何精纠正▼控制点坐标分析:总共采集控制点15个,总体误差0.3971单精度误差最大值0.8113,最小值0.0293控制点的采集符合处理要求,采集点较合适。
▼待正影像的控制点分布:▼卷帘效果Landsat影像数据下载、导入、目视解译与分析1.学会在美国NASA/USGS网上进行Landsat卫星影像的检索和下载,认识影像名编号意义。
2.初步掌握ENVI/IDL影像处理软件的使用方法,熟悉软件的用户界面、功能模块,掌握基本功能的使用。
复习遥感导论课程中的遥感影像目视解译环节,选择自己熟悉的地区,进行Landsat TM/ETM+影像的目视解译。
实习内容1、课前准备:根据自己感兴趣的地区,下载一景Landsat TM/ETM+影像。
2、将单波段分别添加在ENVI中(采用file| open image file工具,注意在添加波段数据时,热红外波段影像应另存为一个文件)3、利用Basic tools | layer stscking视窗选择多波段影像进行添加,在available中选择刚才添加的影像,选择RGB color分别选取4、3、2和7、4、2波段组合进行彩色合成,此时,可以再打开一个视窗(new display)观察地物的色调变化。
列表说明上述地物分别在两种波段组合下的颜色特征。
4、查询并记录影像文件的基本信息、投影信息,以及各个波段直方图信息。
5、将影像缩小、放大、漫游工具识别影像中的土地利用/土地覆盖类型,可能的土地利用/土地覆盖类型包括:(1) 耕地farmland(8) 公路/铁路road or railway(2) 草地grassland(9) 河流stream(3) 裸地barren land(10) 水库reservoir(4) 森林forest(11) 冰雪ice and snow(5) 城镇居民地town(12) 云cloud(6) 农村居民地village(13) 阴影shadow(7) 沙漠desert6、利用load RGB将图像显示,后用视窗中的光谱剖面工具,提取上述地物在不同波段的数值(Digital Number,DN);要求针对影像中的6种地物至少各采集10个样本,取平均值,做光谱剖面图,分析不同地物的灰度值随波段变化的特点。
将光谱剖面数据保存为文本文件。
7、利用步骤6所得数据完成每种地物的样本特征光谱统计表(附后,采用Excel数据分析工具)8、采用Excel打开步骤4生成的数据文件,从TM/ETM+ 1、2、3、4、5、7中任取两个波段的数据,做出散点图,分析上述不同地物在光谱空间中的分布特征(注意要在同一坐标系当中使用不同的符号表示不同的地表覆盖类型)。
9、从可见光/近红外的6个波段影像中选取512×512子区并保存,对影像进行分析,说明子区影像直方图与原影像直方图的差异,并分析子区影像的基本特征。
重新计算原影像统计特征值,注意观察记录结果的变化。
基本原理数据准备一景Landsat 4号星衡阳地区TM/ETM+影像LT41230411989042XXX02lines中分别输入512,获取影像,并且计算影像直方图信息。
结果与分析1、4、3、2和7、4、2波段组合进行彩色合成观察地物的色调变化。
如【图一】、【表一】【图一】【表一】farmland stream town reservoir forest432波段红色淡蓝色淡蓝色墨绿色红色742波段绿色深蓝色紫红色深蓝色绿色2、【直方图】3、【地物光谱剖面图】地物类型:城镇波段号: 1 2 3 4 5 7: 0.45~0.52 0.52~0.60 0.63~0.69 0.76~0.90 1.55~1.75 2.08~2.35 单变量统计均值86.2 35.2 35.8 29.7 41 26.3标准差 2.097618 1.032796 1.549193 2.830391 6.716481 4.547282 方差 4.4 1.066667 2.4 8.01111 45.11111 20.67778 最小值89 37 38 34 52 35最大值83 34 34 26 34 22协方差矩阵1 3.962 1.16 0.963 0.14 0.94 2.162 0.86592 13 0.872175 0.865543 14 0.098903 0.349235 0.325893 15 0.080621 0.396707 0.282051 0.840885 17 -0.28407 0.029531 -0.29394 0.441459 0.678861 1分析:森林采样点集中在(30,60)附近,说明森林在5波段DN比2波段反射大城镇采样点集中在(35,40)附近,城镇在第25波段反射相差不大农村采样点集中在(35,60)附近,说明农村在5波段DN比2波段反射大森林采样点集中在(25,65)附近,森林在5波段DN比2波段反射大河流采样点集中在(37,10)附近,说明河流在5波段DN比2波段反射明显弱阴影采样点集中在(20,20)附近,阴影在第25波段反射相差不大4、子图选取。
分析:1、子区平均值与原图像基本相同。
2、子区图像DN值波峰与原图在曲折变化上有一些差别,原图像地物类型较多。
3、子区与原图像对应区域在显示上颜色有变化。
存在问题与解决办法1、实习过程中影像的分辨率比较低,有些地物的辨别有点小困难,有些地物之间不好辨别。
* *。