python - Sentinel-2:如何从纬度/经度地理坐标中了解 S2 瓦片?
问题描述
在 Sentinel-2 中,如何从经纬度地理坐标中获知 S2 瓦片?
解决方案
您必须登录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 结果是来自给定坐标的所有图块,包括重叠的图块。
推荐阅读
- windows - 在whatsapp桌面上分享链接预览问题
- python - 简单的加密算法 Python
- python - Matplotlib savefig 在笔记本后端工作时崩溃
- python - Python 错误:xlsxwriter.exceptions.FileCreateError:[Errno 2] 没有这样的文件或目录:''
- java - 设计说明:在 RabbitMQ 中获取 Oracle 表插入
- ubuntu - 使用 nmcli 从静态切换到 DHCP
- python - 一次打印一定数量的多维数组行
- c# - 为什么有时 async 关键字不与带有 async 后缀的方法一起使用?
- servicenow - 如何在 servicenow 中创建和显示层次结构列表?
- php - 如何在laravel中执行两个where条件,其中order_id相同但product_ids在数组中