通过上篇文章的介绍,大家应该已经下载到自己感兴趣区的数据了。本节会介绍拿到原始数据后,如何对其进行预处理操作,从而得到可以进行后续计算的L2A级遥感数据(经过辐射定标和大气校正)~
欧空局针对哨兵2号数据的预处理专门开发了Sen2Cor插件,下载它在cmd命令行中简单输入几行代码就可以完成预处理工作,非常方便快捷~
1.下载Sen2Cor插件,网址如下
Sen2Cor | STEPstep.esa.int下载最新版本,文件是个zip文件。
2.插件安装
将zip文件下载到用户文件夹中,在cmd中加载该文件夹。
注意:文件夹名中别出现空格,否则cmd会无法识别,后续会比较麻烦。
双击文件夹中L2A_Process.bat文件
在cmd中输入L2A_Process.bat --help,若出现以下内容,则为配置成功
3.预处理操作
在cmd中输入L2A_Process +原始图像绝对路径,如下图所示,如果出现以下内容则表明预处理操作开始执行,耐心等待其完成即可(可能需要二十分钟左右~)。
注意:1.预处理的结果会和原始图像存储在一个文件夹中,大多数情况下是无需更改的。如果想要自定义输出路径,可在该句代码后加上“--output_dir = 你的自定义路径”;2.有些教程中提到要在语句后加--refresh,但在Sen2Cor2.8版本以上不用加该句,加了会报错。
这里是使用最简单的操作语句,是对所有13个波段都进行操作。若想只对某分辨率数据进行操作(如只对10m分辨率的数据进行预处理),则通过如下代码就可以进行
通过输入L2A_Process --help,可以查看帮助文档,进行更为个性化的操作
4.结果检查
打开SNAP软件,加载原始数据和预处理后数据。通过spectrum view工具(在搜索栏中搜索即可调出),获取像元值的光谱曲线。
下图是原始数据森林上某像元的光谱曲线和预处理后该像元的光谱曲线,可以看到有很大差距。通过比对标准植被光谱曲线可知,原始数据植被光谱曲线明显错误,而预处理后光谱曲线正确。
这表明原始数据已经过正确的辐射定标和大气校正。
补充(转换成ENVI格式)
ENVI5.1无法直接打开Sentinel数据,但大家对ENVI操作会更加熟练,这里补充说明下如何用ENVI5.1打开Sentinel数据,主要就是重采样并转换格式。
打开重采样工具
重采样后的结果就可以导入ENVI5.1中显示了,但其存储方式是一个波段一个文件,在ENVI中融合一下(ENVI中搜索layer stacking工具)就可以得到多波段文件了。
到此,大家就已经学会了Sentinel-2数据预处理啦,下期文章会介绍如何使用预处理后数据进行波段运算(植被指数的计算)。