2020globeland30
产品介绍详见http://www.globallandcover.com/Page/sysFrame/dataIntroduce.html?columnID=81&head=product¶=product&type=data
总的来说，主要数据是含有投影的GeoTIFF格式，一共966个，每个tif连同对应的坐标信息文件、分类影像接图表文件、元数据文件压缩成一个zip，然后所有zip再压缩成一个zip。整个压缩包有7.18G大。
目的
将966个tif拼接成全球经纬度数据。
过程
计划用gdal读取tif然后xarray拼接，如果gdal能直接读取压缩包里指定数据当然是最好，可惜不能。那么就只把tif解压出来，“没用的”其他文件就不必浪费心情解压出来了。
解压
1import zipfile
2from os.path import exists, basename
3from os import makedirs
4
5
6globeland30_zippath = 'G:/GEO_DATA/landcover/globeland30/2020LC030.zip'
7globeland30_tifdir = 'G:/GEO_DATA/landcover/globeland30/tif'
8if not exists(globeland30_tifdir):
9 makedirs(globeland30_tifdir)
10globeland30_zipfile = zipfile.ZipFile(globeland30_zippath)
11for zipinfo in globeland30_zipfile.filelist:
12 a_zipfile = zipfile.ZipFile(globeland30_z
