PIEGEE下载三款10米土地利用数据代码
0.背景
GEE已经上线了Esri和ESA的10米土地利用数据,PIE也上线了FROM_GLC的10米土地利用数据。汇总了这三款数据的下载方法。
1.ESA的10米数据下载
第一步是设置参数,并加载ESA的数据。
//显示参数设置
var imageVisParam_ESA = {"opacity":1,"bands":["Map"],"min":10,"max":110,"palette":["8dffda","14ff31","686dff","ff9b54","ff712d","ff66b4"]};
//浏览全球土地利用数据
var ESA = ee.ImageCollection("ESA/WorldCover/v100")
Map.addLayer(ESA,imageVisParam_ESA,"ESA_landcover")
第二步加载待下载区域ROI。
我的代码上传的是四川省_资阳市_乐至县的行政区数据。
//导入研究区
var ROI = ee.FeatureCollection("users/kitmyfaceplease2/lezhixian").geometry()
Map.centerObject(ROI, 8)
Map.addLayer(ROI, {color: 'FF0000', fillColor: '00000000', width: 1}, "ROI")
第三步是筛选影像,并裁剪:
需要注意是ESA数据在GEE中已经被镶嵌为一张全球影像,可以直接clip。
//筛选数据
var ESA_ROI_landcover=ESA.filterBounds(ROI).mosaic().clip(ROI)//研究区土地利用影像
Map.addLayer(ESA_ROI_landcover,imageVisParam_ESA,"ESRI_ROI_landcover")
第三步是导出影像:
//下载
Export.image.toDrive({image:ESA_ROI_landcover,description: "ESA_ROI_landcover",region:ROI,scale:10,maxPixels:1e13
});
2.ESRI的10米数据下载
ESRI数据的下载大同小异,只需注意一个地方:ESRI数据在GEE中并没有拼接,需要mosaic。
第一步是设置参数,并加载ESRI的数据。
//显示参数设置
var imageVisParam_ESRI = {"opacity":1,"bands":["b1"],"min":1,"max":10,"palette":["8dffda","14ff31","686dff","ff9b54","ff712d","ff66b4"]};//浏览全球土地利用数据
var ESRI = ee.ImageCollection("projects/sat-io/open-datasets/landcover/ESRI_Global-LULC_10m")
Map.addLayer(ESRI,imageVisParam_ESRI,"ESRI_landcover")
第二步加载待下载区域ROI。
//导入研究区
var ROI = ee.FeatureCollection("users/kitmyfaceplease2/lezhixian").geometry()
Map.centerObject(ROI, 8)
Map.addLayer(ROI, {color: 'FF0000', fillColor: '00000000', width: 1}, "ROI")
第三步是筛选影像,并裁剪:
ESRI数据需要使用mosaic函数。
//筛选数据
var ESRI_ROI_landcover=ESRI.filterBounds(ROI).mosaic().clip(ROI)//研究区土地利用影像
Map.addLayer(ESRI_ROI_landcover,imageVisParam_ESRI,"ESRI_ROI_landcover")
第三步是导出影像:
//下载
Export.image.toDrive({image:ESRI_ROI_landcover,description: "ESRI_ROI_landcover",region:ROI,scale:10,maxPixels:1e13
});
3.FROM_GLC的10米数据下载
下载FROM_GLC数据需要切换到PIE-Engine的工作环境。
第一步研究区的加载
//研究区加载 LZ(四川省_资阳市_乐至县)
LZ= LZ.first().geometry();
Map.centerObject(LZ,9);
Map.addLayer(LZ, {color: 'FF0000', fillColor: '00000000', width: 1}, "LZ")
第二步数据的筛选
首先加载影像集,筛选时间,选择波段,镶嵌与裁剪。
//FROM_GLC10 土地利用数据下载
//加载土地利用数据 筛选、镶嵌、裁剪
var LZ_landcover_2017=pie.ImageCollection("THU/FROM_GLC10_2017").filterBounds(LZ).filterDate("2017-1-01", "2017-12-31").select(["B1"]).first().clip(LZ);//显示土地利用数据
Map.addLayer(LZ_landcover_2017,visParam,"LZ_landcover_2017")
在这里,需要注意一个地方:
PIE与GEE的裁剪不同,需要指定波段才能进行裁剪。
第三步进行数据的导出:
这里几乎和GEE差不多,唯一的差异是多了个assetId。
//下载
Export.image({image:LZ_landcover_2017,description: "LZ_landcover_2017",assetId: "LZ_landcover_2017",region:LZ,scale:10,maxPixels:1e13
});
4.代码链接
ESA 10米 2020年 土地利用数据 链接:
https://code.earthengine.google.com/6ba5ff6808c532bc20bce44eeaf85c5e
ESRI 10米 2020年 土地利用数据 链接:
https://code.earthengine.google.com/da15793f593b7ce6470fc3f211dff75d
FROM_GLC 10米 2017年 土地利用数据 链接:
https://engine.piesat.cn/engine/home?sourceId=95120018f4ae488ea92c1f8b35b4c3ee
5.结果展示
6.写在最后
我们放大地图到乐至县的城区:
该图为2020年ESA的10米数据。
该图为2020年Esri的10米数据。
左边为2020年的GLC30的数据,右边为2017年FROM_GLC数据。
除非我的眼睛有问题,号称准确率80%以上的Esri数据真比不上GLC的30米。如果不是做数据对比,Esri数据强烈不建议大家使用。
可以前往“地信遥感数据汇”(https://www.gisrsdata.com/)获取更多数据。
PIEGEE下载三款10米土地利用数据代码相关推荐
- 数据大方送之全球10米土地利用数据
一.前言 很多小伙伴苦于寻找不到高分辨率的土地利用数据,我国科研团队制作的GlobeLand30全球地理信息公共产品,但是最近该网站处于维护阶段,也无法下载. 普天同庆Esri公司在其官网分享一个10 ...
- 10米土地利用数据按行政区下载
0.背景 由于工作原因,接触PIE-Engine的机会比较多.这几天用PIE-Engine时,发现了里面有FROM_GLC的10米分辨率数据(2017年). 有这个数据,就可以直接用PIE-Engin ...
- 欧空局10米土地利用数据下载教程
1.网站注册 首先进入官网(https://viewer.esa-worldcover.org/worldcover) ,然后,点击Register: 然后填写真实信息(鼓励大家填写真实信息,以便于数 ...
- 数据大放送之10米土地利用数据(2017-2022年)
一.前言 土地覆盖数据是我们在各项研究中都非常常用的数据,土地覆盖数据的来源也有很多,Esri10米土地覆盖数据的年份为2017--2022年,坐标为WGS_1984_UTM,范围为全球范围,全球划分 ...
- 2021年最新全国30米土地利用数据(中科院二级分类)
中科院全国30米土地利用数据(1980-2018) 数据来源:中科院地理所/地理国情监测云平台 数据名称:中国全国30米高精度土地利用现状土地利用类型遥感监测数据 1980年-2018年中国土地利用现 ...
- 全国30米土地利用数据分享(1980-2020)
中科院全国30米土地利用数据(1980-2020) 数据来源:中科院地理所/地理国情监测云平台 数据名称:中国全国30米高精度土地利用现状土地利用类型遥感监测数据 1980年-2020年中国土地利用现 ...
- 批量下载1985-2019年研究区的30米土地利用数据
0.背景 PIE-Engine可以下载的6种土地利用数据(仅限本人所知),有10m.30m.1000m分辨率. 这两天连发了两篇使用PIE-Engine下载数据的推文,这一篇是终结篇.剩余的1000m ...
- 全球10米遥感数据免费下载及预处理(附软件下载)
看标题估计大家都猜出来了是Sentinel-2(哨兵2),没错,就是Sentinel-2. 之前我们推过关于高清谷歌数据的下载方法:谷歌卫星影像免费下载?!来看这些软件,今天我们就来看看10米分辨率的 ...
- Google Earth Engine(GEE)——基于s2影像的实时全球10米土地利用/土地覆盖(LULC)数据集
全球动态土地分类数据集是一个10米的近实时(NRT)土地利用/土地覆盖(LULC)数据集,包括九个类别的概率和标签信息. 动态世界的预测适用于2015-06-27至今的Sentinel-2 L1C集合 ...
最新文章
- ImportError: No module named pil
- acm公选课笔记 2020.3.31
- ACM 会场安排问题
- BeetleX进程服务管理组件应用
- python计算定积分_python编程通过蒙特卡洛法计算定积分详解
- mac中一一些常用的命令
- chap6_2 Parallax mapping in OGRE
- [工具]Visual Studio
- npp夜光数据介绍 viirs_NPP-VIIRS年度夜间灯光数据的合成方法与验证
- JavaScript计算器(加减乘除完善)
- 阿里云 centos7.2 安装openstack报错
- 套接字Socket的常见面试题及答案
- 获取WI-FI信息以及蜂窝网络信息
- 利用微搭低代码操作微信缓存
- python画长方形的代码_python使用turtle画一个三角形、正方形(矩形或四边形)
- android sim卡槽,一加5怎么装卡/插卡 一加手机5 SIM卡安装图文教程
- 六大维度层层剖析,一篇文章带你快速读懂信息无障碍
- 消息中心构架设计说明书
- 电脑开机自动登录,断网自动重连以及桌面联网小程序
- 面临公司变相裁员该如何应对