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米土地利用数据代码相关推荐

  1. 数据大方送之全球10米土地利用数据

    一.前言 很多小伙伴苦于寻找不到高分辨率的土地利用数据,我国科研团队制作的GlobeLand30全球地理信息公共产品,但是最近该网站处于维护阶段,也无法下载. 普天同庆Esri公司在其官网分享一个10 ...

  2. 10米土地利用数据按行政区下载

    0.背景 由于工作原因,接触PIE-Engine的机会比较多.这几天用PIE-Engine时,发现了里面有FROM_GLC的10米分辨率数据(2017年). 有这个数据,就可以直接用PIE-Engin ...

  3. 欧空局10米土地利用数据下载教程

    1.网站注册 首先进入官网(https://viewer.esa-worldcover.org/worldcover) ,然后,点击Register: 然后填写真实信息(鼓励大家填写真实信息,以便于数 ...

  4. 数据大放送之10米土地利用数据(2017-2022年)

    一.前言 土地覆盖数据是我们在各项研究中都非常常用的数据,土地覆盖数据的来源也有很多,Esri10米土地覆盖数据的年份为2017--2022年,坐标为WGS_1984_UTM,范围为全球范围,全球划分 ...

  5. 2021年最新全国30米土地利用数据(中科院二级分类)

    中科院全国30米土地利用数据(1980-2018) 数据来源:中科院地理所/地理国情监测云平台 数据名称:中国全国30米高精度土地利用现状土地利用类型遥感监测数据 1980年-2018年中国土地利用现 ...

  6. 全国30米土地利用数据分享(1980-2020)

    中科院全国30米土地利用数据(1980-2020) 数据来源:中科院地理所/地理国情监测云平台 数据名称:中国全国30米高精度土地利用现状土地利用类型遥感监测数据 1980年-2020年中国土地利用现 ...

  7. 批量下载1985-2019年研究区的30米土地利用数据

    0.背景 PIE-Engine可以下载的6种土地利用数据(仅限本人所知),有10m.30m.1000m分辨率. 这两天连发了两篇使用PIE-Engine下载数据的推文,这一篇是终结篇.剩余的1000m ...

  8. 全球10米遥感数据免费下载及预处理(附软件下载)

    看标题估计大家都猜出来了是Sentinel-2(哨兵2),没错,就是Sentinel-2. 之前我们推过关于高清谷歌数据的下载方法:谷歌卫星影像免费下载?!来看这些软件,今天我们就来看看10米分辨率的 ...

  9. Google Earth Engine(GEE)——基于s2影像的实时全球10米土地利用/土地覆盖(LULC)数据集

    全球动态土地分类数据集是一个10米的近实时(NRT)土地利用/土地覆盖(LULC)数据集,包括九个类别的概率和标签信息. 动态世界的预测适用于2015-06-27至今的Sentinel-2 L1C集合 ...

最新文章

  1. ImportError: No module named pil
  2. acm公选课笔记 2020.3.31
  3. ACM 会场安排问题
  4. BeetleX进程服务管理组件应用
  5. python计算定积分_python编程通过蒙特卡洛法计算定积分详解
  6. mac中一一些常用的命令
  7. chap6_2 Parallax mapping in OGRE
  8. [工具]Visual Studio
  9. npp夜光数据介绍 viirs_NPP-VIIRS年度夜间灯光数据的合成方法与验证
  10. JavaScript计算器(加减乘除完善)
  11. 阿里云 centos7.2 安装openstack报错
  12. 套接字Socket的常见面试题及答案
  13. 获取WI-FI信息以及蜂窝网络信息
  14. 利用微搭低代码操作微信缓存
  15. python画长方形的代码_python使用turtle画一个三角形、正方形(矩形或四边形)
  16. android sim卡槽,一加5怎么装卡/插卡 一加手机5 SIM卡安装图文教程
  17. 六大维度层层剖析,一篇文章带你快速读懂信息无障碍
  18. 消息中心构架设计说明书
  19. 电脑开机自动登录,断网自动重连以及桌面联网小程序
  20. 面临公司变相裁员该如何应对

热门文章

  1. 运行程序报错:请检查是否存在数组、列表等越界非法访问,内存非法访问等情况
  2. 【JavaScript-进阶】详解数据类型,内存分配,API元素对象获取
  3. Unity结合HTC Vive开发之VRTK手柄Touchpad键圆盘UI显示
  4. 什么叫爬虫python_网络爬虫是什么?
  5. 开发者工具绕过Paused in debugger
  6. 【MemoryCompression内存占用过高】
  7. Lotus Notes中文档查询(转)
  8. 苏宁大造818发烧节,玩得是哪招?
  9. Javaswing+MySQL实现超市收银管理系统(含源码、包运行)
  10. pdf阅读器如何语音朗读文件