资源描述
单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,专题:基于遥感的草原与沙漠化监测,技术支持邮箱:,ENVI-IDL,技术支持热线:,400-819-2881-7,官方技术博客:,、,大气校正,植被,反演,、,植被覆盖变化监测,使用模块和工具:,ENVI,主模块,FLAASH,大气校正扩展模块,IDL,开发的环境小卫星数据读取补丁,波段运算模块,专题概述,处理流程,数据预处理,植被覆盖度反演,植被覆盖变化监测,流程说明,根据环境小卫星,CCD,数据特点及草原植被变化监测的要求,采用,以下处理流程,:,一、,数据预处理:,1.CCD,数据读取,2.,辐射定标,3.,大气校正,4.,研究区裁剪,二、,反演模型建立,1.,归一化植被指数,2.,植被覆盖度,三、植被变化监测,1.,植被覆盖区提取,2.,植被变化检测,四、后期处理与应用,环境小卫星,CCD-1B,数据,在环保部卫星环境应用中心免费下载获取,环境小卫星,CCD,数据下载地址:,TM,数据下载地址为:,FilesExelisENVI52extensions,目录,下,数据读取和定标,主菜单,-File-Open As-HJ-1A1B,2.1,数据读取和定标,由于整景数据范围,比较,大,,所以在做大气校正前,,,先,将,浑善达克,以及周边区域裁剪出来。,2.2,工程区裁剪,做变化监测,两个时相的数据必须互相配准,2.3,图像配准,获得地物真实的光谱信息是遥感反演的前提,大气校正可消除大气影响,还原地物的真实信息,是定量遥感数据预处理中必不可少的环节。,2.4,大气校正,概述,数据准备,FLAASH,对图像文件有以下几个要求:,数据是经过定标后的辐射亮度(辐射率)数据,单位是:(,W,),/,(,cm2*nm*sr,)。,数据带有中心波长(,wavelenth,)值,如果是高光谱还必须有波段宽度(,FWHM,),这两个参数都可以通过编辑头文件信息输入(,Edit Header,)。,数据类型,支持四种数据类型:浮点型(,floating,)、长整型,(long integer),、整型(,integer,)和无符号整型,(unsigned int),。数据存储类型:,ENVI,标准栅格格式文件,且是,BIP,或者,BIL,。,波谱范围:,400,2500nm,2.4,大气校正,FLAASH,大气校正,基本参数设置,2.4,大气校正,FLAASH,大气校正,多光谱参数设置,单击,Multispectral Setting,按钮,在,Filter Function File,导入光谱响应曲线“环境,1A,星,CCD2,光谱响应,.sli,”,单击,OK,;,2.4,大气校正,FLAASH,大气校正,高级参数设置,2.4,大气校正,FLAASH,大气校正,结果检验,大气校正前后典型地物(植被)光谱曲线对比,2.4,大气校正,结果检验,预处理的最后一步是用浑善达克区,矢量数据,对已经做过辐射定标、大气校正的,HJ-CCD,数据进行不规则,裁剪,矢量数据:,“,4-,研究区裁剪,浑善达克矢量数据”,2.5,裁剪浑善达克地区,3,、植被覆盖度反演,归一化植被指数计算,应用被植被强吸收的红光波段(环境星第,3,波段)和被植被强反射的近红外波段(环境星第四波段)计算归一化植被指数,。,植被覆盖度计算,FC=(NDVI-NDVISoil),(NDVIVeg-NDVISoil),植被覆盖度反演,归一化植被指数,植被覆盖度,打开,ENVI 5,工具箱,/Band Ratio/Band Math,工具,3.1,归一化植被指数计算,波段运算,打开,ENVI 5,工具箱,/Band Ratio/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),3.2,植被覆盖度计算,4,、植被变化监测,植被覆盖度反演,2006,土地利用分类图,2006,植被覆盖区,2009,植被覆盖数据,2009,植被覆盖区,波段运算,植被荒漠化监测,环境小卫星数据,2009,年植被覆盖度大于,0.3,的为植被覆盖区,2006,年土地利用分类图,DN,值,=1,、,2,、,3,为植被覆盖区,4.1,两图像植被覆盖区域提取,土地利用分类类型,影像,DN,值,耕地,1,林地,2,草地,3,水域,4,城市、工矿、居民用地,5,未利用地,6,沙地,61,2009,年,8,月植被覆盖区,(b1 le 0.3)*0+(b1 gt 0.3)*1,b1,为,20090811,植被覆盖度,图像,,0.3,为经验值,2006,年,8,月植被覆盖区,(b1 ge 1 and b1 le 3)*1+(b1 lt 1)*0+(b1 gt 3)*0,b1,为,浑善达克,2006,年,8,月土地利用分类图,4.1,两图像植被覆盖区域提取,波段运算,工具箱,/Band Ratio/Band Math,工具,运算表达,式:,b1-b2,b1:2009,年,8,月的植被覆盖区图像,b2:2006,年,8,月的植被覆盖区图像,4.2,变化区域计算,5,、成果后处理及应用,变化图像:,0,未发生变化的区域,-1,植被荒漠化区域,1,植被增加区域,将背景值掩膜为非,0,、,-1,、,1,的值(如:,2,),5.1,植被变化区域图的背景值处理,对,-1,、,0,、,1,值以及背景值分别进行密度分割,5.2,植被变化区域制图,本专题,从原始的,环境小卫星,CCD-1A,图像,开始,,反演植被覆盖图,、并与前一时相的土地利用图进行差值比较,提取植被发生变化的区域,专题涉及,环境小卫星的数据读取、辐射定标、,图像配准、,大气校正、,植被覆盖,遥感反演过程、植被覆盖变化监测等内容。,掌握知识:,国产卫星数据的处理,多时相数据变化监测的条件,数据处理各个环节在,ENVI,下的操作及参数选择,专题总结,
展开阅读全文