1

在 Sentinel-2 中，如何从经纬度地理坐标中获知 S2 瓦片？

4

2 回答 2

1

您必须登录https://scihub.copernicus.eu/dhus才能请求 API 的用户名和密码。

然后，您可以使用此功能：

lat是地理纬度，lon经度和tiles包含网格中的瓦片：

    def S2tile_fromLATLON(float(lat),float(lon)):
        
        # query scenes
        api = SentinelAPI('USER', 'PASSWORD', 'https://scihub.copernicus.eu/dhus')
        
        footprint = 'POINT(%s %s)' % (lon, lat)
        
      
        product = api.query(footprint, 
                        date=('20190101', '20190301'), 
                        platformname='Sentinel-2', 
                        producttype= 'S2MSI1C', 
                        area_relation='Contains',
                        )
        # get tile
        tiles=[]
        for value in product.values():
            tile = value['tileid']
            if len(tiles)==0:
                print(tile)
                tiles.append(tile)
            aux=0
            for j in range(0,len(tiles)):
                if tile==tiles[j]:
                    aux=1
            if aux==0:
                print(tile)
                tiles.append(tile)
                  
        return tiles

您还可以使用该.Kml文件找到图块。但是，如果您需要查找多个坐标的图块或自动化流程，建议使用 API。API 结果是来自给定坐标的所有图块，包括重叠的图块。

于 2021-10-24T19:27:48.860 回答
0

一个迟到的答案，但也许有人会需要它。

如果您不想使用 API（或帐户或连接），则有一个包含图块定义的 KML 文件：https ://sentinel.esa.int/documents/247904/1955685/S2A_OPER_GIP_TILPAR_MPC__20151209T095117_V20150622T000000_21000101T000000

您可以解析它并根据坐标制作多边形并检查您的点是否在多边形内（例如，形状优美）。不确定这是否会比 API 更快。

不要忘记图块有重叠，即一个坐标可能同时存在于多达 4 个图块中。

于 2021-10-23T11:04:41.963 回答