100 GIS practical application cases (34) - splicing 2020globeland30

2020globeland30

See... For product introduction http://www.globallandcover.com/Page/sysFrame/dataIntroduce.html?columnID=81&head=product&para=product&type=data

in general , The main data contains projection GeoTIFF Format , altogether 966 individual , Every tif Together with the corresponding coordinate information file 、 The classified image is connected to the chart file 、 The metadata file is compressed into a zip, And then all zip Then compress it into a zip. The whole package has 7.18G Big .

Purpose

take 966 individual tif Spliced into global latitude and longitude data .

The process

Plan to use gdal Read tif then xarray Splicing , If gdal It is of course best to be able to directly read the specified data in the compressed package , Unfortunately, I can't . Then just tif Unpack ,“ useless ” Other files do not have to waste the mood to unzip .

decompression

 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