李国春 2021 10 11

哨兵2 （Sentinel-2）数据广受欢迎，数据质量好，还免费。人家欧空局有自己的处理软件，也有控制台命令行的批量处理。RSD也来凑凑热闹沾个光，指不定有人喜欢不同的操作风格，愿意使用RSD处理Sentinel2数据呢。 :)

上一篇曾经介绍过RSD处理Sentinel2数据《哨兵2 （Sentinel2）卫星数据处理》，受到大家给与的肯定和鼓励，这里再介绍一下如何用RSD批量处理Sentinel2数据（RSD软件有更新）。

1. 批量处理10m Sentinel2 L 1C数据

首先定义一个框架，要求像元分辨率10m。定义框架参阅RSD操作手册定义框架部分。或者如果你对RSD框架不熟悉，打开一个该区域的数据，打开数据后就会自动产生一个框架。在《高分6 宽幅（GF6-WFV）卫星数据快速批量处理》有一个建立16m框架的例子，将16m改成10m即可。

打开框架后，执行菜单命令 任务管理->添加其它卫星观测数据->添加 Sentinel 数据-> Sentinel2 L1C /批量，出现图1的对话框。

 图1 指定批量加载/选择分辨率

选择10m分辨率数据，并勾选从文件列表加载。点击确定。

在出现图2的对话框，待搜索的文件名为 MTD_MSIL1C.xml，这是Sentinel2 L1C数据的属性文件名。再指定一个开始搜索的路径，凡是在这个路径下的各级子目录中所有的Sentinel2 L1C数据集都将被搜索。点击 开始搜索 。搜到的数据集显示在图2中间的区域。

图2 Sentinel2 L1C数据集搜索结果

在图2点击确定，出现图3的输出参数对话框。

图3 Sentinel2 L1C数据集输出参数

在图3中选择需要导出的文件类型，在下面指定输出路径，点击确定，即开始处理。处理结果见图4。

图4 Sentinel2 L1C数据集处理结果

22个Sentinel2 L1C数据集的10m数据处理用时15分48秒。这个用时偏少不具代表意义，一是框架范围小，再就是10m数据本来处理起来就快一些。Sentinel2 L1C数据集的20m数据处理起来更费时一些。

再查看处理过程中导出的数据（图5）。

图5 Sentinel2 L1C数据集处理后导出的数据

可见，在指定的输出路径下，为每个数据集创建了一个独立的目录，每个目录里面存放了.rsd、.img和.tif三种导出文件。

注意：有用户频繁反映导出的tif文件格式不正确，使用其它软件打开tif文件时删除那个.hdr文件，那不是为tif准备的。其它软件把它当成tif的头了。

2. 批量处理20m、60m Sentinel2 L1C数据

加载20m和60m Sentinel2 L1C数据和加载10m数据类似，在图1中选择相应的分辨率即可。

3. 批量处理10m、20m、6m Sentinel2 L2A数据

同样创建框架，执行菜单命令 任务管理->添加其它卫星观测数据->添加 Sentinel 数据-> Sentinel2 L2A /批量，注意这里是L2A的那行菜单命令。其它与L1C数据处理一样。

还有在图2的对话框中待搜索的文件名为 MTD_MSIL2A.xml，这是Sentinel2 L2A数据的属性文件名。其它输入和输出参数与L1C数据相同。

L2A输出数据与L1C有的不一样。10m数据输出与L1C相同，20m输出2、3、4、5、6、7、8A、11、12共9个波段。60m输出 1、2、3、4、5、6、7、8A、9、11、12共11个波段。其中10波段卷云和8波段与8A重复被删除了。

4. Sentinel2 L1C所有波段提升至10m的批量处理

Sentinel2数据有10m、20m和60m三种分辨率。如果你的项目需要10m经度，那么10m的数据没有问题。但是当需要使用20m或者60m数据时，就需要对这些数据进行重采样，以适应项目10m需要。既然一定需要重采样，这里推荐使用RSD的LiIFT提升分辨率（参阅《一种光谱无损的遥感数据融合方法——LiIFT插值融合（一）》），将Sentinel2所有分辨率数据统一提升至10m。LiIFT提升后无论在保持光谱不变和使用便利等方面都有很大优势，除了增加了一些存储空间。

4.1 L1C提升分辨率的批量处理

执行菜单命令 开始->打开其它对地观测卫星数据->打开Sentinel数据->Sentinel2 L1C Unified 10m 。点击后对话框如图6。

图6 批量处理选择

因为全部处理为1m0，图6中的分辨率选择全部被变灰禁止。下面的 文件列表批量加载 选中表示批量处理，不选时打开单独的数据集进行处理。

点击确定后，后面的操作与“1.批量处理10m Sentinel2 L 1C数据”一样。

处理结束以后，22个数据集同样创建22个数据目录，见图5左侧部分。但是每个目录里面导出的文件不同。还有一个不同是处理时间增加了，22个数据集大概需要处理2个小时。

分辨率全部提升至10m后每个数据集的导出结果见图7。

 图7 分辨率提升至10m后的导出数据

可见有3种数据的.rsd格式、.img格式和.tif格式的导出数据。

4.2 L1C 10m分辨率数据

第1种本来就是10m数据，含2、3、4、8共4个波段。见图8。

图8 Sentinel2 L1C 10m波段图像

4.3 L1C 20m分辨率数据提升至10m数据

第2种是20m提升至10m数据，含 5、6、7、8A、11、12共六个波段，原始图像二次线性重采样和LiIFT分辨率提升比较见图9。

图9 Sentinel2 L1C 20m波段分辨率LiIFT提升（上）与二次线性插值（下）比较

4.4 L1C 60m分辨率数据提升至10m数据

第三种是60m提升至10m数据，含1、9、10共3个波段，原始图像二次线性重采样和LiIFT分辨率提升比较见图10。

图10 Sentinel2 L1C 60m波段分辨率LiIFT提升（上）与二次线性插值（下）比较

5. Sentinel2 L2A所有波段提升至10m的批量处理

执行菜单命令 开始->打开其它对地观测卫星数据->打开Sentinel数据->Sentinel2 L2A Unified 10m 。在出现图2的对话框中待搜索的文件名为 MTD_MSIL2A.xml，这是Sentinel2 L2A数据的属性文件名。其它L1C数据相同。

Sentinel2 L2A所有波段提升至10m后波段与L1C有点差异。10波段卷云大气校正后删除了，8波段和8A重复，也被删除了。

6. 小结

① 批量处理和导出.rsd、.img、.tif格式数据，自动处理提高了工作效率。

② 可以向预定义框架加载数据，实现自动剪切。

③ 推荐L1C/L2A分辨率提升，高质量完全光谱无损的融合效果。

请下载RSD专业版 SetupRsd3.2.6。下载地址：https://pan.baidu.com/s/1T-LBvaD_zVCwJsGf_hCyCg

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