【教程】利用GEE工具进行区域遥感影像的获取与下载

在日常工作中，我们经常需要获取某一区域的遥感影像数据。在GEE工具（全文GEE工具都应改为rgee ）出现之前，要完成这一工作的往往需要经过影像下载、几何校正、去云处理、影像拼接、矢量裁剪等多个步骤，本期分享如何用GEE工具进行区域遥感影像的获取与下载。

1 rgee初始化

(1) 之前rgee初始化命令不同的是，本案例中涉及到将GEE生成的遥感影像下载到本地，为了实现这一步，需要在rgee初始化语句中定义“drive = TRUE”：

library(rgee)
ee_Initialize(drive = TRUE)

(2) 第一次设置需要给Tidyverse API Packages授权可以访问Google Drive的权限，操作步骤如下：

a. 出现图1中“Google Drive credentials”，输入“2”，即没有Google Drive资格证书，会自动打开网页：

b. 在图2中页面选择对应的谷歌账号：

c. 在图3页面中勾选“查看、修改、创建和删除Google云端硬盘文件”，点击“继续”；

d. 将下图重点的代码复制粘贴到R中：

e. 授权成功：

2 影像的筛选、拼接与裁剪

a. 以提取江苏省的Landsat 8影像为例，代码中使用的矢量边界是提前上传到GEE的，也可以使用st_read()读取本地的shp文件，代码如下：

# 获取江苏省的矢量边界
Jiangsu <- ee$FeatureCollection("users/JinhuaWu/shp/Jiangsu")
Map$addLayer(Jiangsu)

# 影像拼接与裁剪
Jiangsu_image <- ee$Algorithms$Landsat$simpleComposite(
ee$ImageCollection('LANDSAT/LC08/C01/T1')$
filterDate('2016-01-01', '2018-12-31')
)$clip(Jiangsu)

# 定义可视化参数
vizParams <- list(
bands = c("B4", "B3", "B2"),
min = 15,
max = 108,
gamma = 1.3
)

# 影像加载
Map$addLayer(Jiangsu_image, vizParams, "Jiangsu Image")

b. 得到结果：

3 影像下载

(1) 通过ee_as_raster()将影像先传输到Google drive后，再下载保存到本地，代码如下：

# 影像下载
ee_as_raster(
Jiangsu_image,
dsn = 'D:/Downloads/Jiangsu_L8.tif',
region = Jiangsu$geometry(),
via = 'drive',
scale = 30,
maxPixels = 1e13,
)

(2) 江苏省的Landsat 8影像数据大小约1个G，下载的数据用ArcGIS正常打开：

-----END-----

社群交流 / 原创投稿 / 商务合作

（请添加下方小助手微信）

来源：生态遥感笔记

推荐阅读

推荐关注

温馨提示：近期，微信公众号信息流改版。每个用户可以设置 常读订阅号，这些订阅号将以大卡片的形式展示。因此，如果不想错过“测绘之家”的文章，你一定要进行以下操作：进入“测绘之家”公众号 → 点击右上角的 ··· 菜单 → 选择「设为星标」

↓↓↓点击下方“阅读原文”查看更多精彩内容...