关于Sentinel-2快速查询图幅号——使用MGRS_100kmSQ_ID_File_Geodatabase快速查询

关于Sentinel-2快速查询图幅号—— 使用MGRS_100kmSQ_ID_File_Geodatabase快速查询

一、前言

相信大家在欧空局查询Sentinel-2的时候，总是摸不清清楚自己的研究区需要几个图幅来覆盖，或者对于某个图幅的具体覆盖范围不是很明确，遇到网络卡顿的时候，页面迟迟刷不出来，下面我们就来通过ArcGIS来快速查询自己研究区的图幅号。

Sentinel-2的幅宽高达290KM，比Landsat的185KM要大得多，基本上一个条带扫下来可以覆盖大半个省：

欧空局为了方便数据管理和大家的使用将每个条带的数据切成了100KM大小的图幅，如下图所示，这是一个条带的覆盖图，这个条带横跨49和50两个分带，每个哨兵分幅数据的名称里面都有红色框里的分幅图幅号的标记，类似于T50SKD，代表这是UTM50带的SKD图幅：

在欧空局的底图上也能发现明显的条带印记：

二、查询

下面我们以河南省为例来举例说明如何使用ArcGIS自助查询Sentinel-2图幅号：

首先打开ArcMap，加载河南省的行政界线：

然后我们根据UTM投影图来确定大概的带号，关于UTM分带的知识大家可以回顾地图学的相关章节或者在网上自行学习，这里我就不多说了，数学不好的同学可能学起来比较难懂。

这是一张全球的UTM分带示意图：

我们放大一些：
河南大概在粉色圈的位置，东西垮了49-50两个分带，南北位于S和R两个分区，那我们就记住49-50的S和R区间。

下面我们加载MGRS_100kmSQ_ID_File_Geodatabase：
使用Ctrl键将49分带的S和R和50分带的S和R都加载进来。


加载进来是不是发现正好覆盖住了河南：

我们可以将图幅号进行标注显示：

这样河南的图幅号是不是都很直观的显示出来了呢：

我们也可以将它们进行合并后空间筛选：

附录：
MGRS_100kmSQ_ID_File_Geodatabase资料
链接: https://pan.baidu.com/s/1nSEgz6z183JgW5HFoh0JvA 提取码: ffbn

三、拼接

关于哨兵2的拼接问题，例如下面是通过SNAP采样后输出的ENVI格式，在ENVI里进行波段组合后的三幅数据：

由于这三幅数据属于同一个日期拍的同一个条带的数据，严格的说他们是属于一大景数据，所以拼接的时候不要选直方图匹配，但是可以生成一个拼接线：


拼接后的结果：