哨兵2号（Sentinel2）卫星数据处理

李国春

哨兵2号数据是一个重要的共享免费数据源，原来我就知道这个数据的应用很受重视，可是没想到现在在国内是如此空前的重视。哨兵数据处理人家ESA有自己的SNAP处理软件，非常好的软件。RSD也来凑凑热闹，也加了处理哨兵2号卫星数据的功能。

介绍RSD处理哨兵数据之前先夸夸这数据，数据质量真好。但是最让我感动的还是人家的元数据，写的那是真好啊，看来值得我们学习的东西还有很多。好了，言归正传。

Sentinel2使用的是JPG2000压缩，基于小波的，一切都好就是解压稍稍慢了一点，RSD加载一个数据集大概需要3分钟。但是同学们反映说还挺快的，那就好，你们说好就好。

一 加载L1C数据

加载L1C数据，从图1打开。

图1

点击后需要打开一个 .xml，Sentinel2的L1C有好几个.xml，解压后的目录进去，找“MTD_MSIL1C.xml”打开。

打开后出来一个对话框（图2），选择你需要打开的3种10m、20m和60m数据的一种打开。

图2

 

二 加载L2A数据

加载L2A数据，从图1的L2A 10m、20、60m的菜单里面选择。点击后选择打开“MTD_MSIL2A.xml”文件。无论你选择哪个分辨率的数据都是打开这个文件，不用担心。

 

三 向现有任务添加L1C、L2A数据

如果已经有了一个正在处理的任务就可以直接添加，如果没有当前任务，可以打开一个.tpl文件或者重新手动创建一个。然后按图3添加数据。

图3

数据的选择和图2是一样的。

向现有任务添加L2A数据和前面添加L1C数据也一样。

 

四 两个数据集拼接

这里举例盘锦市。先打开一个“盘锦市-UTM-10m.tpl”，然后按图3添加Sentinel2 L1C数据。添加2个数据集进去，见图4。

 

图4

可见两个数据集都添加了一部分进去，多余的部分被切掉了。由于云量的不同造成两景数据成像非常不一致，没关系，不影响数据。加载时间都不到3分钟。

 

下面开始拼接。

勾选层1，按图5选择菜单命令。

图5

在对话框(图6)中选在层2。这里的意思时将层1归并到层2里面，交集部分保留层2的数据。见：

《遥感数据处理流程新概念（一） ——RSD的面向任务》

https://blog.csdn.net/gordon3000/article/details/93171572

《遥感数据处理流程新概念（二） ——再谈拼接》

https://blog.csdn.net/gordon3000/article/details/93974798

 

稍等10几秒，拼接就完成了，结果见图6。

图6

现在层3就是一个拼接好的完整的层。

 

五 剪切一下

你要是打开的“盘锦市-UTM-10m.tpl”，里面直接就有盘锦市的蒙板，否则你要打开一个.shp，从这个矢量数据创建一个蒙板。

右击窗口蒙版区域（绝大部分都在蒙版区域内），出现下面的弹出菜单，选“复制当前蒙板区域”。稍等，你会发现多了一个层4。关闭（不显示，或者删除）其它层，只留层4，就是我们希望的结果（图7）。

 

图7

很多同学需要从L1C进行大气校正，人家ESA有工具，或者直接用L2A也可以。其实RSD也可以加上大气校正要是有需求，看吧。

 

现在新版的RSD更新到3.0.8(2019 06 28), 前面版本的RSD处理哨兵2数据有个大BUG，UTM投影的伪东500000后面多加了一个0，我说大家怎么都没用上呢。

 

软件下载地址：https://pan.baidu.com/s/1T-LBvaD_zVCwJsGf_hCyCg

详情加企鹅群136965427，在这里解答和讨论有关遥感数据处理和RSD平台的有关技术问题。