当前位置:网站首页>100 GIS practical application cases (34) - splicing 2020globeland30
100 GIS practical application cases (34) - splicing 2020globeland30
2022-04-23 13:02:00 【Gatulo book】
2020globeland30
See... For product introduction http://www.globallandcover.com/Page/sysFrame/dataIntroduce.html?columnID=81&head=product¶=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
版权声明
本文为[Gatulo book]所创,转载请带上原文链接,感谢
https://yzsam.com/2022/04/202204230612571856.html
边栏推荐
- Object. The disorder of key value array after keys
- async void 导致程序崩溃
- 拥抱机器视觉新蓝海,冀为好望开启数字经济发展新“冀”遇
- mysql8安装
- Important knowledge of transport layer (interview, retest, final)
- Async void provoque l'écrasement du programme
- 22. Bracket generation
- Free and open source agricultural Internet of things cloud platform (version: 3.0.1)
- GIS practical tips (III) - how to add legend in CASS?
- Homomorphic encryption technology learning
猜你喜欢
100 GIS practical application cases (51) - a method for calculating the hourly spatial average of NC files according to the specified range in ArcGIS
Record Alibaba cloud server mining program processing
Use compressorjs to compress pictures, optimize functions, and compress pictures in all formats
51 single chip microcomputer stepping motor control system based on LabVIEW upper computer (upper computer code + lower computer source code + ad schematic + 51 complete development environment)
教你快速开发一个 狼人杀微信小程序(附源码)
The El table horizontal scroll bar is fixed at the bottom of the visual window
CVPR 2022 & ntire 2022 | the first transformer for hyperspectral image reconstruction
进程虚拟地址空间区域划分
Object. The disorder of key value array after keys
拥抱机器视觉新蓝海,冀为好望开启数字经济发展新“冀”遇
随机推荐
Temperature and humidity monitoring + timing alarm system based on 51 single chip microcomputer (C51 source code)
Free and open source intelligent charging pile SaaS cloud platform of Internet of things
Byte jump 2020 autumn recruitment programming question: quickly find your own ranking according to the job number
How to convert opencv pictures to bytes
CVPR 2022&NTIRE 2022|首个用于高光谱图像重建的 Transformer
HQL statement tuning
The use of dcast and melt in R language is simple and easy to understand
Use of Presto date function
Idea的src子文件下无法创建servlet
有趣的IDEA插件推荐,给你的开发工作增添色彩
SSM框架系列——注解开发day2-2
[51 single chip microcomputer traffic light simulation]
GIS practical tips (III) - how to add legend in CASS?
Introduction to servlet listener & filter
XinChaCha Trust SSL Organization Validated
Kernel error: no rule to make target 'Debian / canonical certs pem‘, needed by ‘certs/x509_ certificate_ list‘
Record Alibaba cloud server mining program processing
Recommended website for drawing result map
CVPR 2022 & ntire 2022 | the first transformer for hyperspectral image reconstruction
Unable to create servlet under SRC subfile of idea