【RS】Sentinel-2哨兵二号L1C数据预处理教程（Sen2Cor大气校正、SNAP重采样、ENVI波段组合）

        刚分享过新版欧空局的数据下载教程，正好把哨兵2号预处理的教程也分享一下，主要就是使用官方插件Sen2or对L1C级数据进行大气校正，使用SNAP软件对L2A级数据进行重采样并导出ENVI可以打开的格式，最后使用ENVI对重采样后的数据进行波段组合，以便于后期的定量分析。

原创作者：RS迷途小书童

博客地址：https://blog.csdn.net/m0_56729804?type=blog

1.软件下载

1.1 Sen2Cor下载

用途：L1C级数据大气校正生成L2A级数据

网站：Sen2Cor – STEP (esa.int)

百度云盘：Sen2Cor– RS迷途小书童

1.2 SNAP下载

用途：欧空局官方公布的哨兵数据处理软件，这里主要使用其波段重采样以及波段导出ENVI格式。

网站：SNAP Download – STEP (esa.int)

百度云盘：SNAP– RS迷途小书童

2.Sen2Cor大气校正

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*即可。成功之后也是返回sucessfully。如下：

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级数据。感兴趣的话可以关注一下。