哨兵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平台的有关技术问题。

 

 

 

 

发布了74 篇原创文章 · 获赞 27 · 访问量 5万+

猜你喜欢

转载自blog.csdn.net/gordon3000/article/details/93992366
今日推荐