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格式影像。