哨兵2号(Sentinel2)卫星数据的处理
-
哨兵2的下载渠道有很多,如USGS、欧空局等等。下载下来的数据为Level-1C数据,这是没有经过大气校正和辐射校正的数据。需要我们自己进行处理,得到Level-2A的数据。
-
哨兵Level-1C往往是通过Sen2cor工具完成大气校正和辐射校正,得到2A数据,Sen2cor工具的官方下载地址如下:
- http://step.esa.int/main/third-party-plugins-2/sen2cor/
-
根据自己的电脑系统下载对应的版本,我的电脑是windows64位系统,就以此为例,下载下来的Sen2cor工具是一个压缩包,解压到一个非中文路径的文件夹内即可(防止出现奇奇怪怪的错误),解压后的界面如下图,对遥感影像的处理用到的是L2A_Process.bat工具
---- -
对于该工具的使用需要通过命令行,主要有两种方法,第一种是通过cd命令进入该批处理文件所在的根目录进行调用,第二种是把该文件所在的目录加入电脑的环境变量之中(网上添加环境变量的教程很多,这里不再赘述)。这里以第一种方式为例。
- > 如果输入bash L2A_Process --help
会弹出L2A_Process.py文件帮助,说明Sen2cor工具正常使用,对于单个文件的代码:
L2A_Process F:\Sentinel-2\2020\S2A_MSIL1C_20200503T023551_N0209_R089_T50RPV_20200503T053405.SAFE #文件所在目录
-
只处理对应分辨率波段在后面加上相应参数即可,如只处理10m分辨率的波段(建议全部处理,可能会出错,而且处理时间差不多)
L2A_Process F:\Sentinel-2\2020\S2A_MSIL1C_20200503T023551_N0209_R089_T50RPV_20200503T053405.SAFE --resolution 10 #处理10m分辨率的波段
-
下面是Sen2cor工具处理的效果图,一般一幅影像处理所花费的时间在半个小时左右
-
如果需要对多幅哨兵2影像进行校正,可以利用Bash命令遍历处理代码如下:
for /D %s in (F:\Sentinel-2\2020\S2A_MSIL1C*) do L2A_process %s