欧空局中国区域数据的筛选

0.背景

最近写了很多ESA土地利用数据的推文，待数据公布后，我也一直在下载。

现实很骨感，一直没下载成功。后来有一位大佬下载了，才得以拥有数据,在此表示感谢:

1.数据筛选

将数据下载后，一共2600余景:

根据中国大陆区域的四至经纬度，编写代码筛选出中国的区域:

#该脚本用于移动影像文件到指定目录，并配套分类原则
import shutil
import os
from tqdm import tqdm
#原始文件位置
file_path=r"原始文件位置"

#第一部分要存放的文件夹的路径
china_file_path=r"存放路径"

#文件列表
lists=[]
for lan in range(18,72,3):
    for long in range(54, 135, 3):
        if long<100:
            long="0"+str(long)
            list_LONG_LAN="ESA_WorldCover_10m_2020_v100_N"+str(lan)+"E"+str(long)+"_Map.tif"
        else:
            list_LONG_LAN = "ESA_WorldCover_10m_2020_v100_N" + str(lan) + "E" + str(long) + "_Map.tif"
        lists.append(list_LONG_LAN)
        print(list_LONG_LAN)
# 第二部分，将名称为file的文件复制到名为file_dir的文件夹中
for list in tqdm(lists):
    try:
        file_src_path = os.path.join(file_path, list)
        #筛选出要存放的文件夹的路径
        tif_file_path = china_file_path+ list
        #将指定的文件file移动到file_dir的文件夹里面
        shutil.copy(file_src_path,tif_file_path)
    except:
             pass

如果需要小岛屿的土地利用数据，可以自己去官网下载。这里主要分享大陆架覆盖区域。

经过纬度18°-72°、经度54°到135°的范围筛选，共获得469景土地利用影像:

2.数据获取

前往gzh，回复“ESA2”获取链接