首页 > 解决方案 > Sentinel-2:如何从纬度/经度地理坐标中了解 S2 瓦片?

问题描述

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

标签: pythontilesentinel2

解决方案


您必须登录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 结果是来自给定坐标的所有图块,包括重叠的图块。


推荐阅读