rgee入门之Sentinel-2影像加载与真彩色、假彩色合成

今天尝试一下使用rgee加载Sentinel-2影像。

主要内容

• 创建Image Collection
• 按时间和经纬度点位筛选数据
• 遥感影像的波段选择与显示

实战过程

• 选择数据集
• 可以在Earth Engine Data Catalog中查看数据集介绍，挑选数据
• 使用Earth Engine Snippet中的代码加载引用数据

数据介绍和加载方式

• 使用日期筛选数据
• 指定start/end
• 使用经纬度点位筛选数据
• 经纬度坐标

library(rgee)
ee_Initialize()
col<-ee$ImageCollection('COPERNICUS/S2_SR')
point <- ee$Geometry$Point(119,38)
start <- ee$Date("2019-07-11")
end <- ee$Date("2019-07-20")
filter<-col$filterBounds(point)$filterDate(start,end)
img <- filter$first()

• 设置可视化参数
• 波段
• 值域
• gamma拉伸

# RGB Visible
vPar <- list(bands = c("B4", "B3", "B2"),min = 100,max = 8000,
gamma = c(1.9,1.7,1.7))

• 设置地图显示参数
• 定义地图中心和级别
• 定义地图叠加内容和参数

Map$setCenter(119,39, zoom = 10)
Map$addLayer(img, vPar, "True Color Image")

程序执行完后，就可以从Viewer面板中查看加载的影像了。

执行结果

点击Zoom，可以弹出一个Viewer Zoom窗口，可以放大查看加载的影像。

真彩色合成的影像

调整一下前面波段选择的参数，这样我们就能得到一幅假彩色合成的影像

vPar <- list(bands = c("B8", "B4", "B3"),min = 100,max = 8000,
             gamma = c(1.9,1.7,1.7))
Map$setCenter(119,39, zoom = 10)
Map$addLayer(img, vPar, "假彩色合成")

数据需要从谷歌服务器上面下载，加载会比较慢，如果影像看着模糊就是没加载出来，需要等一会。下图就是加载号的标准假彩色合成的Sentinel-2数据了。

标准假彩色合成影像

更多文章请关注微信公众号：走天涯徐小洋地理数据科学