Sentinel-2数据下载及处理

数据下载网站：欧空局官网（需注册并登录）

https://scihub.copernicus.eu/dhus/#/home

 哨兵2 L1C数据波段信息

哨兵2 L1C数据时间：2015-06-23至now

        由于数据量大，考虑服务器压力，哨兵2号数据直接下载的时间跨度为一年半至今的数据，offline数据需要提交申请才可下载

• 哨兵2数据检索与下载

1.进入网站并登录

2.通过鼠标左键和滚轮将区域缩放/扩大到感兴趣区域

 

 3.输入要检索的时间和传感器类型，进行搜索

 

 Ingestion Date为影像进入数据库的时间，而Sensing Date为影像拍摄的成像时间，建议选择Sensing Date

卫星运行方式的选择（升轨和降轨） ，可不设置选择

period可设置时间段内的筛选 

勾选Sentinel-2

可选择星S2A或者星S2B，不选则默认两颗星都被检索，建议默认

 

数据类型：

1C级数据（S2MSI1C） ：Level-1C 产品是使用数字地面模型(DEM)将图像投射到制图坐标中的结果，为经正射校正和亚像元级几何精校正后的大气表观(TOA)反射率产品。

2A级数据（S2MSI2A）：是经过Sen2Cor进行大气校正和正射校正后的产品，Level-2A 的主要输出是正射影像大气底部(BOA)校正的反射率产品，2017-03-28之后才生产。

可写入轨迹带号

云量控制百分比，必须按照例子的语法写入 

最后设置完成，点击搜索，可得到检索的影像 

4.点击下载选项即可进行下载，或者放入购物车进行集体下载

• 哨兵2数据预处理

数据处理需要使用官方指定的工具和插件，下载网址如下：

SNAP工具：http://step.esa.int/main/download/snap-download​​​​​​/

Sen2Cor插件：https://step.esa.int/main/snap-supported-plugins/sen2cor/sen2cor-v2-11/

 1.文件名称解读

S2A 为哨兵-2A卫星拍摄的影像，MSI表示传感器（MutiSpectal instrument），L1C表示产品的等级，20150716为拍摄图像的时间，N0204处理基线编号，R050为对应的轨道编号（R001-R143），T35UQR为拼接域编号，SAFE为欧洲标准存档格式

1）单文件大气校正处理 

1.1 找到插件压缩包解压后所在的位置，输入cmd进入到指令窗口

1.2  输入处理命令，L2A_Process.bat+“下载的L1C的位置”

L2A_Process.bat D:\SNAP_process\download_sentinel2\S2A_MSIL1C_20150802T090006_N0204_R007_T35UQR_20150802T090006.SAFE

需要等待很久的时间才可处理完成

 结果显示成功

最后生成的L2A文件会和L1C在同一目录下 

2）批量处理

将sentinel2数据放在一个文件夹下，但是S2A和S2B数据要分开处理，即分成两个文件夹处理。

同样地，先找到插件压缩包解压后所在的位置，再输入cmd进入到指令窗口，再按照如下规则输入代码，注意（）中填写L1C数据所在地址。

for /D %s in (D:\SNAP_process\download_sentinel2\S2A_MSIL1C*) do L2A_process %s

注意：也可以通过输入L2A_Process --help，可以查看帮助文档，进行更为个性化的操作。

如：可以指定处理某一种分辨率（10、20、60m）的band影像，需在代码后面继续写入 --resolution 10

对于处理为20m和60m分辨率的2级数据缺失band 8和band 10的问题，官方文件的Chapter 2.4也给出了解释（见http://s2tbx.telespazio-vega.de/sen2cor/sen2cor-sum-2.0.pdf）

• band 8 missed

Consider the similarity of band 8 and band 8a. They are located at the same spectral position, however band 8 is at 10 m resolution and therefore requires a broader spectral halfwidth It makes little sense resampling this band to 20 m or 60 m when a channel with a smaller bandwidth is available for the same spectral position in that spatial resolution.

• band 10 missed

Sentinel-2 Band 10 is not present in the output of Sen2Cor because “No output generated as it does not contain surface information”.