哨兵2A数据的处理--重采样、波段融合及植被指数计算

1、为什么按10米分辨率重采样后，文件缺少第10波段？

借鉴答案：因为band1用来监测气溶胶，band9监测水蒸气，band10用来监测洋流等。这3个波段空间分辨率较低，如果你重采样成10M的，默认为这3个波段用不到。

2、植被指数按着公式计算为什么错误？

例如：
NDVI=（NIR-R）/(NIR+R)
直接做运算，会出现错误结果，数值很大
当图像为辐射值时，应该变为
NDVI=float(（NIR)-R）/(NIR+R)
计算结果才会是在-1-1之间，正确
检验结果查看：
右击NDVI结果数据–选择Quick Stats–查看范围是否在-1-1之间
3、对于多波段的哨兵数据，重采样用10米分辨率，虽然丢失了第10波段，但问题不大
4、在ENVI中合成（Layer stacking）合成多波段图像
利用(Band math)工具计算植被指数

5、改变颜色

NDVI计算结果：
右击数据—change color table–more–选择Green/white LINEAR色度条–revers（反转颜色深度，越深浓度越大）

6、对NDVI结果二值化

设置阈值为0.35
(b1 ge 0.35)*1+(b1 lt 0.35)*0
注：ge是大于等于（greater equal）， lt是小于（less than）。0.35是阈值，这个可以通过比对NDVI值和真彩色图像人为确定。该公式将大于等于0.35的设为1，小于0.35的设为0，从而实现了二值化操作。其中b1为你的NDVI结果数据。
得到结果：

注：白色区域为植被，黑区域为非植被

7、波段运算中遇到开根号，如何解决？

8、计算CTVI植被指数

不确定其范围，部分数据大于1，有待查证
9、差值植被指数DVI数值很大，待验证
10、

借鉴文章：

1、哨兵2号（Sentinel-2）介绍、下载和预处理、批处理
https://blog.csdn.net/qq_41718357/article/details/83536322?depth_1-utm_source=distribute.pc_relevant.none-task&utm_source=distribute.pc_relevant.none-task
2、【Sentinel-2】（3）：波段运算——聊聊植被指数
https://zhuanlan.zhihu.com/p/75880651
3、ENVI下植被指数模型详解
http://blog.sina.com.cn/s/blog_764b1e9d01010n32.html