【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演

 一、专题概述

  • 植被覆盖度是指植被(包括叶、茎、枝)在地面的垂直投影面积占统计区总面积的百分比。
  • 植被覆盖度常用于植被变化、生态环境研究、水土保持、气候等方面。
  • 植被覆盖度数据来源于地理遥感生态网平台。

二、像元二分法模型

  • 像元二分模型是一种简单实用的遥感估算模型,它假设一个像元的地表由有植被覆盖部分地表与无植被覆盖部分地表组成,而遥感传感器观测到的光谱信息也由这2个组分因子线性加权合成,各因子的权重是各自的面积在像元中所占的比率,如其中植被覆盖度可以看作是植被的权重。

VFC = (S - Ssoil)/ ( Sveg - Ssoil) S为遥感信息

改进的像元二分法模型(2004)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第1张图片


三、处理流程图

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第2张图片


四、图像预处理

1.1 数据读取与定标

  • File-> Open As-> Landsat-> GeoTIFF with Metadata,打开_MTL.txt。(数据读取)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第3张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第4张图片

  • 查看影像数据集(五个等级)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第5张图片

打开矢量数据(北京)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第6张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第7张图片

查看源信息(已经做过地理坐标)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第8张图片

  • 辐射定标:Toolbox/Radiometric Correction/ Radiometric Calibration,选择可见光-红外组(7个波段),打开Landsat定标工具

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第9张图片

选择多光谱数据

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第10张图片

  • 裁剪研究区区域

打开空间裁剪(Spatial Subset),利用所示工具2(Subset by Vector)选择研究区(北京市区Shp)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第11张图片

如下图,自动用矢量边框确定范围(可减少数据量)——对于小部分区域不要先裁剪再定标(导致定标参数缺失),在定标时可用空间子集方法

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第12张图片

Apply FLAASH Settings,保存路径即可

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第13张图片

定标完成

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第14张图片

  • 同样方式对另一个幅影像数据进行定标操作

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第15张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第16张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第17张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第18张图片

两幅影像数据完成定标(北京区域)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第19张图片

1.2 图像镶嵌

  • 当研究区超出单幅遥感图像所覆盖的范围时,通常需要将两幅或多幅图像拼接起来形成一幅或一系列覆盖全区的较大的图像。
  • ENVI的镶嵌功能
    • 基于像素镶嵌和无缝镶嵌工具
    • 自动颜色平衡,边缘直方图匹配
    • 自动生成接边线、接边羽化
    • 提供结果预览功能
    • 运算速度快
    • 占用非常少的虚拟内存空间
  • 两景相邻的OLI数据,镶嵌得到包含北京市区的完整范围/Mosaicking/Seamless Mosaic

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第20张图片

选择数据(完成裁剪的北京区域影像)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第21张图片

Mosaic拓展工具介绍

  • Order工具栏中可以调整图像上下顺序(图层粉色区域为当前选中区域)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第22张图片

  • 颜色平衡工具(直方图匹配——会修改DN值)

基于重叠区域基于完整区域)不同工具产生效果不一,应当选项产生较好质量

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第23张图片

  • 接边线设置(如果不选择会根据图像外边框自然接边),可手动自行绘制

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第24张图片

如下图自动生成了接边线(绿色)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第25张图片

可编辑接边线

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第26张图片

如下图,编辑的接边线(红色)需要把绿色线包含起来,然后应用

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第27张图片

如下图,可查看前后对比

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第28张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第29张图片

在这里不应用,只是示例

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第30张图片

输出Mosaic(镶嵌)文件

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第31张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第32张图片


注:扩展工具(可进行大量数据的快速镶嵌

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第33张图片


1.3 图像裁剪

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第34张图片

打开图像裁剪工具,如下图

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第35张图片

选择镶嵌好的栅格影像

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第36张图片

进行设置,如下图

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第37张图片

北京市影像裁剪完成

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第38张图片

1.4 大气校正——概述

  • 获得地物真实的光谱信息是遥感反演的前提,大气校正可消除大气影响,还原地物的真实信息,是定量遥感数据预处理中必不可少的环节。

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第39张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第40张图片

  • /Radiometric Correction/FLAASH Atmospheric Correction(选择大气校正工具)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第41张图片

  • 设置存储面板

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第42张图片

  • 选好的默认设置

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第43张图片

  • 查看影像源信息,对大气校正剩余面板进行设置补充

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第44张图片

  • 多光谱工具栏设置

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第45张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第46张图片

  • 应用,完成大气校正

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第47张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第48张图片


五、植被覆盖度反演

2.1 NDVI计算

  • Spectral/Vegetation/NDVI(打开NDVI工具栏)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第49张图片

选择L8北京大气校正数据

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第50张图片

  • 设置NDVI波段(查看源数据选择)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第51张图片

如下图,完成NDVI计算

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第52张图片

  • 快速统计(查看影像NDVI信息)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第53张图片

NDVI应该在-1—1之间(如下图,存在异常值)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第54张图片

  • 处理异常值(打开波段计算工具栏)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第55张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第56张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第57张图片

  • 如下图,NDVI异常值已经消除

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第58张图片

2.2 掩膜文件制作——模型参数求解

  • 既使同一景影像 ,对于NDVIsoil 和NDVIveg 不能取固定值,两个值会随土壤类型和土地利用类型的变化而变化。
  • 利用土地利用类型生成不同土地覆盖类型对应的最大值和最小值
  • 打开数据(土地利用覆盖影像)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第59张图片

土地分类图像

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第60张图片

  • 制作掩膜文件

选择处理过异常值的NDVI文件

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第61张图片

选择Mask

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第62张图片

以范围做掩膜

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第63张图片

选择土地分类文件

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第64张图片

  • 林地掩膜

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第65张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第66张图片

  • 田地掩膜

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第67张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第68张图片

  • 城市掩膜

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第69张图片

  • 水域掩膜

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第70张图片

  • 其它类型掩膜

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第71张图片

  • 统计掩膜的各类型NDVI的最大小值

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第72张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第73张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第74张图片

  • 以像元个数达到五位数为界限

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第75张图片

0.38(最小值)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第76张图片

0.86(最大值)

同样方法,

农用地 0.16(最小值) 0.77(最大值)

城市 0.04(最小值) 0.47(最大值)

水域 -0.33(最小值) 0.71(最大值)

其它 0

  • 设置参数文件

b1*0.3804+b2*0.1638+b3*0.0405+b4*0+b5*0.098

  • 在波段计算中设置(最小值)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第77张图片

选择波段对应掩膜文件

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第78张图片

  • 最大值,如下

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第79张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第80张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第81张图片

六、植被覆盖度估算

  • 波段运算的过程 表达式为:(b1-b2)/(b3-b2)
  1. b1—NDVI(去除异常值)
  2. b2—NDVIsoil
  3. b3—NDVIveg
  • 植被覆盖度常用于植被变化、生态环境研究、水土保持、沙漠化、气候变化等方面。
  • 进行波段运算,如下

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第82张图片

选择对应波段

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第83张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第84张图片

  • 快速统计(如下图,存在异常值)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第85张图片

  • 处理异常值

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第86张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第87张图片

统计(0—1正常)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第88张图片

  • 最终影像色带调整

分类(为使影像变得美观清晰)——添加色带

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第89张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第90张图片

添加色带成功

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第91张图片

分割(也可对影像进行分割——条带数目与范围可自行设置)

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第92张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第93张图片

【ENVI精讲】处理专题五:基于像元二分模型的植被覆盖度反演_第94张图片

   赶紧三连关注下, 植被覆盖度指数获取途径如下:

你可能感兴趣的:(ENVI遥感处理教程,地理信息数据获取,植被覆盖度,envi教程,遥感反演,NDVI数据,npp数据)