ENVI预处理GF-1号WFV1遥感数据提取NDVI指数

ENVI预处理GF-1号WFV1遥感数据提取NDVI指数

目录：

  1. 遥感数据获取
  2. 遥感数据预处理
    2.1编辑头文件
    2.2辐射定标
    2.3FLAASH大气校正
    2.4正射校正
    2.5计算NDVI
    2.6图像裁剪

提取结果如下：
在这里插入图片描述

遥感数据获取

https://data.cresda.cn/#/home

打开陆地观测卫星数据服务网站，注册登录之后
在这里插入图片描述
选择高分1号16m数据，搜索河北省廊坊市大厂回族自治县，采集时间按实验需要（我用的2023年5月26日，后面需要实地测量数据验证，因此时间上选对应的）
在这里插入图片描述
选择一个能够覆盖研究区域的图像，加入购物车，提交订单等待网站邮件，然后下载。

遥感数据预处理

编辑头文件

在这里插入图片描述
下载好后，先将tiff格式的遥感影像拖动到ENVI中，右键点击选择View Metadata
在这里插入图片描述
打开Spectral一栏，主要修改红圈中的四个传感器参数，点击Edit Metadata
在这里插入图片描述
找到对应的几个修改参数一栏，如果没有的话，在上方Add中添加。高分1号的Wavelength和FWHM、Offset可以参考上面我图中的参数，第三列的Radiance Gains可以参考下表。
在这里插入图片描述
选择2022年的WFV1传感器参数。具体见下面网站

https://blog.csdn.net/weixin_56529683/article/details/129276880

辐射定标

流程也可以参考如下网站

https://www.cnblogs.com/enviidl/p/16553338.html

在这里插入图片描述
在右侧工具栏中，找到Radiometric Calibration
在这里插入图片描述
在这里插入图片描述

FLAASH大气校正

具体过程可以参考下面网站

https://www.cnblogs.com/enviidl/p/16553338.html

在这里插入图片描述
工具栏中打开
在这里插入图片描述
输入文件找到刚刚辐射定标处理后的影像
在这里插入图片描述
打开后点击下面的选项，默认为1就好
在这里插入图片描述
然后输出路径和文件夹也设置一下，这里有个注意事项，就是输入输出的路径名当中一定不要包含中文，路径最好也不要太长，可以在D盘下面专门放一个文件夹来保存ENVI处理的过程结果
在这里插入图片描述
一些参数设置
在这里插入图片描述
打开Multispectral Settings设置波谱响应函数，选好后点Ok
在这里插入图片描述
点击Advanced Settings，上面的参数设置为100，然后点击Ok，退回到FLAASH参数设置界面，点击左下角的Apply选项。开始大气校正，电脑过程中可能会无响应，等待一段时间即可。

正射校正

具体参考下面网站

https://zhuanlan.zhihu.com/p/400902530

在这里插入图片描述

计算NDVI

过程参考网站

ENVI提取NDVI与植被覆盖度估算 - 知乎 (zhihu.com)

在这里插入图片描述
打开Band Math,输入公式(float(b5)-b4)/(b5+b4)
在这里插入图片描述
点击Add to List
在这里插入图片描述
选中，并点击Ok
在这里插入图片描述
B4选择0.6650的红光波段，B5选择0.8的近红外波段,设置好输出路径后，点击Ok
在这里插入图片描述
得到如下NDVI影像。

图像裁剪

首先需要有自己研究区域的shape文件，该文件如何制作参考我的文章
在这里插入图片描述
拖动大厂回族自治县的shape文件到ENVI中，自动生成红线内的区域

在这里插入图片描述
工具栏中，找到Subset工具
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
点击Cursor Value工具，某点的Data值就是该点的NDVI计算结果
在这里插入图片描述
在这里插入图片描述
最后将该处理结果保存为TIFF格式影像。