哨兵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

发布了21 篇原创文章 · 获赞 8 · 访问量 1313
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 数字20 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览