欧空局中国区域数据分享

0.背景

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

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

1.数据筛选

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

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
#该脚本用于移动影像文件到指定目录,并配套分类原则
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.数据获取

个人主页


文章作者: 锐多宝
版权声明: 本博客所有文章除特別声明外,均采用 CC BY 4.0 许可协议。转载请注明来源 锐多宝 !
评论
 上一篇
2022-03-28 锐多宝
下一篇 
2022-03-28 锐多宝
  目录