欧空局SNAP应用批量处理Sentinel-2数据

zhuzhu1213812138

已于 2024-08-26 13:25:30 修改

阅读量847

点赞数 11

文章标签： sentinel

于 2024-08-26 12:34:52 首次发布

本文链接：https://blog.csdn.net/zhuzhu1213812138/article/details/141553397

版权

因为对python的不熟练，所以学习使用欧空局发布的SNAP应用程序来批处理下载好的sentinel-2的S2A级别数据

1、批量导入.SAFE文件

首先下载好sentinel-2影像，详情可以参考我这篇文章：http://t.csdnimg.cn/jn9XW

打开SNAP软件，选择product library

在repository里面选择AllLocalFolders，点击右边的加号，添加数据所在的本地文件夹。

弹出来一个对话框，我一般选第一个，我的理解就是会挨个搜索该文件夹里的所有子文件夹

然后SNAP就会开始scan文件夹里所有文件。

注意：如果你之前加入过别的文件夹，那么左边folders的下拉框会出现之前添加过的文件夹，等会导入数据的的时候也会顺带着把这里的数据导进去。如果不想混乱，那就在添加新文件夹之前，按一下之前加号旁边的减号，把之前选择的文件夹剔除，再按照上面的步骤加入新文件夹。

数据加载好之后就会这样：

右键点击其中一个影像，然后选择select all，再点击右键选择open（这个操作对每一页的影像都要进行一下，不然你打开的只有一页的20景影像）

打开完成和打开中：（等到所有影像都显示opened了就ok了）

打开完成即可在product explorer看到打开的影像

2、制作数据处理流程脚本（？不知道怎么描述，就是存成一个流程）

>> 1 创建流程

点开上面的graph builder

打开之后长这样：

在上面空白的地方右键（这里演示裁剪、取B2、3、4波段保存成ENVI文件的操作流程）

裁剪：右键-add-Raster-subset

重采样：右键-add-Raster-Geometric-Resample

按照上面点击后，界面会变成下面这样（可以看到加入了subset和resample，read和write是原来就有的）

空白处再点击右键，选择connect graph，界面就会连起来

>> 2 调整每个流程步骤的参数

* read：通常read不用改

* subset：选择下面选项卡的subset，① 选中想要的波段，这里选择了B2/3/4，多选使用ctrl或者shift键，就跟选文件一样。② 裁剪的部分我这里选择使用地理坐标进行裁剪

地理坐标裁剪：① 选择参考波段（因为B1波段分辨率和我输出的波段分辨率不一样，所以我选择了输出的B2波段作为参考波段）； ② 输入裁剪的多边形坐标：输入的格式以及文本获取方式下面会说；③ 点击Update，就可以在地图中看到自己标注的polygon。

输入裁剪的多边形坐标获取方法（本人使用的办法）：

首先，在SNAP中以rgb形式打开真彩影像（下面是局部示意图）：

选择这个画polygon工具

画出你想要裁剪的区域（双击完成绘制）

鼠标选中这个蓝色polygon，右键，选择WKT from Geometry

就可以得到这样的一个框，里面的文本就是这个polygon所有顶点的经纬度，复制这个文本就可以粘贴到上面subset选项中的文本框内。

文字格式如下：

POLYGON ((69.43542995004725 22.80897335988656, 69.43296450961799 22.80293215583949, 69.4377323649236 22.80032652767707, 69.44098455420091 22.79740200939749, 69.44389371780063 22.796528982441913, 69.44656804743549 22.79517326897838, 69.45090364475439 22.80293663748971, 69.44246929551252 22.806307169262155, 69.44246929551252 22.806307169262155, 69.43542995004725 22.80897335988656))