【RS】Sentinel-2哨兵二号L1C数据预处理教程(Sen2Cor大气校正、SNAP重采样、ENVI波段组合)
刚分享过新版欧空局的数据下载教程,正好把哨兵2号预处理的教程也分享一下,主要就是使用官方插件Sen2or对L1C级数据进行大气校正,使用SNAP软件对L2A级数据进行重采样并导出ENVI可以打开的格式,最后使用ENVI对重采样后的数据进行波段组合,以便于后期的定量分析。
原创作者:
博客地址:
1.软件下载
1.1 Sen2Cor下载
用途:L1C级数据大气校正生成L2A级数据
网站:
百度云盘:
1.2 SNAP下载
用途:欧空局官方公布的哨兵数据处理软件,这里主要使用其波段重采样以及波段导出ENVI格式。
网站:
百度云盘:
2.Sen2or大气校正
2.1 单文件处理
使用命令提示符,Win+R,输入cmd(命令提示符)。使用命令cd /d 路径,将目录切换至Sen2Cor的目录。如下:
输入命令L2A_Process.bat 影像路径,注意这里的影像路径必须子目录就是文件,而不是还有一层,解压的时候可能会生成两级目录,切记!!!还有就是不能有中文路径!!!如下:
等待处理成功,返回Successfully即可。
2.2 批量处理
使用命令提示符,Win+R,输入cmd(命令提示符)。使用命令cd /d 路径,将目录切换至Sen2Cor的目录。(与单文件一样,目的是调用L2A_Process插件)
输入命令for /D %s in (G:\1\S2A_MSIL1C*) do L2A_process %s,其中G:\1\是影像的路径,S2A_MSIL1C*是指路径中所有以S2A_MSIL1C开头的文件,如果是S2B数据,就改下S2B_MSIL1C*即可。如下:
3.SNAP重采样
3.1 打开影像
在SNAP中,点击左上角的File,使用Open product打开L2A数据文件夹中的MTD_MSIL2A.xml文件即可。
3.2 重采样
打开Optical>Geometric>S2 Resampling Processor工具,设置导出格式为ENVI,修改导出路径,修改重采样的分辨率为10m,run即可。
4.ENVI波段组合
4.1 加载波段
我们将重采样完成后文件夹中的.img波段文件拖入ENVI打开。
4.2 波段组合
ENVI5.6以上版本使用Build Layer Stack工具,ENVI5.3使用Layer Stacking工具,将波段按顺序依次倒入,坐标系通过影像读取,输出即可!
5.总结
以上即为Sentinel2数据的预处理流程,如果是从官网下载的L2A级数据,就直接重采样+波段组合即可,无需进行大气校正!此外ENVI5.6过后支持打开Sentinel2数据,所以SNAP导出ENVI格式主要是用来ENVI5.3的,我们这里使用SNAP是用来重采样。
有什么问题直接留言沟通即可,后续我可能会更新使用Python一键批量预处理L1C级数据。感兴趣的话可以关注一下。